Численное интегрирование

uses crt;
label start;
var a,b,h,x,s0,s1,e: real;
    i: integer;
function f(x: real): real;
begin
 f:=; {сюда вставить функцию из варианта}
end;
begin
 clrscr;
 write('a/b/e: '); {a - нижний предел интегрирования, b - верхний}
 readln(a,b,e); {e - точность, например, e=0.01}
 h:=(b-a)/2;
 i:=1;
 x:=a;
 s1:=0;
 s0:=0;
 start:
 while x<b do
  begin
   s1:=s1+h*f(x);
   x:=x+h;
  end;
 if abs(s1-s0)<e then
  begin
   writeln('Площадь = ',s1:5:5,' Кол-во итераций = ',i);
   writeln('Шаг разбиения в последней итерации = ',h:10:5);
  end
   else
    begin
     h:=h/2;
     s0:=s1;
     inc(i);
     x:=a;
     s1:=0;
     goto start;
    end;
 readln;
end.

uses crt;
var a,b,S,h,integ: real;
    i,n: integer;
 
function f(x: real): real;
begin
 f:=; {сюда вставить функцию из варианта}
end;
 
begin
 clrscr;
 write('a/b/n: '); {a - нижний предел интегрирования, b - верхний}
 readln(a,b,n);  {n - число разбиений, обычно подхоит n=10} 
 h:=(b-a)/n;
 for i:=1 to n-1 do
  S:=s+f(a + h*i);
 integ:=h*((f(a)+f(b))/2+S);
 writeln('S=',integ);
 readln;
end.

uses crt;
label xx2,xy2;
var i1,i2,x1,x2,x0,s,y0,y1,y2,a,b,h,x,s0,s1,n,e: real;
    i: integer;
function f(x: real): real;
begin
 f:=; {сюда вставить функцию из варианта}
end;
begin
 clrscr;
 write('a/b/e: '); {a - нижний предел интегрирования, b - верхний}
 readln(a,b,e); {e - точность, например, e=0.01}
 i1:=0;
 n:=4;
 h:=(b-a)/4 ;
 xy2:
 i2:=0;
 i:=2;
 xx2:
 x2:=a+i*h;
 x1:=x2-h;
 x0:=x1-h;
 y0:=f(x0);
 y1:=f(x1);
 y2:=f(x2);
 s:=y0+4*y1+y2;
 i2:=i2+s;
 if i<n then
  begin
  i:=i+2;
  goto xx2;
  end
  else
   begin
    i2:=i2*h/3;
    if abs(i1-i2)<e then
    writeln('Result: ',i2:3:2)
    else
     begin
      i1:=i2;
      n:=n*2;
      h:=h/2;
      goto xy2;
     end;
    end;
 readln;
end.