N = \int_{0}^{\infty}g(E)f(E)dE = \frac{2}{\Gamma(\frac{n}{2})}\left(\frac{L}{\hbar}\sqrt{\frac{m_e}{2\pi}}\right)^n\int_{0}^{\infty}\frac{(E)^{\frac{n}{2}-1}}{e^{\frac{E-\mu}{k_BT}}+1}dE
$$
...
...
@@ -140,7 +137,7 @@ Check the source code written in python for solving integral using midpoint rule
3.$g(\epsilon)$ vs $\epsilon$ is a linear plot. Here, the region marked by $-k_B T$ is a triangle whose area gives the number of electrons that can be excited:
3.$g(\epsilon)$ vs $\epsilon$ is a linear plot. Here, the region marked by $-k_B T$ is a triangle whose area gives the number of electrons that can be excited: