【数理考古学】とある円周率への挑戦？

* (2021年2月10日)初投稿
* (2021年4月23日)Tex導入


古来**数学者達**は、それがx=1の時に$\frac{π}{2}$で収束して**円周率**に迫るので**アークタンジェント**(ArcTangent)**近似式**の高速化に血道を上げてきたのです。
[円周率の公式と計算法](http://www.kurims.kyoto-u.ac.jp/~ooura/pi04.pdf)

#ビエトの無限乗積式(<font color=blue>Vieteian Infinite Product Formula</font>)

そうした**円周率近似演算**の原風景として「(初めて既知数を記号化し記号代数の原理と方法を確立して当時の代数学を体系化した事から)**代数学の父**」と呼ばれる**フランソワ・ビエト**(François Viète、1540年~1603年）の**無限乗積式**(Infinite Product Formula)があり、それ自体は「最初のk項までの積=**半径**1の円に内接する正2k+1角形の面積πへの収束」あるいは「半円に内接する半多角形の周長への収束)」を求める内容だったのです。

```R
#π=2×2/sqrt(2)×2/sqrt(2+sqrt(2))×2/sqrt(2+sqrt(2+sqrt(2)))…
> 2*2/sqrt(2)
[1] 2.828427
> 2*2/sqrt(2)*2/sqrt(2+sqrt(2))
[1] 3.061467
> 2*2/sqrt(2)*2/sqrt(2+sqrt(2))*2/sqrt(2+sqrt(2+sqrt(2)))
[1] 3.121445
>2*2/sqrt(2)*2/sqrt(2+sqrt(2))*2/sqrt(2+sqrt(2+sqrt(2)))*2/sqrt(2+sqrt(2+sqrt(2+sqrt(2))))
[1] 3.136548
>2*2/sqrt(2)*2/sqrt(2+sqrt(2))*2/sqrt(2+sqrt(2+sqrt(2)))*2/sqrt(2+sqrt(2+sqrt(2+sqrt(2))))*2/sqrt(2+sqrt(2+sqrt(2+sqrt(2+sqrt(2)))))
[1] 3.140331
#ところで√2/2(0.7071068)は√2/(√2^2)。
#なので分子分母共に√2乗すると1/√2。
#すなわち2×2/sqrt(2)(2.828427)=2×sqrt(2)となる。
[1] 2.828427
> sqrt(2)/2
[1] 0.7071068
> 1/sqrt(2)
[1] 0.7071068
> 2*sqrt(2)
[1] 2.828427
> 2*2/sqrt(2)
[1] 2.828427
#当然2/sqrt(2)=sqrt(2)である。
> 2/sqrt(2)
[1] 1.414214
> sqrt(2)
[1] 1.414214
#すなわち
> 2*sqrt(2)
[1] 2.828427
> 2*sqrt(2)*2/sqrt(2+sqrt(2))
[1] 3.061467
#つまりこの計算は半円に内接する半多角形の周長を求めているのに等しい。
```
半円に内接する半多角形の周長計算と考えた場合の「**ビエトの円近似**(<font color=blue>Vieteian Circle Approximation</font>)」
![Rplot23.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/ba95c71d-093f-e1f1-eb0d-035ceee8dbe4.png)

```R
f0<-function(x) sqrt(1-(x-1)^2)
plot(f0,asp=1,xlim=c(0,2),ylim=c(0,1),main="Vieteian Circle Approximation",xlab="Real",ylab="Imaginal")
segments(0,0,2,0,col=rgb(1,0,0),lwd=2)
#ここでおもむろに「ガウスの巡回群による三角関数演算」概念を導入。
par(new=T)#上書き指定
c0=seq(0,pi,length=3)
cx0=cos(c0)+1
cy0=sin(c0)
plot(cx0,cy0,type="l",asp=1,xlim=c(0,2),ylim=c(0,1),main="",xlab="",ylab="",col=rgb(0,0,1),lwd=2)

par(new=T)#上書き指定
c1=seq(0,pi,length=5)
cx1=cos(c1)+1
cy1=sin(c1)
plot(cx1,cy1,type="l",asp=1,xlim=c(0,2),ylim=c(0,1),main="",xlab="",ylab="",col=rgb(0,1,0),lwd=2)

legend("bottomright", legend=c("2","2×2/√2=2.83","2×2/√2×2/√(2+√2))=3.06"),lty=c(1,1,1,1),col=c(rgb(1,0,0),rgb(0,0,1),rgb(0,1,0)),cex=3/4)
```

* 「**往復**(<font color=blue>1周</font>)」でなく「**片道**(<font color=blue>半周</font>)」しか見ないので、私の投稿では<font color=green>1の冪乗根</font>**を求める演算**<font color=purple>z^n=1</font>の結果集合たる**ガウスの巡回群**を用いて計算を省力化した「(<font color=blue>円を挟む内接多角形と外接多角形の辺数をひたすら増やし続けながら辺長計の差を求める</font>)**挟み撃ち定理**（<font color=blue>Squeeze Theorem</font>）」の半分の計算範囲。<br>[【初心者向け】挟み撃ち定理による円周率πの近似](https://qiita.com/ochimusha01/items/1caff7d8af4aab167aaa)<br>![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fqiita-image-store.s3.ap-northeast-1.amazonaws.com%2F0%2F597772%2F91fd2787-e634-9a4c-64a8-5616bbd12d65.gif?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&w=1400&fit=max&s=c20efc397400fba9c7dfb6c62a4fb04b)
* 一方**ユークリッド計量**(<font color=purple>Euclidean metric</font>)概念を用いた表現では「**xy軸上の直線距離**{<font color=green>0,0</font>}-{<font color=green>2,0</font>}<font color=purple>の</font>**中心**(<font color=green>0,1</font>)に</font>(<font color=purple>ピタゴラスの定理$y=\sqrt{1-x^2}$に従って</font>)<font color=purple>等距離の経過点を置いて行く</font>」イメージとなるが、これを一般化した**計量関数**(<font color=blue>Metric Function</font>)は通常**xy軸上の直線距離**{<font color=green>1,0</font>}{<font color=green>0,1</font>}(<font color=blue>0→π/2=90度の範囲</font>)しか見ないのでその倍の計算範囲。<br>[【初心者向け】方形描画関数②距離空間との関係。](https://qiita.com/ochimusha01/items/06e49a35dffbd12b35aa)<br>![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/155f93ea-e61a-3510-8480-04523fc598d8.png)

こうした**視座**(<font color=blue>Perspective</font>)の概念は、以下の各次元上の制約に対応します。

* **点**=**同心円空間**(<font color=blue>Dot=Concentric Circle Space</font>)<font color=green>0-0.5</font>**次元**の広がりを有する<font color=green>$2^0$</font>**座標系**の場合は**全周範囲**(<font color=blue>2π Range</font>)。<br>[【Rで球面幾何学】等差数列(算術数列)②数直線概念から同心円集合概念へ](https://qiita.com/ochimusha01/items/2a16c3cba3f88986c767)<br>![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fcdn-ak.f.st-hatena.com%2Fimages%2Ffotolife%2Fo%2Fochimusha01%2F20201111%2F20201111061314.gif?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&w=1400&fit=max&s=0ec5ad367a9238d99fae4325c3a95b29)<br>![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fqiita-user-contents.imgix.net%2Fhttps%253A%252F%252Fcdn-ak.f.st-hatena.com%252Fimages%252Ffotolife%252Fo%252Fochimusha01%252F20201202%252F20201202093758.gif%3Fixlib%3Drb-1.2.2%26auto%3Dformat%26gif-q%3D60%26q%3D75%26s%3D167e614307e1813ff07f89cc7e325745?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&w=1400&fit=max&s=eebf5744bdd7ab1d7b98a49c133d54c4)
* **線**=**円柱**/**トーラス空間**(<font color=blue>line=Cylinder/Torus Space</font>)<font color=green>1-1.5</font>**次元**の広がりを有する<font color=green>$2^1$</font>**座標系**の場合は**半周範囲**(<font color=blue>π Range</font>)。<br>[【初心者向け】群論概念①基本定義](https://qiita.com/ochimusha01/items/208dc0e68db64805e4f7)<br>![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fqiita-user-contents.imgix.net%2Fhttps%253A%252F%252Fcdn-ak.f.st-hatena.com%252Fimages%252Ffotolife%252Fo%252Fochimusha01%252F20201202%252F20201202110929.gif%3Fixlib%3Drb-1.2.2%26auto%3Dformat%26gif-q%3D60%26q%3D75%26w%3D1400%26fit%3Dmax%26s%3D2707333ff7b779922e55b3d65cdd454c?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&w=1400&fit=max&s=ddabd89bda0b4aebe5548144265f242f)<br>![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fcdn-ak.f.st-hatena.com%2Fimages%2Ffotolife%2Fo%2Fochimusha01%2F20201202%2F20201202095654.gif?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&w=1400&fit=max&s=b16b4ae7d0364d94bc8506e883ad1132)<br>[【無限遠点を巡る数理】等差数列と等比数列②基底関数概念の導入](https://qiita.com/ochimusha01/items/df8c006439794da89684)<br>![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fqiita-user-contents.imgix.net%2Fhttps%253A%252F%252Fcdn-ak.f.st-hatena.com%252Fimages%252Ffotolife%252Fo%252Fochimusha01%252F20201130%252F20201130011306.gif%3Fixlib%3Drb-1.2.2%26auto%3Dformat%26gif-q%3D60%26q%3D75%26w%3D1400%26fit%3Dmax%26s%3De8c7db6fc0efc1bff4f238b59b3aec23?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&w=1400&fit=max&s=98fd000c088d9e919cbf9d3eb0c1230b)<br>![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/9727379b-88b6-d336-8506-f211ff3a5e46.png)
* **デカルト座標系**=**円錐空間**(<font color=blue>Cartesian Coordinate System=Cone Space</font>)<font color=green>2-3.5</font>**次元**の広がりを有する<font color=green>$2^2$</font>/<font color=green>$2^3$</font>**座標系**の場合は<font color=green>$\frac{1}{4}$</font>**周範囲**(<font color=blue>$\frac{π}{2} $Range</font>)/<font color=green>$\frac{1}{8}$</font>**周範囲**(<font color=blue>$\frac{π}{4}$Range</font>)。<br>[【初心者向け】線形関数や絶対値関数](https://qiita.com/ochimusha01/items/8b7050e65b599f9e15c2)<br>![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fcdn-ak.f.st-hatena.com%2Fimages%2Ffotolife%2Fo%2Fochimusha01%2F20201120%2F20201120124656.gif?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&w=1400&fit=max&s=63f478803a9bfd7a0cd59d7e4d79c0a7)<br>![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fcdn-ak.f.st-hatena.com%2Fimages%2Ffotolife%2Fo%2Fochimusha01%2F20201120%2F20201120125126.gif?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&w=1400&fit=max&s=b4c10d7a3a0a42144baa87b2651b11ae)<br>![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fcdn-ak.f.st-hatena.com%2Fimages%2Ffotolife%2Fo%2Fochimusha01%2F20201120%2F20201120124110.gif?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&w=1400&fit=max&s=3837216a67faea20b81bdf603119c5e8)
* ちなみに**四元数**(<font color=blue>Quaternion</font>)は<font color=green>$2^3$</font>**座標系**とは<font color=green>$2^1$</font>**座標系**からの拡張のされ方が異なる<font color=green>$3^2$</font>**座標系**に分類される(<font color=blue>それぞれが複素数平面を備える二つの回転軸の組み合わせにより全周範囲視座を回復</font>)。<br>[【Rで球面幾何学】ハミルトンの四元数は何を表しているのか？](https://qiita.com/ochimusha01/items/aff8d52f20af01d13bed)<br>[【オイラーの多面体定理と正多面体】正方形における平方対角線と立方体における立方対角線の関係について。](https://qiita.com/ochimusha01/items/9d899779ff8e814beb41)<br>[【Rで球面幾何学】そもそも「内積」とは「外積」とは何か？](https://qiita.com/ochimusha01/items/b34f73d13eef785365cd)<br>![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fcdn-ak.f.st-hatena.com%2Fimages%2Ffotolife%2Fo%2Fochimusha01%2F20201121%2F20201121224728.gif?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&w=1400&fit=max&s=d3218ea17a11411eadc10c33c2ab44f2)<br>![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fcdn-ak.f.st-hatena.com%2Fimages%2Ffotolife%2Fo%2Fochimusha01%2F20201121%2F20201121225859.gif?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&w=1400&fit=max&s=c449f12eea39f9baefef69b3f61866e0)<br>![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fcdn-ak.f.st-hatena.com%2Fimages%2Ffotolife%2Fo%2Fochimusha01%2F20201121%2F20201121230355.gif?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&w=1400&fit=max&s=01d7259857a7d9a3cf58393bea0311a1)<br>人間にとって<font color=green>$2^3$</font>**座標系**より<font color=green>$3^2$</font>**座標系**の方が難しく感じられるのは、**生物が進化過程で獲得してきた視覚の補完機能**に関わってくる話である。<br>[【無限遠点を巡る数理】「ピタゴラスの定理」と「デカルト座標系」の狭間](https://qiita.com/ochimusha01/items/e27e9cd6aad824043a59)<br>こういう話は具体例で考えた方が話が早い。例えば以下のアニメーションが右回りにも左回りにも見えるのが「**視覚的補完**」…<br>[【初心者向け】複素共役（Complex Conjugate）のアニメーション表示について。](https://qiita.com/ochimusha01/items/e0bb2789a420d3317cf1)<br>![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fqiita-user-contents.imgix.net%2Fhttps%253A%252F%252Fcdn-ak.f.st-hatena.com%252Fimages%252Ffotolife%252Fo%252Fochimusha01%252F20201116%252F20201116230848.gif%3Fixlib%3Drb-1.2.2%26auto%3Dformat%26gif-q%3D60%26q%3D75%26w%3D1400%26fit%3Dmax%26s%3D0489492d3b778bcf1b6c8b3886af31c4?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&w=1400&fit=max&s=b1273a9440db7db91f80a947bffa0204)
<br>![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fcdn-ak.f.st-hatena.com%2Fimages%2Ffotolife%2Fo%2Fochimusha01%2F20201113%2F20201113204552.gif?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&w=1400&fit=max&s=b46be23e211ab7df85a44ed1887614b9)
<br>![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fcdn-ak.f.st-hatena.com%2Fimages%2Ffotolife%2Fo%2Fochimusha01%2F20201113%2F20201113205243.gif?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&w=1400&fit=max&s=89f7c0b4a4df01f04f8e790faa13df98)<br>以下の図の「**球面の裏側まで映り込んでいる周辺部**」にどうしても焦点を合わせられないのが「**視覚的切り捨て**」となる。<br>[【初心者向け】誤差関数(ERF)と相補誤差関数 (ERFC)。](https://qiita.com/ochimusha01/items/ee8b32e4b2cc7ce4acb0)<br>![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fcdn-ak.f.st-hatena.com%2Fimages%2Ffotolife%2Fo%2Fochimusha01%2F20201113%2F20201113111114.gif?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&w=1400&fit=max&s=1be74c11a381d9d345d61b6c2be72199)
* いずれにせよ<font color=green>$2^3$</font>**座標系**も<font color=green>3^2</font>**座標系**も**素数族**<font color=green>2^n3^m</font>の部分集合に過ぎず、**平面上**の**正六角形**(<font color=blue>Regular Hexagon</font>)による**平面充填**(<font color=blue>Tiling</font>)をその極限とする。<br>![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/ac398d35-a735-e614-1367-ede5e1aa6b7c.png)

しばしば「<font color=purple>収束が遅く実用的でない</font>」なる一言で片付けられてしまう**ビエトの円近似法**ですが**挟み撃ち定理の特殊例**に過ぎないとはいえ「(<font color=blue>sqrt(2)を含む</font>)**素数族**<font color=green>2^n</font>のみ導入された環境下でも<font color=green>π</font>(<font color=purple>3.141593</font>)へは到達可能である事を示す意味合いにおいて重要です(<font color=blue>sqrt(2)の逆数たる2^2=4すら登場しない</font>)。

```R
#ビエトの円近似には外接多角形の半周の計算に対応する逆算もありそう？
#現段階では未達…

#正方形
>f0(4)
[1] 4
> 2*(2*(2-1))
[1] 4

#正八角形
> f0(8)
[1] 3.313708
> 3.313708/4
[1] 0.828427
> 0.828427/2
[1] 0.4142135
> sqrt(2)-1
[1] 0.4142136
> 2*(2*(2-1))*(2*(sqrt(2)-1))
[1] 3.313708

#正16角形
> f0(16)
[1] 3.182598
> 3.182598/3.313708
[1] 0.9604341
> 0.9604341/2
[1] 0.480217
sqrt(x)-1
> 1.480217^2
[1] 2.191042
[1] 2.191042
> 2*(2*(2-1))*(2*(sqrt(2)-1))*(2*(sqrt(2.191042)-1))
[1] 3.182597

#現段階では、これが精一杯…
```

* ちなみに**ビエトの無限乗積式**(<font color=blue>Vieteian Infinite Product Formula</font>)に拠る演算結果集合は、<font color=green>0</font>あるいは<font color=green>2</font>を**単位元**と定め、<font color=green>π</font>に収束する**元**に<font color=green>-1</font>を掛け、<font color=green>2</font>に収束する演算結果集合を**逆元**と設定すると**結合条件**も成立し**群定義**の少なくとも一部を満たす((<font color=blue>興味深い事に単位元を0に定めた時は2,2に定めた時は0が極限値となるので両者が単位元として並び立つ事はない</font>)<br>[【初心者向け】群論概念①基本定義](https://qiita.com/ochimusha01/items/208dc0e68db64805e4f7)
* 問題は演算結果が閉じない事だが(<font color=blue>そもそも近似過程で現れる全ての数字を含めた元をイメージする事自体がナンセンス?</font>)、同じ問題なら**式形**<font color=purple>(1+1/N)^N</font>の時に**ネイピア数**<font color=green>exp(1)</font>(<font color=blue>2.718282</font>)、**式形**<font color=purple>(1-1/N)^N</font>の時に**ネイピア数の逆数**<font color=purple>exp(-1)</font>(<font color=blue>0.3678794</font>)、**式形**<font color=purple>(1±θi/N)^N</font>の時に**単位元**<font color=green>0</font>/**極限**<font color=green>-1</font>(<font color=blue>±π</font>)へと収束する**自然指数関数**(<font color=blue>Natural Exponential Function</font>)<font color=purple>e^x</font>も抱えており、こういう場合には**実数全体**を元に採用するしかない(<font color=blue>ちなみに私の投稿では実数を「10進方法の演算結果集合」と規定し、無理数もとりあえず計算可能範囲まで計算した結果をその演算結果集合に「仮に」加えている。コンピューター時代ならではの対「ゴルディアスの結び目(Gordian Knot)」解決法?</font>)。<br>[【初心者向け】指数・対数関数の発見とそれ以降の発展について。](https://qiita.com/ochimusha01/items/3e28f4f06abd8cef0af1)

ところで、ここまで**ビエトの無限乗積式**への関心を掻き立てておいて何なのですが、残念なお知らせがあります。実はこの演算、**オイラーの公式**(<font color=blue>Eulerian Formula,奇しくも著名なe^θi=Cos(θ)+Sin(θ)iと同名</font>)<font color=purple>Sin(x)/x</font>(<font color=blue>Sinc関数</font>)=<font color=purple>cos(x/2)×cos(x/4)×cos(x/6)×cos(x/8)×cos(x/16)…</font>に(<font color=blue>(0±1i)^2x→e^πi変換の様に</font>)<font color=green>x=π/2</font>を代入した特殊ケースに過ぎなかったのです。
[ヴィエトの無限積の公式](https://mathtrain.jp/viete)
[sinx/xの極限は？x→0とx→∞の場合を証明付きで東大医学部生が教えます！](https://todai-counseling.com/?p=461)
![image.gif](https://cdn-ak.f.st-hatena.com/images/fotolife/o/ochimusha01/20210205/20210205013804.gif)

```R
Macrolin_expansion_exp<-function(x){

#グラフのスケール決定
Graph_scale_x<-c(-40,40)
Graph_scale_y<-c(-1,1)
#関数と凡例の決定
switch(x,
"0"= f0<-function(x) {cos(x/2)},
"1"= f0<-function(x) {cos(x/2)*cos(x/4)},
"2"= f0<-function(x) {cos(x/2)*cos(x/4)*cos(x/8)},
"3"= f0<-function(x) {cos(x/2)*cos(x/4)*cos(x/8)*cos(x/16)},
"4"= f0<-function(x) {cos(x/2)*cos(x/4)*cos(x/8)*cos(x/16)*cos(x/32)},
"5"= f0<-function(x) {cos(x/2)*cos(x/4)*cos(x/8)*cos(x/16)*cos(x/32)*cos(x/64)}
)

switch(x,
"0"=text<-c("cos(x/2)"),
"1"=text<-c("cos(x/2)*cos(x/4)"),
"2"=text<-c("cos(x/2)*cos(x/4)*cos(x/8)"),
"3"=text<-c("cos(x/2)*cos(x/4)*cos(x/8)*cos(x/16)"),
"4"=text<-c("cos(x/2)*cos(x/4)*cos(x/8)*cos(x/16)*cos(x/32)"),
"5"=text<-c("cos(x/2)*cos(x/4)*cos(x/8)*cos(x/16)*cos(x/32)*cos(x/64)"),
)

f1=function(x){sin(x)/x}

plot(f1,xlim=Graph_scale_x,ylim=Graph_scale_y,type="l",col=rgb(0,0,0),main="Macrolin expansion", xlab="X", ylab="Y")
par(new=T)#上書き指定
plot(f0,xlim=Graph_scale_x,ylim=Graph_scale_y,type="l",col=rgb(1,0,0),main="", xlab="", ylab="")

legend("bottomleft", legend=c("e^x",text),lty=c(1,1),col=c(rgb(0,0,0),rgb(1,0,0)))
}

library("animation")
Time_Code=c("0","0","0","1","1","1","2","2","2","3","3","3","4","4","4","5","5","5")
saveGIF({
for (i in Time_Code){
　　Macrolin_expansion_exp(i)
}
}, interval = 0.1, movie.name = "sync01.gif") 
```
それってどういう事? また歴史的価値しかない式と思いきや、**ロジスティック方程式**/**写像**との絡みがあって大学受験数学の出題範囲にもなってる模様…
[ロジスティック写像と漸化式](https://mathtrain.jp/logistic)
[【初心者向け】ロジスティック方程式とその関連範囲](https://qiita.com/ochimusha01/items/2f0118ecda2e17d366ba)
[ロジスティック写像 - Wikipedia](https://ja.wikipedia.org/wiki/%E3%83%AD%E3%82%B8%E3%82%B9%E3%83%86%E3%82%A3%E3%83%83%E3%82%AF%E5%86%99%E5%83%8F)
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/1ad00292-fe0a-a00e-366d-4681ebf39d6c.png)

#オイラー変換(<font color=blue>Eulerian Transformation</font>)

**オイラー**は**アークタンジェント**(<font color=blue>Arctangent</font>)の求め方の高速化も図っています。具体的には<font color=green>1755年</font>に<font color=green>1671年</font>にスコットランドの**ジェームス・グレゴリー**が発見した**グレゴリー級数**(<font color=blue>Gregory Series</font>)のアルゴリズムを大幅に改善した**オイラー変換**(<font color=blue>Eulerian Transformation</font>)を発表。
[円周率の公式と計算法](http://www.kurims.kyoto-u.ac.jp/~ooura/pi04.pdf)
[グレゴリー級数](https://mathsuke.jp/gregory-series/)
[【Rで球面幾何学】オイラーの公式を導出したマクローリン級数の限界？](https://qiita.com/ochimusha01/items/4a484d5759f96abb0ab3)
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/c8a5c37a-ab3a-58d2-fdba-f5bdc9cc743e.png)

#ライプニッツ級数(<font color=blue>Leibniz series,1674年</font>)とグレゴリー級数(<font color=blue>Gregory series, 1671年</font>)とテイラー級数(<font color=blue>Taylor Series,1715年</font>)

[ライプニッツ級数(Leibniz series,1674年)](https://mathsuke.jp/leibniz-series/)

```math
\sum_{n}^{\infty}\frac{-1^n}{2n+1}=\frac{π}{4}
```
[グレゴリー級数(Gregory series, 1671年)](https://mathsuke.jp/gregory-series/)

```math
\sum_{n}^{\infty}\frac{-1^n}{2n+1}x^{2n+1}=\frac{π}{4}
```
x=1の時、**ライプニッツ級数**と一致する。

[テイラー展開 - Wikipedia](https://ja.wikipedia.org/wiki/%E3%83%86%E3%82%A4%E3%83%A9%E3%83%BC%E5%B1%95%E9%96%8B)

>数学において、**テイラー級数**(<font color=blue>Taylor series</font>)は関数のある一点での導関数たちの値から計算される項の無限和として関数を表したものである。そのような級数を得ることを**テイラー展開**という。
>
>**テイラー級数**の概念はスコットランドの数学者**ジェームズ・グレゴリー**により定式化され、フォーマルにはイギリスの数学者**ブルック・テイラー**によって<font color=green>1715年</font>に導入された。<font color=green>0</font>を中心としたテイラー級数は、**マクローリン級数**(<font color=blue>Maclaurin Series</font>)とも呼ばれる。これはスコットランドの数学者**コリン・マクローリン**にちなんでおり、彼は<font color=green>18世紀</font>に**テイラー級数**のこの特別な場合を積極的に活用した。
>
>関数はその**テイラー級数**の有限個の項を用いて近似することができる。**テイラーの定理**はそのような近似による誤差の定量的な評価を与える。テイラー級数の最初のいくつかの項として得られる多項式は**テイラー多項式**と呼ばれる。関数のテイラー級数は、その関数のテイラー多項式で次数を増やした極限が存在すればその極限である。関数はそのテイラー級数がすべての点で収束するときでさえもテイラー級数に等しいとは限らない。開区間（<font color=blue>あるいは複素平面の開円板</font>）でテイラー級数に等しい関数はその区間上の解析関数と呼ばれる。

なるほど、以下で私が悩んだ「<font color=purple>ベクトルの世界と関数列の世界をどう統合するか</font>」問題の起源はまさにその大源流にあった訳です(<font color=blue>ライプニッツの数列をグレゴリーが関数列化し、これにニュートンやライプニッツの手になる微積分概念が合流して「**関数列**」テイラー級数となった</font>)。
[【無限遠点を巡る数理】等差数列と等比数列③基底関数の線形結合?](https://qiita.com/ochimusha01/items/0c169537670f56591662)

そんな感じで以下続報…