[comicjk]

Probably a symplectic Verlet method instead of Euler or Runge-Kutta, since you want energy and momentum to be conserved in this simulation.

[gus_massa]

I was going to link to the article about sympletic integration, but it's toooo technical. This other article explains the problem https://en.wikipedia.org/wiki/Energy_drift