|
PROG
|
(Magma)
prec := 80;
S := SymmetricMatrix([4, 2, 4, 0, -2, 4, 0, -2, 0, 4, 0, 0, -2, 0, 4, -2, -2, 0, 0, 0, 4, -2, -1, 1, 0, 0, 0, 4, -2, -1, 0, -1, 1, 2, 2, 4, -2, -2, 0, 1, 1, 2, 2, 2, 4, -2, 0, -2, 0, 1, 1, 0, 0, 0, 4, 1, 1, 0, 0, 0, -2, 0, -1, -1, -2, 4, -2, -1, 0, 0, 0, 1, 1, 1, 1, 1, -2, 4, 0, -1, 1, 1, 0, -1, 1, 0, 0, -1, 1, -1, 4, 0, 0, 0, 0, 0, 0, 1, 0, 1, -1, 1, -1, 1, 4]);
L := LatticeWithGram(S);
T<q> := ThetaSeries(L, 42);
M := ThetaSeriesModularFormSpace(L);
B := Basis(M, prec);
Coefficients(&+[Coefficients(T)[2*i-1]*B[i] :i in [1..22]]); // Andy Huchala, May 07 2023
|