Численные методы, Квадратурные формулы Опции

[nyusha]

Здравствуйте, помогите, пожалуйста, доделать программу, в общем случае я написала, а вот обыграть точность в программе не получается. Пытаюсь использовать метод двойного пересчета, вычисляю значение интегралов при шаге h, затем уменьшаю шаг вдвое и вычисляю снова, далее необходимо оценить разность и продолжать уменьшение шага.... Help
Задание.
Найти приближенное значение определенного интеграла по формулам прямоугольников, трапеций и Симпсона с точностью е=0.001
Попытки решить задание
s1-прямоугольников, s2-трапеций, s3-симпсона.

Program lab3;
var
i,n: integer;
a,b,h,s1,s2,s3,x,x1,x2: real;
function f(x: real): real;
   Begin
       f:=(14*exp(x))/((2+x*x)*(2+x*x));
   end;
Begin
a:=0; b:=2;
   n:=3;
   h:=(b-a)/n;
   x1:=a+h/2;
   x2:=a;
      s1:=f(x1);
      s2:=(f(x2)+f(b))/2;
   for i:=1 to n-1 do
            Begin
                 x1:=x1+h; x2:=x2+h;
                 s1:=s1+f(x1); s2:=s2+f(x2);
            end;
   s1:=h*s1; s2:=h*s2; s3:=(2*s1+s2)/3;
   writeln ('s1 =', s1:6:3, '   s2 =' , s2:6:3, '   s3 =', s3:6:3);
end.