— Loadable Function: y = fermi_dirac_int (n, x)
— Loadable Function: [y, err] = fermi_dirac_int (...)

These routines compute 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 y.

This function is from the GNU Scientific Library, see http://www.gnu.org/software/gsl/ for documentation.