I am working in a RK4 orbital propagator, and I have followed the basic structure of: k1=dt×f(tn,yn) k2=dt×f(tn
I wish to set-up an initially circular (e=0) system of planetary rings which I can later perturb over time and see how the eccentricity changes. However, my cal
I am trying to use the velocity Verlet to simulate the 3-body problem. I am 95% sure that I gravitational force formula and the velocity Verlet algorithm is cor