Mathematics Stack Exchange is a question and answer site for people studying math at any level and professionals in related fields. Join them; it only takes a minute:

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

Over at stack overflow I asked a question about how to plot the Riemann zeta zero spectrum from the von Mangoldt function. Then I asked a question about calculating the Riemann zeta function at the critical line with the Fourier transform of a exponential sawtooth function by building upon Heike's answer. After experimenting with anon's answer I believe I plotted the Riemann zeta function at the critical line here at my blog times/divided by a scale factor. By then going back to the beginning I modified the transform as in anon's answer and applied it to the von Mangoldt function and got a zeta zero spectrum plot which appears to have spikes of equal amplitude.

The Mathematica program for the plot looks like this:

scale = 1000000;
f = Range[scale];

f[[1]] = N@MangoldtLambda[1];
  f[[i]] = N@MangoldtLambda[i] + f[[i - 1]], {i, 2, scale}], i]

xres = .002;
xlist = Exp[Range[0, Log[scale], xres]];
tmax = 60;
tres = .015;
s = 1/2;
Monitor[errList1 = 
   Table[((xlist^(-s + I t).(f[[Floor[xlist]]] - xlist)))*(s + 
       I t), {t, Range[0, 60, tres]}];, t]

Print["Variant of the Fourier transform of the von Mangoldt function"]
g1 = ListLinePlot[Re[errList1]/Length[xlist], DataRange -> {0, 60}, 
   PlotRange -> {-.3, 1.3}, Axes -> True, Filling -> Axis];
g2 = Graphics[Line[{{0, 1}, {60, 1}}]];
Show[g1, g2, ImageSize -> Large]

And the output looks like this:

Riemann zeta zero spectrum from von Mangoldt function equal amplitudes

My question is: Are the amplitudes of the frequency spikes equal to $1$ (the upper black line) when the real part of the complex number $s$ is equal to $1/2$?

The matrix formulation of this transform, which is a variant of the Fourier transform, usually gives an a bit un-even output so I guess the answer could be positive.

Edit 23.7.2012:

The program I believe does something like this:

$$\displaystyle X_t = \frac{1}{length(x)} \sum\limits_{x=0}^{x=\log(scale-1)} (e^{x})^{-1/2 + i t} \cdot (f(\lfloor e^{x}\rfloor)-e^{x})) \cdot (-1/2 + i t)$$

This second plot is the same as from the program above except that here the multiplicand in the program was $(-1/2 + i t)$ instead of $(+1/2 + i t)$. This gives a different beginning of the plot.

enter image description here

where $\displaystyle X_t $ should be the graph in the plot.

Edit 23.7.2012:

In the formula above $f$ should be:

$$\displaystyle f = \sum\limits_{n=1}^{n=k} \Lambda(n)$$

where $\Lambda(n)$ is the von Mangoldt function:

Log(1), Log(2), Log(3), Log(2), Log(5), Log(1) ...

share|cite|improve this question
Could you write down an equation for what's being graphed? – anon Jul 14 '12 at 3:19
@anon: I tried to write down in latex what the program does. – Mats Granvik Jul 23 '12 at 14:23

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.

Browse other questions tagged or ask your own question.