
clear all
clc
F=3000;
B=[0;F;0;F;0;0;F;0;0;0;0;0;0;0];
l=2;
h=3;

sinPhi1=h/sqrt((l/2)^2+h^2);
cosPhi1=0.5*l/sqrt((l/2)^2+h^2);
sinPhi2=h/sqrt(l^2+h^2);
cosPhi2=l/sqrt(l^2+h^2);


A=[cosPhi1 1 0 0 0 0 0 0 0 0 0 0 0 0;... %S1
sinPhi1 0 0 0 0 0 0 0 0 0 0 0 0 0;... %S1
-cosPhi1 0 cosPhi1 1 0 0 0 0 0 0 0 0 0 0;... %S2
-sinPhi1 0 -sinPhi1 0 0 0 0 0 0 0 0 0 0 0;... %S2
0 -1 -cosPhi1 0 cosPhi1 0 1 0 0 0 0 0 0 0;... %S3
0 0 sinPhi1 0 sinPhi1 0 0 -1 0 0 0 0 0 0;... %S3
0 0 0 -1 -cosPhi1 cosPhi1 0 0 0 0 0 0 0 0;... %S4
0 0 0 0 -sinPhi1 -sinPhi1 0 0 0 0 0 0 0 0;... %S4
0 0 0 0 0 -cosPhi1 -1 0 cosPhi2 0 0 0 0 0;... %S5
0 0 0 0 0 sinPhi1 0 0 -sinPhi2 -1 0 0 0 0;... %S5
0 0 0 0 0 0 0 0 cosPhi2 0 1 0 0 0;... %S6
0 0 0 0 0 0 0 1 sinPhi2 0 0 0 0 1;... %S6
0 0 0 0 0 0 0 0 0 0 -1 -1 0 0;...% S7
0 0 0 0 0 0 0 0 0 1 0 0 1 0];% S7
X=linsolve(A,B);