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