辛普森公式求math.exp(-x**2)积分，并求估计误差的python代码

下面是求解math.exp(-x**2)积分并求估计误差的Python代码：

import math

# 定义被积函数
def f(x):
    return math.exp(-x**2)

# 定义辛普森公式
def simpson(a, b, n):
    h = (b - a) / n
    sum1 = 0
    sum2 = 0
    for i in range(1, n, 2):
        sum1 += f(a + i * h)
    for i in range(2, n, 2):
        sum2 += f(a + i * h)
    return (b - a) / 3 / n * (f(a) + f(b) + 4 * sum1 + 2 * sum2)

# 计算积分和估计误差
a = 0
b = 1
n = 100
integral = simpson(a, b, n)
error = abs((1 / math.sqrt(math.pi) - integral) / 1 / math.sqrt(math.pi))

# 输出结果
print('积分结果：', integral)
print('估计误差：', error)

在上面的代码中，我们首先定义了被积函数f(x)，然后定义了辛普森公式simpson(a, b, n)。最后，我们使用定义好的函数计算积分和估计误差，并输出结果。