Inverse Cosine

function ArcCos(x : float) : float;
begin
  Return(ATan2(Sqrt(1 - X * X), X));
end;

const  HalfPI     : float = PI/2.0;


function ATan2( y : float; x : float ) : float;
var arctan2: float;
begin
  if X = 0.0 then
    if y = 0.0 then
      // (* Error! Send Error Message *)
    elsif y > 0.0 then
      arctan2 := halfPi;
    else
      arctan2 := - halfpi;
    endif;
  else
    if X > 0.0 then
      arctan2:= ArcTan( y / x );
    elsif X < 0.0 then
      if y >= 0.0 then
        arctan2:= ArcTan( y / X ) + pi;
      else
        arctan2:= ArcTan( y / X ) - pi;
      endif;
    endif;
  endif;
  Return(arctan2);
end;