if x > 0 then atan2 := arctan (y/x) else if x > 0 then atan2 := arctan (y/x) + pi else atan2 := pi/2 * sgn (y);end;
{ Я сделал это некоторое время назад при портировании кода из FORTRAN на Pascal. К сожалению, это было так давно, что я не помню, тестировал я это, или нет. }
