Algorithm 779: Fermi-Dirac functions of order -1/2, 1/2, 3/2, 5/2
Allan J. MacLeod · ACM Transactions on Mathematical Software · 1998
The computation of Fermi-Dirac integrals *** is discussed for the values *** = -1, 1/2, 3/2, 5/2. We derive Chebyshev polynomial expansions which allow the computation of these functions to double precision IEEE accuracy.