3
大学数学基礎解説
文献あり

ガウス日記を追ってみた

320
0

近世数学史談という本ご存じでしょうか.かの高木貞治先生が書かれ,幾人もの人を数学の道へと呼び込んだと言われる数学史に関する本ですが,今回はこの本の 6 つめの話『レム二スケート関数の発見』にある,ガウスの行った計算を追ってみようと思います.まずは本文中からその部分を書き抜いてみましょう(u21とある部分は,原文ではx21ですが,uで良いと思うので書き換えています).

u=0x11x2dx から

           u=x+1215x5+132419x9+135246113x13+13572468117x17+

それから

           x=u110u5+1120u91115600u13+2113536000u171607318240000u21+

級数の逆函数を求めることはラグランジュから学んで,計算は達者である.

まあ,「計算が達者」という次元の話なのかはさておき,ここでは「級数の逆函数を求めることはラグランジュから学んで」の部分から確認していきましょう.これは杉浦光夫『解析入門 Ⅱ』に載っています(Ⅸ 章定理 8.10)(わかりやすさのため多少書き換えています).

ラグランジュの公式

fが領域Dで正則で,一点aにおいてf(a)0をみたすとき,適当なρ>0が存在して,円板D(a,ρ)の正の向きの周をCとすれば,任意のwf(D(a,ρ))に対し

                     f1(w)=12πiCζf(ζ)f(ζ)wdζ

と表される.またf(a)=bを中心とするf1(w)のテーラー展開は,

          f1(w)=n=0bn(wb)n,bn=1nResz=a1(f(z)b)n=1n![dn1dzn1(zaf(z)b)n]z=0

で表される.

これを用いて計算していきます.u=f(z)とおくとf(0)=10なので定理の仮定をみたすので,とりあえずは安心です.f(0)=0周りで,f1の展開を考えます.

b0 の計算

b0=0 は明らかにそうです. 

b1 の計算

b1=[zf(z)]z=0=1です.まだ微分しなくていいので楽ですね.

b2からb4の計算

サボります.サボる理由はb5の計算を見ればわかります ちなみに2のときだけ計算して0になりました.

b5 の 計 算

まずいのは (d4/dz4)(z5/f(z)5)の計算です.簡単のために,D=d/dzとおきます.ライプニッツ則を使おうと思うので,Dk(f(z)5)を計算します.

  D(f(z)5)=5f(z)6f(z)
  D2(f(z)5)=30f(z)7f(z)25f(z)6f(z)
  D3(f(z)5)=210f(z)8f(z)3+90f(z)7f(z)f(z)5f(z)6f(z)
  D4(f(z)5)=1680f(z)9f(z)41260f(z)8f(z)2f(z)+90f(z)7f(z)2+120f(z)7f(z)f(z)5f(z)6f(4)(z)

多分あってるはずです.
ということで,

  D4(z5f(z)5)=k=04(4k)Dk(f(z)5)D4k(z5)

となります.このとき,分母にくるf(z)のうち次数最大は9ですので,全てf(z)9で括ることを考えます.

        =f(z)9{z5(1680f(z)41260f(z)f(z)2f(z)+90f(z)2f(z)2+120f(z)2f(z)f(z)5f(z)3f(4)(z))
         +20z4(210f(z)f(z)3+90f(z)2f(z)f(z)5f(z)3f(z))+120z3(30f(z)2f(z)25f(z)3f(z))
         +240z2f(z)3f(z)+120zf(z)4}

うむ.
ここで,カッコ内を展開したときに z9+といった形になると嬉しいわけですが,どうやらz5の項があるようなので,まずはその係数を見ます.z5の係数がでてくるのは下線部を引いた項のみで,その和は

        16804200+36001200+120=0 ()

となります.余談:この計算は一発で合いました.
さていよいよ本番です.90f(z)2f(z)2の項以外からはすべてz9の項がでてきます.ヒエエ.これを実際に計算すると,

        1680×21260×2+120×612×5+20×(210)×85+20×90×220×5×6
         +120×30×6512001200×45+48
       =12

となります!やった!あとはこれを 5!=120 で割ることで,u5の係数が110になることが実際に確認できました.もっと簡単なやり方があるかもですが,とりあえず数値が合ったので大満足です.

最後に問題を出しておきます.ここまで読んだ方への演習問題です.

u9の係数が1120になることを確かめよ.

参考文献

[1]
高木貞治, 近世数学史談, 岩波文庫, p,37
[2]
杉浦光夫, 解析入門Ⅱ, 東京大学出版会, p.312
投稿日:2022625
更新日:2023114
OptHub AI Competition

この記事を高評価した人

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

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

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

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

投稿者

とと
とと
17
2290
Mathlogを始めてみました。ところでMathlogってどんな対数なんですか?

コメント

他の人のコメント

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