Correct solution with Dirac equation and use of RPA kernel

Background rho error corrected and RPA kernel added to deal with
all RPA contributions. N and alpha included in calculation. Currently
convergent to promising screened potential without N and alpha, but
divergent so far with them.