e1 = [1,0;0,-1] e2 = [0,1;0,0] e3 = [0,0;1,0] m=[a,b;c,(1+bc)/a] m*e1*m^(-1)