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

Program integrali;
const eps=0.0001;                       {h- visota}
var   h,I1,I2,k,i,y,x:real;
      i,n:integer;
begin
     readln(x);
     n:=0;
     I1:=0;
     repeat
     I2:=I1;            {I1 tekushchee znachenie, I2 predidushchee znachenie}
     n:=n+1;                 {pri absolyoutnoi tochnosti}
     I1:=0;
     h:=1/n;               {iz formuli H=(b-a)/n}
     for i:=0 to n do
        begin
        x:=i/n;              {Xi=A+i*H=0+i*1/n}
        y:=exp(2*x)-exp(x)-ln(abs(x)+1)*cos(x);
        if i=0 or i=n then
        I1:=I1+y
         else if i mod 2=0
         then I1:=I1+2*y
         else I1:=I1+4*y;
        end;
     I1:=I1*h/3;
     writeln (I1);
     k:=abs(I1-I2);    {k- uslovia vihoda}
     if I2=0 then k:=1;
     untiln (k<=e) and (k>=0);
end.

Program integral;
const eps=0.0001;                       {h- visota}
var   h,I1,I2,k,y,x:real;
      i,n:integer;
begin
     write('Vvedite x');
     readln(x);
     n:=0;
     I1:=0;
     repeat
     I2:=I1;            {I1 tekushchee znachenie, I2 predidushchee znachenie}
     n:=n+1;                 {pri absolyoutnoi tochnosti}
     I1:=0;
     h:=1/n;               {iz formuli H=(b-a)/n}
     for i:=0 to n do
        begin
        x:=i/n;              {Xi=A+i*H=0+i*1/n}
        y:=exp(2*x)-exp(x)-ln(abs(x)+1)*cos(x);
        if (i=0) or (i=n) then
        I1:=I1+y
         else if i mod 2=0
         then I1:=I1+2*y
         else I1:=I1+4*y;
        end;
     I1:=I1*h/3;
     writeln (I1);
     k:=abs(I1-I2);    {k- uslovia vihoda}
     if I2=0 then k:=1;
     until (k<=eps) and (k>=0);
end.

uses crt;
const
  a = 0;{нижний предел}
  b = 1;{верхний предел}
  e=0.0001;{точность}
function F(x:real):real;
begin
F:=exp(2*x)-exp(x)-ln(abs(x)+1)*cos(x);
end;
 
var n,i:integer;
    h,k,s1,s2: real;
begin
clrscr;
write('Количество разбиений n=');
readln(n);//начальное число
s1:=0;//итоговый интеграл
h:=(b-a)/n;//шаг
s2:=h*F(a);//текущий интеграл
repeat
s1:=s2;//запомним текущий
h:=(b-a)/n;
k:=a;//левый край
for i:=1 to n-1 do
 begin
  k:=k+h;//наращиваем
  if i mod 2 = 0 then s2:=s2+2*f(k)//считаем новый интеграл
  else s2:=s2+4*f(k);
 end;
s2:=(s2+f(a)+f(b))*h/3;//конечное значение
n:=n*2;//если точность не достигнута удваиваем количество разбиений
until abs(s1-s2)<=e;
write('Интеграл=',s2:0:6);
readln
end.