In this project, you are to simulate the trajectory of several satellites as they maneuver into their new orbits. At initial time, the satellite trajectories are produced by the gravitational force from Earth. An engine will be turned on for a short period of time to produce a thrust needed for the satellites to move to higher orbits. You will simulate and analyze the trajectories of the satellites before, during and after the maneuver. From Newton's second law, the orbital motion of the satellites can be described by the following differential equations: Th GMe GM GMe Thy Ot Th2 Y2+ oX Ot aY where t is time (in seconds); X, Y, and Z are position (in meters) of the satellites relative to the center of Earth in rectilinear coordinate. U, V, and W are the velocity components (in meters per second) in the X, Y and Z directions, respectively. Th". ThV and Th* are three components of thrust (in N) produced by the engine. Earth is assumed to be stationary in this project The parameters in the equation 1 above are given as follows: Radius of Earth: Re 6.37 x106 (m) . Mass of the Earth: Me 5.97 x 1024 (kg) . Gravitational constant: G 6.67408 10-11 (m3 kg",-2) Mass of the satellites: m 1500 (kg) Furthermore, the following quantities will be of interest while analyzing the satellite orbits Altitude of the satellites (m): h = + Y2 + Z2-Re Speed of the satellites (m 8-1): Vmag-v U2 + V2 + W2 Acceleration of the satellites (m 8-2): Acc = d(Vinagydt