Premium
The accurate numerical evaluation of half‐order Fermi‐Dirac Integrals
Author(s) -
Mohankumar N.,
Natarajan A.
Publication year - 1995
Publication title -
physica status solidi (b)
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.51
H-Index - 109
eISSN - 1521-3951
pISSN - 0370-1972
DOI - 10.1002/pssb.2221880206
Subject(s) - quadrature (astronomy) , mathematics , scheme (mathematics) , range (aeronautics) , trapezoidal rule , dirac (video compression format) , order (exchange) , mathematical analysis , numerical integration , physics , quantum mechanics , optics , materials science , finance , neutrino , economics , composite material
Abstract An efficient algorithm for the direct evaluation of Fermi‐Dirac integrals of half orders is presented. The computational scheme consists of a modified trapezoidal rule which exactly corrects for the loss of accuracy caused by the presence of poles of the integrand. The analysis of the scheme leads to useful error bounds which indicate the quadrature scheme to be rapidly convergent. As much as 16 steps of trapezoidal summation along with at most 9 steps of pole correction can provide an accuracy of one part in 10 14 for the usual range of values of the argument.