A1=[3 0;   4 5;   0 4];            b1=[5; 6; 7];
A2=[3 0 1; 4 5 1; 0 4 3; 2 2 3];   b2=[5; 6; 7; 8];
%A3=[1 2 3 4; 4 6 8 7; 3 7 5 4; 8 15 16 15]; b3=[10; 25; 19; 54];
