Sicpy求解定积分

Sicpy求解定积分：

\[\int_1^2 { (x+1)} {\rm d}x \]

可以利用Scipy模块下的子模块integrate里的quad函数来求计算值。

from scipy import integrate# 导入支持库Scipy

def f(x): # 被积函数
    return x + 1

v, err = integrate.quad(f, 1, 2)# 计算定积分

print v # 输入计算结果

结果：2.5