As an aside, I tried to implement this on my laptop at the weekend, but
it was hopelessly unstable. Today it seems very stable indeed. I can
only assume this is to do with replacing Euler with 4th-order
Runge-Kutta. I'm surprised it makes quite this much of a difference
though...
