Exponential integral E_n
Returns the exponential integral for integer n and non-negative x and n:
integral(exp(-x*t) / t**n, t=1..inf).
scipy.special.binom
scipy.special.exp1