Two body problem i(t==u(t)/r(t r(t)=()2+v(t v(t)=-v(t)/r(t) () y(t) u(t)/r(t) v(/r()3」 MATLAB function ydot=twobody(t,y) =sqrt(y(1)2+y(2)^2) ydot=y(3)y(4);-y(1)/r^3;-y(2)r^3] %ydot[y(3:4;y(1:2)norm(y(1:2)^3];Two Body Problem 3 3 ( ) ( ) / ( ) ( ) ( ) / ( ) v t v t r t u t u t r t 2 2 r(t) u(t) v(t) ( ) ( ) ( ) ( ) ( ) v t u t v t u t y t 3 3 ( ) / ( ) ( ) / ( ) ( ) ( ) ( ) v t r t u t r t v t u t y t function ydot =twobody(t,y) r=sqrt(y(1)^2+y(2)^2); ydot=[y(3);y(4);-y(1)/r^3;-y(2)/r^3]; %ydot=[y(3:4;-y(1:2)/norm(y(1:2))^3];