Табулирования функции. в pascal в файл tabl.txt вывести таблицу значений для функции y=(2х+1)/(2х-1) (записано дробью)

var

 FD: text;

 y, step, first, last: real;

begin

 write('Введите начальное значение, конечное значение и шаг:');

 read(first, last, step);

 assign(FD, 'tabl.txt');

 rewrite(FD);

 writeln(FD, 'x':4, 'f(x)':10);

 while (first <= last) do

 begin

   if (first.Between(0.41, 0.58)) then  

     writeln(FD, first:4:2, 'infinity':10)

   else  

   begin

     y := (2 * first + 1) / (2 * first - 1);

     writeln(FD, first:4:2, y:10:4);

   end;

   first := first + step;

 end;

 close(FD);

end.