Please answer only question #4 and write an accurate Matlab code.

1. A bar with uniform cross-sectional area A, elastic modulus E, and length L is loaded with distributed loads qx). The bar is fixed on the left end and attach to the spring of stiffhess K, which is relaxed when the bar is undeformed. A as .02m2, the Young's modulus E as 210 GPa (Steel ASTM A36), K = 1 000kN/m Spring constant K afx) F(spring) Cross sectional Area: A(m2) L (m) d2u(x) =-q(x):0