Posts

Python-Scipy进行数值积分

Python的Scipy模块中拥有大量的数值计算函数，方便我们快速进行数值计算。

Scipy中的integrate模块提供了几种数值积分算法，导入方式为：

from scipy import integrate 

使用integrate时，需要先将要进行积分的方程定义为函数。求取一至三重积分的函数分别为：

integrate.quad(func,a,b,args,full_output) integrate.dblquad(func,a,b,gfun,hfun,args,epsabs,epsrel) integrate.tplquad(func,a,b,gfun,hfun,qfun,rfun,args,epsabs,epsrel) 

以三重积分为例。func为运算对象函数，形式为func(z,y,x)。a,b对应变量x的积分区域，gfun,hfun对应变量y的积分区域，依次类推。

注意gfun,hfun等的形式应为函数，其中gfun,hfun是自变量为x的函数，qfun,rfun是自变量为x,y的函数。这些函数可以使用lambda函数进行定义，形式通常为：

lambda x,y:x*y 

如果是常函数，则定义为：