Найти сумму ряда (Pascal)

[pitak]

Здравствуйте. Прошу помочь с задачей.
Задание: найти сумму ряда.
Ряд: cos x +...+(cos(2n-1)*x)/(2n-1)^2
Диапазон изменения аргумента: pi/10<=x<=9pi/10
n=40
Функция Y=(pi^2/8)-(pi/4)*|x|
Программа выдает неверный результат. Подскажите в чем ошибка и как исправить?

Program summa;
Var x, S, h, K, Y, a, b: real;
N1, N2, n: integer;
begin
  a:=Pi/10;
  b:=9*Pi/10;
  N1:=1;
  N2:=40;
  h:=((b-a)/10);
  x:=a;
  S:=1;
  n:=1;
  K:=1;
  Y:=3;
    write('X':2); 
    write('S':4);
    writeln('Y':5);
    While x<= b do
    begin
    S:=0;
    K:=1;
    For n:=N1 to N2 do
    begin
    K:=(cos(2*n-1)*x)/((2*n-1)*(2*n-1));
    S:=S+K;
    end;
    Y:=((Pi*Pi)/8)-(Pi/4)*(ABS(x));
    write(x:4:2); 
    write(S:5:2);
    writeln(Y:5:2);
     x:=x+h;
    end;
    end.

[pitak]

так и не исправил((