Класс комплексных чисел в показательной форме

//---------------------------------------------------//
//-  Класс комплексого числа в символической форме  -//
//---------------------------------------------------//
//-  double A - длина вектора                       -//
//-  double F - угол поворота                       -//
//---------------------------------------------------//
class COMPLEX
{
public:
    double A;
    double F;
    COMPLEX()
    {
        A = 0.0;
        F = 0.0;
    }
    COMPLEX(double fst,double snd)
    {
        A = fst;
        F = snd;
    }
    COMPLEX &operator=(const COMPLEX orig)
    {
        A = orig.A;
        F = orig.F;
    }
    COMPLEX operator*(const COMPLEX &second)const
    {
        COMPLEX result;
        result.A = abs(A * second.A);
        if(result.A == 0) result.F = 0;
        else result.F = F + second.F;
        return result;
    }
    COMPLEX operator/(const COMPLEX &second)const
    {
        COMPLEX result;
        result.A = abs(A / second.A);
        if(result.A == 0) result.F = 0;
        else result.F = F - second.F;
        return result;
    }
    COMPLEX operator+(const COMPLEX &second)const
    {
        COMPLEX result;
        double a = abs(A), b = abs(second.A);
        result.A = sqrt(a*a + b*b + 2*a*b*cos(F-second.F));
        
        if(result.A == 0) result.F = 0; //если вектор равен нулю, обнуляем угол
        else if (F==0 && second.F==0) result.F = 0; //если %операция% два нулевых угла получаем 0
        else
        {
            double x = a*sin(F)+b*sin(second.F);
            double y = a*cos(F)+b*cos(second.F);
            //округление числа для избавления от погрешности вычислений
            x = (double)((int)(x*10000))/10000;
            y = (double)((int)(y*10000))/10000;
            result.F = ArcTan(y,x);
        }  
        return result;
    }
    COMPLEX operator-(const COMPLEX &second)const
    {
        COMPLEX result;
        double a = abs(A), b = abs(second.A);
        result.A = sqrt(a*a + b*b - 2*a*b*cos(F-second.F));
        
        if(result.A == 0) result.F = 0; //если вектор равен нулю, обнуляем угол
        else if (F==0 && second.F==0) result.F = 0; //если %операция% два нулевых угла получаем 0
        else
        {
            double x = a*sin(F)-b*sin(second.F);
            double y = a*cos(F)-b*cos(second.F);
            //округление числа для избавления от погрешности вычислений
            x = (double)((int)(x*10000))/10000;
            y = (double)((int)(y*10000))/10000;
            result.F = ArcTan(y,x);
        }        
        return result;
    }
};

double ArcTan(double y, double x)
{
    if (x==0) return M_PI/2;
    else if (y==0) return 0;
    else
    {        
        //определение четверти
        if (x>0 && y>0) return atan(y/x); //1
        else if (x<0 && y>0) return atan(y/x)+M_PI; //2
        else if (x<0 && y<0) return atan(y/x)-M_PI; //3
        else if (x>0 && y<0) return atan(y/x); //4
    }
}