How to use the shooting method in matlab

The governing differential equation for the deflection of a cantilever beam subjected to a point load at its free end (Fig. 1) is given by: EI_zz d^2 y/dx^2 = -Px where E is classic modulus, is beam moment of inertia, y is beam deflection, P is the point load, and x is the distance along the beam measured from the free end. The boundary conditions arc the deflection y (L) is zero and the slope (dy/dr) at x=L is zero, where L is beam length. Solve for the deflection of the beam using the shooting method. Use 1=2 m, /77= 0.0005 m4, P= 10 kN, and E= 200,000 MPa. Provide a plot comparing your results against the analytical solution given by: y(x) = - P/6 EI_zz (L - x)^2 (2L + x)