Exponential integrals

Syntax: y = EXPINT(x)
y = EI(x)
y = EXPN(x,n)

where         and    

  Elliptic integrals
  Fermi-Dirac distribution