Вычислить интеграл с заданной точностью

type s1=record;
end;
var s:s1; d:file of s1; d1:text;
var
 n,i: integer;
 a,b,itog,h,M,x,sp: real;
function F(x:real):real;
begin
  F:=sin(sqrt(abs(ln(x)))*sqrt(abs(ln(x))));
end;
 
Begin
  assign(d,'1.dat');assign(d1,'Интеграл.txt');
  rewrite(d);rewrite(d1);
  write('Введите начало отрезка = '); readln(a);
  write('Введите конец  отрезка = '); readln(b);
    sp:=0;h:=(b-a)/20;
  i:=1;
while x<=b do begin
   x:=a+h*((2*i-1)/2);
   sp:=sp+F(x)*h;
   writeln(d1,sp:0:5);
   writeln(sp:0:5);   i:=i+1;
  end;
  close(d);close(d1);
END.

procedure pram(a,b,n)
var sum,integ,i: float;
begin
 
    integ=0;
        i=a;
    sum=(b-a)/n;
    while i<b do
    begin
        integ=integ+F(i)*sum;
                i=i+sum;
    end;
end;
 
BEGIN
 
while
 (abs(integ1-integ2)>eps)
do
 begin
    integ1=integ2;
    n=2;
    integ2=pram(a,b,n); 
                       {
                          pram - это процедура подсчета интеграла методом прямоугольников, 
                          с значениями а-начало отрезка, b-конец, n-количество разбиений
                        }
 end;
 
END.