1239.sagews

a = [[(i/10).n(), (pari(1/2 + i*I/10).zeta().abs().precision(1))] for i in range(0,800)]

p = line(a)

p.save("zetacritical.png")

b = [[(pari(1/2 + i*I/10).zeta().real()).precision(1),(pari(1/2 + i*I/10).zeta().imag()).precision(1)] for i in range (0,800)]

q = line(b)

q.save("zetaparamplot.png")