Simulation now works perfectly.

My implementation of RK4 was wrong, which resulted in a spiralling
trajectory of the planets.

The orbits are still spiralling, but at a much, much slower rate. This
is probably due to the instability of the RK4 method.