procedure forms(N)
   assert N in {105, 210, 420};
   eps := BaseExtend(KroneckerCharacter(105),Rationals()); 
   eps := DirichletGroup(N)!eps;
   S := NewSubspace(CuspidalSubspace(ModularForms([eps],2)));
   N := Newforms(S);
   for i in [1..#N] do
      f := N[i][1];
      CoefficientField(f);
      PowerSeries(f,100);
   end for;
end procedure;

print "Level 105";
forms(105);
print "Level 210";
forms(210);
print "Level 420";
forms(420);
