So it’s basically using the midpoint method? I’m not sure if RK4 is particularly expensive but it’s a bit more of a pain to implement than these explicit methods. Implicit methods are nice though because they handle numerically stiff problems without exploding, and can usually do much larger step sizes.
In any case, you guys did a pretty good job with it =D