please help me with these coding question. use a different coding method to solve it. this answer at the botom of the second photo. thank you so much!!!!

This computing language is called Mathematica software

Question 2. Solve the same Euler Bernoulli beam using the virtual work method with a fourth order polynomial defined over the whole domain. Compare the fourth order polynomial obtained with the one obtained using the RR method. In[109] = Clear(EI, L, P, yexact, yapprox, a4, a4s ] yapprox = a + a1*x + a2 + x^2 + a3 + x^3 + a4 *x^4; EI = 1; L = 1; P = 1; th =D [yapprox, x); y1 = yapprox /. x 0; y2 = yapprox/. X L ; th1 = th /. x 0; th2 = th /. xL; a = Solve[(y1 = 0, y2 = 0, th1 = 0, th2 - 0), (ao, ai, a2, a3}); yapprox = yapprox /. a[[1]); ystar = yapprox /. (a4 a4s}; EVW = FullSimplify[-Pystar /. X L /2]; IVW = Integrate[EID[yapprox, {x, 2}] D[ystar, fx, 2}}, {x, 0, L)]; Eq1 = Coefficient [IVW, a4s) - Coefficient (EVW, a4s); S = Solve(Eq1 - 0, a4]; Print("Approximate displacement:"] yapprox = yapprox /. s[1]] Approximate displacement: 5x 5x 5x Out[125]=-- 64 32 64 Question 2. Solve the same Euler Bernoulli beam using the virtual work method with a fourth order polynomial defined over the whole domain. Compare the fourth order polynomial obtained with the one obtained using the RR method. In[109] = Clear(EI, L, P, yexact, yapprox, a4, a4s ] yapprox = a + a1*x + a2 + x^2 + a3 + x^3 + a4 *x^4; EI = 1; L = 1; P = 1; th =D [yapprox, x); y1 = yapprox /. x 0; y2 = yapprox/. X L ; th1 = th /. x 0; th2 = th /. xL; a = Solve[(y1 = 0, y2 = 0, th1 = 0, th2 - 0), (ao, ai, a2, a3}); yapprox = yapprox /. a[[1]); ystar = yapprox /. (a4 a4s}; EVW = FullSimplify[-Pystar /. X L /2]; IVW = Integrate[EID[yapprox, {x, 2}] D[ystar, fx, 2}}, {x, 0, L)]; Eq1 = Coefficient [IVW, a4s) - Coefficient (EVW, a4s); S = Solve(Eq1 - 0, a4]; Print("Approximate displacement:"] yapprox = yapprox /. s[1]] Approximate displacement: 5x 5x 5x Out[125]=-- 64 32 64