scipy.special.exp1

scipy.special.exp1(x[, out])

y=exp1(z) returns the exponential integral (n=1) of complex argument z: integral(exp(-z*t)/t,t=1..inf).