Исправить ошибки в программе и дописать

[Anton94.by]

Исправъте ошибки в программе.

Задание для программы: Составить программу в Паскале, подсчитывающую с точностью значение определенного интеграла по формуле. Для достижения заданной точности использовать метод двойного пересчета. Результат сравнить со значением интеграла, подсчитанным с помощью первообразной F(x) по формуле Ньютона-Лейбница: .
Метод должен быть запрограммирован с помощью процедуры. Отрезок интегрирования, точность вычисления вводятся с клавиатуры, подынтегральная функция описывается с помощью функции пользователя.
Распечатать: отрезок интегрирования, точность вычисления, приближенное значение определенного интеграла, значение интеграла, подсчитанное по формуле Ньютона-Лейбница.
Метод: центральных прямоугольников:
\int_{0}^{1}{}2^{3}dx, e=10^{-3}, F(x)=\frac{1}{3ln2}(2^{3x}-1)
( если непонятна эта формула, то посмотрите внизу вложение т.к. я незнаю как вставить здесь картинку).

Program integral1; 
Uses Crt; 
Const eps=0.01;
Var
a,b,eps,integral :Real;
nn :Integer;
Function F(z:Real):Real;
Begin
f:=1/(z-1);
End;
Procedure Trap(a1,b1:Real; n:Integer; 
eps1:Real; var tk:Real);
Var
i :Integer;
x,t0,s,h :Real;
Begin
t0:=f(a1)+f(b1)/2; 
While True do 
begin
h:=(b1-a1)/n; {шаг}
s:=f(a1)+f(b1)/2;
For i:=1 to n-1 do
begin
x:=a1+i*h;
s:=s+f(x);
end;
tk:=s*h;
If abs(tk-t0)<eps then Exit;
t0:=tk;
n:=2*n; 
end;
End;
Begin
Clrscr;
Write('Введите пределы интегрирования a,b’);
ReadLn (a,b);
nn:=6; 
Trap(a,b,nn,eps,integral);
WriteLn('Значение интеграла ',integral:8:4);
Repeat Until KeyPressed;
End.

Я смотрел подобное задание и по подомному делал это, у меня получилась эта программа.

Если код правильный, то напишите что код правильный.