Implement Fermi-dirac integrals
Implement Fermi-dirac integrals