Algorithm 779: Fermi-Dirac Functions of Order –1/2, 1/2, 3/2, 5/2: The computation of Fermi-Dirac integrals ℱ k is discussed for the values k=-1/2,1/2,3/2,5/2. We derive Chebyshev polynomial expansions which allow the computation of these functions to double precision IEEE accuracy. (Source: http://dl.acm.org/)

This software is also peer reviewed by journal TOMS.