値を計算する方法として、現代人が素朴に思いつくのは(1)いくつかの角度に於ける値を事前に計算しておき、一般の場合は、それを補間した値を使う(2)Taylor展開の有限項近似の二つの方法だと思う。Taylor展開を使う場合、角度をラジアン単位に変換する必要があるので、円周率を、ある程度の精度で知っていないといけない。コンピュータ用に、もう少し凝ったアルゴリズムが使われることもある/あったらしいけど、今のコンピュータでは、(2)の方法が使われることが多い。libc)では、アーキテクチャ独立な実装は、単純なTaylor展開を利用するものになっている。https://android.googlesource.com/platform/bionic/+/refs/heads/master/libm/upstream-freebsd/lib/msun/src/k_sin.c【追記】コメントで指摘をいただいた。小さいけど、上ソースコードで使われてる係数だと、(コード中のコメントによれば)誤差は、もう少し小さい。計算機イプシロンより小さくなるように、多項式の次数が決められているのかもしれないけど、そうすると、13次までのTaylor展開では、少し精度が足りない人間が手で計算していた時代には、(1)の方法で計算されていた。事前に計算された値は、三角関数表として利用可能になっている必要がある。級数に頼ることなく、計算しやすい三角関数の値は、等分点(360度の有理数倍)に於けるもので、それは、代数方程式の解となるから、(2)の方法に頼ることなく、数値的に計算できる。この場合、正弦値の計算は、本質的には、円に内接する正多角形の一辺の長さを計算するのと同等でもある。19世紀に、チャールズ・バベッジが計画した階差機関は、人手で作られた数表(三角関数表や対数関数表)の誤りをなくしたいという動機があったとされている。階差機関は、多項式を自動で計算(し印字)できる機械として設計されていて、三角関数なども、(2)に基づいて計算することを想定していたのだろう。現存する写本はないものの、紀元前2世紀のニカイア(当時は、ビテュニア王国という国に属していたはず)出身のヒッパルコスが、三角関数表の一種として、「弦の長さの表」を作っていたらしいことは、他の本などの記述から確実とされ、最古の三角関数表だったと考えられている。ヒッパルコスの動機は天文学にあり、それ以後も、三角関数表は、多くの場合、天文学とセットで出てくる。測量術や航海術でも使われるようになったのは、大航海時代より以降のことらしい。新しい土地へ行くことが増えたり、遠洋航海が増えたためだろう。地理学や地図作成技術自体は、非常に古くからあり、経緯度を定義したのはヒッパルコスだとされ、地図の起源は、もっと古い。三角関数表のあるところに天文学があるというのは正しいけど、天文学があっても、三角関数表がないことはある。中国では、(おそらく紀元前から)円周率の正多角形近似による計算が行われていたにも関わらず、一般の角度に於ける三角関数の値を計算するという意識はなかったのか、三角関数表のようなものは殆ど見られない。中国で三角法が発展しなかったことは、正弦や余弦という用語が作られたのが、17世紀であることからも察せられる。中国でも、夏(BC1900~1600)の時代に、日食の予報ができず死刑になった人の話とかあるし、それは伝説としても、暦を作るのは重要な仕事だったので、天文学がなかったわけではない。天文学に動機づけられないと、"三角関数"のような考え方が便利だとは、なかなか思わないものなのかもしれない。考えてみると、角度の概念は日常にもありふれてる気がするのに、一回転360度というのさえ、天文学の影響が見られる。中国では、長い間、一回転は、"365.25度"だったらしく、天文学の影響が、より明らか。円周率が3.16だった時代(?)三角関数表を作るのに必要な計算技術は、円周率の多角形近似と同じで、三角関数表と違って、円周率の多角形近似は、複数の地域で追求がされた。歴史上使われた円周率の近似値を見ていると、3や3.14...以外に、3.16に近い値が、いくつかあるのに気付く。外接正24角形を使って近似計算した結果、出てきた数値なのかもしれない。仮にそうだとして、外接正24角形が頻出する理由が何かあったのか気になる。内接正多角形より外接正多角形を好む理由として、"円周と直径の比=円の面積と半径の二乗の比"という事実が、容易に分かるということが考えられる。外接正多角形での近似を認めれば、"円周と直径の比"及び"円の面積と半径の二乗の比"は、同じ数列の極限となって、等しいことは自明となる。内接正多角形での近似を使った場合は、2つの量は、異なる数列の極限になるので、極限値が等しいことは、自明じゃない不思議なことに、ユークリッドの『原論』には、円の面積と直径の二乗が比例することは書かれている(12巻)のに、"円周と直径の比=円の面積と半径の二乗の比"であることは書かれてないらしく、一般的には、これは、アルキメデスの成果とされている。『原論』は計算が殆ど出てこない特殊な本だし、また、小アジア〜地中海圏以外の地域の人にまで、アルキメデスの議論が伝播して理解されたと考えるのは、無理があるように思える。文献上では、内接正多角形を使った近似値だけでなく、外接正12角形や外接正48角形を使ったと思われる円周率の近似値も、あまり見られないっぽい。外接正96角形による近似計算→で、外接正96角形は、アルキメデスが、円周率の上界の評価に用いたものでもあり、得られる値は、円周率の最良近似分数の一つ22/7に近い。外接正24角形は、やっぱ特別らしい。内接正12角形の一辺の長さ→などとなって、一辺の長さが、正有理数の平方根の和の形になる場合に限れば、外接正24角形で近似精度が一番高くなるのが理由かもしれない。7289)があるらしい。methodと呼ばれるアルゴリズム(Newton法の特殊形)が存在しているけど、バビロニアで、このアルゴリズムが使われてた直接的な証拠はないよう。Euclidのアルゴリズムより古い、最古のアルゴリズムかもしれない。エジプトでは、どうだったか調べてないけど、平方根の計算が、非常に古くから行われてた可能性はあるだろう。アルゴリズムとしては、無理数の平方根も整数の平方根も同じだけど、手計算しかなかった時代に、無理数の平方根を計算する労苦を厭う気持ちは理解出来る。methodは収束が速い有能なアルゴリズムなので、例えば、1を初期値として僅か3ステップで、√2と√3の近似値として、577/408と97/56を得る。この近似値から外接正24角形の周長によって、円周率を概算するととなる。整数解になっていて、最良近似分数を与える。古代インドのシュルバ・スートラ(BC800~BC200)には、√2の近似値として、17/12と577/408が記載されているらしい。methodを使ったなら、1→3/2→17/12→577/408→...が順に求まる。古代のインド人が、どういう方法を使ったか詳細は知られてないようである。Circle"の中で、説明なしに、√3は、265/153と1351/780の間にあると述べたらしい。methodでは得られない。収束する分数変換は無数にある。methodより収束は速いはず。x=1に繰り返し適用すれば、265/153と1351/780を得られる。アルキメデスは、これに類した変換を用いた可能性が考えられている。平方根の計算のたびに、分数変換を探してると面倒なので、系統的な方法が望まれるけど、そのような一般論を知ってたかどうかは定かでない外接正24角形による近似が本当に多用されたのかは分からないけど、以上のようなシナリオはありえないこととも思えないので、宇宙人や異世界人も、初期文明では、3.16に近い近似値を採用する可能性があるエジプトやバビロニア、インド、古代ギリシャで、平方根の計算が、どれくらい一般性を持って理解されていたかは、よく分からないが、ある程度は計算する手段を持っていたようである古代オリエントの三角関数表ヒッパルコス(BC190?~120?)以前に、三角関数表に類するものがなかったかは分からない。2017年に、バビロン第一王朝(紀元前19〜16世紀)時代のバビロニアに、"三角関数表"があったという論文が出ていた。形の数値で、1列目の大きさ順でソートされている。いう解釈もできるらしいけど、どっちを選んでも、大きな影響はない。(p,q)は正整数の組で、一行目から順に、(12,5),(64,27),(75,32),(125,54),(9,4),(20,9),(54,25),(32,15),(25,12),(81,40),(60,30),(48,25),(15,8),(50,27),(9,5)とすると、表が再構成される。表は、60進表記で、各桁が、60の何乗に対応するのか不明らしい。11行目の2,3列目は、(45,75)かもしれないけど、それぞれを60倍した(2700,4500)と解釈されている。いうことは以前から分かっていたけど、これらの組を選んだ基準とか、この表の目的については、共通見解がない。Priceの1964年の提案通り、38行目まで続きがあったとして解釈をしている。ピタゴラス数なので、表は、直角三角形のリストで、鋭角の大きさに従って、ソートされているとも解釈できる。角度の分布は、残っている15行目までだと、45度弱〜32度程度だけど、38行目まで補うと、最小角は2.34度弱となる。表には、12709と18541のような大きな数が含まれているので、バビロニア人が、原始ピタゴラス数(12709,13500,18541)を作り出す方法を知ってたことは、間違いないと思われる論文には、表の使い方とかについても書いてあるけど、仮に、主張が完全に正しくても、通常の意味での三角関数表とは違って、三辺の長さが整数比という縛りがあるので、等分点に於ける弦長や正弦値を計算しているようなものではない。"三角関数表"と呼ぶのは、語弊が大きい。ピタゴラス数やペル方程式は、初等整数論に於ける主要なテーマだけど、元々の動機は、"応用数学"にあると言えるのかもしれない。数論には最近まで何の応用もなかったと言われることもある気がするけど、単に人類が、その起源を忘れてるだけだろう。バビロニアは、バビロン第一王朝以後も、長い間、少なくとも天文学では、先進地域であり続けた。当時、天文学と占星術は結びついていて、ヘレニズム時代になると、バビロニア占星術が伝わった結果、ホロスコープ占星術として、アレキサンドリアで発展したのは、その表れと言える。バビロニアの天文記録は、現存する一番新しいものでは、西暦74〜75年のもので、楔形文字記録の最後のものでもある。Bなるルールがあるらしいけど、これは、バビロニア由来のものと思われている。紀元前4世紀の人であるKidinnuが考えたという説があり、ギリシャに紹介したのはヒッパルコスであるかもしれない。ヒッパルコスの知識の内、どれがバビロニア由来のものかは明らかではないヒッパルコスの表は、7.5度刻みだったとされている。これは、内接正48角形の一辺の長さを計算できれば十分。計算したものなので、3.75度刻みで正弦値を計算したのと同等。3.75度は、30度の1/8なので、30度の正弦値から始めて、15度、7.5度、3.75度の正弦値を計算していけばいい。これくらいの刻み幅でも線形補間した時の精度は、それほど悪いものではない。3.75度と7.5度の正弦値から5度に於ける値を線形補間で計算した場合、0.08711081689...となり、5度の正弦値は、0.087155742...で、誤差は0.1%より小さい少し後の時代になって、AD2世紀に、アレキサンドリアのプトレマイオス(83年頃〜168年頃)が書いたとされる天文学書『アルマゲスト』に掲載されている「弦長の表」は、現代まで写本が伝わっていて、0.5度刻みになっているらしい。chordshttps://en.wikipedia.org/wiki/Ptolemy%27s_table_of_chordsアルマゲストの表を作るには、0.25度における正弦値が必要だけど、30度や90度に於ける正弦値から始めて、二等分していくだけでは、この値を知ることはできない。計算方法は、色んなところで説明されているけど、基本的な方針は、既知の値から始めて線形補間が十分な精度となるまで二等分を繰り返し、最後に線形補間によって、0.5度に於ける弦長(0.25度の正弦値の2倍)を決定したらしい。出発点には、12度に於ける弦長を使ったらしい。正五角形の一辺の長さは、72度に於ける弦長であり、これと、60度に於ける弦長と合わせて、加法公式から、12度に於ける弦長が決まる。後は、二等分を繰り返して得られる、12/16度と12/32度に於ける弦長を線形補間して、0.5度の弦長を決定できる。プトレマイオスは、この表を使って線形補間するように書いているらしい。ヘレニズム文明圏の書写材料はパピルスで、紀元前2世紀頃、アナトリアのペルガモン王国で、羊皮紙が量産されるようになったと伝えられる。プトレマイオスは、占星術書『テトラビブロス』も書き、これが後の西洋占星術の基礎になったらしい古代オリエント→インドメソポタミア地域とインドは、シュメール人とインダス文明の時代には既に交易があったと考えられている。ヒッパルコス以降のある時期に、インドに、三角法が伝播したと推測されている。紀元前後、仏教とギリシャ文化は、パルティア、バクトリア、ガンダーラ地方や中央アジア方面で共存してたし、紀元1世紀にはインド洋で季節風貿易が行われてたので、どこで、どういう形で伝わっても不思議はないけど、正確な時期や伝来経路は、よく分かってない。Samratによって、アラビア語版をベースとして、なされたらしい)。文献レベルでは、三角法に関する記述が見られるようになる時期と、ギリシャ占星術に関する記述や影響が見られるようになる時期は、かなり近いっぽい4世紀後半〜5世紀初期に書かれたと推測されている天文学書『スーリヤ・シッダーンタ』には、三角法と天文学への応用が書かれていたようである。この本は、ヘレニズム数学で見られなかったsinとcosと等価な概念を使っているらしい。Sûrya-Siddhântahttps://books.google.co.jp/books?id=jpE7AAAAcAAJ紀元500年頃、インドの数学者・天文学者アーリヤバタは、『アーリヤバティーヤ』という本を書いて、これの内容は、現在まで残ってるらしい。アーリヤバタは、パータリプトラにあったナーランダー僧院の学長を務め、ナーランダー僧院には、大きな天体観測所があったとWikipediaには書いてある。ナーランダー僧院は、当時最大の仏教の学院で、少し後の時代には、玄奘三蔵も、636〜641年まで、ここに滞在した。tablehttps://en.wikipedia.org/wiki/%C4%80ryabha%E1%B9%ADa%27s_sine_table3.75度刻みなので、多分、ヒッパルコスと同じような方法で計算したと推測される。プトレマイオスは、円周率を377/120=3.1416666...としたらしいけど、アーリヤバタは、62832/20000=3.1416を採用したと書いてある。外接正384角形の周長から評価した場合、3.141662747...となり、内接正384角形の周長を使うと、3.1415576079...となり、アーリヤバタは、後者の数値の近似として、3.1416を採用したらしい。円周率はともかくとして、基本的には、この時代のインドに於ける三角法は、『アルマゲスト』より、大きく優れてる面は、ないっぽい。紀元7世紀には、紀元前から存在する都市ウッジャインの天文台の天文台長だったらしいブラフマグプタ(598〜668)はPell方程式を調べていて、恒等式を発見したとされている。これは、今日、Brahmaguptaの恒等式と呼ばれる。これは、加法定理でもあるけど、ブラフマグプタが三角関数と結びつけて考えたかは分からない。ブラフマグプタは、他にもブラフマグプタの補間公式という二次の補間法で知られる。ブラフマグプタの頃には、0を用いた十進記数法が完成していたようである。北インドでは、550年にグプタ朝が滅び、中心都市であったパータリプトラは衰退したが、滅んだわけではなく、ナーランダー僧院も1193年まで残っていたようである。古代インドの書写材料には、貝葉という、ヤシの葉を加工したものを使っていたらしい。中国では、3世紀頃には、正192角形や正3072角形の周長から、円周率を評価していたらしいから、正弦の特殊値の計算として見れば、プトレマイオスのものより細かい計算ができていたことになる。概念的な面では、三角法には繋がっていないっぽい。"(『周官』大司徒に「夏至の日の南中時、8尺の晷針を立て、その影が1尺5寸になる地を「地中」と謂う」と。(鄭玄の)説に「そこは、(夏至の南中時に)太陽が真上にくる地点より1万5000里離れている」)という記述があるらしい。(道具)https://ja.wikipedia.org/wiki/%E8%A1%A8_(%E9%81%93%E5%85%B7)のことだろう。"8尺の晷針を立て、その影が1尺5寸になる地"は、北緯34度あたりで洛陽や長安のある緯度に近い。夏至の南中時に太陽が真上にくる地点とは、今で言う北回帰線(北緯23.4度)で、南越国とかあったあたりまで南下すれば、確認できたはずだけど、実際に行って測定したかは定かでない。方法で"太陽の高さ"が既知(上の記述では、夏至の南中時には、高さ8万里にあるという想定らしい)になれば、夏至の南中時に影の長さを測定するだけで、他の地点でも、北回帰線との距離が決定できる。この計算は、ある意味では、三角法と言えなくもないけど、三角関数表は全く必要ない。エラトステネスと違って、結果的には、古代中国では、誤った仮定に基づいて計算していたわけで、太陽と地球の距離に関しては大きな差が出たけれど、地上での距離計測に関して言えば、古代中国人の生活圏では、十分実用的で、実際に二地点間の距離を測定して、自分たちの計算の正しさへの確信を深めた可能性さえある。「地中」と北回帰線の距離は、1200km弱で、これが一万五千里だと言ってるから、一里80mくらいの換算になる。普通は、一里4〜500mくらいと考えられているけど、こうした計算があることは、魏志倭人伝の解釈で使われることがある短里説(一里75〜90m)の根拠とされている。紀元718年に、インド人を祖先に持つとされる瞿曇悉達は、唐で『開元占経』を編纂し、これの104巻には、位取り記数法の説明とかが書いてあり、この中にアーリヤバタの正弦表も含まれている。本とされ、冒頭は、古代に於ける宇宙論のような内容を含む。『開元占経』は、写本が沢山残っているらしく、検索すると、いくつか見つかる。殘三卷http://kanji.zinbun.kyoto-u.ac.jp/db-machine/toho/html/C004menu.htmlにある。右天竺算法,用上件九個字,乘除其字,皆一舉禮而成。..."とか書いてある部分(と、その続き)が、位取り記数法の説明っぽい(?)。第一段,二百二十五第二段,二百二十四,並四百四十九。"(以下略)のあたり。3.75度ごとの正弦値が書かれている。...と書かれている。この時に、日本にも、正接表が持ち込まれた可能性はある。測影鉄尺は詳細がよく分からないけど、冬至の影を測る鉄の棒などと解釈されてることが多い。圭表と同じものなのかもしれない。成立して、これらの技術は、江戸時代まで国家の管理下にあり、秘匿されていた。世界最初のTAN表の制作方法http://www.osaka-kyoiku.ac.jp/~jochi/j18.htmでは、一回転を364"度"として、実測によって、表を作ったのではないかと考察をしている。インド式として、一周360度が紹介されてるのに、中国人が、頑なに360を採用しなかった理由は謎である。一行の正接関数表(724AD)http://www.wasan.earth.linkclub.com/sugakusipdf/153.pdfには、宣明暦でも、一行の表が踏襲されたと書いてある。宣明暦は、江戸時代に改暦されるまで、日本で、長期に渡って用いられていた暦。これ以後、清代に、ヨーロッパ人がやってくるまで、三角関数表が中国で積極的に活用された例は、知られていない。圭表は、これ以後の中国でも使われていたようなので、正接表は使われていても、不思議はないけれどインド→アラビアローマ帝国末期には、各都市で、非キリスト教系の施設が閉鎖・破壊され、異教徒は殺されたり追放された(アレクサンドリア図書館は、紀元400年頃、破壊され、アテネのアカデメイアも、紀元529年に閉鎖された)ため、多くの学者が、亡命先としてササン朝へ移動したと考えられている。ササン朝の学術中心は、ジュンディーシャープールだったらしい。ササン朝では、ビザンツ帝国やインド、唐など各国の書物を集め、パフラヴィー語に翻訳したとされる。ササン朝は651年にイスラム教徒に征服されて滅亡した。皇族は、唐に亡命したらしいけど、首都クテシフォンの人口は、維持されたようだし、ジュンディーシャープールの学院は、832年にバグダードで設立された知恵の館と競合したと言われるので、その後も、200年以上、学術の中心地であり続けたようだ。して、ササン朝の技術と科学が、後のイスラム科学の基礎となったと考えられている。歴史書を見ても、その詳細は分からない8世紀には、アッバース朝のアル・ファザーリが、インドの数学書、天文学書を、アラビア語翻訳して、『シンドヒンド』と呼ばれる本に、まとめた。中国と異なり、アラビアの方では、インドの数学や天文学を継承して発展させた。al-Hindī)を書いて、この中には、インド式の十進記数法や三角法が含まれるらしい。当時のインドでは、まだ球面三角法はなかったけど、ヘレニズム数学の成果を継承したイスラム世界では、実用的な球面三角法が発展した(アレキサンドリアのメネラウスが、球面幾何学の創始者とされている)。8〜12世紀のイスラム世界で起こったことの詳細は、様々な本で書かれているので省略。第4代モンゴル帝国皇帝モンケの時代に、弟のクビライは中国地域の征服を担当し、三弟のフレグは西アジアの征服を担当して、それぞれ、後に、元とイルハン朝の創始者となった。フレグは、1258年に、バクダードを征服して、この時、大量の本が失われたと伝えられている。イスラム世界の中心地は、カイロに移る一方で、モンゴル帝国も、イスラムの科学者を連行し、その知識を継承・発展させようと試みている。フレグは、1259年には、中央アジアのマラーゲ(現在のイラン東部)に、マラーゲ天文台を建設させ、イスラムの学者を連れてきて、研究させたようである。マラーゲ天文台で働いた人々は、マラーゲ学派と呼ばれる元では、1271年、北京に回回天文台(イスラム天文台)が作られ、初代天文台長には、ジャマールッディーンという人が任命された。ジャマールッディーンも、それ以前はフレグに仕えていた人とされる。この頃の中国には、イスラムの数学・科学が、かなり持ち込まれた。元代に、イスラム圏の学術書が中国語翻訳されたという話は、聞かないけど、郭守敬は授時暦を作るに当たって、イスラム天文学の観測機器を利用したと言われる(具体的に、何を使ったのかは知らないけど)。クビライは、ヨーロッパの科学に精通した人も集めようとしていたらしい。ヨーロッパは遠いので、ヨーロッパ人は、あんまり行かなかったようだけど。マルコ・ポーロは有名人になった。後に、(ヨーロッパではないけど、中国までの距離は同じ)モロッコから、はるばる泉州、北京までやってきたイブン・バトゥータ(1304~1368)のような人もいる少し後になって、中央アジアには、モンゴル帝国の継承政権の一つであるティムール朝が興り、首都サマルカンドには、1420年代に、ウルグ・ベク天文台が建設された。ウルグ・ベクはティムール朝第四代君主であるけど、当人も、天文学や暦学の研究を行ったらしい。アル=カーシー(1380〜1429)は、ウルグ・ベクに招かれて、サマルカンドを中心に活動したとされる。彼の本には、1度の正弦値を求めるのに、3度の正弦値から、三次方程式を使って計算する方法が書かれているらしい。本当かどうか知らないけど、フランスでは、余弦定理を「アル・カーシーの定理」と呼ぶと、Wikipediaに書いてある。アル=カーシーは多角形近似の計算を推し進め、1424年の『円周論』で、円周率を高い精度で計算したらしい。アル=カーシーは60進数表記を使っているが、10進表記だと、3.1415926535897932...くらいまで合ってる計算になるよう。アル=カーシーは、別に趣味で、この計算をしたわけではなく、高精度の円周率が天文学で必要だと考えたのだそうだインド数学と三角関数のべき級数展開インドでは、12世紀に解析学の始まりが見られる。バースカラ2世の『シッダーンタ・シローマニ』(1150)に、三角関数の微分を含む解析学の萌芽が見られるとされる。ヘレニズム数学にはなかったsinとcosを導入しておいたことが、役立った。バースカラ2世は、インドで球面三角法を導入した人でもあるバースカラ2世は、ウッジャインの天文台の天文台長だったらしい。1350~1425)という人が現れ、ケーララ学派と呼ばれる数学と天文学の学派の創始者と伝えられている。ケーララ州は、インド南西部の一地方で、それ以前のインド数学の中心は、もっと北部にあったけど、多分、イスラム勢力が南下してきたため、Hindu数学の中心地は、南インドの方に移動したのだろうケーララ州には、コーリコード(カリカット)を中心とする王国や、コーチン王国などがあり、彼らは、この地方のみで使われるマラヤーラム語を話していたらしい。この時代の南インドには、Hindu系の大国ヴィジャヤナガル王国(1336〜1649)があり、関係がよく分からない。属国の類かもしれないけど、コーチン王国には、1409年に、明の鄭和が訪れ、明の保護国となっている。1500年に、ポルトガルの第2回インド遠征隊は、コーチン王国で歓迎され、1503年から1663年まで、ポルトガルの同盟国となったある人口の推定では、首都ヴィジャヤナガルは、1400〜1550年頃、全てのヨーロッパ都市とカイロを超える人口を擁する大都市であった。ヴィジャヤナガルの最盛期と、ケーララ学派が活躍した時期は、一致してるので、無関係ではないだろうと思うマーダヴァは、三角関数のテイラー展開や、円周率のライプニッツの公式を導いていたとされていて、また、Taylor級数のtrunctionを使って、三角関数表を作成したと信じられている。tablehttps://en.wikipedia.org/wiki/Madhava%27s_sine_table1530年頃までに、ケーララ学派が得た成果の多くは、『ユクティバーサ』という本に(証明付きで)まとめられ、世界初の微積分学の教科書と評される。実際には、微分については、大きな発展はなかったようで、速度が位置の時間微分であるという程度の認識はあったようだけど、微積分学の基本定理には到達しなかったっぽい。この本は、サンスクリットではなく、マラヤーラム語で書かれているらしい。どっちにしろ、読めないことには変わりないので、解説を見る。http://www.ms.uky.edu/~sohum/ma330/sp12/india_calculus.pdf解説も長いし、幾何学的議論は追うのが怠いけど、現代の視点から見ると、ケーララ学派は、積分を長方形近似で評価して極限を取ることで、多くの解析的結果を得たようだ。基本的なところでは、x^nの積分を、によって計算したらしい。後に、ヨーロッパでも、同様の計算は、行われたようである(FermatやFaulhaberなど)π/4の級数展開は、1/(1+x^2)の級数展開を項別積分するのと等価な式を示すことによって得た(上記解説中の式(124)など)sin(x)の級数展開については、上記解説中の式(169)(170)がキーになる式だけど、これらの式は、以下の2つの積分に対応するものと考えることができる。後者の積分は、両辺にsin関数が入っているので、積分を繰り返すことにより、冪級数展開の高次の項が決定される。やってることは、Picardの逐次近似に近いけど、微積分学の基本定理を持ってなかったケーララ学派に、何らかの一般論があったわけではないと思われる。ケーララ学派では、上記2つの式を、(現代的に解釈すれば)、積分の長方形近似で、直接計算して得たようである。極限を取れば、求める積分を遂行できる。cos関数について、同様の計算をやれば後者の"sin関数に対する積分方程式"を証明することができるこうして得られた冪級数を使って三角関数表を作るためには、円周率を、良い精度で計算しておく必要がある。マーダヴァはを証明したと考えられているけど、そのまま計算すると収束は遅い。上の解説によれば、『ユクティバーサ』には、初等的変形で得られる、もう少し効率のいい級数が書いてあるらしい。具体的にはなので、とおけばという級数に変形できる。選び方によっては、級数を有限項で打ち切った時の誤差を小さくすることができる。適用例としてを使ってなどが得られる。『ユクティバーサ』には、よりよい補正項としてと同値な式が書かれているらしい。誤差項の連分数展開を途中で打ち切って得られる。Euler数とする時、形式的べき級数間の等式として、以下が成立し、とおけば、この形式的べき級数と、上記誤差項との関係は明らか。Leibniz級数の誤差項は、良い漸近展開を持つ。この漸近展開は、例えば、以下の論文に載っている。Revisitedhttps://doi.org/10.1080/00029890.1997.11990646誤差項は、のような形でも評価でき、実際上は、どれを使っても、大差なさそうに思う。この級数の証明はに対して、部分積分によってを示せばいい。級数の収束も証明は難しくない。a=0の時を得る(この級数は、少なくともEulerによって知られていたようである)。この級数の最初の10項を評価すると、1.5702...を得るので、円周率が、3.14より大きいことが分かる。a=2として、最初の6項を足したものをライプニッツ級数の初項のみの評価値である1から引くと、0.78570...で、円周率に対する評価としては、22/7より少しいい。外接正96角形による近似計算に勝つには、7項まで評価しないといけないけど、計算量は大分少なくて済む。ケーララ学派では、誤差項の式については、経験的に得られただけのものだった可能性が高く、また、マーダヴァが計算した円周率は、同時代のアル=カーシーの計算に精度では劣っているものの、多角形近似によらない点で、方法としては新しい。did."と書かれている。astronomershttp://adsabs.harvard.edu/full/1998BASI...26...11R予測精度の観点からすると、地動説か天動説かより、円軌道か楕円軌道かの方が重要ではある。ケーララ学派について書かれた資料を読んでいても、観測の話が出てこないけど、彼らが、観測を重んじなかったのかはわからない。EUROPEhttp://ckraju.net/Joseph/PA-3-Manchester-2007-paper.pdfで、色々と考察がされているけど、決定的な証拠はない。アラビア→ヨーロッパヨーロッパでは、カスティーリャ王国が、1085年に、トレドを征服し、トレド翻訳学派と呼ばれる集団が形成され、アラビア語文献のラテン語翻訳(後には、カスティーリャ語翻訳)が行われるようになった。1000年頃のイスラム世界の中心地は、コルドバ、バクダード、カイロにあった。とある推計によれば、この三都市の中で、1200年頃までに人口を伸ばしたのはカイロのみで、コルドバは、1236年に、カスティーリャ王国に征服された。神聖ローマ皇帝フリードリヒ2世(1194〜1250)は、プトレマイオスの「アルマゲスト」のラテン語訳を支援し、フリードリヒ2世の宮殿には、イスラムの学者も多く滞在していたと言われる。この頃のヨーロッパの詳細は、それなりに分かっているようだけど、基礎科学に関して、他の地域で見られなかった発展を開始するのは、概ね、1600年頃からと言える。バグダードも、1258年に、モンゴル帝国に征服され、中央アジア〜東アジアでも、モンゴルの継承政権の元で、イスラム科学は暫く存続していたようだけど、サマルカンドのウルグ・ベク天文台が、1449年に破壊されて以降は、東アジア寄りの地域でのイスラム科学は消滅していった。Qushjiは、ウルグ・ベクの死後、1470年には、イスタンブールに移り、1474年に、そこで死んだそうである。100年近く後のイスタンブールでは、天文台が建設されたものの数年で破壊されたとなっている。al-Dinと綴られてることもある。ティコ・ブラーエ(1546〜1601)の同時代人なので、よく比較されるらしい。この頃でも、イスラム科学は、まだヨーロッパ科学の後塵を拝していたわけではない。ad-jinの足跡を見ると、当時のイスラム世界の学術中心は、イスタンブール、カイロ、ダマスカスなどにあったと推測される。13〜15世紀頃のマラーゲ学派やウルグ・ベク天文台で得られた成果も継承されており、それは、ヨーロッパにも伝わったようである。この頃、コペルニクスの弟子であるレティクス(1514-1574)は、10秒(1/360度)刻みの三角関数表を作ろうと試みたらしい。Ceulenという人が、正多角形近似に基づいて、1590年代に、3.14159265358979323846...まで正しく計算した。libcのmath.hlibc/include/math.hhttps://android.googlesource.com/platform/bionic/+/refs/heads/master/libc/include/math.hなどでは、M_PIが、3.14159265368979323846で定義されている。この精度の近似値が歴史上初めて得られた。Ceulenは、計算を進めて35桁目まで評価を行ったそうまた、1600年前後のヨーロッパでは、天文学以外に、航海術や測量術のような分野での三角法の利用も一般的になってきた。17世紀後半には、ヨーロッパで様々な関数の冪級数展開が与えられた。Infinitorumには、x^nの積分や、項別積分の考え方が書かれていて、実質的に√(1-x^2)のべき級数展開も知っていたようだ。人類に知られていた最古のべき級数は、1/(1+x)のものだろうけど、項別積分すると、log(1+x)のべき級数展開を得られる。Mercatorが"Logarithmo-technia"で書いたとされている。1/xの積分が、対数であることは、メルカトル以前に知られていた。※)メルカトル図法の"メルカトル"とは別人。Mercator(1512~1594)という人に由来する。Briggsが常用対数を提案した。Bürgiが対数の発見者と言われるので、その成果かもしれない。場合が、未解決として残った。Saint-Vincentは、正数a,b,c,dに対して、現代の記号で書くとを示したらしい。今では、積分の変数変換をするだけで、すぐ分かる。なって、対数と同じ性質を持つことを認識したらしい。当時は、1/xの積分によって自然対数が定義され、幾何級数に項別積分を適用することで、log(1+x)の級数展開を得たのだろう。このべき級数を得た人は、他にもいたかもしれないけど、現在では、メルカトル級数と呼ばれることがある。Newtonは、Wallisの結果に触発されて、一般二項級数から始めて、逆三角関数の級数展開に到達したと信じられている。形式的に逆関数を計算して、三角関数のべき級数展開を得たようだ。infinitas"を書き、この手稿は当時の学者の間で回覧されたらしい。infinitashttp://www.newtonproject.ox.ac.uk/catalogue/record/NATP00204が、その内容っぽい。ラテン語だけど、それほど長くはない。Google翻訳で訳すと、それなりの英語になるけど、読むのが面倒なので、書かれている係数だけ雰囲気で見ると、1/(1+x)や√(1+x)の級数展開らしきものが前半にあり、1/(1+x^2)と、その積分なども級数の形で書かれている。後半の方には、arcsinやsinの級数展開らしきものも、書かれている。関係を、どのくらい理解していたのか、イマイチ分からない。Gregoryは、メルカトル級数の結果を見て、1/(1+x^2)の積分から、arctan(x)の級数展開を発見したとされ、現在ではグレゴリー級数と呼ばれているけど、この級数自体は、Newtonの手稿にも見られる。Newtonが、この級数が、arctan(x)に等しいことを理解してなくて、Gregoryが、それに気付いたということなら、グレゴリー級数という名称は妥当だと思われる。Oldenburgを通じて、Leibnitzにも伝わったとされる。Sequence)"と呼ばれる二重数列が収束することを証明したらしい。収束先がarctan(x)を使って書けることは分かっていただろうし、この二重数列が満たす漸化式を使って、arctan(x)のべき級数展開を見つけることも不可能ではなかったはずだけど、Gregoryは、そういう道は、辿らなかったようだ。いずれにせよ、ヨーロッパでは、1670年前後に、三角関数や逆三角関数や対数関数のべき級数展開が知られるようになった。現代の標準的な微積分学のカリキュラムでは、微分→積分の順で習うけど、インドでもヨーロッパでも、微分法が発展する以前に、まず積分と無限級数の計算が先に発展したようである【補足】メルカトル級数が知られる以前に、自然対数の数値計算に興味を持つ人がいたかどうかは分からないけど、常用対数表が存在してたわけだから、x=1以外のどこか一点で、自然対数の値を決めれば、底の変換によって、自然対数の計算を行うことはできたはず。自然対数の値を計算する方法として、(1)数値積分(2)交代調和級数の利用(3)アルキメデスの漸化式などが考えられる。Wrightが、緯線距離を計算するために、正割関数の積分を、数値的に計算したという例がある。Gregoryが公表したようなので、Wrightの計算結果を流用したことはないだろうけど、数値積分の発想自体は存在したようである。Mengoliが発見したらしい。elementaという本を出版しているので、これに含まれるのかもしれないけど、400ページ以上あるので、中を探す気は起きない。Godementの解析学の教科書には、Mengoliは、積分の長方形近似から得られる式と、容易に得られる等式から、交代調和級数の値に至ったと書いてある。この初等的な議論は、1650年代のヨーロッパでも可能だっただろうし、また、Mengoliが、メルカトル級数を発見しなかったことも納得がいく。hyperbolahttps://royalsocietypublishing.org/doi/10.1098/rstl.1668.0009で、交代調和級数の初等変形で得られる級数を3つ書いている。Brounckerの細かい議論までは見てないけど、級数と収束値はとなる。最後の級数は、比較的少ない計算量で、それなりの評価を与える。Iを見ると、おそらく、積分から、この式を得たのだろう痕跡がある)があって、類似の等式を導いている。この続きには、同じような公式が、いくつか記載されている(Brounckerの級数は載ってない)。Ramanujanも、最初は、簡単なことから始めたということが分かる。VIIIのN13や14)。これは、対数関数やarctanを含むとはいえ、有理関数の積分だから、今となっては別に面白いものではない(n,mを固定すれば、Mathematicaでも計算できる)けど、Gaussの日記でも触れられている(Entry54)ので、賢い人たちも、色々な一般化を試行錯誤したらしいことが窺える(Gaussは、単に、Eulerの著作で、この積分を見つけたのかもしれない)(3)17世紀でも、自然対数が、"アルキメデスの漸化式"で計算できることに気付ける可能性はあったと思うけど、多分、誰もやってないっぽい。Arctangentsという論文が書かれている(調和平均の代わりに算術平均を使ってるけど、同じこと)。17世紀だと一般項が書けなくて、(M(a,b)で、初期値(a,b)に対するアルキメデスの漸化式の収束値を表すとする時)、漸化式から導ける関係式のみだと、収束先を定数倍を除いて決められず、"x=1以外に、最低でも、もう一点の値を決めれば..."という問題に戻ってしまう気もする。ヨーロッパ→東アジア中国では、イエズス会の宣教師らによって、1631〜34年頃に、『崇禎暦書』が編纂され、これにも、割円八線表という名前で、三角関数表が含まれていたらしい(sin,cos,tan以外に、1/sin,1/cos,cot,versine,coversineの合計8つの三角関数があったことが、八線表という名前の由来のよう)。前者は、授時暦を研究した成果、後者は、回回暦を研究した成果で、清を宗主国としていた時は、独自の暦と清の暦(時憲暦)を併用したらしい。(朱世傑、1299年)が教科書として使われて残っていた。この本には、天元術と呼ばれる初等代数が含まれ、初期の和算家が受容・発展させたとされる。秀吉が、李氏朝鮮を侵略した時に、技術者も連行したらしいから、初期の和算家が、彼らから、数学を学んだことはありそうに思うけど、そういう話は聞かない。医者の岡野井玄貞は、1643年の朝鮮通信使で来日した朴安期(容螺山)に暦学に関する問答を行ったと伝えられている。岡野井玄貞は、日本で最初に作られた暦(貞享暦)の作成者である渋川春海が暦学を学んだ一人でもある現代の韓国は、キリスト教徒が多いことで知られ、これは李氏朝鮮の時代に、中国経由で伝来したものだけど、宣教師の入国などは許可されていなかった。一方で、1631年に、マテオ・リッチのキリスト教に関する著作『天主実義』が輸入されたらしいから、書籍の輸入に関しては、それほど制限されてなかったのかもしれない。中国の暦書などを介して、ヨーロッパの数学や天文学が伝わることもあったかもしれないけど、詳細は分からない江戸時代前半の日本では、鎖国と共に、ヨーロッパの書物は、中国で書かれたものも含めて、全面的に輸入が禁止されていたので、日本には、1723年に清で刊行された『曆算全書』が、1726年に輸入された時に、三角法が伝わったとされている。三角関数表の伝来と二つの経路https://repository.kulib.kyoto-u.ac.jp/dspace/handle/2433/232876また、江戸時代初期には、イエズス会士カルロ・スピノラやウィリアム・アダムスのような人が滞在していて、数学や天文学を教えたと伝えられているから、その時にも、ヨーロッパ数学伝達の機会は存在した。いえ、17世紀の日本に三角関数表が存在していたという物証はないっぽい。成立年不明の『算暦雑考』という本は、内容などから、建部賢弘(1664〜1739)の著作と推測されていて、日本最初の三角関数表を含むとされている。成立年は1727年頃とされている。成立年不明なので、どっちが先か分からない無限級数は、18世紀前半に、何人かの和算家が使うようになった。(1740):1/(1-x),√(1-x)の級数展開arcsin(x)^2のべき級数は、1737年にEulerが証明したとされている。建部賢弘は、証明や成立理由を知っていたわけでなく、数値計算の結果から予想しただけらしい。Jartoux(1669〜1720)が中国に、いくつかのべき級数を伝えたとされ、それが日本にも伝わった可能性がある。Jartouxは、Leibnizと文通していたらしいから、高い学術能力を持った人だったようである。Mingda)のような人たちが、無限級数の研究を継続していたらしい。項明達は、楕円の周長と長径の比を、離心率の無限級数で書けることを発見していたようだ。同様の結果は、ヨーロッパでは、1742年に、マクローリンが公表したとされている。中国は、楕円の周長計算によって、微妙に、楕円積分に踏み込んでいたようだけど、一方、18世紀ヨーロッパで、楕円積分の研究が盛り上がった理由は、平面弾性曲線にあったっぽい。historyhttps://www2.eecs.berkeley.edu/Pubs/TechRpts/2008/EECS-2008-103.htmlによれば、弧長を独立変数にして弾性曲線の微分方程式を書き、振り子の運動方程式と等価になる(つまり、平面弾性曲線が可積分系で記述される)ことに気付いた人は、Kirchhoffらしい(1859)。acceptihttps://archive.org/details/methodusinvenie00euleだそうだけど、300ページ以上あるので、確認する気にはならない。topなども発見されたのだろうし、これらの解を記述しようと思って、楕円関数に辿り着いた人がいないのは不思議に思える。cf)弾性曲線と可積分系の関係については、以下の解説が詳しい。天動説であれ、地動説であれ、昔の人が最初に考えたモデルは、等速円運動だったらしい。等速円運動モデルが観測と合わないことは、割とすぐ認識されて、色々と複雑なモデルが作られることになった。等速円運動モデルは、今から見ると二次元調和振動子の特殊解になっている。Keplerが楕円運動を提案した時、ガリレイ(やデカルト?)の反対があったらしいので、円運動への拘りは天動説よりも強かったのだろう。corpshttps://projecteuclid.org/euclid.acta/1485887161Levi-Civitaの目的は三体問題だったので、特に明記されてないけど、この変数を使うと、Kepler系は調和振動子に帰着する。思うところあって少し調べた。»はてなブログをはじめよう!m-a-oさんは、はてなブログを使っています。ブログを報告するスターをつけました引用をストックしましたストック一覧を見る閉じる引用するにはまずログインしてくださいログイン閉じる引用をストックできませんでした。再度お試しください閉じる限定公開記事のため引用できません。