Вычисление определенного интеграла методом трапеций

[Alexander886]

Здравствуйте, помогите найти ошибку. Необходимо вычислить определенный интегралИнтеграл.pngметодом трапеций с точностью 0,0001. В коде должно быть два цикла repeat. Программа считает неверно.

procedure TForm1.Button2Click(Sender: TObject);
var a,b,Eps,S,h,x,S0:real;
begin
 a:=StrToFloat(Edit1.Text);
 b:=StrToFloat(Edit2.Text);
 Eps:=StrToFloat(Edit3.Text);
 h:=b-a;
    repeat
      begin
        S:=(f(a)+f(b))/2;
        x:=a;
          repeat
            x:=x+h;
            S:=S+f(x);
          until x>b;
        S0:=S;
        h:=h/2;
      end;
    until abs(S-S0)<Eps;
 Label6.Caption:=FloatToStrf(S,fffixed,6,4);
end;
end.