Excel VBA 数学教室ではアフィリエイトプログラムを利用して商品を紹介しています。

πが無理数であることの証明

円周率 πsinπ=0,cosπ=1 を満たすような定数であると定義して、円周率が無理数であることを証明します。

円周率が無理数であることの証明

まず、次のような関数を定義します。
 f(x)=pnn!xn(πx)n
定義から明らかなように、f(x)
 f(x)=f(πx)
という周期性があり、m 階導関数について
 f(m)(π)=(1)mf(m)(0)
という関係をみたします。この f(x)sinx をかけて 0 からπ まで積分した値を I とおきます。
 I=0πf(x)sinxdx
この積分値 I について2通りの方法で評価します。

0<x<π より
 0<x(πx)<π(πx)<π2
n 乗して n! で割ります。
 0<pnn!xn(πx)n<pnn!π2n0<f(x)<pnn!π2n
0<x<π の範囲では 0<sinx<1 なので
 0<f(x)sinx<pnn!π2n
0 から π まで積分すると
 0<I<(pπ2)nn!π
ここで
 (1)limnαnn!=0
が成り立つので、適当に大きな n をとれば
 0<I<1
となります。

② 部分積分によって実際に I を計算します。ここで
 (1)f(x)sinxdx=k=0(1)k[f(2k+1)(x)sinxf(2k)(x)cosx]
という積分公式を用いると、
 I=2k=0n(1)kf(2k)(0)
が得られます。ここで二項展開の公式
 (a+b)n=k=0nnCkakbnk
を用いて f(x) を展開すると
 f(x)=pnn!k=0n(1)knCkπnkxn+k
xn の係数を簡単に aN で表すと、
 f(x)=anxn+an+1xn+1++an+kxn+k+
となります。これを n+k 回微分すると、xn+k 以降を微分した項だけが残り、また xn+k を微分した項のみ定数となります。よって、x=0 を代入すると xn+k を微分した項以外は全て消えます。よって
 fn+k(0)=pnn!(1)knCk(n+k)!πnk
となります。ここで円周率が自然数 a, b によって
 π=ba
と書ける(つまり有理数である)と仮定します(背理法)。f(x)p は任意ですから pa で置き換えてみると
 pnπnk=anπnk=an(ba)nk=akbnk
これは整数です。つまり I は整数ということになります。しかし①で I0 から 1 の間にある数だと評価しました。これは矛盾します。よって円周率無理数です。

エクセルや数学に関するコメントをお寄せください