N := 2^6*67;
bound := 4000;


   F<alp> := GF(25);
   G<a2,b2,c>:=DirichletGroup(N, F);
   assert Order(a2) eq 2;
   assert Order(b2) mod 2 eq 0;
   assert Order(c) mod 2 eq 0;
   eps := b2^(Order(b2) div 2)*(c^(Order(c) div 2));

   time M := ModularSymbols(eps,5,+1);
   "Created M = ", M;

   R<x>:= PolynomialRing(F);
   I := [<19,x>, <23,x>, <47,x>, <3, x+2>, <7, x+alp^8>];
   time V := Kernel(I, M);
   "Created V = ", V;
   save "4288.space";

   if Dimension(V) ne 2 then
      error "Dimension of V is not 2.";
   end if;
   
printf "Computing T_2, T_3, T_5, ..., T_%o.\n", PreviousPrime(bound+1);
time T := [DualHeckeOperator(V,p) : p in [2..bound] | IsPrime(p)];

file := Open(IntegerToString(N) cat ".Tp", "w");
fprintf file, "F<alp> := GF(%o);\n", #F;
fprintf file, "T%o := \n %m;\n", N, T;

print "All done!";
