E = EllipticCurve('389a')
f = CDF['x'](E.anlist(100))
h0(x) = exp(2*pi*I*abs(f(x)))
h0 = fast_callable(h0, CDF)
h = lambda x : h0(x) if x.abs() < 1 else 1e999
complex_plot(h, (-1,1), (-1, 1), plot_points=600, axes=False, figsize=[8,8])