|
t=15；
y[n_]：=（1+总和[下标[p，k]/n^k，{k，1，t}]）；
mul=1；开始=9；如果[t>9，mul=n^（t-9）；start=t]；
w=分开[Expand[mul*Simplify[
y[n]*n*（n-1）*（n-2）*（n-3）*
-（（3*n-30）*y[n-11]
+（6*n-45）*y[n-10]*（n-10）
+（5*n+18）*y[n-9]*（n-9）*（n-10）
-（8*n-139）*y[n-8]*（n-8）*（n-9）*（n-10）
-（26*n-204）*y[n-7]*（n-7）*（n-8）*（n-9）*（n-10）
-（4*n-30）*y[n-6]*（n-6）*（n-7）*
+（26*n-148）*y[n-5]*（n-5）*（n-6）*
+（8*n-74）*y[n-4]*（n-4）*
-（9*n-18）*y[n-3]*（n-3）*（n-4）*
-（2*n-15）*y[n-2]*（n-2）*（n-3）*
+（n+2）*y[n-1]*（n-1）*（n-2）*；
sol=求解[表[系数[w，n，j]==0，{j，start，start-t+1，-1}]]；
asympt=y[n]/.sol[[1]；
表[分母[系数[asympt，n，-j]]，{j，0，t}](*瓦茨拉夫·科特索维奇2012年4月6日*）
| 