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

Программа, вычисляющая синус как сумму ряда вида 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));
}