Вычислить и вывести на экран в виде таблицы значения функции

var n: integer;
    xnac, xkon, dx, x, a, s, eps: real;
begin
  write('    Начальное значение аргумента: ');
  readln(xnac);
  write('     Конечное значение аргумента: ');
  readln(xkon);
  if xnac > xkon
    then begin
      writeln('Начальное больше конечного, значения обменяны местами.');
      x := xnac;
      xnac := xkon;
      xkon := x
    end;
  repeat
    write('Шаг изменения аргумента (dx > 0): ');
    readln(dx)
  until dx > 0;
  if xkon - xnac < dx then writeln('Шаг слишком большой, будет выведено только одно значение функции.');
  repeat
    write('   Точность вычислений (eps > 0): ');
    readln(eps)
  until eps > 0;
  writeln;
  writeln('Таблица значений ряда Тейлора для функции f(x)=ln((x+1)/(x-1))');
  writeln('на интервале от ', xnac, ' до ', xkon);
  writeln('с шагом ', dx, ' с точностью ', eps);
  writeln;
  writeln('+------------------------+------------------------+------------------------+');
  writeln('|           x            |          S(x)          |  Просуммировано членов |');
  writeln('+------------------------+------------------------+------------------------+');
  x := xnac;
  xkon := xkon + dx / 2;
  eps := eps / 2;
  while x < xkon do
    begin
      write('| ', x:22, ' | ');
      if abs(x) <= 1
        then writeln('не определено':22, ' | ', 'не определено':22, ' |')
        else begin
          a := 1 / x;
          s := a;
          n := 1;
          while abs(a) >= eps do
            begin
              a := a * (n - 0.5) / (n + 0.5) / x / x;
              s := s + a;
              inc(n)
            end;
          writeln(2 * s:22, ' | ', n:22, ' |')
        end;
      x := x + dx
    end;
  writeln('+------------------------+------------------------+------------------------+');
  readln
end.