%time B = 25 i = CDF(I) def f(z): return -i*zeta(z) save(complex_plot(f, (-B,B),(-B,B), frame=False, axes=False, plot_points=2400), '7.png', figsize=[8,8])