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

[Solidera]

Пожалуйста, помогите с заданием:
Разработать алгоритм и составить программу вычисления значений суммы S заданного ряда при различных значениях аргумента Х: во внутреннем цикле для фиксированного значения Х поочередно суммировать члены ряда, пока их абсолютная величина превышает заданную точность eps. Вычисление очередного член ряда осуществлять по рекуррентной формуле. Во внешнем цикле аргумент Х изменяется от 0,5 до 0,75 с шагом 0,05.

Где напортачила?

float Summ(float x, float eps)
{
	float elem,summ;
	int n=1;
	summ=elem=(x*(3+x))/6;
	while (x*(2*n+1+x)*(2*n+2)*(2*n+3)/(2*n+x)<eps)
	{
		elem*=-(x*(2*n+1+x)*(2*n+2)*(2*n+3))/(2*n+x);
		summ+=elem;
		n++;
	}
	return summ;
}

[Stilet]

Где напортачила?

Не пойму где считается факториал? Да и возведение в степень не вижу.

[ViktorR]

1. Неверна рекурентная формула. У меня получилось так:
an = x*(3+x)/(2*3);
k = -pow(x,2)*(2*n+3+x)/((2*n+1+x)*(2*n+2)*(2*n+3));
2. По условию надо вычисления закончить тогда, когда элемент станет меньше eps.
Таким образом получаем:

float Summ(float x, float eps)
{
   float elem,summ, k;
   int n = 1;
   summ = elem = x*(3+x)/(2*3);
   while (fabsf(elem) < eps)
   {
     k = -pow(x,2) * (2*n+3+x) / ((2 * n+1+x) * (2 * n+2)  * (2 * n + 3));
     elem* = k;
     summ+ = elem;
     n++;
/* Полезно, в этом месте вставить ограничение по n, например,
   if n > 100
    break;
*/
   }
return summ;
}

fabsf(elem) - получить абсолютное значение числа (модуль). Эту функцию нашел на одном из форумов, но как это на С++ уверенно не скажу.


Как-то так, ...

[Rin]

У вас же в знаменателе (2n+1)! - факториал то есть, а Вы в формуле ставите (2n+x). Ещё не могу найти X в степени 2n-1. Так же не могу найти у вас пресловутую -1 в степени n-1. Может вы как-то преобразовали формулу, я в числовых рядах не очень. Возможно, косяк кроется где-то тут.

[Solidera]

Спасибо всем, сейчас исправлю и посмотрю. У меня обычно косяки в мелочах, стандарт - перепутать плюс с минусом и т. д.

[Solidera]

Сообщение от ViktorR

fabsf(elem)

На плюсах просто abs.

Огромное Вам спасибо!
А ограничитель по n там не просто полезен, а жизненно необходим.