Найти значение интеграла (методом прямоугольника)

const a=1; b=2; //интервал интегрирования
      tn=2;tk=12; //интервал вычислений
      ht=2; //шаг вычислений
      e=0.01; //точность вычислений
//подинтегральная функция
function f(x,t:real):real;
begin
f:=exp(ln(x)*t*exp(ln(ln(x))*t));
end;
//функция вычисления интеграла методом средних прямоугольников
//при числе разбиений=n
function IntRect(n,t:integer):real;
var h,x,s:real;
    i:integer;
begin
h:=(b-a)/n;
x:=a+h/2;
s:=0.0;
for i:=1 to n do
 begin
  s:=s+f(x,t);
  x:=x+h;
 end;
IntRect:=s*h;
end;
var h,x,s1,s2:real;
    t,n,i:integer;
begin
t:=tn;
while t<=tk do
 begin
  n:=10;
  s1:=IntRect(n,t);
  repeat
   n:=n*2;
   s2:=s1;
   s1:=IntRect(n,t);
  until abs(s1-s2)<=e;
  writeln('t=',t:2,' S=',s1:5:2);
  t:=t+ht;
 end;
end.