2

バーゼル問題の比較的分かりやすい解法

716
0

導入

今回は、バーゼル問題についてです。

バーゼル問題(バーゼルもんだい、英: Basel problem)は、級数の問題の一つで、平方数の逆数全ての和はいくつかという問題である。ヤコブ・ベルヌーイやレオンハルト・オイラーなどバーゼル出身の数学者がこの問題に取り組んだことからこの名前で呼ばれる。

(引用:Wikipedia)
要するに、
n=11n2=1+122+132+
の収束値を問うているわけです。驚くべきことに、この答えはπ26であることが、かの大数学者オイラーにより導かれました。
その求め方として、sinxの因数分解とマクローリン展開の係数比較をする、という手法が有名です。

しかし、マクローリン展開はある程度納得できるにしても、sinxが因数分解できることはそこまで自明でないように感じます。
他にも、フーリエ級数展開を用いた手法や、重積分を用いた手法などがありますが、大学程度の予備知識が必要で、(個人的には)
まだ分かりにくいと感じていました。

今回私が見つけた手法は、逆三角関数、広義積分と、ライプニッツの積分法則(ファインマンの手法)さえ知っていれば高校生でも十分理解可能なものです。
この記事では、その手法でバーゼル問題を解いていきます。
(※数学的な厳密性は、分かりやすさのため欠いています。一応、私自身各過程で成り立つことは証明しているはずですが、微妙な部分が含まれているかもしれません。)

証明

級数の変形

まず、求める級数を、偶数部分と奇数部分に分けます。
n=11n2=n=11(2n1)2+n=11(2n)2
左辺第2項を移行して、
34n=11n2=n=11(2n1)2
今、n=11n2=π26であることを示すには、n=11(2n1)2=34π26=π28であることを示せばよいです。

そのため、ここから下は、
n=11(2n1)2=π28
を示します。

関数I(a)とその導関数

まず、次のような関数を用意します。
I(a)=0tan1(ax)1+x2dx

ここで、tan1(x)は逆正接関数で、その導関数が11+x2になることは、大学受験を経験した方ならご存じだと思います。

x=tan(y)   (π2<y<π2)とする。両辺をyで微分すると、dxdy=1cos2(y)=1+tan2(y)=1+x2となる。
一方で、yの定義域内でxとyは一対一対応するから、tan1(x)=tan1(tan(y))=yすなわちy=tan1(x)である。
よって、逆関数の微分公式より、
ddx(tan1(x))=dydx=1dxdy=11+x2 

次に、I(a)の導関数を求めます。実は、次のことが成り立ちます。

dIda=ln(a)a21ただし、a=1での値は、左辺は極限lima1ln(a)a21で定義する。

これが成り立つことを確認するために、ライプニッツの積分法則(の特殊形)を用います。

ライプニッツの積分法則

二変数関数f(x,t)および実数a,bにたいして、
ddx(abf(x,t)dt)=abxf(x,t)dt

厳密な証明ではない、これが成り立つ説明は、例えば次のようなものが考えられます。
g(x)=abf(x,t)dtとすると、微分の定義から、
g(x)=limh0g(x+h)g(x)h
です。一方で、
limh0g(x+h)g(x)h=limh01h(abf(x+h,t)dtabf(x,t)dt)=limh0abf(x+h,t)f(x,t)hdt
です。
h0の極限と積分が交換可能ならば、
limh0abf(x+h,t)f(x,t)hdt=ablimh0f(x+h,t)f(x,t)hdt=abxf(x,t)dt
で、確かに成り立ちそうです。

ただし、このままでは、元のI(a)の積分範囲の上端がであることから、
直ちには用いることはできません。そこで少し工夫をします。
まず、積分範囲を分けます。
0tan1(ax)1+x2dx=01tan1(ax)1+x2dx+1tan1(ax)1+x2dx
右辺第2項について、x=1tの変数変換をすることで、
1tan1(ax)1+x2dx=01tan1(ax)1+x2dx
であることが分かります。
よって、
I(a)=01tan1(ax)1+x2dx+01tan1(ax)1+x2dx
と、積分範囲を01にできました。この下で先ほどの補題1および逆正接関数の微分を用いると、
dIda=01a(tan1(ax)1+x2)dx+01a(tan1(ax)1+x2)dx
=0111+x2x1+(ax)2dx+0111+x211+(ax)21xdx
=01x(x2+a2+1+x2a2)(1+x2)(1+(ax)2)(x2+a2)dx
=01x(1+a2)(1+(ax)2)(x2+a2)dx
=1a2101(a2x1+(ax)2xx2+a2)dx
=1a21[12ln(1+a2x2)12ln(x2+a2)]01
=1a21[12ln(1+a2x2x2+a2)]01=1a21(12ln(1+a212+a2)12ln(1+a20202+a2))=ln(a)a21

これで、dIda=ln(a)a21が示せました。
(※a=1の値については、部分分数分解を行わず、x=tanθの置換をすることで、実際の微分係数と、lnを用いた表示での極限が一致していることが分かります。)

I(a)の特殊値を求値

I(0)およびI(1)を、定義に従って求めていきます。
I(0)=001+x2dx=0
I(1)=0tan1(x)1+x2dx=[12(tan1(x))2]0=12(π2)2=π28

定積分の値

ここで、
J=01ln(a)a21daという定積分を考えます。dIda=ln(a)a21およびI(0),I(1)の値から、
J=01dIdada=I(1)I(0)=π28
です。
一方で、幾何級数公式11x=n=0xn(|x|<1)を用いると、
J=01ln(a)1a2=01ln(a)n=0a2nda
このとき、積分と極限の交換について、ルベーグの優収束定理を用いれば交換可能であることが分かります。ここでは交換可能であることを認めると、
01ln(a)n=0a2nda=n=001ln(a)a2nda
ここで、01ln(a)a2ndaについて、部分積分をして、
01ln(a)a2nda=[ln(a)a2n+12n+1]01+12n+1011aa2n+1da
各項について、
[ln(a)a2n+12n+1]01=0+limϵ0+ln(ϵ)ϵ2n+12n+1=0
12n+1011aa2n+1da=1(2n+1)2[a2n+1]01=1(2n+1)2
よって、
01ln(a)a2nda=1(2n+1)2
となります。これを先ほどの級数に代入して、
J=n=01(2n+1)2=n=11(2n1)2
J=π28であったので、
n=11(2n1)2=π28
であることが示せました。
よって、
n=11n2=π26
となります。

終わりに

いかがだったでしょうか。煩雑な式処理が随所にあるものの、難解な部分はほとんどなく、高校生でも十分理解可能であると思います。
まあ、簡単な証明かどうかはそこまで重要ではなく、自分自身で新たな解法を得られた、という点が非常に嬉しいことです。
つらい浪人生活も、「まあ、私はバーゼル問題の新解法を思いついたし」と思うと、不思議とやっていける、そんな気がします。
以上です。

投稿日:2024831
OptHub AI Competition

この記事を高評価した人

高評価したユーザはいません

この記事に送られたバッジ

バッジはありません。
バッチを贈って投稿者を応援しよう

バッチを贈ると投稿者に現金やAmazonのギフトカードが還元されます。

投稿者

n=1 帰納法の失敗

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中