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

These routines compute the incomplete Fermi-Dirac integral with an index of zero, F_0(x,b) = \\ln(1 + e^{b-x}) - (b-x).

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.