\[\left(-\frac{\hbar^{2}}{2m}\frac{{\mathrm{d}}^{2}}{{\mathrm{d}r}^{2}}+\frac{%
\hbar^{2}l(l+1)}{2mr^{2}}+V(r)\right)\psi_{n,l}=\epsilon_{n,l}\psi_{n,l}(r),\]