Разложение функции в ряд Паскаль

[HuntMary]

Для х изменяющегося от a до b с шагом (b-a)/k, где (k=10), вычислить функцию f(x), используя ее разложение в степенной ряд в двух случаях:
а) для заданного n;
б) для заданной точности e (e=0.0001).
Для сравнения найти точное значение функции.

y=Xarctgx - lnкорень(1+x^2)
диапазон изменения 0,1<=x<=0,8
n=10
sum =

[newerow1989]

function f(x:real):real;
begin
   Result:=x*arctan(x)-ln(Sqrt(1+x*x));
end;

function f1(x:real):real;
const e=0.0001;
var n:integer;
    z:-1..1;
    xn,s1:real;
begin
   Result:=0;
   n:=0;
   z:=-1;
   xn:=1;
   Repeat
      n:=n+2;
      z:=-z;
      xn:=xn*x*x;
      s1:=z*xn/n/(n-1);
      Result:=Result+s1;
   Until Abs(s1)<e;
end;

const k=10;
      a=0;
      b=1;
      dx=(b-a)/k;
var x,y:real;
begin
   x:=a-dx;
   Repeat
      x:=x+dx;
      y:=f1(x);
      Write(x:5:2,'  ',y:5:5,'  ');
      y:=f(x);
      Writeln(y:5:5);
   Until x>=b;
   Readln;
end.