Using Matlab to solve for the satellite motion in a central gravitational field
F=-GmM r/r2.
Motion is in a plane.Variables are position(x,y),velocity(vx,vy) and accelration (ax,ay).
Apply the leapfrog method.I need the theoretical solution alone to input into Matlab.
Thanks.