Accurate Evaluation of Fermi-Dirac Integrals and Their Derivatives for Arbitrary Degeneracy and Relativity

The equation of state of an ideal Fermi gas is expressed in terms of Fermi-Dirac integrals. We give formulae for evaluation the Fermi-Dirac integrals of orders 1/2, 3/2, and 5/2 and their derivatives in various limits of non- and extreme degeneracy and relativity. We provide tables and a Fortran subroutine for numerical evaluation of the integrals and derivatives when a limit does not apply. The functions can be evaluated to better than 1% accuracy for any temperature and density using these methods.