求助：python中如何写出指数积分的表达式。

崔东山

                               
登录/注册后可看大图

楼主有一个函数的表达式，其中温度变化值及对应的时间已经测得了，Q=q/c，q和r暂时是定为已知的，楼主想通过非线性拟合得出参数c和阿尔法，但是现在卡着就这个函数表达式写不出来，Ei是指数积分（楼主卡在这python中指数积分如何写的）。
求助各位朋友能提供一个思路，感谢了。

崔东山

from scipy.special import *
'''
expn(n,x):
#n:阶数
#x:指数上参数的系数exp(xt)
integral(exp(-x*t) / t**n, t=1..inf).
'''
expn(1, 3)
#>>>0.01304838109419704
from scipy.integrate import quad,dblquad,nquad
quad(lambda t:np.exp(-3*t)/t**1,1,np.inf)
#>>>(0.013048381094197026, 1.1074269344388145e-10)
def expn1(n, x):
    return quad(lambda t:math.exp(-x*t)/t**n,1,inf)
在csdn上面找到一点，抄进去发现不得行，我理解为第6行代码是那个指数积分的表达式写法.....