I find it surprising that some configurations are as stable as they are with the simple integration method. At least (if I'm reading this correctly) you calculate positions and velocities at different times.

I really like you simulation.