Introduction

On s’intéressera ici à la fonction zêta de Riemann, définie sur $\left] {1; + \infty } \right[$ par :
$$\zeta \left( x \right) =  \sum _{n = 1}^{ + \infty } \frac{1}{n^x}$$
ainsi qu’à plusieurs autres fonctions apparentées.

Ces dernières sont d’abord les deux fonctions $\alpha $ et $\beta $ directement liées à $\zeta $ par simples regroupements de termes, respectivement définies sur $\left] {1; + \infty } \right[$ par :

$$\alpha ( x ) = \sum _{m = 1}^{ + \infty } \frac{1}{(2m)^x}\;\;\;\;\;;\;\;\;\;\;\beta ( x ) =  \sum _{m = 0}^{ + \infty } \frac{1}{( 2m + 1 )^x}$$

Les fonctions apparentées considérées sont aussi celles correspondant aux séries alternées qui sont naturellement associées aux trois séries précédentes, respectivement définies sur $\left] {0; + \infty } \right[$ par :

$$\xi ( x ) =  \sum _{n = 1}^{ + \infty } \frac{(  - 1 )^{n - 1}}{n^x}\; ; \;\varphi ( x ) =  \sum _{m = 1}^{ + \infty } \frac{( - 1)^{m - 1}}{( 2m )^x}\; ; \;\psi ( x ) =  \sum _{m = 0}^{ + \infty } \frac{( - 1)^m}{( 2m + 1)^x}$$

Le problème considéré ici est le calcul des images des entiers naturels par ces six fonctions.

Il a suscité une vaste littérature depuis les travaux initiés à son sujet au XVIIIe siècle. La connaissance de $\psi \left( 1 \right)$, égal à $\frac{\pi}{4}$, est même très antérieure puisqu’on peut non sans justifications la considérer comme datant au plus tard du début du XVe siècle : on la trouve sous une forme ne faisant certes pas l’usage de l’infini actuel dans les travaux de l’astronome indien Mādhava de Sangamagrama1.

Le mathématicien et philosophe allemand Gottfried W. Leibniz fut le premier Européen à l’établir, en 1673. Le problème fameux de la valeur de $\zeta \left( 2 \right)$, connu sous le nom de « problème de Bâle », avait été posé dès 1644 par le mathématicien italien Pietro Mengoli. Il fallut attendre 1735 pour le voir résolu par le mathématicien suisse Leonhard Euler : il trouva la valeur $\frac{\pi ^2}{6}$ par une méthode ingénieuse —mais logiquement incorrecte— consistant à appliquer au développement en série entière de $\frac{\sin \left( x \right)}{x}$ un argument sur les relations entre racines et coefficients des polynômes, et détermina les valeurs de $\zeta \left( 4 \right)$, $\zeta \left( 6 \right)$, $\zeta \left( 8 \right)$, $\zeta \left( {10} \right)$ et $\zeta \left( {12} \right)$ par la même méthode2.

Les progrès de l’analyse, surtout la théorie des séries de Fourier, ont permis au XIXe siècle de trouver les valeurs des cinq premières fonctions en les entiers pairs, et celles de la dernière aux entiers impairs. Mais les valeurs des cinq premières fonctions aux entiers impairs et celles de la sixième aux entiers pairs demeurent éminemment problématiques. Bien que plusieurs découvertes importantes à ce sujet aient été faites au cours des dernières décennies3, le problème que posent ces valeurs demeure largement ouvert et fait toujours l’objet d’actives recherches.

Une approche relativement élémentaire du problème est présentée ici, qui le met en relation étroite avec quatre thèmes de parité : les notions de parité d’une fonction et de parité du degré d’un polynôme rejoignent ici les distinctions de parité concernant tant l’argument entier des six fonctions considérées que les entiers dont certaines puissances inverses sont sommées.

La méthode adoptée a essentiellement pour but de familiariser avec le problème considéré par une démarche simple.

Prérequis : les deux théorèmes de Dirichlet sur les séries de Fourier

Tous les résultats de cet article découlent de deux propositions dont les démonstrations reposent sur l’utilisation de deux théorèmes qui seront ici dits « de Dirichlet », même si historiquement, seul le premier peut être attribué au mathématicien allemand Johann P.G. Lejeune Dirichlet, qui l’énonça et le démontra partiellement en 1829.

Ils seront utilisés dans la version suivante :

Images des entiers naturels pairs non nuls par $\zeta$, $\xi $ et $\varphi$.

Dans l’approche du problème proposée ici, le calcul des images des entiers naturels pairs non nuls par $\zeta$, $\xi $ et $\varphi $ repose sur la proposition suivante :

Démonstration :

Soit $F$ la fonction $2\pi $-périodique et paire qui coïncide avec $f$ sur $\left[ {0;\pi } \right]$. $F$ est continue et de classe ${C^1}$ par morceaux. D’après le théorème de Dirichlet, la série de Fourier de $F$ converge donc uniformément vers $F$. On en déduit que pour tout $t \in \left[ {0;\pi } \right]$ :

$$f\left( t \right) = \frac{1}{\pi } \int _0^\pi  f\left( t \right)\;dt + \mathop \sum \limits_{n = 1}^{ + \infty } {a_n}{\rm{cos}}\left( {nt} \right).$$

où :

$$\forall n\in {\mathbb N}^*, \; a_n=\frac{2}{\pi}\int_0^\pi f(t)\cos(nt)dt.$$

Par conséquent :

$$f\left( 0 \right) = \frac{1}{\pi }\int _0^\pi  f\left( t \right)\;dt + \mathop \sum \limits_{n = 1}^{ + \infty } {a_n}.$$

Or, pour tout $n\in {\mathbb N}^*$ fixé, on établit par récurrence sur $j\in {\mathbb N}^*$ que :

$$\begin{split} {a_n} = & \frac{2}{\pi } \sum _{k = 1}^j ( - 1)^{k - 1} \frac{( - 1 )^n {f^{(2k - 1 )}}{ } {( \pi )} - f^{(2k - 1 )} (0)}{n^{2k}} \\ & + \frac{2( - 1 )^j}{\pi n^{2j}}\int _0^\pi  f^{( 2j )}( t )\cos ( nt )dt.\end{split}$$

En effet, deux intégrations par parties successives permettent d’obtenir :

$$\begin{split} {a_n} = & - \frac{2}{{\pi n}} \int _0^\pi  f'\left( t \right)\sin \left( {nt} \right)dt = \frac{2}{{\pi {n^2}}}\left( {{{\left( { - 1} \right)}^n}f'\left( \pi  \right) - f'\left( 0 \right)} \right) \\ & - \frac{2}{{\pi {n^2}}} \int _0^\pi  f''\left( t \right)\cos \left( {nt} \right)dt.\end{split}$$

Ce qui montre que l’identité est satisfaite pour $j = 1$. Supposons que l’identité est satisfaite pour un certain $j\in {\mathbb N}^*$. On obtient en intégrant par parties deux fois que :

$$ \int _0^\pi  {f^{\left( {2j} \right)}}\left( t \right)\cos \left( {nt} \right)dt   =  - \frac{1}{n} \int _0^\pi  {f^{\left( {2j + 1} \right)}}\left( t \right)\sin \left( {nt} \right)dt $$ $$  = \frac{1}{n^2}\left( {{{\left( { - 1} \right)}^n}{f^{\left( {2j + 1} \right)}}\left( \pi  \right) - {f^{\left( {2j + 1} \right)}}\left( 0 \right)} \right) - \frac{1}{{{n^2}}} \int _0^\pi  {f^{\left( {2j + 2} \right)}}\left( t \right)\cos \left( {nt} \right)dt.$$

En utilisant l’hypothèse de récurrence, on en déduit que l’identité est satisfaite au rang $j + 1$. L’identité annoncée, vraie par récurrence pour tout $j\in {\mathbb N}^*$, est donc vraie en particulier en $j = p$ ; c’est-à-dire que :

$$\begin{split} {a_n} = & \frac{2}{\pi }\mathop \sum \limits_{k = 1}^p {\left( { - 1} \right)^{k - 1}}\frac{{{{\left( { - 1} \right)}^n}{f^{\left( {2k - 1} \right)}}\left( \pi  \right) - {f^{\left( {2k - 1} \right)}}\left( 0 \right)}}{{{n^{2k}}}} \\ & + \frac{{2{{\left( { - 1} \right)}^p}}}{{\pi {n^{2p}}}} \int _0^\pi  {f^{\left( {2p} \right)}}\left( t \right)\cos \left( {nt} \right)dt.\end{split}$$

Or, $f$ étant polynomiale de degré inférieur ou égal à $2p$, la fonction ${f^{\left( {2p} \right)}}$ est constante sur $\left[ {0;\pi } \right]$. On en déduit que :

$$ \int _0^\pi  {f^{\left( {2p} \right)}}\left( t \right)\cos \left( {nt} \right)dt = 0.$$

D’où résulte :

$${a_n} = \frac{2}{\pi }\mathop \sum \limits_{k = 1}^p {\left( { - 1} \right)^{k - 1}}\frac{{{{\left( { - 1} \right)}^n}{f^{\left( {2k - 1} \right)}}\left( \pi  \right) - {f^{\left( {2k - 1} \right)}}\left( 0 \right)}}{{{n^{2k}}}}.$$

On obtient donc finalement :

$$ \begin{split} f\left( 0 \right)  &=  \frac{1}{\pi } \int _0^\pi  f\left( t \right)\;dt \\ & + \mathop \sum \limits_{n = 1}^{ + \infty } \frac{2}{\pi }\mathop \sum \limits_{k = 1}^p {\left( { - 1} \right)^{k - 1}}\frac{{{{\left( { - 1} \right)}^n}{f^{\left( {2k - 1} \right)}}\left( \pi  \right) - {f^{\left( {2k - 1} \right)}}\left( 0 \right)}}{{{n^{2k}}}} \\ & =  \frac{1}{\pi }\int _0^\pi  f\left( t \right)\;dt \\ & - \frac{2}{\pi }\mathop \sum \limits_{k = 1}^p {\left( { - 1} \right)^{k - 1}}\mathop \sum \limits_{n = 1}^{ + \infty } \left( {\frac{{{{\left( { - 1} \right)}^{n - 1}}}}{{{n^{2k}}}}{f^{\left( {2k - 1} \right)}}\left( \pi  \right) + \frac{1}{{{n^{2k}}}}{f^{\left( {2k - 1} \right)}}\left( 0 \right)} \right).\end{split}$$

Et cette identité n’est autre que celle annoncée dans la proposition 1.$\square$

Pour obtenir la valeur de $\xi \left( {2p} \right)$ et $\zeta \left( {2p} \right)$ pour tout $p\in \mathbb N^*$, il suffit alors de choisir convenablement la fonction $f$ considérée dans la proposition 1.

Démonstration :

Soit $p\in\mathbb N^*$. On applique la proposition 1 avec $f\;:t \mapsto {t^{2p}}$. On obtient alors :

$$0 = \frac{1}{\pi }\frac{{{\pi ^{2p + 1}}}}{{2p + 1}} - \frac{2}{\pi }\mathop \sum \limits_{k = 1}^p {\left( { - 1} \right)^{k - 1}}\frac{{\left( {2p} \right)!}}{{\left( {2p + 1 - 2k} \right)!}}{\pi ^{\left( {2p + 1} \right) - 2k}}\;\xi \left( {2k} \right).$$

Soit aussi :

$$\mathop \sum \limits_{k = 1}^p {\left( { - 1} \right)^{k - 1}}\frac{{\left( {2p} \right)!}}{{\left( {2p + 1 - 2k} \right)!}}\;\frac{1}{{{\pi ^{2k}}}}\;\xi \left( {2k} \right) = \frac{1}{{2\left( {2p + 1} \right)}}$$

En posant ${A_q} = \frac{\xi \left( {2q} \right)}{\pi ^{2q}}$ pour tout $q\in {\mathbb N}^*$, on obtient la relation exprimant ${A_p}$ lorsque $p \ge 2$.

Dans le cas où $p = 1$, on obtient de plus avec la relation précédemment obtenue que $\xi \left( 2 \right) = \frac{{\pi ^2}}{12}$, donc que $\frac{\xi (2)}{{\pi^2}}=A_1=\frac{1}{12} \in\mathbb Q$.

Supposons alors que, pour un $p\in {\mathbb N}^*$ fixé, $A_j \in \mathbb Q$ pour tout $j \in [\![ 1;p ]\!] $. Il est clair que

$$\begin{split} {{A}_{p+1}}& =\frac{\xi \left( 2p+2 \right)}{{{\pi }^{2p+2}}} \\ & = \frac{{{\left( -1 \right)}^{p}}}{\left( 2p+2 \right)!}\left( \frac{1}{2\left( 2p+3 \right)}-\underset{k=1}{\overset{p}{\mathop \sum }}{{\left( -1 \right)}^{k-1}}\frac{\left( 2p+2 \right)!}{\left( 2p+3-2k \right)!}{{A}_{k}} \right)\\ & \in \mathbb{Q} \end{split}$$

Par récurrence forte sur $p$, on en déduit que $A_p\in\mathbb Q$ pour tout $p\in {\mathbb N}^*$. $\square$

Démonstration :

Il suffit d’observer que :

$$\varphi \left( {2p} \right) = \mathop \sum \limits_{m = 1}^{ + \infty } \frac{{{{\left( { - 1} \right)}^{m - 1}}}}{{{{\left( {2m} \right)}^{2p}}}} = \frac{1}{{{4^p}}}\mathop \sum \limits_{m = 1}^{ + \infty } \frac{{{{\left( { - 1} \right)}^{m - 1}}}}{{{m^{2p}}}} = \frac{1}{{{4^p}}}\xi \left( {2p} \right). \; {\color{red}\square}$$

Démonstration :

Soit $p\in {\mathbb N}^*$. On applique encore la proposition 1, mais avec cette fois-ci la fonction $f$ définie sur $\left[ {0;\pi } \right]$ par : $f\;\left( t \right) = {\left( {t - \pi } \right)^{2p}}$. On obtient alors :

$$\begin{split} {\pi ^{2p}} = & \frac{1}{\pi } \frac{{{\pi ^{2p + 1}}}}{{2p + 1}} \\ & - \frac{2}{\pi }\mathop \sum \limits_{k = 1}^p {\left( { - 1} \right)^{k - 1}}\frac{{\left( {2p} \right)!}}{{\left( {2p + 1 - 2k} \right)!}} \; ( - 1 )^{2p - 2k + 1} \; \; \pi ^{2p + 1 - 2k} \; { }\;  \zeta \left( {2k} \right).\end{split}$$

Soit aussi :

$$\frac{p}{{2p + 1}} = \mathop \sum \limits_{k = 1}^p {\left( { - 1} \right)^{k - 1}}\frac{{\left( {2p} \right)!}}{{\left( {2p + 1 - 2k} \right)!}}\;\frac{1}{{{\pi ^{2k}}}}\;\zeta \left( {2k} \right).$$

En posant ${B_q} = \frac{\zeta \left( {2q} \right)}{{\pi ^{2q}}}$ pour tout $q\in {\mathbb N}^*$, on obtient la relation exprimant ${B_p}$ lorsque $p \ge 2$.

Dans le cas où $p = 1$, on obtient de plus avec la relation précédemment obtenue que $\zeta \left( 2 \right) = \frac{\pi ^2}{6}$, donc que $\frac{\zeta (2)}{\pi^2}=B_1=\frac{1}{6}\in \mathbb Q$.

Supposons alors que, pour un $p\in {\mathbb N}^*$ fixé,  $B_j\in \mathbb Q$ pour tout $j \in [\![ 1;p ]\!]$. Il est clair que

$$\begin{split} {{B}_{p+1}}& =\frac{\zeta \left( 2p+2 \right)}{{{\pi }^{2p+2}}} \\ & =\frac{{{\left( -1 \right)}^{p}}}{\left( 2p+2 \right)!}\left( \frac{p+1}{2p+3}-\underset{k=1}{\overset{p}{\mathop \sum }}\,{{\left( -1 \right)}^{k-1}}\frac{\left( 2p+2 \right)!}{\left( 2p+3-2k \right)!}~{{B}_{k}} \right)\in \mathbb{Q}\end{split}$$

Par récurrence forte sur $p$, on en déduit que $B_p\in \mathbb Q$ pour tout $p\in {\mathbb N}^*$.

La dernière relation annoncée, qui lie$\;{B_p}$ à ${A_p}$, sera établie plus loin (proposition 5). $\square$

Images des entiers naturels non nuls par $\alpha $ et $\beta $

Le calcul des images des entiers naturels non nuls par les fonctions $\alpha $ et $\beta $ est élémentaire à partir de celui effectué avec les fonctions $\zeta $ et $\xi $. Leurs valeurs aux entiers pairs non nuls sont donc connues à partir des propositions 2 et 4 :

Démonstration :

Pour tout $p\in {\mathbb N}^*$ :

$$\zeta \left( p \right) =  \sum _{n = 1}^{ + \infty } \frac{1}{n^p} =  \sum _{m = 1}^{ + \infty } \frac{1}{( 2m )^p} +  \sum _{m = 0}^{ + \infty } \frac{1}{( 2m + 1 )^p} = \alpha \left( p \right) + \;\beta \left( p \right)$$

et

$$\begin{split} \xi \left( p \right) & =  \sum _{n = 1}^{ + \infty } \frac{(  - 1 )^{n - 1}}{n^p} \\ & =  \sum _{m = 1}^{ + \infty } \frac{( - 1 )^{2m - 1}}{(2m )^p} + \sum _{m = 0}^{ + \infty } \frac{( - 1 )^{2m}}{( 2m + 1 )^p} =  - \alpha \left( p \right) + \;\beta \left( p \right).\end{split}$$

D’où les deux identités générales annoncées, puis leurs applications aux entiers pairs. Pour ces dernières, on utilise aussi les observations que :

$$\alpha \left( {2p} \right) = \sum _{m = 1}^{ + \infty } \frac{1}{( 2m )^{2p}} = \frac{1}{4^p}\zeta ( {2p} )$$

$$\zeta ( 2p ) =  \sum _{m = 1}^{ + \infty } \frac{1}{(2m )^{2p}} +  \sum _{m = 0}^{ + \infty } \frac{1}{( 2m + 1 )^{2p}} = \frac{1}{4^p}\zeta ( 2p ) + \beta ( 2p ).$$

dont résulte :

$$\beta \left( {2p} \right) = \frac{{{4^p} - 1}}{{{4^p}}}\zeta \left( {2p} \right).$$

C’est ce résultat qui justifie l’égalité annoncée à la fin de la proposition 3. En effet, on obtient par identification des deux expressions connues de $\beta \left( {2p} \right)$ la nouvelle identité :

$$\frac{1}{2}\left( {{B_p} + {A_p}} \right) = \frac{{{4^p} - 1}}{{{4^p}}}{B_p}$$

dont résulte :

$${A_p} = \frac{{{4^p} - 2}}{{{4^p}}}{B_p}.\; {\color{red} \square}$$

Images des entiers naturels impairs par $\psi $

Une approche similaire à celle employée au 3. permet le calcul des images des entiers naturels impairs par $\psi $ :

Démonstration :

Soit $F$ la fonction $2\pi $-périodique qui coïncide avec $f$ sur $\left] { - \pi ;\pi } \right[$ et est telle que $F\left( {n\pi } \right) = 0$ pour tout $n\in {\mathbb N}$. $F$ est impaire et de classe ${C^1}$ par morceaux. La série de Fourier de $F$ converge donc simplement vers la fonction régularisée de $F$ d’après le théorème de Dirichlet. On en déduit que pour tout $t \in \left] { - \pi ;\pi } \right[$ :

$$f\left( t \right) = \mathop \sum \limits_{n = 1}^{ + \infty } {b_n}\sin \left( {nt} \right)$$

où :

$$\forall n\in\mathbb N ^* , b_n=\frac{2}{\pi}\int_0^\pi f(t)\sin (nt) dt$$

Par conséquent :

$$f\left( {\frac{\pi }{2}} \right) = \mathop \sum \limits_{m = 0}^{ + \infty } {\left( { - 1} \right)^m}{b_{2m + 1}}$$

Or, pour tout $n\in {\mathbb N}^*$ fixé, on établit par récurrence sur $j\in {\mathbb N}$ que :

$${b_n} = \frac{2}{\pi } \sum _{k = 0}^j ( - 1 )^{k - 1}\frac{( - 1 )^nf^{(2k )}( \pi )}{n^{2k + 1}} + \frac{2( - 1 )^j}{\pi n^{2j + 1}} \int _0^\pi  f^{( 2j + 1 )}( t )\cos ( nt )dt.$$

En effet, on obtient en intégrant trois fois par parties et en tenant compte du fait que l’imparité de $f$ implique celle de toutes ses dérivées d’ordre pair (et donc leur nullité en 0) :

$$\begin{split} {b_n} & =  - \frac{2}{{\pi n}}{\left( { - 1} \right)^n}f\left( \pi  \right) + \frac{2}{{\pi n}} \int _0^\pi  f'\left( t \right)\cos \left( {nt} \right)dt \\ & =  - \frac{2}{{\pi n}}{\left( { - 1} \right)^n}f\left( \pi  \right) + \frac{2}{{\pi {n^3}}}{\left( { - 1} \right)^n}f''\left( \pi  \right) - \frac{2}{{\pi {n^3}}} \int _0^\pi  {f^{\left( 3 \right)}}\left( t \right)\cos \left( {nt} \right)dt.\end{split}$$

Ce qui montre que l’identité est satisfaite successivement pour $j = 0$ et $j = 1$. Supposons que l’identité annoncée vaut pour un certain $j\in {\mathbb N}$. On obtient, en intégrant deux fois par parties :

$$  \int _0^\pi  {f^{\left( {2j + 1} \right)}}\left( t \right)\cos \left( {nt} \right)dt  =  - \frac{1}{n} \int _0^\pi  {f^{\left( {2j + 2} \right)}}\left( t \right)\sin \left( {nt} \right)dt$$

$$\;\;\; = \frac{1}{{{n^2}}}{\left( { - 1} \right)^n}{f^{\left( {2j + 2} \right)}}\left( \pi  \right) - \frac{1}{{{n^2}}} \int _0^\pi  {f^{\left( {2j + 3} \right)}}\left( t \right)\cos \left( {nt} \right)dt.$$

En utilisant l’hypothèse de récurrence, on en déduit que l’identité est satisfaite au rang $j + 1$. L’identité annoncée, vraie par récurrence pour tout $j\in {\mathbb N}$, est donc vraie en particulier en $j = p$ ; c’est-à-dire que :

$${b_n} = \frac{2}{\pi }\mathop \sum \limits_{k = 0}^p {\left( { - 1} \right)^{k - 1}}\frac{{{{\left( { - 1} \right)}^n}{f^{\left( {2k} \right)}}\left( \pi  \right)}}{{{n^{2k + 1}}}} + \frac{{2{{\left( { - 1} \right)}^p}}}{{\pi {n^{2p + 1}}}} \int _0^\pi  {f^{\left( {2p + 1} \right)}}\left( t \right)\cos \left( {nt} \right)dt.$$

Or, $f$ étant polynomiale de degré inférieur ou égal à $2p + 1$ sur $\left[ {0;\pi } \right]$, la fonction $f^{( 2p + 1 )}$ y est constante. On en déduit :

$$\int _0^\pi  {f^{\left( {2p + 1} \right)}}\left( t \right)\cos \left( {nt} \right)dt = 0.$$

D’où résulte :

$${b_n} = \frac{2}{\pi }\mathop \sum \limits_{k = 0}^p {\left( { - 1} \right)^{k - 1}}\frac{{{{\left( { - 1} \right)}^n}{f^{\left( {2k} \right)}}\left( \pi  \right)}}{{{n^{2k + 1}}}}.$$

On obtient donc finalement :

$$\begin{split} f\left( {\frac{\pi }{2}} \right) & = \mathop \sum \limits_{m = 0}^{ + \infty } \frac{{2{{\left( { - 1} \right)}^m}}}{\pi }\mathop \sum \limits_{k = 0}^p {\left( { - 1} \right)^{k - 1}}\frac{{{{\left( { - 1} \right)}^{2m + 1}}{f^{\left( {2k} \right)}}\left( \pi  \right)}}{{{{\left( {2m + 1} \right)}^{2k + 1}}}}\\ &  = \frac{2}{\pi }\mathop \sum \limits_{k = 0}^p {\left( { - 1} \right)^k}\left( {\mathop \sum \limits_{m = 0}^{ + \infty } \frac{{{{\left( { - 1} \right)}^m}}}{{{{\left( {2m + 1} \right)}^{2k + 1}}}}{f^{\left( {2k} \right)}}\left( \pi  \right)} \right).\end{split}$$

Et cette identité n’est autre que celle annoncée dans la proposition 6. $\square$

Démonstration :

Soit $p\in {\mathbb N}$. On applique la proposition 6 avec $f\;:t \mapsto {t^{2p + 1}}$. On obtient alors :

$$\frac{{{\pi ^{2p + 1}}}}{{{2^{2p + 1}}}} = \frac{2}{\pi } \sum _{k = 0}^p {\left( { - 1} \right)^k}\frac{{\left( {2p + 1} \right)!}}{{\left( {2p + 1 - 2k} \right)!}} \; {\pi ^{\left( {2p + 1} \right) - 2k}}\; \; \psi \left( {2k + 1} \right).$$

Soit aussi :

$$\mathop \sum \limits_{k = 0}^p {\left( { - 1} \right)^k}\frac{{\left( {2p + 1} \right)!}}{{\left( {2p + 1 - 2k} \right)!}}\;\frac{1}{{{\pi ^{2k + 1}}}}\;\psi \left( {2k + 1} \right) = \frac{1}{{{4^{p + 1}}}}.$$

En posant ${C_q} = \frac{\psi \left( {2q + 1} \right)}{\pi ^{2q + 1}}$ pour tout $q\in {\mathbb N}$, on obtient la relation exprimant ${C_p}$ lorsque $p \ge 1$.

Dans le cas où $p = 0$, on obtient de plus avec la relation précédemment obtenue que $\psi \left( 1 \right) = \frac{\pi}{4}$, donc que $\frac{\psi (1)}{\pi} = C_0 = \frac{1}{4}\in \mathbb Q$.

Supposons alors que, pour un $p\in {\mathbb N}$ fixé, $C_j\in\mathbb Q$  pour tout $j \in [\![ 0;p ]\!]$. Il est clair que

$$\begin{split} {{C}_{p+1}}& =\frac{\psi \left( 2p+3 \right)}{{{\pi }^{2p+3}}} \\ & =\frac{{{\left( -1 \right)}^{p+1}}}{\left( 2p+3 \right)!}\left( \frac{1}{{{4}^{p+2}}}-\underset{k=0}{\overset{p}{\mathop \sum }}\,{{\left( -1 \right)}^{k}}\frac{\left( 2p+3 \right)!}{\left( 2p+3-2k \right)!}{{C}_{k}} \right)\in \mathbb{Q}\end{split}$$

Par récurrence forte sur $p$, on en déduit que $C_p\in\mathbb Q$ pour tout $p\in {\mathbb N}$. $\square$

Typologie liant la résolution du problème à quatre « niveaux » de parité

Tous les calculs possibles permettant de résoudre les problèmes initialement posés au moyen de l’approche utilisée ont été effectués. Concernant les fonctions envisagées au 3. et au 5., leur évaluation en d’autres réels de $\left[ {0;\pi } \right[$ que 0 et $\frac{\pi}{2}$ ne permet pas de progresser – le lecteur pourra notamment vérifier que l’application à $\frac{\pi}{2}$ d’une fonction $f$ satisfaisant les hypothèses de la proposition 1 ne fournit aucun résultat nouveau. On retrouve ici exactement les cas problématiques actuellement encore non réglés, qui restent soumis aux investigations contemporaines ; c’est-à-dire celui du calcul des images des entiers impairs par $\zeta ,\;\alpha ,\;\beta ,\;\xi $ et $\varphi $, et du calcul des images des entiers pairs par $\psi $. On ignore en fait toujours si, pour $p$ donné, le rapport à ${\pi ^{2p + 1}}$ des images de $2p + 1$ par les cinq premières fonctions est un nombre rationnel, et la même ignorance vaut pour le rapport à ${\pi ^{2p}}$ de $\psi \left( {2p} \right)$, pour tout $p$. En particulier, on ignore même si, hormis le cas de $\zeta \left( 3 \right)$, démontré irrationnel en 1978 par le mathématicien français Roger Apery, l’un de ces nombres arbitrairement choisi est bien irrationnel.

Le cas des images de 1 par les fonctions $\zeta ,\;\alpha ,\;\beta ,\;\xi $ et $\varphi $ mérite enfin d’être relevé car sa spécificité interroge au regard de tout ce qui précède. Il est clair que $\zeta \left( 1 \right)$, $\alpha \left( 1 \right)$ et $\beta \left( 1 \right)$ n’existent pas. Il n’en va par contre pas de même de $\xi \left( 1 \right)$ et de $\varphi \left( 1 \right)$, mais il se trouve que leurs valeurs ne sont pas, contrairement à la logique qui se dégage des cas résolus, un multiple rationnel de $\pi $. En effet :

$$\xi \left( 1 \right) = \ln \left( 2 \right){\rm{\;\;\;et\;\;\;}}\varphi \left( 1 \right) = \frac{1}{2}\ln (2).$$

ll est d’ailleurs remarquable que la détermination de ces deux valeurs ne s’obtient pas par le même type de procédure fondée sur les séries de Fourier. La seconde se déduit sans difficulté de la première. La valeur de $\xi (1)$ se déduit du développement en série entière de $\ln (1+x)$ et de l’argument de la convergence de la série en 1 par le critère spécial de convergence des séries alternées. Elle peut aussi être établie de manière moins sophistiquée comme suit :

$$\forall t \in {\mathbb R}, \forall n \in {\mathbb N}, \sum_{k=0}^n (-1)^k t^k = \frac{1}{1+t}+(-1)^n\frac{t^{n+1}}{1+t}$$

D'où résulte par intégration sur $[0;1]$ :

$$\sum_{k=0}^n \frac{(-1)^k}{k+1} = \ln 2 + (-1)^n\int_0^1\frac{t^{n+1}}{1+t}dt$$

La valeur de $\xi \left( 1 \right)$ s’en déduit par passage à la limite lorsque $n$ tend vers $ + \infty $, compte tenu de :

$$\forall ~n\in \mathbb{N}, \; 0\le \int_0^1 \frac{{{t}^{n+1}}}{1+t}~dt\le \int_0^1{{t}^{n+1}}dt=\frac{1}{n+2}$$

qui implique

$$\mathop {\lim }\limits_{n \to  + \infty } \int _0^1 \frac{{{t^{n + 1}}}}{{1 + t}}\;dt = 0.$$

L’approche utilisée dans cet article a le mérite de révéler de manière relativement élémentaire deux aspects intimement liés auxquels se rapporte la résolution du problème considéré : (1) le déphasage des fonctions sinus et cosinus, dont la conséquence est manifeste dans l’inégalité du nombre de cas résolus par les propositions 1 et 6 ; (2) les liens profonds existant entre les résolutions cherchées et quatre situations différentes de parité : (a) la parité des fonctions $F$ utilisées ; (b) la parité du degré des restrictions polynomiales $f$ de ces fonctions à $\left[ {0;\pi } \right[$ ; (c) la parité des entiers formant les arguments des six fonctions sur lesquels porte la question initialement posée ; (d) la parité des entiers dont les puissances inverses interviennent dans les sommes définissant les quatre fonctions $\alpha ,\;\beta ,\;\varphi $ et $\psi $.

Les deux tableaux suivants synthétisent une typologie des cas de résolution et des cas de non résolution relativement à ces quatre « niveaux » de parité :