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

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