... function EF=forcesInElementGrid(numberElements, elementNodes,xx,yy,E,I,G,J,displacements)% forces in elementsEF=zeros(6,numberElements); for e=1:numberElements;% elementDof: element degrees ... stiffness=formStiffnessGrid(GDof, numberElements,elementNodes,xx,yy,E,I,G,J)% function to form global stiffness for grid element stiffness=zeros(GDof); for e=1:numberElements;% elementDof: element ... transformation matrixa=[1 0 0; 0 C S;0 -S C];R=[a zeros(3);zeros(3) a];stiffness(elementDof,elementDof)= stiffness(elementDof,elementDof)+R’*k*R;endThe code for computing the forces in elements...