Prototype: simple orbital model, 1 satellite orbiting 1 body
v1: Euler step for planets and satellites
v2: Planets follow circular path, satellite has integrated solution
v2.1: stops integrating at terminal event (inside planet radius)
v2.2: Minor changes to varilable names