function[Ke]=Kelement(A,E,L,teta) l=cos(teta); m=sin(teta); Ke=A*E/L*[ l^2 l*m -l^2 -l*m l*m m^2 -l*m -m^2 -l^2 -l*m l^2 l*m -l*m -m^2 l*m m^2];