| 数据
|
1, 1, 7, 87, 1577, 37921, 1143991, 41734167, 1793837945, 89100737537, 5038278258759, 320488252355991, 22712229678267017, 1778818548078114337, 152926844472960316055, 14348332105800041202903, 1461880180517958608890585, 161034066043430013259095681, 19105043857756090069661974951, 2432865068875486088572762200535, 331511875063241457659846364208233, 48205214775404458968179455649349921, 7461345443274460130807423699070922103
|
| 黄体脂酮素
|
(PARI)/*A(x)=x+(p+q)*G((p*A(x)+q*x)/(p+q));G(x)=积分A(x)dx:*/
{a(n,p=2,q=-1)=my(a=x,G)
对于(n=1,30,打印1(a(n,2,-1),“,”)
(PARI)/*A(x-积分p*A(x)dx)=x+积分q*A(x)dx:*/
{a(n,p=2,q=-1)=my(a=[1],F=x);对于(i=1,n,a=concat(a,0);F=x*Ser(a);G=intformal(F);a[#a]=-polcoeff(subst(F,x,x-p*G)-q*G,#a));n!*a[n]}
对于(n=1,30,打印1(a(n,2,-1),“,”)
(PARI)/*生成前N个术语的非正式代码:*/
{N=20;p=2;q=-1;A=x;对于(i=1,N,G=整数(A+x*O(x^N));A=x+(p+q)*子集(G,x,(p*A+q*x)/(p+q));Vec(serlaplace(A))}
|