improved plotting, pow() should, apparently, not take fractional arguments, though changing 1.5 -> 3/2 really appears to improve integrator performance