matlab is being very difficult for me in this question

20 points] The following nonlinear differential equations describe the motion of a body in orbit around two much heavier bodies. An example would be an Apollo capsule in an Earth-moon orbit. The three bodies determine a two-dimensional Cartesian plane in space The origin is at the center of mass of the two heavy bodies, the x-axis is the line through these two bodies, and the distance between their centers is taken as the unit. Thus, if is the ratio of the mass of the moon to that of Earth, then the centers of the moon and the earth are located at coordinates (1-, 0) and (-u, 0), respectively, and the coordinate system moves as the moon rotates about Earth The third body, the Apollo spacecraft, is assumed to have a mass that is negligible compared to the other two bodies, and its position as a function of time is (r(t), y(t)). The equations are derived from Newton's law of motion and the inverse square law of gravitation; the first-order derivatives in the equations come from the rotating coordinate system. The governing equations for motion of the spacecraft are where =82.45 and the initial conditions are z(0)1.2, ' (0) 0, y(0) 0, and y(0)-1.04935751 The solution will be periodic with a period T = 6.19216933. This means that Apollo starts on the far side of the moon with an altitude about 0.2 times the earth-moon distance and a given initial velocity. The resulting orbit brings Apollo in close proximity to earth, out in a big loop in the opposite side the the earth from the moon, back in close to the earth again, and finally back to its original position and velocity on the far side of the moon (a) Solve the equations of motion for the given initial conditions in the domain 0 t 2T using the MATLAB ode45 function. Use the following ode45 call in your program optionsodeset ('RelTol' 1e-6); [TIME, Z] = ode45 ('filename ,, [02*T], [initial conditions], options); and plot Apollo's orbit y(x) for the above conditions. In your plot, clearly mark the position of the moon and the earth, and set the axis properties using axis([-1.5 1.5-1.5 1.5]); axis equal, grid on (b) Try half and double the initial velocity and plot the resulting orbits for 0 t