femi qua Thu Apr1816:46:022002 d2u/dx2: F [kn2,kn3)=K(kn2,kn3}-didx·(-4/3}; t quadratic elenents Kixn3, kn3) w K(kn3, kn31-dxi dx'[ 7/6 a use caussian quadrature to evaluate forcing term integ 1/61·xx,“3-(1/12).xx,4-1/12}·xx Gals point b set number of Gauss point.a for forcing cern NGEa N3 Inxs.10, G rand at Causs point for weight functions at nodes x:2: riaht-hand sid N=3 f-age(rn x=linspace 10. 1.2.Ne+:> t Zero stiffness a Set dirichlet conditions at xa Loop over all el △ nd calculate stiffness and residua】s g Set Dirichlet conditions ac x=l tor 11=1: Ne kn2·1+2 n1z(1,2 Anl, kn1 ution x3-1 A plot it and compare K(knl, knl)=K4knul, knl I -dxicx(7/5 Plot =20; t Nunber ot points per element to plot KIkn2, knl= K(kn2, knll-dxidx(-4/2 kn1=1+2*[4=1)