『無限解析入門』 における誤差について 徳島大学工学部 高橋浩樹 1.『無限解析入門』の誤差とは? 2.誤差に注目する理由 3.ある誤差リストの解釈 4.仮説の全体像 1 1.『無限解析入門』の誤差とは? Leonhard Euler(1707-1783) 生誕300周年 「オイラーは人が息をするように、 また鷲が空を舞い遊ぶように、 見た目には何の苦労もなく計算した」 天文学者 フランソア・アラゴ 2 オイラーの解析三部作 『無限解析入門』 『微分計算教程』 『積分計算教程』 著述 1745 1748 1763 出版 1748 1755 1768~1770 ・『入門』には,数値データに60個を超える間違いがある. ・単なる誤植や数値の表現(四捨五入・切捨て)といった 理由ではないようである. 疑問 オイラーは数値計算が不得意だったのか? 3 log (1+x)/(1-x)のマクローリン展開を用いた計算 53527 16963 04904 4 sin (m/n π/2)の マクローリン展開 プログラム 16 38 94 39 08 85 77 32 56 74 60 05 51 39 50 5 cos (m/n π/2)の マクローリン展開 00 49 63 72 62 30 70 10 36 23 86 15 30 84 67 28 6 tan (m/n π/2) cot (m/n π/2)の マクローリン展開 35 15 373 69 33 53 60 27 58 7 (ゼータ値の近似値) ×(1-1/2n) 548818 8 (ゼータ値の近似値) ×1/2n 10 75 04 70 96 73 15 271503 08 53 9255 55 13 86 84 08 85 58 10 77 45 86 71 68 9 素数ベキ和の近似値 ゼータ値の近似値を組 み合わせて求められる 1065 246 37 18 74437 35 25 19 02 24 46 84 55 33 26 30 07 51 10 2.誤差に注目する理由 ゼータ関数に関わる2つの疑問 A.ゼータ関数の「美しい関係」 なぜオイラーは,関数等式を美しいと 形容したのか? B.ゼータ関数の特殊値 なぜオイラーは,特殊値を数多く計算 したのか? 11 A. 『美しい関係』(E352) 美しい? 定義 12 注)1749年に著述 太陽と月の記号 13 1748年の金冠日食に関する文献 E117 [著述 1748, 出版 1749] Reflexions sur la derniere eclipse du Soliel du 25 julliet a. 1748. 金冠日食の観測データ E142 [著述 1748, 出版 1750以降?] Sur l'atmosphere de la lune prouvee par la derniere eclipse annulaire du soleil. 月の大気に関する推測 14 解釈 上の式・・・太陽が月の背後に位置する. 下の式・・・太陽の光の環がなおも輝く. Z(n-1)/Z(-n) Γ(n) cos(nπ/2) πn 15 B. ゼータ関数の特殊値のリスト E41 [著述 1735, 出版 1740.] バーゼル問題の解決 12 までのゼータ値 E101 [著述 1745, 出版 1748.] 『無限解析入門』 ・・・「いくつか書き添えておく」 26までのゼータ値 E212 [著述 1748, 出版 1755.] 『微分計算教程』 30までのゼータ値 E352 [著述 1749, 出版 1768.] 美しい関係 ・・・ 「計算した限りを示す」 34 までのゼータ値a b c 16 ゼータ関数の特殊値 ベルヌーイ数との関係(オイラーの発見) 17 予測できる素数たち 18 周期p-1 予測が困難な素数たち 19 分子に現れない素数 ⇔B0,B2,B4,…..,Bp-3のどの分子もpで割れない。 正則素数 2、3、5、7、11、13、17、19、 23、29、31、41、・・・・・・ 分子に現れる素数 ⇔B0,B2,B4,…..,Bp-3のどれかの分子がpで割れる。 非正則素数 37、59、67、101、103、131、・・・、 283、・・・、593、・・・、617、・・・、 683、691、・・・・・・ 20 疑問 ・691という素数がζ(12)に突如現れたことに 興味を抱かなかったのか? ・34までのゼータ値を求めたのに,最小の非正則 素数37に気づかなかったのか? 仮説 『無限解析入門』のゼータ値の近似値の誤差に, オイラーは非正則素数を書き記した. 基本的には,誤差の素因数として書き記した. 1998=54・37 1998-1526=472= 8・59 804=12・67=37+59+67+101+103+131+149+157 公表しなかった理由は後述 21 60を超える誤差の解釈 『無限オイラー解析』 現代数学社 この本はオイラーの問題集. 未発見の解答が残っている. 「理系への数学」(現代数学社) 9月号より詳細な探究の連載 『オイラー数学の源流』 第1回 巨人オイラー 第2回 超越への助走 22 3.ある誤差リストの解釈 sin と cos のマクローリン展開の係数の誤差リスト 01 -0.0000000000000000000000000000 03 +0.0000000000000000000000000002 05 -0.0000000000000000000000000006 07 +0.0000000000000000000000000007 09 -0.0000000000000000000000000003 11 +0.0000000000000000000000000005 13 -0.0000000000000000000000000006 15 +0.0000000000000000000000000008 17 -0.0000000000000000000000000006 19 +0.0000000000000000000000000004 21 -0.0000000000000000000000000004 23 +0.0000000000000000000000000002 25 -0.0000000000000000000000000001 27 +0.0000000000000000000000000000 29 -0.0000000000000000000000000001 00 -0.0000000000000000000000000000 02 +0.0000000000000000000000000004 04 -0.0000000000000000000000000004 06 +0.0000000000000000000000000008 08 -0.0000000000000000000000000004 10 +0.0000000000000000000000000004 12 -0.0000000000000000000000000005 14 +0.0000000000000000000000000002 16 -0.0000000000000000000000000006 18 +0.0000000000000000000000000003 20 -0.0000000000000000000000000005 22 +0.0000000000000000000000000003 24 -0.0000000000000000000000000003 26 -0.0000000000000000000000000001 28 -0.0000000000000000000000000002 30 +0.0000000000000000000000000002 23 気づくこと ・31個中28個もの誤差がある. ・数値に対する誤差の割合は急激に膨張している. ・誤差は最終一桁の範囲に収まっている. ・ひとつのデータのみ絶対値が正値より大きい. →奇妙な数値データに思える. 偶然の間違いではなく,意図的なものではないか? sin(正弦),cos(余弦) -楽譜であれば,面白いだろう. 24 楽譜への変換 -1-シ 0-ド 1-レ 2-ミ 3-ファ 4-ソ 5-ラ 6-シ 7-ド 8-レ 演奏0 どうすれば,曲になるのだろうか? オイラーの公式 eix=cos x+i sin x 係数を組み合わせて,元の指数関数の展開の順番で曲にする! 演奏1 25 リズムの調整 賛美歌を参考 26 楽譜 27 オイラーの曲であると推測する理由 28 29 4.仮説の全体像 『無限解析入門』がオイラーの問題集だとすると, 出題者の義務として解答集も残すだろう. 解答集であると考えられる著作 E352 『美しい関係』 E343 『ドイツ王女への手紙』 30 非正則素数と指数 (1000まで) 31 3年前に描いた非正則素数たちの宇宙 (p cos (2πk/(p-1)), psin (2πk/(p-1))) 32 E343 ドイツ王女への手紙 33 仮説の続き(公表しなかった理由) 非正則素数に関して,オイラーは神秘的な事実を見出した. けれども,神秘の排除という学界の流れの中では, その事実の公表は立場上差し控えるべきだった. 他方,オイラーは彼自身の驚くべき発見を何らかの形では 残しておきたかった. (ゼータ関数⇔天体⇔ ・・・) 1741年7月25日 ロシアからベルリンに到着 ↓ 1748年7月25日 金冠日食 1748年 『無限解析入門』出版 問題集 ↓ 1768年 『美しい関係』出版 解答集1 1768年 『ドイツ王女への手紙』出版 解答集2 34 仮説の反証可能性 ・誤差の生成理由を説明する. ・このような数値データは稀ではないことを示す. (誤差が曲になり,しかも生成理由が説明できな いデータの例示) 35 参考 log105の計算 0 0 7 3 5 5 1 5 log(ab)1/2=(log a+log b)/2 0 0 0 0 0 0 0 0 4 5 7 6 6 1 7 4 3 6 3 5 7 5 0 0 1 3 7 9 5 7 9 8 8 36
© Copyright 2024 ExpyDoc