n(\varepsilon)=\frac{1}{e^{(\varepsilon-\mu)/k_B T}+1}