A program for computing the Fermi-Dirac functions - Bañuelos, A.; Depine, R. A.; Mancini, R. C. (Comput. Phys. Comm., 1981)
Double-precision Fortran, maximum accuracy 10D.