In a two-dimensional intrinsic semiconductor the density of states is defined as
$$
\begin{gathered}
D_{\mathrm{c}}(\varepsilon)=\frac{4 \pi m_{\mathrm{c}}^{*}}{h^{2}}, \quad \varepsilon>E_{\mathrm{c}} \\
D(\varepsilon)=0, \quad E_{\mathrm{v}}<\varepsilon<E_{\mathrm{c}} \\
D_{\mathrm{v}}(\varepsilon)=\frac{4 \pi m_{\mathrm{b}}^{*}}{h^{2}}, \quad \varepsilon<E_{\mathrm{v}}
\end{gathered}
$$
Use Fermi-Dirac statistics without making any approximation to show that
$$
\varepsilon_{\mathrm{F}}=E_{\mathrm{c}}-\frac{1}{2} E_{\mathrm{g}}+\frac{k_{\mathrm{B}} T}{2} \ln \frac{8}{3}-k_{\mathrm{B}} T \ln \cos \frac{\phi}{3}
$$
with $\phi=\tan ^{-1}\left[\left(\frac{32}{27}\right) \exp \left(-E_{\mathrm{g}} / k_{\mathrm{B}} T\right)-1\right]^{1 / 2}$