\[\epsilon_{p,l}=\frac{-Z^{2}}{2\left(p+l+1\right)^{2}}.\]