out ss;
%off nat;
%load "algint";

k:=mat((2,-1),(-1,2));
%k:=mat((2,-1),(-2,2));


%procedure fip(ii,x); rrr:=sp(ii)*e**(sss*x);
%procedure fim(ii,x); rrr:=sm(ii)*e**x;
%procedure fip(ii,x); rrr:=x**-n;
%procedure fim(ii,x); rrr:=x**-n;
%procedure fip(ii,x); rrr:=n;
%procedure fim(ii,x); rrr:=n;

procedure fip(ii,x); rrr:=(x**7)*e**x;
procedure fim(ii,x); rrr:=(x**7)*e**x;


%procedure integ(f,x,a,b); rrr:=sub(x=b,int(f,x)) - sub(x=a,int(f,x));
procedure integ(f,x,a,b); rrr:=sub(x=b,int(f,x)) ;


for all ii let ap(ii) = integ(fip(ii,t1),t1,c,zp);
for all ii let am(ii) = integ(fim(ii,t1),t1,c,zm);

for all ii,jj let ap(ii,jj) = integ(fip(jj,t2)*integ(fip(ii,t1),t1,c,t2),t2,c,zp);
for all ii,jj let am(ii,jj) = integ(fim(jj,t2)*integ(fim(ii,t1),t1,c,t2),t2,c,zm);


gg1:=1 - ap(1)*am(1) - k(1,2)*ap(1,2)*am(2,1);
gg2:=1 - ap(2)*am(2) - k(2,1)*ap(2,1)*am(1,2);


ro(1):=log( - df(log(gg1),zp,zm) );
ro(2):=log( - df(log(gg2),zp,zm) );


df(ro(1),zp,zm) - k(1,1)*e**ro(1) - k(1,2)*e**ro(2);
df(ro(2),zp,zm) - k(2,1)*e**ro(1) - k(2,2)*e**ro(2);


ap(1,2);
ap(1);
ap(2);


shut ss;
bye;
end;