If V is the initial speed of the ball at angle of elevation a, if the origin of the coordinate system is chosen to coincide with the propulsion device, then the initial conditions (t = 0) may be assumed to be
.
Use the Matlab ODE45 solver to obtain the numerical solution of the above differential equations, taking the data for V, a (convert to radians), m, d and c given on the data sheet, g = 9.81 ms^{-2}, and p=1.255kg/m^{3}. Plot the trajectory of the ball (y against x) and compare it with the corresponding trajectory when the drag is zero. Include a printout of the Matlab code you have used.