out aa2;

%depend j1, x, y, p;
%depend j2, x, y, p;

%j2:=atan(y/x) + p;

procedure jj1(x,y,p); rrr:=x*cos(p) - y*sin(p);
procedure jj2(x,y,p); rrr:=x*sin(p) + y*cos(p);
procedure ss(x,y,p); rrr:=s;

procedure xx(j1,j2,s); rrr:=j1*cos(s) + j2*sin(s);
procedure yy(j1,j2,s); rrr:=-j1*sin(s) + j2*cos(s);
procedure pp(j1,j2,s); rrr:=p;



xx(jj1(x,y,p), jj2(x,y,p), ss(x,y,p));
yy(jj1(x,y,p), jj2(x,y,p), ss(x,y,p));
pp(jj1(x,y,p), jj2(x,y,p), ss(x,y,p));


%j1,j2,s);
%yy(j1,j2,s);
%pp(j1,j2,s);


bye;
end;



%df(j2,p);
%df(j2,x);
%df(j2,y);

df(j1,p)+ y*df(j1,x) -x*df(j1,y);
df(j2,p)+ y*df(j2,x) -x*df(j2,y);
%---------------------------------------

let x=jj1*cos(s) + jj2*sin(s);
let y=-jj1*sin(s) + jj2*cos(s);


j1;
j2;




shut aa2;
bye;
end;