Работа с функциями Pascal

[Harrr]

в программе
function integral(a:real; b: real; h: real):real;
var
x, y:real;
begin
x:=a;
y:=0;
while x<b do begin
y:=y+x*x*h;
x:=x+h;
end;
integral:=y;
end;

label l1;
var a,b,h,x,y, eps, y1 : real;
n: integer;
begin
writeln('Vvedite nijniy predel a ');
readln(a);
writeln('Vvedite verhniy predel b ');
readln(b);
writeln('Vvedite tochnost eps ');
readln(eps);
writeln('Vvedite kolichestvo otrezkov n ');
readln(n);
h:=(b-a)/n;
y1:=integral(a,b,h);
l1: n:=n*2;
h:=(b-a)/n;
y:=integral(a,b,h);
if (abs(y-y1)>eps) then begin goto l1; y1:=y; end;
Writeln('Integral= ', y:6:2);
end.
выполнение зависает на строке x:=x+h;
есть подозрение, что в функцию integral не передаются значения некоторых переменных. помогите, пожалуйста!

[Z1000000]

if (abs(y-y1)>eps) then begin goto l1; y1:=y; end;

В этой строке никогда не выполняется операция y1:=y;
y1:= в программе не пересчитывается. И погрешность вычисления не уменьшается. Надо так.

if (abs(y-y1)>eps) then begin  y1:=y; goto l1; end;