2015-05-17-104830-sato-tate-walk.sagews

E = EllipticCurve('11a')

%time
B = 100000
v = E.aplist(B)
w = prime_range(B)
s = [float(v[i]/sqrt(w[i])) for i in range(len(v))]
t = stats.TimeSeries(s)
t.plot_histogram(bins=100)

CPU time: 5.12 s, Wall time: 5.16 s

t.sums().plot()

t[:200].sums().plot()

t[:1000].sums().plot()