A = [  1    -3    -4    -3
      -1     6     4     1
      -2     6    10     7
      -3     3    14    15];
b = [ -5;     16;      9;     -9];

