MetodeElement Hinggga

k1=400;
k2=100;
k3=300;
k4=700;
a1=[k1 -k1 0 0 0; -k1 k1 0 0 0; 0 0 0 0 0; 0 0 0 0 0; 0 0 0 0 0]
a2=[0 0 0 0 0; 0 k2 -k2 0 0; 0 -k2 k2 0 0; 0 0 0 0 0; 0 0 0 0 0]
a3=[0 0 0 0 0; 0 0 0 0 0; 0 0 k3 -k3 0; 0 0 -k3 k3 0; 0 0 0 0 0]
a4=[0 0 0 0 0; 0 0 0 0 0; 0 0 0 0 0; 0 0 0 k4 -k4; 0 0 0 -k4 k4]
k=a1+a2+a3+a4
L=[k(2,2) k(2,3) k(2,4) k(2,5); k(3,2) k(3,3) k(3,4) k(3,5); k(4,2) k(4,3) k(4,4) k(4,5); k(5,2) k(5,3) k(5,4) k(5,5)]
t=[0; 0; 0; 10]
d2x_d5x=inv(L)*t
dx=[0;d2x_d5x(1,1);d2x_d5x(2,1);d2x_d5x(3,1);d2x_d5x(4,1)]
Fx=k*dx