out aa4;
k(1):=0.3$
k(2):=0.6$
k(3):=0.1$

f(1):=1;
f(2):=2;
f(3):=3;

mk:=mat((k(2), k(1), 0),
       (k(3), k(2), k(1)),
       (0   , k(3), k(2))   );

mk**-1;

ma:=mat((a(1), 0   , 0),
       (0   , a(2), 0),
       (0   , 0   , a(3))   );



%mn:=mat((n(1), 0   , 0),
%       (0   , n(2), 0),
%       (0   , 0   , n(3))   );

matrix mn(3,3);
for ii:=1:3 do for jj:=1:3 do mn(ii,jj):=n(ii,jj);
mn;


matrix mf(3,3);
%for ii:=1:3 do for jj:=1:3 do mf(ii,jj):=f(ii,jj);
for ii:=1:3 do mf(ii,ii):=f(ii);

mf;


ma*mf;

ma*mf - mf*ma;


mg:=mk*(mf**-1)*(mk**-1);


mg*mk*mf-mk;

%---------------------------------
%mr:=mat((r(1), 0   , 0),(0   , r(2), 0), (0   , 0   , r(3))   );


mr:=mk*ma;

mr-mg*mk*ma*mf;

mr*mf-mf*mr;
%---------------------------------

matrix mrn(3,3);
for ii:=1:3 do mrn(ii,ii):=rn(ii);



mrn:=mr + mn;

mg*mrn*mf - mrn;

mg*mn*mf - mn;


%eqqn1-eqqn;


shut aa4;
bye;
end;