$D_c(\epsilon)=\frac{1}{2\pi^2}(\frac{2m_n}{\hbar^2})^{3/2}(\epsilon-E_c)^{1/2}$ and
$D_v(\epsilon)=\frac{1}{2\pi^2}(\frac{2m_p}{\hbar^2})^{3/2}(E_v-\epsilon)^{1/2}$.
Simplify the Fermi function before calculating the integral and use the substitutions $x=(\epsilon - E_{\text{c}})/k_{\text{B}}T$ and $x=(E_{\text{v}}-\epsilon)/k_{\text{B}}T$.
$D_c(\epsilon)=\frac{1}{2\pi^2}(\frac{2m_n}{\hbar^2})^{3/2}(\epsilon-E_c)^{1/2}$ and
$D_v(\epsilon)=\frac{1}{2\pi^2}(\frac{2m_p}{\hbar^2})^{3/2}(E_v-\epsilon)^{1/2}$.
Simplify the Fermi function before calculating the integral and use the substitutions $x=(\epsilon - E_{\text{c}})/k_{\text{B}}T$ and $x=(E_{\text{v}}-\epsilon)/k_{\text{B}}T$.