Formule de Riemann-Siegel

En mathématiques, et plus précisément en analyse, la formule de Riemann-Siegel est une estimation asymptotique de l'erreur de l'équation fonctionnelle d'approximation de la fonction zêta de Riemann, c'est-à-dire une approximation de la fonction zêta par la somme de séries de Dirichlet finies.

Histoire

Cette formule a été découverte par Carl Siegel en 1932 alors qu'il analysait les notes manuscrites non publiées de Bernhard Riemann, lesquelles dataient des années 1850.

Siegel l'a obtenue en partant de l' « intégrale de Riemann–Siegel », une expression de la fonction zêta qui fait appel à une intégrale curviligne. Au début du XXI^e siècle, cette formule est régulièrement utilisée pour calculer les valeurs de la fonction zêta, dans le but de rechercher la position de ses zéros (voir hypothèse de Riemann). Elle est parfois combinée à l'algorithme de Odlyzko–Schönhage, ce qui rend le calcul plus rapide.

Description

Si M et N sont des entiers strictement positifs, alors la fonction zêta peut être définie par

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

où

$\displaystyle \gamma (s)=\pi ^{1/2-s}\Gamma (s/2)/\Gamma ((1-s)/2)$

est le facteur qui apparaît dans l'équation fonctionnelle $\displaystyle \zeta (s)=\gamma (s)\zeta (1-s)$ et où

$R(s)={\frac {-\Gamma (1-s)}{2\pi i}}\oint {\frac {(-x)^{s-1}e^{-Nx}dx}{e^{x}-1}}$

est une intégrale curviligne sur une ligne commençant et se terminant en $\displaystyle \infty$ et encerclant les pôles de module inférieur ou égal à $2\pi M$.

L'équation fonctionnelle d'approximation produit une estimation de la taille de l'erreur. Siegel en 1932 et Edwards en 1974 ont obtenu la formule de Riemann-Siegel en appliquant la méthode du point col à cette intégrale pour en déduire un développement asymptotique de l'erreur $R(s)$ en série de puissances négatives de $\scriptstyle \Im (s)$.

Dans la pratique, s est habituellement sur la ligne critique et les entiers positifs M et N ont comme valeur $\scriptstyle 2\pi \Im (s)^{1/2}$. En 1979, Gabcke a découvert de bonnes limites pour l'erreur de la formule de Riemann–Siegel.

Formule intégrale de Riemann

Riemann a démontré que

$\int _{0\searrow 1}{\frac {e^{-i\pi u^{2}+2\pi ipu}}{e^{\pi iu}-e^{-\pi iu}}}~\mathrm {d} u={\frac {e^{i\pi p^{2}}-e^{i\pi p}}{e^{i\pi p}-e^{-i\pi p}}}$

où la courbe d'intégration est une droite de pente −1 qui passe entre 0 et 1.

Il a utilisé cette équation pour donner une formule intégrale de la fonction zêta :

$\displaystyle \pi ^{-s/2}\Gamma (s/2)\zeta (s)=$

$\pi ^{-s/2}\Gamma (s/2)\int _{0\swarrow 1}{\frac {x^{-s}e^{\pi ix^{2}}}{e^{\pi ix}-e^{-\pi ix}}}~\mathrm {d} x+\pi ^{-(1-s)/2}\Gamma ((1-s)/2)\int _{0\searrow 1}{\frac {x^{s-1}e^{-\pi ix^{2}}}{e^{\pi ix}-e^{-\pi ix}}}~\mathrm {d} x.$

Notes et références

(en) Cet article est partiellement ou en totalité issu de l’article de Wikipédia en anglais intitulé « Riemann–Siegel formula » (voir la liste des auteurs).

Annexes

Bibliographie

• (en) Michael V. Berry, « The Riemann–Siegel expansion for the zeta function : high orders and remainders », Proceedings of the Royal Society. London. Series A. Mathematical, Physical and Engineering Sciences, vol. 450,‎ 1995, p. 439–462 (ISSN 0962-8444, DOI 10.1098/rspa.1995.0093)
• (en) John Derbyshire (en), Prime Obsession (en) : Bernhard Riemann and the Greatest Unsolved Problem in Mathematics, Plume, 2004, 448 p. (ISBN 978-0-452-28525-5)
• (en) H. M. Edwards, Riemann's Zeta Function, New York, Dover Publications, 1974, 315 p. (ISBN 978-0-486-41740-0, lire en ligne)
• (de) Wolfgang Gabcke, Neue Herleitung und Explizite Restabschätzung der Riemann-Siegel-Formel, Georg-August-Universität Göttingen, 1979 (lire en ligne)
• (de) Carl Ludwig Siegel, Über Riemanns Nachlaß zur analytischen Zahlentheorie, 1932, 80 p. (lire en ligne), p. 45–80
  Republié in Gesammelte Abhandlungen, vol. 1., Berlin: Springer-Verlag, 1966

Liens externes

• (en) X. Gourdon, « Numerical evaluation of the Riemann Zeta-function »
• (en) Eric W. Weisstein, « Riemann–Siegel Formula », sur MathWorld