out ss;
%off nat;
load "algint";
procedure fi(ii,x); rrr:=d(ii)*e**(c(ii)*x);
procedure integ(f,x,a,b); rrr:=sub(x=b,int(f,x)) - sub(x=a,int(f,x));





fi(1,s);

int(fi(2,t1),t1);

int((x**2/2)*e**x,x);

sub(x=3,x**2);

integ(x,x,1,2);

integ(fi(1,x),x,c,s);

integ(integ(fi(1,x),x,c,s),s,c,zp);

%-----------------------------------------
XP1:=mat((0,1,0),(0,0,0),(0,0,0));
XP2:=mat((0,0,0),(0,0,1),(0,0,0));
XP3:=mat((0,0,1),(0,0,0),(0,0,0));

XM1:=TP(XP1);
XM2:=TP(XP2);
XM3:=TP(XP3);
%-----------------------------------------
ONE:=mat((1,0,0),(0,1,0),(0,0,1));
DEPEND AP1,ZP;
DEPEND AP2,ZP;
DEPEND AP12,ZP;
DEPEND AM1,ZM;
DEPEND AM2,ZM;
DEPEND AM12,ZM;

MP:=ONE + AP1*XP1 + AP2*XP2 +AP12*XP3;
MM:=ONE + AM1*XM1 + AM2*XM2 +AM12*XM3;
MP**-1;
KK:=(MP**-1)*MM;


v1:=mat((1),(0),(0));
v2:=mat((0),(1),(0));
GG1:=TP(V1)*KK*V1;
GG2:=TP(V2)*KK*V2;

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


GG1:=TP(V1)*KK*V1;
GG2:=TP(V2)*KK*V2;


RO(1):=LOG(-DF(LOG(GG1(1,1)),ZP,ZM));
RO(2):=LOG(-DF(LOG(GG2(1,1)),ZP,ZM));

%DF(RO(1),ZP,ZM) -2*EXP(RO(1)) + EXP(RO(2));


shut ss;
bye;
end;