$$ f(E) = \frac{1}{1+e^{\frac{E-E_f}{K_bT}}} $$