Here is a heavily used method in astronomy, this involves a higher order polynomial expansion than RK4.

This method has been extended a few times, my code uses a variation of it, and I know of several other projects which are also descended from it.