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