I need help with this matlab problem.

The deflection $y$ of a sailboat mast subjected to a wind force can be modeled with the following differential equation:

$\frac{d^{2}y}{d{z}^2}=\frac{f(z)}{2EI}(L-z{)}^2$

where $z$ is the position along the mast starting from $0$ at the base. The function $f(z)$ is the force of the wind transferred to the mast by the sail and can be modeled as

$f(z)=\frac{200z}{5+z}{e}^{\frac{-2z}{30}}$

parameter $L$ is the length of the mast.

Write a MATLAB function P$2\_$odesolve that will solve for deflection $y(z)$ and the slope $\frac{dy}{dz}$ given that the deflection and the slope are both equal to $0$ at the base of the mast with the following outputs $($in given order$)$:

Solutions using Euler's method as a $2-$column matrix $($store solution for deflection $y(z)$ in the first column and the solution for slope $\frac{dy}{dz}$ in the second column$)$

Solutions using Heun's method $($without corrector iteration$)$ as a $2-$column matrix $($store solution for deflection $y(z)$ in the first column and the solution for slope $\frac{dy}{dz}$ in the second column$)$

Solutions using Midpoint method as a $2-$column matrix $($store solution for deflection $y(z)$ in the first column and the solution for slope $\frac{dy}{dz}$ in the second column$)$

with the following input

vector defining the span of the mast with intial and final position values $(12$ row vector$)$

initial value vector for the deflection and the slope $(12$ row vector$)$

solver step size $h($scalar$)$ Function $($

$1$ function $[yE,yH,yM]=P2$ odesolve$($zspan$,y0,h)$

$2\%$type your code below, do not change the name of the function and the order of the variables$\%$

Code to call your function $(8)$