リーマンゼータ関数 $\zeta (s) = \sum_{n=1}^\infty \frac{1}{n^s}$は，その零点となる$s$が負の偶数と，実部が1/2の複素数に限られるというリーマン予想に結びつく重要な関数であり，数理物理学の難しい議論だけでなく，統計力学の入り口のあたりにも少しだけ登場する。
自由電子気体やデバイ模型のあたりの積分の計算に必要となるのが，$\zeta(2)$や$\zeta(4)$の値だ。ただし，時間がもったいないので，結果だけ与えて終りたくなる。簡単な導出法は，$\sin x $ の級数展開によるもので，バーゼル問題を最初に解いたオイラーの方法だ。
$\sin x = x \Bigl( 1 - \dfrac{x^2}{3!} + \dfrac{x^4}{5!} -\dfrac{x^6}{7!} \cdots \Bigr) = x \ \Pi_{n=1}^\infty \Bigl( 1- \dfrac{x^2}{n^2 \pi^2} \Bigr)$
ここで，両辺の$\ x^{2n}\ $の係数を比較すると，$\zeta(2n)\ $の値が求まる。例えば，$n=1\ $の場合，$-\frac{1}{6} = - \frac{1}{\pi^2} \Sigma_{n=1}^\infty \frac{1}{n^2}\ $，すなわち，$\zeta(2) = \frac{\pi^2}{6}$であり，これがバーゼル問題の解であった。
$a_n=\frac{1}{\pi^2 n^2}$とおいて，$\zeta(4)$と$\zeta(6)$の場合を考えてみたい。$\theta_{ij}=1\ {\rm if}\ i<j \quad {\rm otherwise}\ =0$，$\theta_{ijk} = 1 \ {\rm if}\ i<j<k \quad {\rm otherwise}\ =0$として，対称和の一般的な分解方法を考える。
n=2の場合：
$\Sigma_{ij} a_i a_j = \Sigma_{ij} a_i a_j (\delta_{ij} + \theta_{ij}+\theta_{ji}) = \Sigma_{ij} a_i a_j (\delta_{ij} + 2 \theta_{ij}) $
$+\frac{1}{5!} = \Sigma_{i<j} a_i a_j = \Sigma_{ij} a_i a_j \theta_{ij}= \Sigma_{ij} a_i a_j (1-\delta_{ij}) /2 = ( ( \Sigma a_i) ^2 - \Sigma a_i^2 ) /2 $
$\therefore \zeta(4)/\pi^4 = \Sigma a_i^2 = (\frac{1}{6})^2 - \frac{1}{60} = \frac{1}{90}$
n=3の場合：
$\Sigma_{ijk} a_i a_j a_k= \Sigma_{ij} a_i a_j a_k (\delta_{ijk} + 6 \theta_{ijk}+ 3 \theta_{ij}\delta_{jk} + 3 \delta_{ij}\theta_{jk}) $
$-\frac{1}{7!} = -\Sigma_{i<j<k} a_i a_j a_k = -\Sigma_{ijk} a_i a_j a_k \theta_{ijk} = -\Sigma_{ijk} a_i a_j a_k \bigl( 1-\delta_{ijk}-3 \theta_{ij}\delta_{jk} - 3 \delta_{ij}\theta_{jk} \bigr ) /6 \\ = -\Bigl ( (\Sigma a_i)^3 -\Sigma a_i^3 - 3 \bigl\{ (\Sigma a_i)(\Sigma a_i^2)-\Sigma a_i^3 \bigr \} \Bigr )/6 $
$\therefore \zeta(6)/\pi^6 = \Sigma a_i^3 = \Bigl ( \frac{6}{7!} - (\frac{1}{6})^3 + 3 \frac{1}{6 \cdot 90} \Bigr ) / 2= \frac{1}{945}$
付録：
$X=(\Sigma a_i) (\Sigma a_i^2) = \Sigma_{ij} a_i a_j^2 = \Sigma_{ij} a_i a_j^2 (\delta_{ij} + \theta_{ij} + \theta_{ji}) = \Sigma a_i^2 + \Sigma_{i<j} a_i a_j^2 + \Sigma_{i<j} a_i^2 a_j $
$Y=\Sigma a_i a_j a_k (\theta_{ij}\delta_{jk} + \delta_{ij} \theta_{jk}) = \Sigma_{i<j} a_i a_j^2 + \Sigma_{i<j} a_i^2 a_j$
$\therefore X - \Sigma a_i^2 =Y$
