Проверьте код (интегралы на паскали)

[Айдар]

Проверьте плиз код программы нахождения интеграла методом трапеции..
Суть такая ищем интеграл из фун-ции (х*х+3), с заданным кол-вом шагов, нашли, теперь ищем еще раз с увеличенным кол-вом шагов в 2 раза, сравниваем результат если разность интегралов меньше e=0,001 то выводим результат (любой из двух), если нет, то снова ищем интеграл n:=2*n;

Код:

var
n,i:integer;
a,b,e,shag,sum,itog,itog2:real;
function F(x:real):real;
begin
f:=x*x+3;
end;

procedure int1 (var a,b,e:real;n:integer);
begin
shag:=(b-a)/n;
sum:=0;
for i:=1 to n-1 do
begin
sum:=sum+f(a+i*shag);
end;
sum:=sum+(f(a)+f(b))/2;
itog:=shag*sum;
writeln ('Ответ при  ', n, ' Повторах',itog:10:5);
{until itog2-itog<e;}
readln;
end;

procedure int2 (var a,b,e:real;n:integer;var itog:real);
begin
repeat
n:=2*n;
shag:=(b-a)/n;
sum:=0;
for i:=1 to n-1 do
begin
sum:=sum+f(a+i*shag);
end;
sum:=sum+(f(a)+f(b))/2;
itog2:=shag*sum;
writeln ('Ответ при  ', n, ' Повторах',itog2:10:5);
until itog2-itog<e;
readln;
end;

begin
write ('a=');
readln (a);
write ('b=');
readln (b);
write ('N=');
readln (n);
e:=0.001;
writeln ('e=0.001');
int1 (a,b,e,n);
int2 (a,b,e,n,itog);
readln;
end.

PS кто может проконсультировать по ведюхам?