out aee;
%off nat;
on rounded;

procedure as(s,a,r); rrr:=(1/k)*sqrt(2*s + r**2) + a - r/k;
procedure rs(s,a,r); rrr:=sqrt(2*s + r**2);


as(0,a,r);
rs(0,a,r);

a1:=as(s1,a,r);
r1:=rs(s1,a,r);

a2:=as(-s1,a1,r1);
r2:=rs(-s1,a1,r1);


ax:=as(-s1,as(s1,a,r),rs(s1,a,r));
rx:=rs(-s1,as(s1,a,r),rs(s1,a,r));


ax:=as(-s1,as(s1,a,r+n),rs(s1,a,r+n));
rx:=rs(-s1,as(s1,a,r+n),rs(s1,a,r+n));


shut aee;
bye;
end;