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

[Anton94.by]

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

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

1) \int_{0}^{1}{}2^{3^{x}}dx, e=10^{-3}, F(x)=\frac{1}{3ln2}(2^{3x}-1)
Если непонятно задание по этим формулам то посмотрите вложение, т.к. я не умею здесь выставлять картинки.

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

function I(fn:TIntFunc;a,b,e:real):real;
var
s0,s,w,xi:real;
j,n:integer;
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

{Zadacha 1}

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=',FF1(b)-FF1(a):15:5);
writeln;

writeln('I=',I(f2a,0,3,1E-3):15:5);
writeln('I=',I(f2b,3,4,1E-2):15:5);

end.