（PARI）
N=13；F=7；gmax（n）=n\2；
Q=矩阵（N+1，N+1）；
Qget（n，g）={如果（g<0||g>n/2，0，Q[n+1，g+1）}；
Qset（n，g，v）={Q[n+1，g+1]=v}；
二次（{x=1}）={
Qset（0，0，x）；
对于（n=1，长度（Q）-1，对于（g=0，gmax（n），
my（t1=（1+x）*（2*n-1）/3*Qget（n-1，g），
t2=（2*n-3）*（2*n-2）*（2*n-1）/12*Qget（n-2，g-1），
t3=1/2*总和（k=1，n-1，总和（i=0，g，
（2*k-1）*（2*（n-k）-1）*Qget（k-1，i）*Qget（n-k-1，g-i）））；
Qset（n，g，（t1+t2+t3）*6/（n+1）））；
};
二次（'x+O（'x^（F+1）））；
v=向量（N+2-F，N，向量（1+gmax（N-1），g，polcoeff（Qget（N+F-2，g-1），F））；
凹面（v）
