Вычисляем значение exp ( -x * x) ряд Тейлора

[kuzlar]

Дана функция exp( - x * x) имеет вот такой ряд
Замечаем что n+1 член ряда вычисляется после n-ого, поэтому программа получиться более эффективной, если находить член ряда не "с нуля", а умножением предыдущего члена на некоторую величину. Найдем эту величину ( Для этого запишем формулу для n+1, и разделим её на формулу для нахождения n) и получим -x / n + 1.
Теперь пишем код

program ex;
const Maxlter = 500;         {max колво интераций}
var x, eps : double;          {аргумент и точность}
   c,y  : double;               {член ряда и его сумма}
   n    : integer;               {номер члена ряда}
   done  : boolean;   
 
begin
writeln('BBeduTe argymeHT u To4HocTb:');
readln(x, eps);
done := true;
c := 1 - (x * x);   {первый член ряда}
y := c;                 {задаем начальное значение суммы}
n := 1;
while abs(c) > eps do begin
c := -x /(n + 1);   {очередной член}
y := y + c;              {плюсуем к сумме}
inc(n);
if n <= Maxlter then continue;
writeln('P9D PaccxoduTbc9!');
done := false; break
end;
if done then
writeln('ArgumeHT: ', x:10:6, #13#10,
'3Ha4EHue func:    ', y:10:6, #13#10,
'eps ', eps:8:6, ' 3A ',   n, 'uHTepacuu',
'EXP(-x*x) = ', exp( - x * x):10:6;                                {*прямое вычисление}
);
readln;
end.

Почему прямое вычисление не совпадает с вычислением через ряд Тейлора? ПОМОГИТЕ пожалуйста!

[Serebro]

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

Дана функция exp( - x * x) имеет вот такой ряд
Замечаем что n+1 член ряда вычисляется после n-ого, поэтому программа получиться более эффективной, если находить член ряда не "с нуля", а умножением предыдущего члена на некоторую величину. Найдем эту величину ( Для этого запишем формулу для n+1, и разделим её на формулу для нахождения n) и получим -x / n + 1.
Теперь пишем код

program ex;
const Maxlter = 500;         {max колво интераций}
var x, eps : double;          {аргумент и точность}
   c,y  : double;               {член ряда и его сумма}
   n    : integer;               {номер члена ряда}
   done  : boolean;   
 
begin
writeln('BBeduTe argymeHT u To4HocTb:');
readln(x, eps);
done := true;
c := 1 - (x * x);   {первый член ряда}
y := c;                 {задаем начальное значение суммы}
n := 1;
while abs(c) > eps do begin
c := -x /(n + 1);   {очередной член}
y := y + c;              {плюсуем к сумме}
inc(n);
if n <= Maxlter then continue;
writeln('P9D PaccxoduTbc9!');
done := false; break
end;
if done then
writeln('ArgumeHT: ', x:10:6, #13#10,
'3Ha4EHue func:    ', y:10:6, #13#10,
'eps ', eps:8:6, ' 3A ',   n, 'uHTepacuu',
'EXP(-x*x) = ', exp( - x * x):10:6;                                {*прямое вычисление}
);
readln;
end.

Почему прямое вычисление не совпадает с вычислением через ряд Тейлора? ПОМОГИТЕ пожалуйста!

Наверное здесь

c := -x /(n + 1);   {очередной член}

ошибка.
Должно быть
c:=c*(-x/(n+1))

[kuzlar]

перерешал, очередной член равен - х / n, но все равно не правильно! =( помогите решить!