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;

r:=tp( mat(( d(1),d(2) )) );


matrix SM(n,n);
for tt:=1:n do for pp:=1:n do sm(tt,pp):=s(tt,pp);



qq:=tp( mat(( q(1),q(2) )) );

matrix QM(n,n);
for tt:=1:n do for pp:=1:n do qm(tt,pp):=q(tt,pp);

SM*R - R;

eqq:=SM*(R+Qq) - R;

solve({eqq(1,1), eqq(2,1)}, {s(1,2),s(2,1)});



shut bb;
bye;
end;