用python计算定积分

Python可以使用scipy.integrate库中的quad函数计算定积分。quad函数的参数是被积函数和积分区间的下界和上界。以下是一个计算定积分的示例代码：

from scipy.integrate import quad

def integrand(x):
    return x ** 2 + 2 * x + 1

result, error = quad(integrand, 0, 1)
print("定积分结果：", result)
print("误差估计：", error)

这个例子计算了函数 $x^2 + 2x + 1$ 在区间 $[0,1]$ 上的定积分。函数quad返回两个值，第一个是计算得到的定积分值，第二个是误差的估计值。输出结果如下：

定积分结果： 1.6666666666666667
误差估计： 1.8503717077085944e-14