Please include script code and output (Use MATLAB, do NOT do by hand)

Using MATLAB R2017b:

?Please write the program according to the instructions!

9. Vertical motion and energy. If an object is shot upward with an initial vertical velocity Vo, its height at any time t can be calculated by the equation: h(t) = ho + vot - where he is the initial height of the object, and g=9.8 m/s, the acceleration due to gravity. The object's potential energy PE (its energy due to its position above the ground). kinetic energy KE (the energy associated with its movement) and total mechanical energy TE (potential energy + kinetic energy) can also be calculated at each time t, using the following equations: PE = mgh KE = 1/2mv TE = PE + KE where m is the mass of the object, and vis the vertical velocity of the object. The vertical velocity of the object is the derivative of the height with respect to time, dh / dt. Write a program, verticalMotion.m, for a 58 g tennis ball that is hit upward from an initial height of 0.8 m at a speed of 60 mph (R17.8 m/s), which computes for Nt=500 timesteps and Tfinal = 4s. a. The height, h, of the ball. b. The ball's vertical velocity, dh/ dt. c. Its potential energy (PE), kinetic energy (KE) and total energy (TE). Using the subplot command (use doc subplot to get information on the subplot command) to plot and animate (with Nstride=20) the following three graphs: a. The height, h, of the tennis ball. b. The vertical velocity, dhfdt, of the ball vs. time. c. The PE, KE, and TE of the ball vs. time. Use the legend command to create a legend that distinguishes the graphs from one another. Each plot should be titled and each axis should be labeled. Note: After the plot is complete, the legend can be dragged using the cursor to get a better view of the energy plot