Метод парабол (Симпсона)

From

Eduard Vatutin (2:5035/43.25)

To

ALl

Date

2002-07-26T01:37:59Z

Area

RU.ALGORITHMS

 _Приветик_ ALl!

При реализации сабжевого алгоритма возникли определенные трудности. При четном числе итераций погрешность составляет 1е-14, при нечетном - 1е-1. Так и должно быть или я где-то ошибся?

S := 0;
for I := 0 to N do
  S := S + Func(A+I*h)*q(I);

N - число итераций.
[A, B] - отрезок интегрирования.
h - шаг.
q(I) вычисляется следующим образом:

if (I=0) or (I=N) then
  Result := 1
else if Odd(I) then
  Result := 4
else
  Result := 2;

--- _/Пока, ALl/_
 * Origin: Хороший хирург поможет плохому танцору! (2:5035/43.25)