\qquad G_{n,lj}(r) = \left(1-\frac{E}{mc^2}\right)e^{-\frac{\rho}{2}}(F_1(\rho)-F_2(\rho)),