plot(im(zeta (x))) from -42 to -40