Приближенное вычисление значения функции cos(x)

Uses crt;
const n=28;
function CosX(x:real):extended;
var p:real;
    s,ch,zn:extended;
i:integer;
begin
ch:=1;
i:=1;
s:=0;
zn:=1;
p:=1;
repeat
s:=s+p*ch/zn;
i:=i+1;
ch:=ch*x*x;
zn:=zn*(2*i-2)*(2*i-3);
p:=p*(-1);
until abs(ch/zn)<0.0001;
CosX:=s;
end;
var x,h,yp,yt,pa,po:extended;
    i:integer;
begin
clrscr;
h:=2*pi/(n-1);
writeln('Результаты вычислений:':50);
writeln('-------------------------------------------------------------');
writeln('|   |        |        |        |         Погрешность        |');
writeln('| № |   x    | cos(x) |  CosX  |-----------------------------');
writeln('|   |        |        |        | абсолютная | относительная |');
writeln('-------------------------------------------------------------');
x:=-pi;i:=0;
while x<=pi do
 begin
  i:=i+1;
  yp:=CosX(x);
  yt:=cos(x);
  pa:=abs(yp-yt);
  po:=abs(pa/yp);
  gotoXY(1,whereY);
  write('|',i:2);
  gotoXY(4,whereY);
  write(' |',x:7:4);
  gotoXY(14,whereY);
  write('| ',yt:7:4);
  gotoXY(23,whereY);
  write('| ',yp:7:4);
  gotoXY(32,whereY);
  write('| ',pa:10:7);
  gotoXY(45,whereY);
  write('| ',po:10:7);
  gotoXY(61,whereY);
  write('|');
  writeln;
  if i=18 then readln;
  x:=x+h;
 end;
writeln('-------------------------------------------------------------');
readln
end.