def xy_to_zeta_size(x,y):

return abs(zeta(N(x+I*y)))

cvals = [e^i for i in srange(-7,1,.25)]

cp = contour_plot(xy_to_zeta_size,(-6,3),(-3,3),contours = cvals, fill=False, plot_points = 201)

rzeta(z) = zeta(z)/norm(zeta(z))^(.25)

rzf = fast_callable(rzeta,domain=CDF)

cparg = complex_plot(rzf,(-6,3),(-3,3))

show(cparg+cp,figsize=[18,12])

## Monday, March 21, 2011

### Plotting the zeta function

I lectured a tiny bit on the Riemann zeta function for the first time in my complex analysis course, which inspired me to make the following plot. Colors are the argument of the Riemann zeta function, brightness is proportional to magnitude. The default brightness map gives very low contrast, so I modified the magnitudes. To help with seeing the magnitudes a contour map with exponentially spaced contours is overlaid.

Subscribe to:
Post Comments (Atom)

## 1 comment:

Nice picture!

I hope you don't mind, but I stole it and turned it into a stackoverflow question. http://stackoverflow.com/q/5385103/421225

Post a Comment