■多角形連鎖(その3)

 オイラー級数を

  Hn=1/1^2+1/2^2+1/3^2+1/4^2+・・・+1/n^2

と定義します.オイラー級数は次第に減少する項からなりますが,nを無限大にしたとき,無限級数

  H∞=1/1^2+1/2^2+1/3^2+1/4^2+・・・=π^2/6

に収束します.

 1734年,オイラーが当時懸案の大問題(バーゼル問題)を解いたのですが,今回のコラムではこの漸近挙動を示す直角三角形連鎖を取り上げてみることにします.まずその前に発散する直角三角形連鎖について調べてみることにしましょう.

===================================

【Q1】平方根の螺線

 原点をO,点P1を(1,0)とします.点P1において長さ1の垂線P1P2⊥OP1を立て点P2(1,1)とします.すると,OP2=√2となります.さらに点P2において長さ1の垂線P2P3⊥OP2を立てます.OP3=√3となります.これを繰り返せば,1,√2,√3,・・・,√(n−1),√nが得られ,点P1,P2,・・・は螺線のような図形を作ります.

 動径ベクトルOPkの長さは√kになりますが,その偏角がπ,2πを越すときのkの値は?

[A1]平方根の螺線はピタゴラスの定理に基づく有名な図形です.

  θ=arctan(1/√(k−1))

としてΣθを求めると,

  π→ k=7

  2π→k=18

===================================

【Q2】オイラー級数の螺線

 原点をO,点P1を(1,0)とします.点P1において長さ1/2の垂線P1P2⊥OP1を立て点P2(1,1/2)とします.すると,OP2=√(1/1^2+1/2^2)となります.さらに点P2において長さ1/3の垂線P2P3⊥OP2を立てます.OP3=√(1/1^2+1/2^2+1/3^2)となります.これを繰り返せば,1,√(1/1^2+1/2^2),√(1/1^2+1/2^2+1/3^2),・・・,√(Σ1/k^2)が得られ,点P1,P2,・・・は螺線のような図形を作ります.

 動径ベクトルOPkの長さは√(Σ1/k^2)になりますから,k→∞のとき,  |OPk|→π/√6.また,

  |OP1|+|P1P2|+・・・+|Pn-1Pn|=1+1/2+1/3+1/4+・・・→∞.

それでは動径ベクトルOPkの偏角がπ,2πを越すときのkの値は?

[A2]オイラー級数の収束速度は非常に遅く,その結果,点Pkはゆっくり外向きの螺線を描きながら半径π/√6の極限円に限りなく近づきます.

 単精度で計算したので誤差があると思われるのですが,

  θ=arctan(1/k√Hk-1)

としてΣθを求めると,

  π→ k=69

  2π→k=3858

===================================

【3】雑感

 ζ(1)は発散し,オイラーはζ(2)=π^2/6を証明した.すべての偶数sに対しζ(s)の値は無理数であるが,アペリは1979年にζ(3)が無理数であることを証明した.その後,2000年にリボールが無限個の奇数sに対しζ(s)が無理数であることを証明した.2001年にリボールはこの結果を精密化し,ζ(5)からζ(21)までの奇数sのうち少なくとも1つのsについて無理数であることを証明した.同年,ズディリンはこの範囲をζ(5)からζ(11)までに狭めることに成功した.

−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−

(Q)y=x^xとy=x^1/xを微分せよ.

(A)それぞれy’=(logx+1)x^xとy’=(1−logx)x^1/x-2.x^1/xはx=eのとき最大値1.4446・・・をとる.

(Q)x^x=exp(xlogx)の級数展開を使うことにより

  ∫(0,1)x^xdx=1−1/2^2+1/3^3−1/4^4+・・・

を導け.

−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−

 (その2)で

  an→Π1/cos(π/i)→8.70

に収束すると書いたことが気にかかっている.

  an→Π1/cos(π/i)→およそ12

と記載した本もあり,収束値が大きく食い違っているからである.

 調和級数を

  Sn=1/1+1/2+1/3+1/4+・・・+1/n

と定義する.調和級数Snは次第に減少する項からなるが,nを無限大にしたとき,無限調和級数

  S∞=1/1+1/2+1/3+1/4+・・・

は発散する.

 しかし,実際に計算してみると数列{Sn}は収束の遅い数列であることはわかるが,発散することは確認できないであろう.同様に,数列{an}も収束の遅い数列であると思われるが,このように収束の遅い数列の場合,収束値に関しては何ともいえないのである.それではどのように計算すれば精密な収束値を求めることができるのだろうか?

===================================