|
q:=x-4+I*（（x+4）*（16-x））^（1/2）；
f:=x*（q/8）^4；
s：=（（q-2）/（8*I-6））^（1/4）；
y1:=超几何（[1/8，1/8]，[3/4]，f）*s/x^（1/8）；
r:=2/（x*（（x+4）*（16-x））^（1/2）*y1^2）；
h:=表皮（[1/2，1/2，1/2，1/2]，[1，1，1]，256*x）；
u：=（15*（223*x+72）*x^2*diff（h，x，x）+（14579*x+3226）*x*diff
+（9969*x+1002）*diff（h，x）+320*h）/（16*（16-x）*（x+4）*x^2）；
ogf：=y1^2*整数（r*（1+整数（r*Int（u/（r*y1）^2，x），x））；
#通过计算级数展开来检查o.g.f.：
SER:=程序（a，x）系列（a，x，20）结束：
INT:=进程（a，x）INT（SER（a，x），x）结束：
SER（评估（ogf，Int=Int），x）#马克·范·霍伊2013年4月2日
| 