Нужно посчитать сумму ряда. В программе ошибка. Пожалуйста помогите найти

[Рита Белова]

Даны действительные числа x и eps ( х не равен 0, eps>0). Вычислить с точностью до eps сумму ряда
(-1)^k * x^(4k+3)/ ((2k+1)!*(4k+3))
Программу написала, но ответ, просчитаный мной не совпадает с ответом написаной програмы. Скажите что не правильно:
procedure TForm1.BitBtn2Click(Sender: TObject);
Var a,u,S,eps,x:Real;
i:Integer;
begin
x:=StrToFloat(Edit1.Text);
eps:=StrToFloat(Edit2.Text);
S:=0;
i:=0;
a:=x*x*x;
u:=x*x*x/3;
While abs(u)>=eps do
begin
S:=S+u;
i:=i+1;
a:=a*(-x*x*x*x/2*i/(2*i+1));
u:=a/(4*i+3);
end;
Edit3.Text:=FloatToStr(S);
end;

[Greblin]

Ошибка вот тут:

a := a * (-x*x*x*x/2*i/(2*i+1));

2 * i забыли в скобки заключить

[Voody]

Рита, вам не кажется, что вы не правильно считаете члены ряда.
Ваша формула:

a:=a*(-x*x*x*x/2*i/(2*i+1));
u:=a/(4*i+3);

идеально подходит для элемента i=1, но уже дальше считает неправильно.
Вам лучше составить функцию вычисления факториала.

И второе (говорю как учившийся на матмехе, может вам по-другому сказали): точность вычисляется разностью, т.е. (S от i+1) минус (S от i). Если разность меньше eps, значит указанная точность достигнута.

[Greblin]

Voody, да вроде всё правильно
a0 = x^3; u0 = x^3 / 3
a1 = a0 * (-x^4) / ((2*1) * (2*1+1)) = - x^7 / 3!; u1 = -x^7 / (3! * 7)
a2 = a1 * (-x^4) / ((2*2) * (2*2+1)) = x ^ 11 / 5!; u2 = x^11 / (5! * 11)
и т.д.
А так как S := S + u; то S(i+1) - S(i) = u, что и стоит в условии

[Voody]

о, тогда прошу прощения.))

[Рита Белова]

Greblin, спасибо. Но все равно что-то не то получается, хотя ответ уже ближе к истине))