Fermi-Dirac functions
Fermi-Dirac functions as per the Gnu Scientific Library, reference manual section 7.18. These functions are declared in header file gsl_sf_fermi_dirac.h
fermi_dirac_m1(x, give=FALSE, strict=TRUE) fermi_dirac_0(x, give=FALSE, strict=TRUE) fermi_dirac_1(x, give=FALSE, strict=TRUE) fermi_dirac_2(x, give=FALSE, strict=TRUE) fermi_dirac_int(j, x, give=FALSE, strict=TRUE) fermi_dirac_mhalf(x, give=FALSE, strict=TRUE) fermi_dirac_half(x, give=FALSE, strict=TRUE) fermi_dirac_3half(x, give=FALSE, strict=TRUE) fermi_dirac_inc_0(x, b, give=FALSE, strict=TRUE)
x,j,b
: input: real values
give
: Boolean with TRUE
meaning to return a list of three items: the value, an estimate of the error, and a status number
strict
: Boolean, with TRUE
meaning to return NaN
if status is an error
https://www.gnu.org/software/gsl/
Robin K. S. Hankin
x <- seq(from=0,to=2,by=0.01) fermi_dirac_m1(x) #table 7.5 of Ab and St