Abstract
A general numerical method for calculation of the Fermi integrals based on the approximation of a dilogarithm by a series of Chebyshev polynomials is developed. The accuracy of the method proposed is analysed in detail. As an example, the Fermi integral of the model electron density for iron at the Curie temperature is considered.