Вычислить значения заданной функции - Turbo Pascal - Ответ 620433

uses crt;
var xn,xk,e,dx,f,s,si,t:real;
    i,k:integer;
begin
clrscr;
repeat
write('Введите конец интервала xk<-1 xk=');
readln(xk);
until xk<-1;
repeat
write('Начало интервала xn<',xk:0:2,' xn=');
readln(xn);
until xn<xk;
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);
write('Press Enter...');
readln;
clrscr;
writeln('      Таблица табулирования функции arctg(x)');
writeln('на интервале ',xn:0:2,'  ',xk:0:2,' с шагом ',dx:0:2,' с точностью ',e:0:6);
writeln('----------------------------------');
writeln('|   x   | вычислено |  arctg(x)  |');
writeln('----------------------------------');
k:=0;
while xn<=xk+e do
 begin
  s:=0;//внутренний цикл вычисления суммы
  t:=-1/xn;
  si:=-1/xn;
  i:=1;
  while abs(s-si)>e do//считаем пока модуль соседних сумм больше точности
   begin
    s:=si;
    i:=i+1;
    t:=-t/(xn*xn);
    si:=si+t/(2*i-1);
   end;
  f:=-pi/2+s;//добавляем -pi/2
  writeln('|',xn:6:2,' |',f:10:6,' |',arctan(xn):10:6,'  |');
  k:=k+1;
  if k mod 18=0 then{если больше 18 итераций, все на экран не войдет}
   begin
    write('Press Enter...');
    readln;{делаем задержку}
   end;
  xn:=xn+dx;
 end;
writeln('----------------------------------');
readln
end.