python 中常用的定积分求解方法

【例1】

解：本例题使用 Scipy 科学计算库的 quad 函数，它的一般形式是 scipy.integrate.quad(f,a,b)，其中 f 是积分的函数名称，a和b分别是下线和上线。

【代码如下】：

import numpy as np
from scipy.integrate import quad
func = lambda x:np.cos(np.exp(x)) **2
solution = quad(func,0,3)

【结果说明】：

输出结果：(1.296467785724373, 1.3977971853986262e-09)。前一个是积分值，后一个是误差。

【例2】

其中D={(x,y) | 0 <=x <=10, 0<= y <=10}。

【代码如下】：

import numpy as np
from scipy.integrate import dblquad
def integrand(x,y):return np.exp(-x**2-y**2)
x_a=0
x_b=10
y_a=0
y_b=10
solution,abserr= dblquad(integrand,x_a,x_b,lambda x :y_a,lambda x:y_b)
print(solution,abserr)

【结果说明】：

输出结果：(0.7853981633974476, 1.3753098510218537e-08)。前一个是积分值，后一个是误差。