программа в паскале

[Anton94.by]

Отредактировать программу чтобы работала, я её написал но выдаёт ошибку, если можете, то помогите. Спасибо.
Задание:
1. Составить программу, подсчитывающую с точностью значение определенного интеграла по формуле. Для достижения заданной точности использовать метод двойного пересчета. Результат сравнить со значением интеграла, подсчитанным с помощью первообразной F(x) по формуле Ньютона-Лейбница: .
Метод должен быть запрограммирован с помощью процедуры. Отрезок интегрирования, точность вычисления вводятся с клавиатуры, подынтегральная функция описывается с помощью функции пользователя.
Распечатать: отрезок интегрирования, точность вычисления, приближенное значение определенного интеграла, значение интеграла, подсчитанное по формуле Ньютона-Лейбница.

Метод: центральных прямоугольников:

1) \int_{0}^{1}{}2^{3^{x}}dx, e=10^{-3}, F(x)=\frac{1}{3ln2}(2^{3x}-1)

Внизу вложен конец задания если непонятно по этой формуле.

Исправьте ошибки в программе, если я что-то написал не в том порядке то исправьте. Спасибо...

program NV; 
uses crt; 
type TIntFunc=function ( x:real):real; 
function F(fn: TIntFunc; a, b, e:real):real; 
var 
s0, s, i, w, xi:real; 
j, n:integer; 
const eps=0.001; 
begin 
n:=1000; 
if a>b then 
begin 
s:=a; 
a:=b; 
b:=s; 
end; 
repeat 
s0:=s; 
s:=0; 
w:=(b-a)/n; 
for j:=0 to n-1 do 
begin 
xi:=a+(2*j+1)*(b-a)/2/n; 
s:=s+w*fn(xi); 
end; 
inc(n,10); 
until (n>12) and (abs(s-s0)<=e); 
I:=s; 
end; 
var 
a, b, e:real; 
begin 
function f1(x:real):real; 
far; 
begin 
f1:=exp(3*x*ln(2)); 
end; 
function f2(x:real):real; 
far; 
begin 
f2:=((1/3*ln*2*(exp(3*x-1))); 
end; 
write('Input A B E (separated by space):'); 
readln(a, b, e); 
writeln('i=',i(f1, a, b, e):15:5); 
writeln(' Checking:'); 
writeln('i=',f2(b)-f2(a):15:5); 
writeln 
end.

Ошибка там где я отметил курсивом.