Converted everything to scipy's integrator. No more reference to my old
"numerics" package. It's faster now! However, blast off to infinity
still occurs. Maybe I can set an adaptive timestep, or something of the sort?