Протабулировать функцию arcctg(x), заданную в виде ряда Тейлора

const e=0.001;
function S(x:real;var n:integer):real;
var t:real;
begin
n:=0;
t:=-x;
result:=t;
while abs(t/(2*n+1))>e do
 begin
  n:=n+1;
  t:=-t*x*x;
  result:=result+t/(2*n+1);
 end;
result:=pi/2+result;
end;
function f(x:real):real;
begin
result:=pi/2-arctan(x)
end;
var xn,xk,dx,x,y:real;
    n:integer;
begin
repeat
 writeln('Введите интервал -1<xn<xk<1');
 readln(xn,xk);
until(xn>-1)and(xk>xn)and(xk<1);
repeat
write('Введите шаг dx<',xk-xn:0:2,' dx=');
readln(dx);
until dx<xk-xn;
clrscr;
writeln('            Таблица значений функции arcctg(x)');
writeln('представленной разложением в ряд Тейлора на интервале [',xn:0:2,',',xk:0:2,']');
writeln('        с шагом ',dx:0:2,' с точностью 0.001');
writeln('---------------------------');
writeln('|  x  |  S(x) | n |  f(x) |');
writeln('---------------------------');
x:=xn;
while x<xk+dx/2 do
 begin
  y:=S(x,n);
  writeln('|',x:5:2,'|',y:7:3,'|',n:3,'|',f(x):7:3,'|');
  x:=x+dx;
 end;
writeln('---------------------------');
end.

repeat
 writeln('Введите интервал -1<xn<xk<1');//например -0.9 и 0.9
 readln(xn,xk);
until(xn>-1)and(xk>xn)and(xk<1);
repeat
write('Введите шаг dx<',xk-xn:0:2,' dx=');//например 0.2
readln(dx);
until dx<xk-xn;