Sampling the Fermi-Dirac density

A method is given for sampling the nonrelativistic Fermi-Dirac electron energy density for all values of the degeneracy parameter eta on the range - infinity < eta less than or equal to 50. The efficiency of the various rejection techniques employed is never less than 30%, and drops below 50% only for a short range of eta values around eta = 2. The range can certainly be extended beyond eta = 50, the efficiency there being 71%, and decreasing very slowly. 5 references.