tmp.sagews

%load modular_symbol_map.pyx
def ms(E, sign=1):
    g = E.modular_symbol(sign=sign)
    h = ModularSymbolMap(g)
    d = float(h.denom)  # otherwise get int division!
    return lambda a,b: h._eval1(a,b)[0]/d

Compiling ./modular_symbol_map.pyx...

M = 1000000
for lbl in ['11a', '37a', '389a', '5077a']:
    name = 'plots/walk-%s-Z%s.pdf'%(lbl,M)
    if os.path.exists(name): continue
    print name
    s = ms(EllipticCurve(lbl))
    v = [s(a, M) for a in range(M)]
    stats.TimeSeries(v).sums().plot().save(name)

ms

Error in lines 1-1 Traceback (most recent call last): File "/projects/c6b42aa2-a189-445d-8bb7-880848c5b633/.sagemathcloud/sage_server.py", line 879, in execute exec compile(block+'\n', '', 'single') in namespace, locals File "", line 1, in <module> NameError: name 'ms' is not defined

s = ms(EllipticCurve('5077a'))