Рекуррентная формула

[Damning]

Помогите с заданием. Не получается правильно составить выражение.
Вычислить значение некоторой функции в заданном диапазоне изменения аргумента Х с заданной точностью.

 double factorial(int N)
{
    if(N < 0) 
        return 0; 
    if (N == 0) 
        return 1; 
    else 
        return N *  factorial(N - 1); 
}

void func::formula(float X)
{
	x = X;
	float s = 1;
	k = 1;
	sum = s;
	y = 2 * (powf(cos(x), 2) - 1);
	do 
	{
		k++;
		s *= powf((-1),k) * ( (powf((2*x),(2*k))) / (2*factorial(k)));
		sum += s;
	} while (fabs(s) > eps);
}

[min@y™]

функция factorial() здесь не нужна: в каждой итерации достаточно домножать знаменатель на k.
Знак очередного слагаемого не надо вычислять, возводя –1 в степень:

sign = (k & 1) ? –1 : +1;

[Damning]

Все равно значение функции Y и суммы не сходятся. Что я делаю не так?

void func::formula(float X)
{
	x = X;
	float s = 1;
	int sign;
	k = 1;
	sum = s;
	y = 2 * (powf(cos(x), 2) - 1);
	do 
	{
		k++;
		s *= (powf((2*x),(2*k))) / (2*k);
		sign = (k&1) ? -1 : +1;
		sum += s*sign;
		
	} while (fabs(s) > eps);
}

[min@y™]

я ж сказал

в каждой итерации достаточно домножать знаменатель на k

а ты чо творишь? вынеси знаменатель в отдельную переменную z = 1, затем в каждой итерации перед вычислением члена ряда:

z *= k;

[Damning]

Все равно не получается. Не могли бы вы исправить предыдущий код? Очень прошу.

[Аватар]

Рекуррентная формула

Не вижу вообще. Рекуррентная - это когда кадый последующий член ряда выводится из предыдущего. Здесь:

R(1)=-x*x
R(k)=-R(k-1)*4*x*x/k

и получается алгебраическим делением k-го члена на (k-1)-ый

[Damning]

-------------