please write a matlab code for this

Figure P5.15a shows a uniform beam subject to a linearly increasing distributed load. The equation for the resulting elastic curve is (see Fig. P5.15b) y=120EIw0(x3+2L2x3L4x)(P5.15) Use bisection to determine the point of maximum deflection (i.e., the value of x where dy/dx=0 ). Then substitute this value into Eq. (P5.15) to determine the value of the maximum deflection. Use the following parameter values in your computation: L=600cm,E=50,000kN/cm2,l=30,000cm4, and n0=2.5kN/cm. FIGURE P5.15