R_k.sagews

%load code.sage

R = RiemannPiApproximation(10, 100); R

Riemann explicit formula for pi(x) for x <= 100 using R_k for k <= 10

plot(lambda x : R.Tk(x,1), (10,100), plot_points=100)

R = RiemannPiApproximation(100, 500, 50)

%time g = [R.plot_Rk(k, 3, 100, plot_points=100, ymin=0, ymax=25) for k in [1,..,50]]

CPU time: 17.22 s, Wall time: 17.45 s

%time
animate(g)

CPU time: 14.06 s, Wall time: 20.05 s

%time g = [R.plot_Rk(k, 3, 100, plot_points=200, ymin=0, ymax=25) for k in [1,..,100]]

CPU time: 38.13 s, Wall time: 38.91 s

%time
animate(g)

CPU time: 27.02 s, Wall time: 39.10 s