clear
disp('Solution of the bridge problem')
W = -20000
theta = pi/6
pause
c = cos(theta);
s = sin(theta);
A = [-1 0 -c 0 0 0 -1  0;
    0 -1 -s 0 0 0  0 0;
    0 0 -c -c 0 0  0 0;
    0 0 -s s 0 0 0 -1;
    0 0 0 -c 0 -1  0 0;
    0 0 0 s -1 0 0 0;
    0 0 0 0 0 1 -1 0;
    0 0 0 0 0 0 0 1]
pause
d = zeros(8,1);
d(8) = W
pause
disp('Use A\d')
sol1 = A\d
pause
disp('Use inv(A)*d')
disp('inv(A) =')
inv(A)
pause
disp('A*inv(a) =')
A*ans
pause
sol2 = inv(A)*d