# Riemann–Siegel formula

In mathematics, the Riemann–Siegel formula is an asymptotic formula for the error of the approximate functional equation of the Riemann zeta function, an approximation of the zeta function by a sum of two finite Dirichlet series. It was found by Siegel (1932) in unpublished manuscripts of Bernhard Riemann dating from the 1850s. Siegel derived it from the Riemann–Siegel integral formula, an expression for the zeta function involving contour integrals. It is often used to compute values of the Riemann–Siegel formula, sometimes in combination with the Odlyzko–Schönhage algorithm which speeds it up considerably. When used along the critical line, it is often useful to use it in a form where it becomes a formula for the Z function.

If M and N are non-negative integers, then the zeta function is equal to

$\zeta (s)=\sum _{n=1}^{N}{\frac {1}{n^{s}}}+\gamma (1-s)\sum _{n=1}^{M}{\frac {1}{n^{1-s}}}+R(s)$ where

$\gamma (s)=\pi ^{{\tfrac {1}{2}}-s}{\frac {\Gamma \left({\tfrac {s}{2}}\right)}{\Gamma \left({\tfrac {1}{2}}(1-s)\right)}}$ is the factor appearing in the functional equation ζ(s) = γ(1 − s) ζ(1 − s), and

$R(s)={\frac {-\Gamma (1-s)}{2\pi i}}\int {\frac {(-x)^{s-1}e^{-Nx}}{e^{x}-1}}dx$ is a contour integral whose contour starts and ends at +∞ and circles the singularities of absolute value at most 2πM. The approximate functional equation gives an estimate for the size of the error term. Siegel (1932) and Edwards (1974) derive the Riemann–Siegel formula from this by applying the method of steepest descent to this integral to give an asymptotic expansion for the error term R(s) as a series of negative powers of Im(s). In applications s is usually on the critical line, and the positive integers M and N are chosen to be about (2πIm(s))1/2. Gabcke (1979) found good bounds for the error of the Riemann–Siegel formula.

## Riemann's integral formula

Riemann showed that

$\int _{0\searrow 1}{\frac {e^{-i\pi u^{2}+2\pi ipu}}{e^{\pi iu}-e^{-\pi iu}}}\,du={\frac {e^{i\pi p^{2}}-e^{i\pi p}}{e^{i\pi p}-e^{-i\pi p}}}$ where the contour of integration is a line of slope −1 passing between 0 and 1 (Edwards 1974, 7.9).

He used this to give the following integral formula for the zeta function:

$\pi ^{-{\tfrac {s}{2}}}\Gamma \left({\tfrac {s}{2}}\right)\zeta (s)=\pi ^{-{\tfrac {s}{2}}}\Gamma \left({\tfrac {s}{2}}\right)\int _{0\swarrow 1}{\frac {x^{-s}e^{\pi ix^{2}}}{e^{\pi ix}-e^{-\pi ix}}}\,dx+\pi ^{-{\frac {1-s}{2}}}\Gamma \left({\tfrac {1-s}{2}}\right)\int _{0\searrow 1}{\frac {x^{s-1}e^{-\pi ix^{2}}}{e^{\pi ix}-e^{-\pi ix}}}\,dx$ 