Write and test MatLAB code implementing the mathematical models
of the boundary value problem to evaluate the elastic deflection of
the beam based on Euler-Bernoulli and Timoshenko theories with
finite difference discretisation (for numerical integration and
differentiation). The results must be plotted on a graph with
labelled local maxima and minima