Computation of the integral Fermi-Dirac function

The Fermi-Dirac functions appear in many areas of physics. The problem of their computation has been studied for a long time. In high energy physics it is also required to compute more complex special functions that use Fermi-Dirac functions. A numerical method is suggested to allow precise computation of one of such functions which is used to calculate the exchange correction to electron energy for arbitrary temperature.