Not the problem, it was the pow(x,1/3), where x < 0
A test that nonpolar works with multipolar
Added multipolar terms; doesn't compile on mac anymore
Appears to be working, additional testing required. Generalizing to handle G&V and G&T correlation integrals remains