Continued debugging attempts, with experimentation on safer branch
still seing rapoval divergence, but improved initial rguess in root finder
Fixed psi2 error in integrand (hopefully)
Added mathcalG and f, still no resolution for g, divergent integrals in f
Fixed ggood (problem with drhodr)
Merged nukergdebug with master