Нахождение интеграла

[Baburek]

Использую метод трапеции. циклится. Помогите пожалуйста разобраться. Прикладываю проект.

integral.rar

unit zadach2;

interface

uses
  Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
  Dialogs, StdCtrls;

type
  TForm1 = class(TForm)
    Button3: TButton;
    edt1: TEdit;
    edt2: TEdit;
    edt3: TEdit;
    Label1: TLabel;
    lbl1: TLabel;
    ComboBox: TComboBox;
    mmo1: TMemo;
    procedure Button3Click(Sender: TObject);
  private
    { Private declarations }
  public
    { Public declarations }
  end;

var
  Form1: TForm1;

implementation

{$R *.dfm}

procedure TForm1.Button3Click(Sender: TObject);
var
x,k,e,h,a,b,sum,q,Newton,nulx,nulk:real;
N,dvaN,i:integer;

Function F(x:real):real;
begin
F:=sqrt(2*x*x+3);
end;

Function P(k:real):real;
begin
P:=sqrt(2*k*k+3);
end;


begin
a:=StrToFloat(edt1.Text); // пределы
b:=StrToFloat(edt2.Text);
e:=StrToFloat(edt3.Text); // искомая точность
N:=2;
dvaN:=N*2;
Repeat
Inc(N);
sum:=0;
Newton:=0;
h:=(b-a)/N;
For i:=1 to N-1 do begin // вычисление интеграла по N разбиениям
x:=a+h*i;
nulx:=a+h*(i-1);
end;
h:=(b-a)/(2*N);
For i:=1 to dvaN do begin // вычисление интеграла по 2N разбиениям
nulk:=a+h*(i-1);
k:=a+h*i;
end;
sum:=sum+(h*((F(nulx)+F(x))/2));
Newton:=Newton+(h*((P(nulk)+P(k))/2));
q:=Abs(Newton-sum);
mmo1.Lines.Add(FloatToStr(q));

Until Abs(Newton-sum)<e;
Label1.Caption:=FloatToStr(sum);
lbl1.Caption:=FloatToStr(h);
end;
end.

[Greek9000]

не совсем понятно, почему у Вас цикл repeat ... until, ведь по указанной точности сразу ясно, на сколько интервалов надо разбить заданный диапазон. Поэтому тут нужен цикл for.
Кстати, рекомендую расчёт площади одного интервала вынести в отдельную функцию. Тогда восприятие алгоритма сильно упростится.
Получится что-то типа...

s := 0.0; // Текущее значение площади (интеграла)
x := x1; // Начало диапазона
dlt := dlt1; // Точность (шаг приращения ф-ии)
//... рассчитываем кол-во шагов N по длине интервала и точности
for i:=1 no N do begin
  // TrapSquare - функция расчёта площади трапеции
  s := s + TrapSquare(x); 
  x := x+dlt;
end;