Нужно исправить ошибку

program zadanie;
 
const
  a = 2; { нижний предел }
  b = 4; { верхний предел }
 
function F(x: real): extended;
begin
  F := ln(x) * exp((ln(x) * (1 + x)) * ln(x));
  { после математического преобразования }
end;
 
function simpson(a, b: real; n: integer): real;
var
  s, h, k: real;
  i: integer;
begin
  h := (b - a) / n;
  k := a; { левый край }
  s := F(a) + F(b);
  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;
  simpson := s * h / 3; { Возвращается вычисленное значение }
end;
 
var
  n: integer;
  s1, s2, e: real;
 
begin
  n := 16;
  e := 0.1;
  // writeLn('Количество разбиений n=');
  // readLn(n); { начальное число }
  // writeLn('Точность e=');
  // readLn(e);
  s1 := simpson(a, b, n); { Начальное значение }
  repeat
    s2 := s1; { Запоминается предыдущее значение }
    n := n * 2; { Увеличивается количество шагов }
    s1 := simpson(a, b, n); { Рассчитывается новое значение }
  until abs(s1 - s2) <= e; { Повторяется до достижения необходимой
    точности }
  writeLn('Интеграл=', s1:0:6);
  readLn;
 
end.