h:=(b-a)/n;

y1:=LevPr.integralf1(a, b, h);

l1: n:=n*2;

h:=(b-a)/n;

y:=LevPr.integralf1(a, b, h);

if (abs(y-y1)>eps) then begin y1:=y; goto l1; end;

levprf1:=y1;

end;

function LevPr.levprf2(a:real; b:real; eps: real; n:integer):real;

label l1;

var

i: integer;

x, y, h, y1:real;

begin

h:=(b-a)/n;

y1:=LevPr.integralf2(a, b, h);

l1: n:=n*2;

h:=(b-a)/n;

y:=LevPr.integralf2(a, b, h);

if (abs(y-y1)>eps) then begin y1:=y; goto l1; end;

levprf2:=y1;

end;

Обработка измерений

(процедура, подсчитывающая коэффициенты)

begin

label1.visible:=true;

label4.visible:=true;

n:=strtoint(edit1.Text);

k:=0; l:=0; m:=0;

v:=0; w:=0; r:=0;