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

Program integr;
Uses Crt;
Const eps=0.001;
Var i:integer;
h,a1,b1,s2,d1,d2,s1,SUM,plosh:real; raznica,n:real;
Function Ingr (x: real) : real; 
begin
Ingr:=1/(x*ln(x)); 
end;
Procedure Go(a,b:real); var i1:integer; 
begin
Write('Vvedite nomer uchastka integrirovania n= ');
Read(n); 
h:=(b-a)/n; 
i1:=1; 
s1:=0;
while i1<=(n-1) 
do
begin s1:=s1+Ingr(a+i1*h); 
i1:=i1+2; 
end; 
s2:=0; 
i1:=2;
While i1<=(n-2) 
do
begin
s2:=s2+Ingr(a1+i1*h); 
i1:=i1+2; 
end ;
d1:=Ingr(a); 
d2:=Ingr(b);
SUM:=h*(d1+d2+4*s1+2*s2)/3;
end;
begin
ClrScr;
write('Vvedite niznuju granizzu a=');
Read(a1);
Write('Vvedite verhnuju granizzu b=');
Read(b1);
WriteLn('Integral  Newleib   Eps    shag');
repeat
Go(a1,b1);
plosh:=(ln(ln(b1))-ln(ln(2)))*ln(10)-(ln(ln(a1))-ln(ln(2)))*ln(10); 
raznica:=(abs(plosh-sum));
WriteLn(sum:6:3,' ',plosh:8:5,' ',raznica:8:5,' ',h:7:4);
until(abs(raznica)<eps);
readln;
end.

Function Ingr (x: real) : real;
begin
Ingr:=ln(10)/(x*ln(x));
end;

uses crt;
function F(x:real):real;
begin
F:=ln(10)/x/ln(x);
end;
function P(x:real):real;
begin
P:=ln(10)*(ln(ln(x))-ln(ln(2)));
end;
var a,b,h,k,s:real;
    n,i:integer;
begin
clrscr;
a:=2;
b:=4;
n:=36;
h:=(b-a)/n;
s:=0;
k:=a;
for i:=1 to n-1 do
 begin
  k:=k+h;
  if i mod 2=0 then s:=s+2*f(k)
  else s:=s+4*f(k);
 end;
s:=(s+f(a)+f(b))*h/3;
writeln('Интеграл=',s:0:5);
writeln('Проверка=',P(b)-P(a):0:5);
readln
end.

uses crt;
function F(x:real):real;
begin
F:=ln(10)/x/ln(x);
end;
function P(x:real):real;
begin
P:=ln(10)*(ln(ln(x))-ln(ln(2)));
end;
 
var x,a,b,h,s:real;
    n:integer;
begin
clrscr;
a:=2;
b:=4;
n:=36;
h:=(b-a)/n;
s:=0; x:=a+h;
while x<b do
begin
s:=s+4*f(x);
x:=x+h;
s:=s+2*f(x);
x:=x+h;
end;
s:=h/3*(s+f(a)-f(b));
writeln;
writeln('Интеграл = ',s:0:5);
writeln('Проверка = ',p(b)-p(a):0:5);
readkey;
end.