Sharedwww / tables / artin / 3556.mOpen in CoCalc
Author: William A. Stein
1
N := 2^2*7*127;
2
bound := 4000;
3
4
5
F<alp> := GF(25);
6
G<a,b,c>:=DirichletGroup(N, F);
7
assert Order(a) eq 2;
8
assert Order(b) mod 2 eq 0;
9
assert Order(c) mod 3 eq 0;
10
eps := b^(Order(b) div 2) * (c^(Order(c) div 3));
11
assert Order(eps) eq 6;
12
13
time M := ModularSymbols(eps,5,+1);
14
"Created M = ", M;
15
16
R<x>:= PolynomialRing(F);
17
I := [<19,x>, <29,x>, <89,x>, <3, x+alp^4>];
18
time V := Kernel(I, M);
19
"Created V = ", V;
20
save "3556.space";
21
22
if Dimension(V) ne 2 then
23
error "Dimension of V is not 2.";
24
end if;
25
26
printf "Computing T_2, T_3, T_5, ..., T_%o.\n", PreviousPrime(bound+1);
27
time T := [DualHeckeOperator(V,p) : p in [2..bound] | IsPrime(p)];
28
29
file := Open(IntegerToString(N) cat ".Tp", "w");
30
fprintf file, "F<alp> := GF(%o);\n", #F;
31
fprintf file, "T%o := \n %m;\n", N, T;
32
33
print "All done!";
34