for N in [24..100]: S = CuspForms(N,2) print N print S.newforms('a')
S = CuspForms(41,2) D = S.newforms('a') D
f = D[0]
K = f.hecke_eigenvalue_field(); K
K.signature()
proof.number_field(False) %time U = CyclotomicField(29).unit_group()