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
如果是常函数,则定义为:
lambda x:0
lambda x,y:1
args可选,为传递给func的格外参数;
full_output可选,非零则返回积分信息的dictionary。如果非零,则还会禁止显示警告消息,并将消息附加到输出元组。
epsabs可选,绝对容差直接传递到内部1-D正交积分。默认值为1.49e-8。
epsrel可选,内