Are you stimulating gravity from all the planets, or do you use only the nearest one?
Are you using F = m * M / r^2 with and tweaked the value of M until it look right?
Or you changed the 2 in the exponent?
Some parts of the trajectories in the demo look weird for plain mM/r^2