Почему зацикливается следующая программа и как её можно преобразовать к более удобной форме???

program LabWork4;
var
 x,x0,xn,y,y0,yn,hx,f:real;
 n:integer;
begin
 f:=2*x+3;
 writeln('График функции ограничен прямоугольником со сторонами x=a, x=b, y=c, y=d. Введите значения этих сторон.');
 write('x=');
 readln(x0);
 write('x=');
 readln(xn);
 write('y=');
 readln(y0);
 write('y=');
 readln(yn);
 writeln;writeln('Введите значения шага аргумента');
 readln(hx);
 writeln;writeln('Исходные данные: График функции ограничен прямоугольником со сторонами x=',x0,', x=',xn,', y=',y0,', y=',yn,'. Шаг значений аргумента hx=',hx);
 {Вывод таблицы значений функции с помощью оператора while}
 writeln;writeln('Таблица значений функции (оператор while):');
 x:=x0;
  while x<=xn do begin
   y:=f;
   if y>=y0 then
    if y<=yn then
    writeln ('x=',x,' y=',y)
   else x:=x+hx;
  end;
 {Вывод таблицы значений функции с помощью оператора repeat}
 writeln;writeln('Таблица значений функции (оператор repeat):');
 x:=x0;
  repeat
   y:=f;
   if y>=y0 then
    if y<=yn then
    writeln ('x=',x,' y=',y)
   else x:=x+hx;
  until x>xn;
 {Цикл исследования функции f(x)}
 x:=x0;
 n:=0;
  repeat
  y:=f;
   if y>=y0 then
    if y<=yn then
    n:=n+1
    else
    x:=x+1
  until x>xn;
 {Вывод результатов исследования функции}
 writeln;writeln('Количество точек графика функции лежащих в прямоугольнике со сторонами x=',x0,', x=',xn,', y=',y0,', y=',yn,' равно',n);
end. {Program LabWork4}

program LabWork4;
function F(x:real):real;
begin
F:=2*x+3;
end;
var
 x,x0,xn,y,y0,yn,hx:real;
 n:integer;
begin
 writeln('Введите значения координат прямоугольника.');
 write('x0='); read(x0);
 write('y0='); read(y0);
 write('xn='); read(xn);
 write('yn='); read(yn);
 write('Введите значения шага аргумента hx='); read(hx);
 writeln('Исходные данные:');
 writeln('График функции ограничен прямоугольником со координатами x0=',x0,',y0=',y0,',xn=',xn,',yn=',yn);
 writeln('Шаг значений аргумента hx=',hx);
 {Вывод таблицы значений функции с помощью оператора while}
 writeln('Таблица значений функции (оператор while):');
 x:=x0;
  while x<=xn do
   begin
    y:=F(x);
    writeln ('x=',x:4:1,' y=',y:5:2);
    x:=x+hx;
  end;
 {Вывод таблицы значений функции с помощью оператора repeat}
 writeln;
 writeln('Таблица значений функции (оператор repeat):');
 x:=x0;
  repeat
   y:=F(x);
   writeln ('x=',x:4:1,' y=',y:5:2);
   x:=x+hx;
  until x>xn;
 {Цикл исследования функции f(x)}
 x:=x0;
 n:=0;
 repeat
  y:=F(x);
  if(x>=x0)and(x<=xn)and(y>=y0)and(y<=yn) then n:=n+1;
  x:=x+hx;
 until x>xn;
{Вывод результатов исследования функции}
writeln('Количество точек графика функции лежащих в прямоугольнике=',n);
end. {Program LabWork4}