scipy.special.expn

scipy.special.expn(x1, x2[, out])

y=expn(n,x) returns the exponential integral for integer n and non-negative x and n: integral(exp(-x*t) / t**n, t=1..inf).