【Python演算処理】特異点(Singularity)を巡る数理の自分なりのまとめ。

(2021年4月03日)投稿。
(2021年4月12日)回転方向を反時計回りで統一。

[【数理考古学】冪乗算(Exponentiation)の微積分](https://qiita.com/ochimusha01/items/0bd18370cfa96e83347a)
冪乗算(Exponentiation)の微積分定理は以下となります。

* $\frac{1}{β+1}x^{β+1}$を微分すると$x^β$となる。
* $x^β$を積分すると$\frac{1}{β+1}x^{β+1}+C$となる(<font color=blue>不定積分</font>)。
* ただしβ=-1の時は分母が0となり計算出来ない。

**等差数列**(Arithmetic Sequence=算術数列)や**等比数列**(Geometric Sequence=幾何数列)との関係からこの現象を追っていくと以下となります。
[【数理考古学】とある実数列の規定例②等比数列から乗法群へ](https://qiita.com/ochimusha01/items/58afced6ea8bebf391c7)

* 初項1、公差1の時に添字と演算結果集合が一致する等差数列は、公比1,公差0の時に公比1,公差1の等比数列と重なる(Y=1)。**ラップ級数**(Rap Series)$N_n(n=-\infty→0→)=(1,…,1,1,1,…,1)$を生成する演算たる1の冪乗$y=1^n$同様、Xの値によらずYの値が1に定まるのである。<br>[【Python演算処理】単位円筒を巡る数理](https://qiita.com/ochimusha01/items/fd8e51940b5e0e1a28c5)
* 等比数列は理論上初項0,公比0の時にy=0に到達するがコンピューターは負数の場合それをちゃんと計算しない($n\geqq0$の時は0、$n<0$の時はNaNを返す)。一方冪乗算$y=0^n$は正数$(n>0)$に対して0,0$(n=0)$に対して1、負数$(n<0)$の時に無限大$\infty$を返す。
* この反応のズレは$n\geqq0$で等比数列と冪乗算の添字が一つズレる事と密接な関係があると目される。<br>![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/c955e98d-de08-0b92-733b-e666dea09479.png)

この様な形で現れる到達不可能地点を数学の世界では**特異点**(Singularity)と呼びます。
[特異点 - Wikipedia](https://ja.wikipedia.org/wiki/%E7%89%B9%E7%95%B0%E7%82%B9)

#正負概念導入に伴う特異点の発生

**ピタゴラスの定理**$y=±\sqrt{1-x^2}$も半円ずつしか描けません($Y=0$の場合が特異点)。$x=±\sqrt{1-y^2}$も同様($x=0$の場合が特異点)。<br>[【Rで球面幾何学】「半円しか描けなかった」世界の思い出？](https://qiita.com/ochimusha01/items/cf74b1fff650edb0a496)<br>![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fqiita-image-store.s3.ap-northeast-1.amazonaws.com%2F0%2F597772%2F0ba98cbe-b449-9b1d-b0f0-4a6b60e55f35.gif?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&w=1400&fit=max&s=dacd205b94963cc117e9977dda71d05d)
<br>![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fqiita-image-store.s3.ap-northeast-1.amazonaws.com%2F0%2F597772%2Fdae51c2d-b2ce-515b-c55e-d8dec59a0760.gif?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&w=1400&fit=max&s=59f7dd6225b5486bdfc394f56e8c396b)

かかる特異点は「**正負の概念導入に伴う円弧概念の二分**」に由来し、**偶奇性**(Parity)や「**対蹠間を結ぶ円弧状に2つ、球表面上には無数に存在する経路**」**複素共益**(Complex Conjugate)の概念と同時発生する訳です。
[【Python演算処理】単位円筒を巡る数理](https://qiita.com/ochimusha01/items/fd8e51940b5e0e1a28c5)
![image.gif](https://cdn-ak.f.st-hatena.com/images/fotolife/o/ochimusha01/20210412/20210412045910.gif)
![image.gif](https://cdn-ak.f.st-hatena.com/images/fotolife/o/ochimusha01/20210412/20210412050306.gif)
![image.gif](https://camo.qiitausercontent.com/9e756e40dceca25c93c9517bc61e134658c8ae2e/68747470733a2f2f63646e2d616b2e662e73742d686174656e612e636f6d2f696d616765732f666f746f6c6966652f6f2f6f6368696d7573686130312f32303230313132362f32303230313132363032303830302e676966)
[【初心者向け】複素共役のアニメーション表示について。](https://qiita.com/ochimusha01/items/e0bb2789a420d3317cf1)
![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fqiita-image-store.s3.ap-northeast-1.amazonaws.com%2F0%2F597772%2F0aeed895-bafb-002c-82f5-2565c3145d80.gif?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&w=1400&fit=max&s=b3ec3e935c84f41e5178fc5290e94539)
![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fqiita-image-store.s3.ap-northeast-1.amazonaws.com%2F0%2F597772%2F98cb3b85-8775-d5ac-4f21-5608fbb7db13.gif?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&w=1400&fit=max&s=67dbcf7af7ae818104e752b345cc2e01)
![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fqiita-image-store.s3.ap-northeast-1.amazonaws.com%2F0%2F597772%2F161ed342-36aa-8db1-c463-f1ea040eb6d9.gif?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&w=1400&fit=max&s=4e34d38444b79004490be7610c262a10)
[【Python演算処理】単位円筒(Unit Cylinder)を巡る数理](https://qiita.com/ochimusha01/items/fd8e51940b5e0e1a28c5)
![image.gif](https://staging-qiita-user-contents.imgix.net/https%3A%2F%2Fqiita-user-contents.imgix.net%2Fhttps%253A%252F%252Fqiita-user-contents.imgix.net%252Fhttps%25253A%25252F%25252Fcdn-ak.f.st-hatena.com%25252Fimages%25252Ffotolife%25252Fo%25252Fochimusha01%25252F20210105%25252F20210105003122.gif%253Fixlib%253Drb-1.2.2%2526auto%253Dformat%2526gif-q%253D60%2526q%253D75%2526w%253D1400%2526fit%253Dmax%2526s%253D1c5e7641bb15c183f387fb9c6c92759a%3Fixlib%3Drb-1.2.2%26auto%3Dformat%26gif-q%3D60%26q%3D75%26s%3D8b3c1ac01fe5a4de036898e40dd0d907?ixlib=rb-4.0.0&auto=format&gif-q=60&q=75&w=1400&fit=max&s=d1618c45e01ce045850f462a826425e9)
![image.gif](https://staging-qiita-user-contents.imgix.net/https%3A%2F%2Fqiita-user-contents.imgix.net%2Fhttps%253A%252F%252Fqiita-user-contents.imgix.net%252Fhttps%25253A%25252F%25252Fcdn-ak.f.st-hatena.com%25252Fimages%25252Ffotolife%25252Fo%25252Fochimusha01%25252F20210105%25252F20210105071334.gif%253Fixlib%253Drb-1.2.2%2526auto%253Dformat%2526gif-q%253D60%2526q%253D75%2526w%253D1400%2526fit%253Dmax%2526s%253D6cc2d838c056f87959b9d35678b0d60a%3Fixlib%3Drb-1.2.2%26auto%3Dformat%26gif-q%3D60%26q%3D75%26w%3D1400%26fit%3Dmax%26s%3D302dce6d5da2947a562c26895875c993?ixlib=rb-4.0.0&auto=format&gif-q=60&q=75&w=1400&fit=max&s=d954c48b7a8ac00c138ead17c14d8273)

まさしく「騙し絵版画家」**エッシャー**（Maurits Cornelis Escher, 1898年～1972年）「**版画画廊**(Print Gallery,1956年)」における中央空白部のイメージ。
[エッシャーはどのようにして「奇想の版画家」になったのか？](https://bijutsutecho.com/magazine/insight/18015)
[『版画の画廊』における「ミーズ・アン・ナビーム」](http://tetsusala.seesaa.net/article/400611052.html)
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/04774cba-47df-a5e5-8bfe-03dc37ab8ebf.png)

エッシャーの騙し絵については以前「**正四面体の空間充填性**(単独では備えず、正八面体と組み合わせる必要がある)」でも触れました。人間の空間認識能力はかくも騙されやすいのです。
[【オイラーの多面体定理と正多面体】正四面体と正六面体と正八面体の連続性と「テトラパックの思わぬ正体」について。](https://qiita.com/ochimusha01/items/3356f7c805917660632c)
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/42b10708-f485-9c20-0046-a06715f8cb08.png)
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/22c9d8ff-ebff-a39b-181e-bd8b940b5ab4.png)
![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fqiita-image-store.s3.ap-northeast-1.amazonaws.com%2F0%2F597772%2Ffe7a0a07-fc3f-c8fc-418b-ef57c8d8438e.gif?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&w=1400&fit=max&s=cd927a335b4e862ee2ce92bc9d572714)
![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fqiita-image-store.s3.ap-northeast-1.amazonaws.com%2F0%2F597772%2Fa0a7c416-ad7b-6c2d-ad72-eb0170710366.gif?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&w=1400&fit=max&s=6d360ce6da74c0b6dbae8390b26590b9)
![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fcdn-ak.f.st-hatena.com%2Fimages%2Ffotolife%2Fo%2Fochimusha01%2F20210125%2F20210125150335.gif?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&w=1400&fit=max&s=1fb9c7cd4d03786cb5aae9078b8e131e)
![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%252F20201129%252F20201129055848.gif%3Fixlib%3Drb-1.2.2%26auto%3Dformat%26gif-q%3D60%26q%3D75%26s%3D8aaababc6b4f244b6086ae1fe5f1bdc1?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&w=1400&fit=max&s=de13e20e0ff195c72ddef7275db87d64)
上掲の複素共益アニメーションが観察者の認識によって進行方向を変える様に(<font color=blue>カンブリア爆発期以降一部の生物が備える様になった視覚とその情報を処理する脊髄に由来する</font>)人間の空間認識能力(<font color=blue>およびその発展形たる数理処理能力</font>)は欠陥だらけ。おそらく**負数**や**虚数**の概念はそれを補う為に発明されたと考えるべきなのです。

#物理学における等速円運動の解釈

[等速円運動についての物理学と数学の立場の違い?](https://qiita.com/ochimusha01/items/d5526764c7d1199671df)
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/5619f730-5987-5919-a26b-7fb42af30280.png)
物理学の世界は等速円運動を円軌道上の直線運動として捉え様としますが、この考え方だけでは「**X方向に1進もうとしてY方向にも1進んでしまう過程において移動距離が伸びる**」主観的体験を説明し尽くす事が出来ません。移動距離に注目するなら等速円運動は1動くのに$\frac{π}{2}$、2動くのにπ進まねばならない運動でもある訳で、それがどういう条件でどれだけズレるか数理的に明らかにするのが**距離関数**(Metric Function)の課題設定。この場合は曲率が0の場合、あるいは＋∞の場合が特異点として追加される事になります。
[【初心者向け】方形描画関数②距離空間との関係。](https://qiita.com/ochimusha01/items/06e49a35dffbd12b35aa)
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/c1c9fb6f-5290-020a-0dc8-2074961e24c2.png)
おや、特異点が増えてしまいました!!　かかる状況を説明する為、座標系における**視界**(Visibility)なる概念を導入したいと思います。

* 実軸が1本も設定されてない時、**視界**は**片側無限数列**(One-sided Infinite  Sequence)で表現可能な範囲で無制限と映る。**平面状極座標系**(Planar Polar Coordinates System=二次元)では…<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>同じく**片側無限数列**に立脚する**球状極座標系**(Spheroidal Polar Coordinates System=三次元)では…<br>[【初心者向け】「単位円筒」から「単位球面」へ](https://qiita.com/ochimusha01/items/08f5c5ad1aba4a59a27d)<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)
<br>ただし**両側無限数列**(Both-sided Infinite Sequence)に立脚する**円柱状極座標系**(Cylindrical Polar Coordinates System=三次元)は、この段階では以下の不完全な形でしか表現し得ない。<br>[【Python演算処理】単位円筒(Unit Cylinder)を巡る数理](https://qiita.com/ochimusha01/items/fd8e51940b5e0e1a28c5)<br>![image.gif](https://cdn-ak.f.st-hatena.com/images/fotolife/o/ochimusha01/20210411/20210411223344.gif)
* 正負の概念を導入すると(実軸を1本だけ導入すると)、今度は**平面状極座標系**上において半円しか描けなくなる。視界$\frac{1}{2}$(x=1→-1,y=0→±1)のy=±$\sqrt{1-x^2}$と同じく視界$\frac{1}{2}$(x=0→±1,y=-1→1)のx=±$\sqrt{1-y^2}$の対称式一対…<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%2F20201111043648.gif?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&w=1400&fit=max&s=30c4597a9f5a932f98335f99b2d0895d)<br>![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fcdn-ak.f.st-hatena.com%2Fimages%2Ffotolife%2Fo%2Fochimusha01%2F20201110%2F20201110234311.gif?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&w=1400&fit=max&s=1e84bc84fde30af53a4af9e9e295350b)<br>かかる制限を撤廃すべく、**複素平面**(Complex Plane)はこれらの対称式をそれぞれ$\cos(θ)+\sin(θ)i$(初期位相0ラジアン,反時計回り)と$\sin(θ)+\cos(θ)i$(初期位相π/2ラジアン,時計回り)と表す様になった訳である。<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%2F20201108%2F20201108040152.gif?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&w=1400&fit=max&s=393ceabfd8aac2d18252a315f9915419)<br>![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fcdn-ak.f.st-hatena.com%2Fimages%2Ffotolife%2Fo%2Fochimusha01%2F20201110%2F20201110223555.gif?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&w=1400&fit=max&s=1ab6aaf4879c0f9ea59a871231550b19)<br>一方、円柱状極座標系はやっと完全な形で表現可能となるが…<br>[【Python演算処理】単位円筒(Unit Cylinder)を巡る数理](https://qiita.com/ochimusha01/items/fd8e51940b5e0e1a28c5)<br>![image.gif](https://cdn-ak.f.st-hatena.com/images/fotolife/o/ochimusha01/20210411/20210411224605.gif)
<br>複素平面上での判断の反映させ方に(分枝切断点を1つだけ置いて)周期を1に取る方法と(分枝切断点を2つ置いて)2に取る方法の2通りがあり、後者を選ぶと偶奇の概念が導入される。<br>[【Python演算処理】単位円筒(Unit Cylinder)を巡る数理](https://qiita.com/ochimusha01/items/fd8e51940b5e0e1a28c5)<br>![image.gif](https://cdn-ak.f.st-hatena.com/images/fotolife/o/ochimusha01/20210412/20210412045131.gif)<br>![image.gif](https://staging-qiita-user-contents.imgix.net/https%3A%2F%2Fcdn-ak.f.st-hatena.com%2Fimages%2Ffotolife%2Fo%2Fochimusha01%2F20210329%2F20210329071623.gif?ixlib=rb-4.0.0&auto=format&gif-q=60&q=75&w=1400&fit=max&s=99b64d9425d977658b2b572bb97b8d65)<br>![image.gif](https://cdn-ak.f.st-hatena.com/images/fotolife/o/ochimusha01/20210412/20210412045910.gif)<br>![image.gif](https://cdn-ak.f.st-hatena.com/images/fotolife/o/ochimusha01/20210412/20210412050306.gif)<br>![image.gif](https://camo.qiitausercontent.com/9e756e40dceca25c93c9517bc61e134658c8ae2e/68747470733a2f2f63646e2d616b2e662e73742d686174656e612e636f6d2f696d616765732f666f746f6c6966652f6f2f6f6368696d7573686130312f32303230313132362f32303230313132363032303830302e676966)

こうした考え方の延長線上において実数軸を2本導入する方法は2通りあります。

* いわゆる**デカルト座標系**(Cartesian Coordinate System)をそのまま単純に導入し、横軸(X軸)の視界$\frac{1}{2}$と縦軸(Y軸)の視界$\frac{1}{2}$を合算すると$\frac{3}{4}$となってしまうので、全体像を視界$\frac{3}{4}(-i→-(i^2)=1→-(-i)=i)$の**実軸**(Real Axis)と**虚軸**(Imaginal Axis)の交代と捉えて「**死角**」を封じたのが**複素数表現**となる。それ自体は実軸を1本導入しただけの一次元表現である点に注意。<br>[【Rで球面幾何学】「半円しか描けなかった」世界の思い出？](https://qiita.com/ochimusha01/items/cf74b1fff650edb0a496)<br>![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fqiita-image-store.s3.ap-northeast-1.amazonaws.com%2F0%2F597772%2F9c314439-c60a-7804-c2ae-26d0f6382a23.gif?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&w=1400&fit=max&s=65069edfaaf08c32d62cbf8e9a02cdea)<br>ところで、実軸を2本導入した時点で$\frac{1}{4}$となる「**死角**」は理論上、導入する実軸を無限に増やし続ける事で限りなく0ラジアンに近づけていける(それに応じて視界は全周=2πラジアンに接近する)が、これを求めるベルヌーイ試行「**N回の試行で確率p/Nの事象が1度も起らない確率**」の計算式$(1-\frac{1}{n})^n$は、実は$e^{-1}$(ネイピア数の逆数)を求める計算式そのものなのである。その事を式で表したのが**オイラーの等式**$e^{πi}=(1±\frac{πi}{N})^N=(0±1i)^{2x}=-1$となり、**共益複素数**(Conjugate Complex Number)が「**1から-1に向かう**(平面上は2つ、球表面上には無数にある)**経路**」として把握される展開を迎える。<br>[ベルヌーイ試行 - Wikipedia](https://ja.wikipedia.org/wiki/%E3%83%99%E3%83%AB%E3%83%8C%E3%83%BC%E3%82%A4%E8%A9%A6%E8%A1%8C)<br>[【初心者向け】指数・対数関数の発見とそれ以降の発展について。](https://qiita.com/ochimusha01/items/3e28f4f06abd8cef0af1)<br>![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fqiita-image-store.s3.ap-northeast-1.amazonaws.com%2F0%2F597772%2F241a9d60-380c-d71e-08e3-085d596ec8fe.gif?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&w=1400&fit=max&s=08fb0b1e6c6390220480c2b75dea42e9)

そして特異点に到達不可能だという事は、**視界**が常に開空間としてのみ存在する(すなわち観測結果集合も開集合としてのみ存在する)事を意味します。
[近傍 (位相空間論) - Wikipedia](https://ja.wikipedia.org/wiki/%E8%BF%91%E5%82%8D_(%E4%BD%8D%E7%9B%B8%E7%A9%BA%E9%96%93%E8%AB%96))

#Γ関数の登場
ところで円の面積を求める公式$πr^2$を半径rで微分すると円周を求める公式$2πr$になりますが、πを変数に見立て、これで微分すると$\frac{\sqrt{π}}{2}r^2$となります。そしてこの見慣れない$\frac{\sqrt{π}}{2}$なる数字、Γ関数に1/2(すなわち1と2の中間値)を与えた場合の返り値としても登場してくるのです。
[【Python演算処理】階乗と順列と組み合わせ](https://qiita.com/ochimusha01/items/90d9ef0a224bdd604721)
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/5e20d3eb-dbfd-bf0b-8455-3a307685f345.png)

```python
import numpy as num
from scipy.special import gamma

# [Γ(1) Γ(1.5) Γ(2)]
a=gamma([1, 1.5, 2])
a=num.round(a, 8)
print("Γ(1,1.5,2)=",a)

#sqrt(π)/2
b=num.sqrt(num.pi)/2
b=num.round(b, 8)
print("sqrt(π)/2=",b)

#出力結果
Γ(1,1.5,2)= [1.         0.88622693 1.        ]
sqrt(π)/2= 0.88622693
```

* 階乗計算を小数点下の添字にも対応させたΓ関数はガウス関数$e^{-x^2}$を実数域全体で広義積分した**ガウス積分**（Gaussian integral）$\int_{-\infty}^{\infty}e^{-x^2}dx=\sqrt{π}$と密接な関係がある。

[ガウス積分 - Wikipedia](https://ja.wikipedia.org/wiki/%E3%82%AC%E3%82%A6%E3%82%B9%E7%A9%8D%E5%88%86)

* 被積分関数が偶関数なので$\int_{-\infty}^{\infty}e^{-x^2}dx=2\int_{0}^{\infty }e^{-x^2}dx$が成立。計算が単純化される。
* ここで変数変換$x=t^{\frac{1}{2}}$を行えばオイラー積分が可能となり、半整数値の階乗が$\sqrt{π}$の有理数倍となる事が示される。 

```math
2\int_{0}^{\infty }e^{-x^2}dx=2\int_{0}^{\infty }\frac{1}{2}e^{-t}t^{\frac{1}{2}}dt=2Γ(\frac{1}{2})=\sqrt{π}
```

* これを一般すれば$\int_{0}^{\infty }e^{-ax^b}dx=\frac{1}{b}a^{\frac{-1}{b}}Γ(\frac{1}{b})$となる。 

このΓ関数の概念を用いると円概念、すなわち半径/直径/円弧/球面/体積それぞれの計算式の往復をスッキリ一つの体系に統合する事が出来ます。
[【Python演算処理】「N次球概念導入による円/球関係の数理の統合」？](https://qiita.com/ochimusha01/items/f7c9077eccfdbd0dd47a)

かかる円概念統合において「**円弧と面積の計算式の往復**」や「**球面と体積の計算式の往復**」と関連してくるので、重積分を用いた証明も示しておきましょう。
[Pythonで数学の勉強：微分積分学](https://qiita.com/tibigame/items/aebbac176d9bbdaf3d15)

* 半球D、すなわち$z=\sqrt{1-x^2-y^2}$と$z=0$に囲まれた部分の体積Vを求める。

```math
v=\iint_{D}\sqrt{1-x^2-y^2}dxdy
```
* Dの第一象限は$y=\sqrt{1-x^2}$なので以下となる。

```math
v=4\int_{0}^{1}\int_{0}^{\sqrt{1-x^2}}\sqrt{1-x^2-y^2}dydx
```
この積分は普通に求まらないから数値積分してみる。

```python
import sympy as sym
sym.init_printing()
Pi = sym.S.Pi # 円周率
E = sym.S.Exp1 # 自然対数の底
I = sym.S.ImaginaryUnit # 虚数単位
oo = sym.oo # 無限大

# 使用する変数の定義(小文字1文字は全てシンボルとする)
(a,b,c,d,e,f,g,h,i,j,k,l,m,n,o,p,q,r,s,t,u,v,w,x,y,z) = sym.symbols('a b c d e f g h i j k l m n o p q r s t u v w x y z')

import numpy as np
import matplotlib.pyplot as plt
%matplotlib inline

from mpl_toolkits.mplot3d import Axes3D
from scipy import integrate

print(2*np.pi/3) # 理論値
integrate.dblquad(lambda x, y: 4*np.sqrt(1-x**2-y**2), 0, 1, lambda x: 0, lambda x: np.sqrt(1-x**2))
```
2.0943951023931953
(2.0943951023931957, 3.5335157022586827𝑒−10)

* ここで$x=r\cosθ,y=r\sinθ$とおく。

```math
V=\int_{0}^{2π}\int_{0}^{1}\sqrt{1−r^2}\begin{vmatrix}
\frac{∂x}{∂r} & \frac{∂x}{∂θ} \\
\frac{∂y}{∂r} & \frac{∂y}{∂θ}
\end{vmatrix} drdθ
```

```python
import sympy as sym
sym.init_printing()
Pi = sym.S.Pi # 円周率
E = sym.S.Exp1 # 自然対数の底
I = sym.S.ImaginaryUnit # 虚数単位
oo = sym.oo # 無限大

# 使用する変数の定義(小文字1文字は全てシンボルとする)
(a,b,c,d,e,f,g,h,i,j,k,l,m,n,o,p,q,r,s,t,u,v,w,x,y,z) = sym.symbols('a b c d e f g h i j k l m n o p q r s t u v w x y z')

import numpy as np
import matplotlib.pyplot as plt
%matplotlib inline

from mpl_toolkits.mplot3d import Axes3D
from scipy import integrate

X = r * sym.cos(t)
Y = r * sym.sin(t)
Z = sym.sqrt(1-x**2-y**2).subs([(x, X), (y, Y)]).trigsimp() # 簡略化は必要(1行ずつやった方がいい)
Jacobian = sym.Matrix([[sym.diff(X, r), sym.diff(X, t)], [sym.diff(Y, r), sym.diff(Y, t)]])
sym.integrate(sym.integrate(Z * Jacobian.det().trigsimp(), (r, 0, 1)), (t, 0, 2*Pi))
```

```math
\frac{2π}{3}
```
見方を変えればVの全領域を1で重積分しても同じ結果となる。

```math
v=\iint_{V}dxdydz
```
これも$x=r\sinθ\cosφ,y=r\sinθ\sinφ,z=r\cosθ$とおく。

```math
V=2\int_{0}^{\frac{π}{2}}\int_{0}^{π}\int_{0}^{1}\begin{vmatrix}
\frac{∂x}{∂r} & \frac{∂x}{∂θ}  & \frac{∂x}{∂φ}\\
\frac{∂y}{∂r} & \frac{∂y}{∂θ}  & \frac{∂y}{∂φ}\\
\frac{∂z}{∂r} & \frac{∂z}{∂θ} & \frac{∂z}{∂φ}
\end{vmatrix} drdθdφ
```

```python
import sympy as sym
sym.init_printing()
Pi = sym.S.Pi # 円周率
E = sym.S.Exp1 # 自然対数の底
I = sym.S.ImaginaryUnit # 虚数単位
oo = sym.oo # 無限大

# 使用する変数の定義(小文字1文字は全てシンボルとする)
(a,b,c,d,e,f,g,h,i,j,k,l,m,n,o,p,q,r,s,t,u,v,w,x,y,z) = sym.symbols('a b c d e f g h i j k l m n o p q r s t u v w x y z')

import numpy as np
import matplotlib.pyplot as plt
%matplotlib inline

from mpl_toolkits.mplot3d import Axes3D
from scipy import integrate

X = r * sym.sin(t) * sym.cos(u)
Y = r * sym.sin(t) * sym.sin(u)
Z = r * sym.cos(t)
Jacobian = sym.Matrix([
    [sym.diff(X, r), sym.diff(X, t), sym.diff(X, u)],
    [sym.diff(Y, r), sym.diff(Y, t), sym.diff(Y, u)],
    [sym.diff(Z, r), sym.diff(Z, t), sym.diff(Z, u)]
])
2*sym.integrate(sym.integrate(sym.integrate(Jacobian.det().trigsimp(), (r, 0, 1)), (t, 0, Pi)), (u, 0, Pi/2))
```

```math
\frac{2π}{3}
```
とりあえず現状では転記が精一杯。ちなみに**線績分**(Line Integral)・**面積分**(Surface Integral)・**体積分**(Volume Integral)を追求する過程で**ヤコビアン**(Jacobian=ヤコビ数列)が登場する過程については以下が詳しい。
[EMANの物理学>物理数学>線積分](https://eman-physics.net/math/calculus02.html)
[EMANの物理学>物理数学>面積分と体積分](https://eman-physics.net/math/calculus03.html)
[EMANの物理学>物理数学>ヤコビアン](https://eman-physics.net/math/calculus04.html)

```math
\int U(x,y,z)dV=\int_{0}^{2π}\int_{0}^{π}\int_{0}^{1}U(r,θ,φ)r^2\sinθdrdθdφ\\
=\int_{0}^{2π}\int_{0}^{π}\int_{0}^{1}\begin{vmatrix}
\frac{∂x}{∂r} & \frac{∂x}{∂θ}  & \frac{∂x}{∂φ}\\
\frac{∂y}{∂r} & \frac{∂y}{∂θ}  & \frac{∂y}{∂φ}\\
\frac{∂z}{∂r} & \frac{∂z}{∂θ} & \frac{∂z}{∂φ}
\end{vmatrix} drdθdφ
```

```python
import sympy as sym
sym.init_printing()
Pi = sym.S.Pi # 円周率
E = sym.S.Exp1 # 自然対数の底
I = sym.S.ImaginaryUnit # 虚数単位
oo = sym.oo # 無限大

# 使用する変数の定義(小文字1文字は全てシンボルとする)
(a,b,c,d,e,f,g,h,i,j,k,l,m,n,o,p,q,r,s,t,u,v,w,x,y,z) = sym.symbols('a b c d e f g h i j k l m n o p q r s t u v w x y z')

import numpy as np
import matplotlib.pyplot as plt
%matplotlib inline

from mpl_toolkits.mplot3d import Axes3D
from scipy import integrate

X = r * sym.sin(t) * sym.cos(u)
Y = r * sym.sin(t) * sym.sin(u)
Z = r * sym.cos(t)
Jacobian = sym.Matrix([
    [sym.diff(X, r), sym.diff(X, t), sym.diff(X, u)],
    [sym.diff(Y, r), sym.diff(Y, t), sym.diff(Y, u)],
    [sym.diff(Z, r), sym.diff(Z, t), sym.diff(Z, u)]
])
sym.integrate(sym.integrate(sym.integrate(Jacobian.det().trigsimp(), (r, 0, 1)), (t, 0, Pi)), (u, 0, 2*Pi))
```

```math
\frac{4π}{3}
```
[ガウス積分 - Wikipedia](https://ja.wikipedia.org/wiki/%E3%82%AC%E3%82%A6%E3%82%B9%E7%A9%8D%E5%88%86)

* 初等関数としての不定積分$\int e^{-x^2}dx$は存在しないが、定積分$\int_{-\infty}^{\infty} e^{-x^2}dx$は評価可能である。
* 極座標を用いてガウス積分を求めるアイデアはポアソン(1781年~1840年)まで遡る。

①平面$\mathbb {R}^2$上の函数$e^{−(x2+y2)} = e^{−r^2}$を考え、これを2通りの方法で計算する。

* 直交座標系に関する二重積分として計算し、その値は求める値の平方になることを確かめる。
* 極座標系に関する二重積分（いわゆるバウムクーヘン積分）として計算し、その値が π となることを確かめる。

②広義積分が現れることに注意して、これら2つの計算を比較して積分の値が求まる。即ち、面積要素$dA$が xy-直交座標系では$dA=dxdy$, rθ-極座標系では$dA=rdrdθ$で与えられることに注意すれば、前者と後者はそれぞれ以下の様に計算される。

```math
\int _{\mathbb {R}^2}e^{-(x^2+y^2)}dA=\int_{-\infty}^{\infty}\int_{-\infty}^{\infty}e^{-(x^2+y^2)}dxdy=(\int_{-\infty }^{\infty }e^{-t^2}dt)^2
```
および

```math
\int _{\mathbb {R}^2}e^{-(x^2+y^2)}dA=\int_{0}^{2π}\int_{0}^{π}e^{-r^2}rdrdθ\\
=2π\int_{0}^{\infty }r e^{-r^2}dr=π\int_{-\infty }^{0}e^sds=π
```
③後者では$s=−r^2$なる置換を行って、$ds=−2rdr$となることを用いた。これらの結果から以下となる。

```math
(\int_{-\infty }^{\infty }e^{-x^2}dx)^2=π
```
④符号を考慮して以下となる。

```math
\int_{-\infty }^{\infty }e^{-x^2}dx=\sqrt{π}
```

* 直交座標を用いてガウス積分を求めるアイディアはラプラス(1812年)にまで遡れる。

①$y=xs,dy=xds$と置くと、yを±∞へ近づけるときsの極限はxの符号で決まる。そして$e^{-x^2}$は偶函数なので実数全体にわたる積分が正の実数全体にわたる積分の2倍になる事から以下の様な形で計算を簡単にする事が出来る。

```math
\int_{-\infty}^{\infty}e^{-x^2}dx=2\int_{0}^{\infty }e^{-x^2}dx
```
②すなわち積分範囲を x ≥ 0 に限れば、変数 y と s とは同じ極限を持ち以下が成立する。

```math
I^2=4\int_{0}^{\infty}\int_{0}^{\infty }e^{-(x^{2}+y^{2})}dydx
```
③そこから以下の展開で所期の$I=\sqrt{π}$を得る。

```math
\frac{I^2}{4}=\int _{0}^{\infty}(\int_{0}^{\infty}e^{-(x^2+y^2)}dy)dx=\int_{0}^{\infty}(\int_{0}^{\infty}e^{-x^2}(1+s^{2})xds)dx\\
=\int_{0}^{\infty}(\int_{0}^{\infty}e^{-x^2}(1+s^{2})xdx)ds=\frac{1}{2}\int_{0}^{\infty}{\frac{ds}{1+s^2}}\\
={\frac{1}{2}}\arctan s|_{0}^{\infty }={\frac  {\pi }{4}}
```
ところで、**ガウス関数**の正体は「**対数尺に放り込む事によって発散の方向をX軸に歪めた放物線**」に他なりません。
[【初心者向け】指数・対数関数の発見とそれ以降の発展について。](https://qiita.com/ochimusha01/items/3e28f4f06abd8cef0af1)
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/92ff2f7d-07e4-3fb6-992d-9b9a11ea498f.png)

この何処に円要素が？　実は**正規分布**(Normal Distribution)として知られる**ベルカーブ**$y=e^{-(x^2+y^2)}$は、3次元的には以下の様に同心円的堆積を構成するのです。
[【初心者向け】正規分布(Normal Distribution)とは何か？](https://qiita.com/ochimusha01/items/fe73ddf6e3e8cb08fbee)

![image.gif](https://cdn-ak.f.st-hatena.com/images/fotolife/o/ochimusha01/20210402/20210402034523.gif)

```python
import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import axes3d
from matplotlib import cm
import matplotlib.animation as animation

fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')

x = y = np.arange(-15, 15, 0.5)
X, Y = np.meshgrid(x, y)

def dg(n):
    plt.cla()
    sigma = 4
    Z = np.exp(-(X**2 + Y**2)/(2*sigma**2)) / (2*np.pi*sigma**2)
    ax.plot_surface(X, Y, Z, rstride=1, cstride=1,cmap="summer")
    # グラフを回転
    ax.view_init(elev=45, azim=Time_code[n])
    
Time_code=num.arange(0,360,6)    

ani = animation.FuncAnimation(fig, dg, interval=50,frames=len(Time_code))
ani.save("output400.gif", writer="pillow")
```
という事はつまり…そんな感じで以下続報。