out tt;
off nat;

depend ap,zp;
depend am,zm;

h:=mat((1,0),(0,-1));
xp:=mat((0,1),(0,0));
xm:=mat((0,0),(1,0));


Mp:=mat((1,ap),(0,1));
Mm:=mat((1,0),(am,1));


MM:=(MP**-1)*MM;
ro:=log( - df(log(MM(1,1)),zp,zm) );
df(ro,zp,zm) - 2*exp(ro);

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

RO := LOG((DF(AP,ZP)*DF(AM,ZM))/(AP**2*AM**2 - 2*AP*AM + 1));
df(ro,zp,zm) - 2*exp(ro);

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

am:=(-zp - zm +1)/ap;
ro:=log( - df(log(MM(1,1)),zp,zm) );
df(ro,zp,zm) - 2*exp(ro);

%-------
df(ro,zp,zp) - 2*exp(ro);
df(ro,zp);

%-----------------------------------------
ap:=e**zp;
am:=e**zm;

ap:=c*e**zp;
am:=e**zm;

ro:=log( - df(log(MM(1,1)),zp,zm) );
df(ro,zp,zm) - 2*exp(ro);
df(ro,zp,zp) - 2*exp(ro);

df(ro,zp);


%-----------------1D------------------------

%zm:=zp;
%RO := LOG((DF(AP,ZP)*DF(AM,ZM))/(AP**2*AM**2 - 2*AP*AM + 1));
%df(ro,zp,zp) - 2*exp(ro);

out t5;
write("out tt5");
write("eqq1:=","ro-",ro);
write("eqq2:=","pro-",df(ro,zp));
%df(ro,zp,zp) - 2*exp(ro);

write("shut tt5;bye;end;");

shut t5;


shut tt;
bye;
end;