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

These routines compute the complete Fermi-Dirac integral F_{-1/2}(x).

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.