|
PROGRAM
|
(PARI) a(n)=if(n<0, 0, prod(k=1, n, if(k%2, k+2, 1-k)))
(PARI) {a(n)=local(A); if(n<0, 0, A=matrix(2, n+1); for(k=0, n, A[2, k+1]=if(k%2, 3, 1); A[1, k+1]=if(k<2, (-1)^k, if(k%2, -(k-2)*(k-1), -k*(k+1)))); contfracpnqn(A)[2, 1])} /* Michael Somos Jul 15 2003 */
|