zeta_circle.sagews

v = zeta_zeros()

len(v)

100000

def rmod(a,b):
    return a -  a//b *b

rmod(7,float(pi))

0.7168146928204138

def f(per):
    per = float(per)
    w = [rmod(a,per) for a in v]
    show(stats.TimeSeries(w).plot_histogram())

f(2*pi)

f(1)

def circ():
    per = float(2*pi)
    w = [rmod(a,per) for a in v]
    g = points([(cos(a),sin(a)) for a in w])
    show(g, aspect_ratio=1)

circ()