Just for fun, I did an n-body simulation of KSP's solar system. The video covers about 80 years starting from the orbital elements given for each planetary body in the KSP wiki, as of version 0.21.1.
The main views on the left show the whole solar system as seen from above and from the side. The center column shows a closeup of the inner solar system (above) and Jool and its moons (below). The right column shows views of the three other planets with moons: Eve, Kerbin and Duna from top to bottom.
Jool very quickly loses a moon, and some years later loses a second one. All the other moons, and the planets themselves, seem relatively stable.
The animation was made using Octave, which gave me fits. I'm open to suggestions for better ways to render the output.
The simulation was performed using an embedded fifth-order Runge Kutta scheme, with the per-step position error limit set at 0.001 m. I've made a version of the program I used to perform this simulation available here: [ Ссылка ]
![](https://s2.save4k.ru/pic/qKp1M4T6z24/maxresdefault.jpg)