... e=1:numberElements;% elementDof: element degrees of freedom (Dof)indice=elementNodes(e,:) ;elementDof=[ (indice(1 )-1 )*3+1 (indice(1 )-1 )*3+2 (indice(1 )-1 )*3+3 (indice(2 )-1 )*3+1 (indice(2 )-1 )*3+2 ... elementNodes,numberNodes,nodeCoordinates,E,A,Iz,Iy,G,J);stiffness=zeros(GDof);% computation of the system stiffness matrixfor e=1:numberElements;% elementDof: element degrees of freedom (Dof)indice=elementNodes(e,:) ;elementDof=[6*indice(1 )-5 6*indice(1 )-4 ... forces in elementsEF=zeros(6,numberElements);for e=1:numberElements;% elementDof: element degrees of freedom (Dof)indice=elementNodes(e,:) ;elementDof= [(indice(1 )-1 )*3+1 (indice(1 )-1 )*3+2...