off nat;


%on rounded;


matrix RM(2,2);


for tt:=1:2 do Rm(tt,tt):=r(tt);




matrix NM(2,2);


for tt:=1:2 do for pp:=1:2 do nm(tt,pp):=n(tt,pp);




matrix RnM(2,2);


for tt:=1:2 do for pp:=1:2 do Rnm(tt,pp):=rn(tt,pp);




matrix sM(2,2);


for tt:=1:2 do sm(tt,tt):=s(tt);



S*rm*s**-1;


MAT((R(1),0),(0,R(2)))$




%-----------------

shut aa7;