Python编写如下代码：用复化simpson自动控制误差方式计算积分。输入是积分区间a,b，精度控制e，函数fx。输出是积分值S

好的，下面是代码实现：

def simpson(f, a, b):
    return (b-a) * (f(a) + 4*f((a+b)/2) + f(b)) / 6

def adapt_simpson(f, a, b, e):
    mid = (a + b) / 2
    S1 = simpson(f, a, b)
    S2 = simpson(f, a, mid) + simpson(f, mid, b)
    if abs(S1 - S2) < 15 * e:
        return S2 + (S2 - S1) / 15
    else:
        return adapt_simpson(f, a, mid, e/2) + adapt_simpson(f, mid, b, e/2)

def integral(f, a, b, e):
    return adapt_simpson(f, a, b, e)

其中，`simpson` 函数是计算 Simpson 积分的函数，`adapt_simpson` 函数是自适应 Simpson 积分的函数，`integral` 函数是对外提供的计算积分的函数。

你可以将你想要计算的函数 `fx` 作为参数传入 `integral` 函数中，例如：

def fx(x):
    return x ** 2

result = integral(fx, 0, 1, 0.0001)
print(result)

这里计算的是 $x^2$ 在 $[0, 1]$ 区间内的积分值，精度控制为 $0.0001$。