Вычислить значение функции ln((x+1)/(x-1)), используя ряд Тейлора

var x,e,t,s:real;
    n:integer;
begin
repeat
write('|x|>1 x=');
readln(x);
until abs(x)>1;
e:=0.0001;
n:=1;
t:=1/x;
s:=t;
while abs(t)/(2*n-1)>e do
 begin
  n:=n+1;
  t:=t/x/x;
  s:=s+t/(2*n-1);
 end;
s:=2*s;
writeln('s=',s:0:4);
write('ln((x+1)/(x-1))=',ln((x+1)/(x-1)):0:4)
end.

const e=0.0001;
function f(x:real;var n:integer):real;
var t,s:real;
begin
n:=1;
t:=1/x;
s:=t;
while abs(t)/(2*n-1)>e do
 begin
  n:=n+1;
  t:=t/x/x;
  s:=s+t/(2*n-1);
 end;
f:=2*s;
end;
var xn,xk,x,s,dx:real;
    n:integer;
begin
repeat
write('Введите начало интервала xn>1 xn=');
readln(xn);
until xn>1;
repeat
write('Введите конец интервала, больше ',xn:0:1,' xk=');
readln(xk);
until xk>xn;
repeat
write('Введите шаг табуляции, меньше ',xk-xn:0:1,' dx=');
readln(dx);
until dx<xk-xn;
x:=xn;
writeln('   Табулирование функции F(x)=ln(((x+1)/(x-1))');
writeln('представленной разложением в степенной ряд');
writeln('   на интервале ',xn:0:1,'..',xk:0:1,' с шагом ',dx:0:1);
writeln;
writeln('        -------------------------');
writeln('        |   X   |    S    |  n  | ');
writeln('        -------------------------');
while x<=xk+dx/2 do
 begin
  s:=f(x,n);
  writeln('        |',x:6:2,' |',s:8:4,' |',n:5,'|');
  x:=x+dx;
 end;
writeln('        -------------------------');
//проверка последнего члена по аналитической формуле
x:=x-dx;
writeln('Проверка последнего члена');
writeln('x=',x:0:2);
writeln('ln((x+1)/(x-1))=',ln((x+1)/(x-1)):0:4);
end.