higher-weight.sagews

for N in range(4, 50):
    M = ModularSymbols(Gamma0(N), 10, sign=-1)
    p = 2
    while N%p == 0: p = next_prime(p)
    print p, M.hecke_matrix(p).charpoly()
    T = M.hecke_matrix(p).change_ring(GF(2))
    if not T.is_invertible(): print N

3 x^3 + 84*x^2 - 46800*x - 5548608 4 2 x^3 + 18*x^2 - 904*x - 7872

Error in lines 1-7 Traceback (most recent call last): File "/projects/sage/sage-6.10/local/lib/python2.7/site-packages/smc_sagews/sage_server.py", line 957, in execute exec compile(block+'\n', '', 'single') in namespace, locals File "", line 6, in <module> File "sage/matrix/matrix_rational_dense.pyx", line 1351, in sage.matrix.matrix_rational_dense.Matrix_rational_dense.change_ring (/projects/sage/sage-6.10/src/build/cythonized/sage/matrix/matrix_rational_dense.c:15168) raise TypeError("matrix denominator not coprime to modulus") TypeError: matrix denominator not coprime to modulus

p = 2^30 - 35
t = magma(True)

for N in prime_range(2, 500):
    M = magma.ModularSymbols(N, 10, -1)
    T = M.HeckeOperator(2)
    Tp = magma.Matrix(GF(p), T)
    for i in [-45..45]: 
        if (Tp-i).IsInvertible() != t:
            print N, i