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)*(n-3)*(n-4)*(n-5)*(n-6)*(n-7)*(n-8)*(n-9)*(n-10))],n],n];
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日*)