(PARI)我的(N=30,x='x+O('x^N));Vec(塞拉普拉斯(cos(2*log(1+x)))
(PARI)a(n)=和(k=0,n \ 2,(-4)^k*斯特林(n,2*k,1));
(PARI)a(n)=(-1)^n*(prod(k=0,n-1,2*I+k)+prod(k=0,n-1,-2*I+k))/2;
(PARI)a_vector(n)=我的(v=向量(n+1));v[1]=1;v[2]=0;对于(i=2,n,v[i+1]=-(2*i-3)*v[i]-(i^2-4*i+8)*v[1]);v;