Extend the spline to cover the r=0 point
move integrate_all switch to each function, and switch it to true
More accurate calculation of the rpa kernel

-- Use integration all way through
-- Set regularization parameter epsilon to a smaller value
-- Make its value r-dependent