function k=Evaluate_Local_Stiffness(anElement)

E=anElement.aMaterial.Eyoung;
A=anElement.aProperty.A;
L=anElement.Length;

k=E*A/L*[1 0 -1 0;
                 0 0 0 0;
               -1 0  1 0;
                 0 0  0 0];

end