atan2

[DrD]

Here is a little atan2(y,x) program especially useful as a subroutine:

Code:


EXPORT atan2(y,x)
BEGIN

// atan2 returns coordinates in correct angle for all four quadrants  
  CASE
    if x==0 AND y==0 then return "undefined"; end;         // x=0, y=0  
    if x>0 then return atan(y/x); end;                     // x>0
    if x<0 AND y>=0 then return atan(y/x)+pi; end;         // x<0, y>=0
    if x==0 AND y>0 then return pi/2; end;                 // x=0, y>0
    if x<0 AND y<=0 then return atan(y/x)-pi; end;         // x<0, y<0
    if x==0 AND y<0 then return -pi/2; end;                // x=0, y<0 
  END;
 
  return;
 
END;