Вычислить значение функции arctg(x), заданной с рядом Тейлора

var x,e,t,s:real;
    n:integer;
begin
repeat
write('Введите x>1 x=');
readln(x);
until x>1;
e:=0.00001;
n:=0;
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:=s+pi/2;
writeln('S=',s:0:4);
write('arctg(x)=',arctan(x):0:4)
end.

function f(x,e:real;var n:integer):real;
var t,s:real;
begin
n:=0;
t:=-1/x;
s:=pi/2+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:=s;
end;
var xn,xk,dx,e,x,y:real;
    n:integer;
begin
repeat
write('xнач>1 xнач=');
readln(xn);
until xn>1;
repeat
write('xкон>',xn:0:2,' xкон=');
readln(xk);
until xk>xn;
repeat
write('xнач<dx<xкон dx=');
readln(dx);
until dx<xk-xn;
repeat
write('0<e<1 e=');
readln(e);
until (e>0)and(e<1);
writeln('   Табуляция функции arctg(x) при х>1');
writeln('представленной разложением в ряд Тейлора');
writeln('на интервале [',xn:0:2,';',xk:0:2,'] с шагом ',dx:0:2);
writeln('       с точностью е=',e:0:6);
writeln('       ------------------------');
writeln('       |  x  |   f(x)   |  n  |');
writeln('       ------------------------');
x:=xn;
while x<xk+dx/2 do
 begin
  y:=f(x,e,n);
  writeln('       |',x:5:2,'|',y:8:4,'  |',n:4,' |');
  x:=x+dx;
 end;
writeln('       ------------------------');
writeln('Проверка последнего значения arctg(',x-dx:0:2,')=',arctan(x-dx):0:4)
end.