Please use MATLAB !!

In this project, we want to solve a differential equation that describes small deflections of a beam that is subjected to lateral loads. The vertical displacement of the beam is represented by a function y (x), where 0 s x s L, along the beam of length L. The displacement y(x) satisfies the following equation Ely"-f(x) Eq. () where E and / are constant along the beam. f(x) is the applied load, including the weight of the beam, in force per unit length. Discretization will convert (1) into a system of linear equations. A reasonable approximation for the fourth derivative is y(x - 2h) -4y(x-h) + 6y(x)-4y(x+h) + y(x + 2h) for small increment h. The discretization error of this approximation is proportional to i?. Our strategy is to consider the beam as the union of many segments of length h, and to apply the discretized version of the differential equation on each segment. For a positive integer n, set h- L. Consider the evenly spaced grid 0 Xo S ?