#include #include /* atan2 discovers what quadrant the angle is in and calls atan. */ double math·atan2(double arg1, double arg2) { if(arg1+arg2 == arg1) { if(arg1 >= 0) return math·PIO2; return -math·PIO2; } arg1 = math·atan(arg1/arg2); if(arg2 < 0) { if(arg1 <= 0) return arg1 + math·PI; return arg1 - math·PI; } return arg1; }