Cannot get RK4 to solve for position of orbiting body in Python
Physics The Newton law gives you a second order ODE u”=F(u) with u=[x,y]. Using v=[x’,y’] you get the first order system u’ = v v’ = F(u) which is 4-dimensional and has to be solved using a 4 dimensional state. The only reduction available is to use the Kepler laws which allow to reduce the … Read more