M. Goano, ALGORITHM 745 - COMPUTATION OF THE COMPLETE AND INCOMPLETE FERMI-DIRAC INTEGRAL, ACM transactions on mathematical software, 21(3), 1995, pp. 221-232
Portable Fortran subroutines computing the Fermi-Dirac integral F-j(x)
and the incomplete Fermi-Dirac integral F-j(x,b) are presented. For t
he first time a set of series expansions is implemented allowing these
special functions to be evaluated efficiently within a prescribed acc
uracy for real j and x;