Вычислить и вывести на экран в виде таблицы значения функции, заданной с помощью ряда Тейлора

uses crt;
{процедура вычисления суммы ряда(s) и числа итераций(i)}
procedure F(x,e:real;var s:real;var i:longint);
var t:real;
begin
s:=1/x;{начальные значения при i=1}
t:=1/x;
i:=1;
while abs(t/(2*i+1))>e do{пока модуль очередного члена ряда больше точности}
 begin
   i:=i+1;{наращиваем счетчик}
   t:=t/(x*x);{домножаем на x^2}
   s:=s+t/(2*i-1);{делим на очередное нечетное и суммируем}
 end;
s:=2*s; {умножаем на 2}
end;
var xn,xk,e,dx,s:real;
    n:longint;
begin
clrscr;
repeat
write('Введите начало интервала x>=-1 или x>=1 xn=');
readln(xn);
until (xn<=-1)or(xn>=1);
repeat
write('Введите конец интервала по модулю больше ',abs(xn):0:1,' xk=');
readln(xk);
until abs(xk)>xn;
repeat
write('Введите точность в интервале (0,1) e=');
readln(e);
until(e>0)and(e<1);
repeat
write('Введите шаг табуляции, положительное число меньше ',xk-xn:0:2,' dx=');
readln(dx);
until (dx<xk-xn)and(dx>0);
clrscr;
writeln('      Таблица табулирования функции ln((x+1)/(x-1))');
writeln('на интервале ',xn:0:2,'  ',xk:0:2,' с шагом ',dx:0:2,' с точностью ',e:0:6);
writeln('----------------------------------------');
writeln('|       |           |     x+1    |     |');
writeln('|   x   | вычислено |  ln-----   |  n  |');
writeln('|       |           |     x-1    |     |');
writeln('----------------------------------------');
while xn<=xk+e do
 begin
  F(xn,e,s,n);
  writeln('|',xn:6:2,' |',s:10:6,' |',ln((xn+1)/(xn-1)):10:6,'  |',n:4,' |');
  xn:=xn+dx;
 end;
writeln('----------------------------------------');
readln
end.