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