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)
t.sums().plot()
t[:200].sums().plot()
t[:1000].sums().plot()