In this paper, we present a simple and most accurate and economical nu
merical method to evaluate the generalized Fermi-Dirac integral, which
is needed in many astrophysical problems. The method is based on trap
ezoidal integration with pole correction after a suitable variable tra
nsformation.