7.sagews

%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])

CPU time: 2532.82 s, Wall time: 2774.58 s