#Maple代码来自N.J.A.斯隆,2020年12月12日。对于相关序列也很有用。
#在Maple中获取标准θ级数。首先设置max degree,maxd。
最大值:=201:
#得到th2,th3,th4=Jacobiθ常数的最大值(参考Conway-Sloane,p.102)
温度0:=trunc(evalf(sqrt(maxd)))+2:
a: =0:对于从-temp0到temp0的i,执行a:=a+q^((i+1/2)^2):od:
th2:=系列(a,q,maxd);
a: =0:对于从-temp0到temp0的i,执行a:=a+q^(i^2):od:
th3:=系列(a,q,maxd);
th4:=系列(subs(q=-q,th3),q,maxd);
#获取Dedekind eta函数
a: =q^(1/24):对于从1到maxd的m,做a:=a*(1-q^m);日期:
eta:=a;
#得到phi0和phi1(参考Conway Sloane,第103页)
phi0:=系列(subs(q=q^2,th2)*subs(q=q^6,th2;
phi1:=系列(subs(q=q^2,th2)*subs(q=q^6,th3)+subs(q=q^2、th3)*sub(q=q ^6,th2),q,maxd);
#得到delta12(参考Conway Sloane,第204页,其中它被称为Delta_6)
增量12:=系列((subs(q=q^3,eta)*eta)^6,q,maxd);
增量12:=系列(子(q=q^2,增量12),q,最大值);
#要获得当前序列:(参考斯隆笔记本页面)
x: =phi0;y: =增量12;
w1:=x^12-72*x^6*y-1728*y^2;w1s:=系列(w1,q,maxd);w2:=子(q=sqrt(t),w1s);w3:=系列(w2,t,101);
|