Вычисление синуса как суммы ряда

Программа, вычисляющая синус как сумму ряда вида x-x^3/3+x^5/5-... с заданной точностью и сравнивающая полученный результат с результатом вычисления синуса стандартной библиотечной функцией С++.

// синус как сумма ряда
double CompareSin(double Eps, double X)
{
   double Fx = 0, Fxold = 0, PoxX = X;
   long i = 1;
   int minus = 1;

   do 
    {
       Fxold = Fx;
       Fx += minus*PoxX/i;
       PoxX *= X*X;
       i+=2;
       minus *= -1;
    } 
   while((Fx-Fxold)>=Eps);

   return (Fx - sin(X));
}