prod sato-tate

948 days ago by wstein

E = EllipticCurve('11a') F = EllipticCurve('37a') 
       
aE = finance.TimeSeries(E.aplist(10^6)) aF = finance.TimeSeries(F.aplist(10^6)) P = prime_range(10^6) 
       
time v = finance.TimeSeries([aE[i]*aF[i]/(4*P[i]) for i in range(len(aE))]) 
       
Time: CPU 0.82 s, Wall: 0.82 s
Time: CPU 0.82 s, Wall: 0.82 s
v.plot_histogram(bins=500)