len=12;s[p]:=(-1)^p*BernoulliB[2p+2]/(2p+1)/(2p+2);做[m[n,1]=0,{n,0,len}];做[m[n,2]=s[n+1]/s[n],{n,0,len-1}];做[m[n,k]=
如果[OddQ[k],m[n+1,k-2]+m[n+1,k-1]-m[n,k-1],
m[n+1,k-2]*m[n+1,k-1]/m[n,k-1]],{k,3,len},{n,0,
len-k+1}];Do[m[n,1]=s[n],{n,0,len}];
表[m[0,k],{k,1,len}]//分母