Maple Programming - Attempting to program the calculation of the B matrix and sigma (Stresses) after finding the forces and displacements of two bars in Maple. Here is what I have so far, but the code is not working. Topic - Finite Element Analysis

KU:=2.333333333108u39.999999999108u35.333333333108u3eqs:={80000.0=9.999999999108u3,Fx1=2.333333333108u3,Fx2=2.333333333108u3,Fx4=5.333333333,vars:={Fx1,Fx2,Fx4,u3}sols:={Fx1=18666.66667,Fx2=18666.66667,Fx4=42666.66667,u3=0.00008000000001} \#Post-processing \#efine Stresses for i from 1 to nel do n_point :=Vector[ column ]([connect[i,1], connect [i,2]]); flocal [i]:=k_local[i][[[U[n_point[1]]],[U[n_point [2]]]]]; end do; \#Version 2 B1:=Matrix([[A11,0],[0,A11]])E12Matrix([[1,L121],[L121,1]])B2:=Matrix([[A21,0],[0,A11]])E12Matrix([[1,L121],[L121,1]])B3:=Matrix([[A31,0],[0,A31]])E3Matrix([[1,L31],[L31,1]];sigmal:=U:=LinearAlgebra[LUDecomposition](K)solve(F);B1:=Matrix([[A11,0],[0,A11]])E12Matrix([[1,L121],[L121,1]])B2:=Matrix([[A21,0],[0,A21]])E12Matrix([[1,L121],[L121,1]])B3:=Matrix([[A31,0],[0,A31]])E3Matrix([[1,L31],[L31,1]]);sigmal:=B1(Matrix([U[1],U[3]])Matrix([U[3],U[4]]));sigma2:=B2(Matrix([U[2],U[3]])Matrix([U[3],U[4]]));sigma3:=B3(Matrix([U[3],U[4]])Matrix([U[1],U[2]])); Figure 1: Problem 1