Altered QM/MM forces routine to also calculate energy terms (what intsol.f does) - first version
-Still need to optimize a lot of things
-Fixed what I think is a small bug in intsol.f