
ep = 1.e-10;

A = [ 1   1   1  1;
      ep  0   0  1;
      0   ep  0  1;
      0   0   ep 1] ;


