Create a matlab code, using ode$45,$ that simulates a bouncing ball.

Using the following:

m$=10$kg

g$=9\mathrm{.}81$

p$=1\mathrm{.}225$

Cd$=0\mathrm{.}285$

Area of the ball$=0\mathrm{.}0314$

Initial V$=202\mathrm{.}22$

Inital Angle $=43$

e$=0\mathrm{.}9$motion of the ball are as follows:

$m{x}^{}=-F_D*\frac{(x^{})}{v},m{y}^{}=-F_D*\frac{(y^{})}{v}-mg,F_D=\frac{1}{2}{v}^2{C}_{D}A,v=\sqrt[\phantom{2}]{x^{}{?}^2+y^{}{?}^2}$

$x(0^{-})=0,y(0^{-})=0,x^{}(0^{-})=v_0*cos(),y^{}(0^{-})=v_0*sin()$

Create two empty figures in MATLAB using figure$()$ command. Each figure would correspond to position and velocity plots respectively.

Implement the ode$45$ solver to solve the differential equations for ball's projectile motion, until the ball hits the ground. You will have to set up an event based ode$45$ termination similar to what was discussed in the lecture.

Generate a plot of the ball's position and velocity in their respective figures. Don't forget to make the axes have the same scaling.

Update the initial time and initial conditions for the projectile motion after the impact with the ground using the equations given below,

$t_i^{-}=t_i,x(t_i^{-})=x(t_i),y(t_i^{-})=y(t_i),x^{}(t_i^{-})=x^{}(t_i),y^{}(t_i^{-})=-e*y^{}(t_i)$

where, $e$ is known as the coefficient of restitution $($COR$)$ which relate the pre$-$impact and post$-$impact velocities. You will use COR of $0\mathrm{.}9$ for this problem.

Go back to step $3$ and repeat the process until

I NEED HELP CREATING THE CODE PLEASE. Thank you