Computes the complete Fermi-Dirac integral with an integer index of j, F_j(x) = (1/\Gamma(j+1)) \int_0^\infty dt (t^j /(\exp(t-x)+1)).
err contains an estimate of the absolute error in the value z.
This function is from the GNU Scientific Library, see http://www.gnu.org/software/gsl/ for documentation.
Package: gsl