python的scipy.integrate.quad具有复杂的集成边界

积分应该在水平半直线上 a 向右。(如果 一 是一个负实数,但无论如何,这是一个分支切割区域)。此半直线由参数化 a+t 其中t是实的,从0到无穷大。所以,整合 exp(-(a+t))/(a+t) 从0到无穷大。

也, quad 从scipy需要实值函数,所以把它分成实部和虚部。别忘了函数必须作为可调用函数传递,比如 lambda t: np.exp(-t) 不仅仅是 exp(-t)

from scipy import integrate
myfun2_re = integrate.quad(lambda t: np.real(np.exp(-(a+t))/(a+t)), 0, np.inf)[0]
myfun2_im = integrate.quad(lambda t: np.imag(np.exp(-(a+t))/(a+t)), 0, np.inf)[0]
myfun2 = myfun2_re + 1j*myfun2_im
print(myfun2)

这张照片

(0.3411120086192922-0.36240971724285814j)

相比于 myfun ,

array([ 0.34111201-0.36240972j])

顺便说一下,如果您只想转换一个数字,就不需要将myfun包装成一个数组: complex(gammainc(...)) 会的。