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