out bb;
%off nat;
on rounded;

procedure delta(ii,jj);
begin 
if ii eq jj then return 1 else return 0; 
end;

procedure addone(x); rrr:=x+1;

n:=2;
m:=2;


matrix KM(n,n);
for tt:=1:n do for pp:=1:m do km(tt,pp):=k(tt,pp);

AA:=tp( mat(( a(1),a(2) )) );
matrix AM(n,n);
for ii:=1:n do am(ii,ii):=a(ii);


matrix RM(n,n);
for tt:=1:n do for pp:=1:m do Rm(tt,pp):=r(tt,pp);

%-----------------------------------------------------
matrix fm(n,n);
%for tt:=1:n do for pp:=1:m do fm(tt,pp):=f(tt,pp);
fm(1,1):=1;fm(2,2):=2;
fm;

matrix gm(n,n);
%for tt:=1:n do for pp:=1:m do gm(tt,pp):=g(tt,pp);
gm:=km*(fm**-1)*(km**-1);

gm*km*fm-km;
gm*rm*fm-rm;

RM;

matrix NOM(n,n);
for tt:=1:n do for pp:=1:m do nom(tt,pp):= NO(tt,pp);

NOM;
matrix RNM(n,n);


%RM:=KM*AM;
RNM:=RM+NOM;

rnm(1,1):=d(1) + 10; rnm(1,2):=d(1) - 10;
rnm(2,1):=d(2) + 5; rnm(2,2):=d(2) - 5;

k(1,1):=4; k(1,2):=2;
k(2,1):=2; k(2,2):=4;
%a(1):=5; a(2):=10;
d(1):=140;
d(2):=60;

Km*am;
gm*rnm*fm;

eqq:=gm*rnm*fm - rnm;
eqqn:=gm*nom*fm - nom;

eqqx:=eqq - eqqn;

eqq(1,1) + eqq(1,2);


eqqlist:={eqqx(1,1), eqqx(1,2), eqqx(2,1), eqqx(2,2)};
arglist:={no(1,1), no(1,2), no(2,1), no(2,2)};



solve(eqqlist, arglist);


shut bb;
bye;
end;