
A=[0.001,1,1; 1,-1,0];
M1=[1000,0;0,1]
M2=[1,0;-1,1];
M3=[1,0;0,-1/1001];
M4=[1,-1000;0,1];

A
M1*A
M2*M1*A
M3*M2*M1*A
M4*M3*M2*M1*A

fprintf(1, '*************\n')

A=[0.001,1,1; 1,-1,0];
M1=[0,1;1,0]
M2=[1,0;-1/1000,1];
M3=[1,1/1.001;0,1];
M4=[1,0;0,1/1.001];

A
M1*A
M2*M1*A
M3*M2*M1*A
M4*M3*M2*M1*A

fprintf(1, '*************\n')

A=[0.001,1,1; 1,-1,0];
M1=[0,1;1,0]
M2=[1,0;-1/1000,1];
M3=[1,1/1.001;0,1];

A
M1*A
M2*M1*A
M3*M2*M1*A
