out cc;
off nat;

%for all xxx let cos(xxx)**2 + sin(xxx)**2 = 1;
%for all xxx let sinh(xxx) = sqrt(cosh(xxx)**2 - 1);
%for all xxx let sinh(xxx)**2 = cosh(xxx)**2 - 1;
for all xxx let sinh(xxx)**2 + 1 = cosh(xxx)**2;
for all xxx let asinh(sinh(xxx)) = xxx;
for all xxx let sinh(asinh(xxx)) = xxx;


depend kaptt,p,s;
depend c1,p,s;
depend c2,p,s;

%c1*att(kaptt);
c2:=sqrt(kaptt);


df(kaptt,p) - c1*df(kaptt,s) - c2*kaptt**2 - df(c2,s,2);

df(c2,s,2);



shut cc;
bye;
end;


