Algorithm 779: Fermi-Dirac functions of order -1/2, 1/2, 3/2, 5/2 - MacLeod, Allan J. (ACM Trans. Math. Software, 1998)
Double-precision Fortran, maximum accuracy 16D.