roll our own RK4 solver

[tmcclintock]

Or RK5. doesn't matter. scipy's solver is slow, so anything will do.

AC:

• let the solver be passable
• demonstrate that our own is faster