We approximate the Fermi-Dirac integral �by means of composite ten-point Gauss-Legendre quadrature, for values of x in the range �For any x, the integral is approximated by composite quadrature on the interval �which is subdivided into a number of subintervals. We achieve a relative error of no more than �with 10 subintervals, and an error of no more than �with 89 subintervals. On our computational platform, the real-time duration of the computation is faster than two milliseconds for any