A3 = [-1      1      0      0      0;...
      1  -2002      1   1000   1000;...
      0      1     -1      0      0;...
      0   1000      0  -1000      0;...
      0   1000      0      0  -1000]
 
detA3 = det(A3)
 
A2 = [-2 1 0 1;...
1 -2 1 0;...
0 1 -2 1;...
1 0 1 -2]
 
detA2 = det(A2)
 
A = zeros(7, 7);
A(1+0,1+0)= -1;
A(1+1,1+0)= 1;
A(1+0,1+1)= 1;
A(1+1,1+1)= -2002;
A(1+2,1+1)= 1;
A(1+3,1+1)= 1000;
A(1+4,1+1)= 1000;
A(1+1,1+2)= 1;
A(1+2,1+2)= -1;
A(1+3,1+3)= -1000;
A(1+1,1+3)= 1000;
A(1+5,1+3)= -1;
A(1+4,1+4)= -1000;
A(1+1,1+4)= 1000;
A(1+6,1+4)= -1;
A(1+3,1+5)= 1;
A(1+4,1+6)= 1;
 
x1 = [0,0,0,0,0,1,0]';
x2 = [0,0,0,0,0,0,1]';
 
A^(-1);
 
A^(-1) * x1
A^(-1) * x2;
 
disp(A);

QTMgPSBbLTEgICAgICAxICAgICAgMCAgICAgIDAgICAgICAwOy4uLgogICAgICAxICAtMjAwMiAgICAgIDEgICAxMDAwICAgMTAwMDsuLi4KICAgICAgMCAgICAgIDEgICAgIC0xICAgICAgMCAgICAgIDA7Li4uCiAgICAgIDAgICAxMDAwICAgICAgMCAgLTEwMDAgICAgICAwOy4uLgogICAgICAwICAgMTAwMCAgICAgIDAgICAgICAwICAtMTAwMF0KICAgICAgCmRldEEzID0gZGV0KEEzKQoKQTIgPSBbLTIgMSAwIDE7Li4uCjEgLTIgMSAwOy4uLgowIDEgLTIgMTsuLi4KMSAwIDEgLTJdCgpkZXRBMiA9IGRldChBMikKCkEgPSB6ZXJvcyg3LCA3KTsKQSgxKzAsMSswKT0gLTE7CkEoMSsxLDErMCk9IDE7CkEoMSswLDErMSk9IDE7CkEoMSsxLDErMSk9IC0yMDAyOwpBKDErMiwxKzEpPSAxOwpBKDErMywxKzEpPSAxMDAwOwpBKDErNCwxKzEpPSAxMDAwOwpBKDErMSwxKzIpPSAxOwpBKDErMiwxKzIpPSAtMTsKQSgxKzMsMSszKT0gLTEwMDA7CkEoMSsxLDErMyk9IDEwMDA7CkEoMSs1LDErMyk9IC0xOwpBKDErNCwxKzQpPSAtMTAwMDsKQSgxKzEsMSs0KT0gMTAwMDsKQSgxKzYsMSs0KT0gLTE7CkEoMSszLDErNSk9IDE7CkEoMSs0LDErNik9IDE7Cgp4MSA9IFswLDAsMCwwLDAsMSwwXSc7CngyID0gWzAsMCwwLDAsMCwwLDFdJzsKCkFeKC0xKTsKCkFeKC0xKSAqIHgxCkFeKC0xKSAqIHgyOwoKZGlzcChBKTs=