function a = coef_cos (input1, input2)
  N = 30;
  n = [0:N-1];
  a = zeros(1, 2*N);
  a(1:2:2*N) = (-1).^n./factorial(2*n);
endfunction
