【Python演算処理】微分の概念を高校数学段階からグレードアップ？

以下の関数の定義から出発します。
[【Python演算処理】冪算と乗除算の関係について。](https://qiita.com/ochimusha01/items/8a66fcab8285a146265f)

>二つの数の集合X,Yがあって、Xの要素xに対してYの要素yがただ一つ定まる時、この対応を**関数**(Function)と呼びy=f(x)などと書き表す(fはfunctionの略だが、複数の関数への言及が不可避となる状況ではそれをy=g(x),y=h(x)などと呼び分ける)。

* この時xを**独立変数**(Independent Variable)、yを**従属変数**(Dependent Variable)と呼ぶ。
* また集合Xを関数fの**定義域**(Domain)、集合Yの部分集合$(y|y=f(x),x \in X)$をこの関数の**値域**(Range)という。

>関数y=fxの変化の様子はxの値を横座標、yの値を縦座標としxを**定義域**内で変化させた時の点P(x,y)の軌跡を描くと一目瞭然となる。これが関数y=f(x)のグラフである。

正直言って高校までの数学においては**微分**(Differential)は事実上以下の公式の暗記と再生に過ぎませんでした。

[微分公式一覧（基礎から発展まで） | 高校数学の美しい物語](https://manabitimes.jp/math/1109)

>微分計算は素直。積分計算は超エキサイティング。**初等関数**（三角関数や指数関数など）の四則演算や合成で表現できる関数は，基本的な公式を組み合わせるだけで必ず微分できます（一方，不定積分は必ずしも初等関数で表せるとは限らない）。

##基本演算

###①線形性
$(kf(x))'=kf'(x)$(ただしkは定数=スカラー倍の概念に対応)
$(af(x)±bg(x))'=af'(x)±bg'(x)$
[高校数学における線形性の8つの例 | 高校数学の美しい物語](https://manabitimes.jp/math/684)
別の記法ではこうなる。

```math
\frac{d}{dt}(kX(t))=k\frac{d}{dt}X(t)\\
\frac{d}{dt}(aX(t)±bY(t))=a\frac{d}{dt}X(t)±b\frac{d}{dt}Y(t) 
```

###②積の微分：
$(fg)'=f'g+fg'$
[積の微分公式とその証明の味わい](https://manabitimes.jp/math/1079)

**【例題】**関数$y=f(x)=x^2$と関数$y=g(x)=\sin(x)$の積の微分。

```math
\frac{dy}{dx}x^2=2x\\
\frac{dy}{dx}\sin(x)=\cos(x)
```
よって**積の微分公式**f'g+fg'に従って

```math
f'g+fg'=2x \sin(x)+x^2 \cos(x)
```
pythonによる実装結果

```python
import sympy as sp
x,y = sp.symbols('x,y')

#原始関数
f=x**2*sp.sin(x)

#微分
g=sp.diff(f,x,1)

#積分
h=sp.integrate(g,x)

#画面表示
sp.init_printing()
display(f) 
print(sp.latex(f))
display(g) 
print(sp.latex(g))
display(h) 
print(sp.latex(h))
```
原始関数

```math
x^{2} \sin{\left(x \right)}
```
微分結果

```math
x^{2} \cos{\left(x \right)} + 2 x \sin{\left(x \right)}
```
積分による復元

```math
x^{2} \sin{\left(x \right)}
```
###③商の微分（<font color=blue>分数の微分</font>）：
$(\frac{g}{f})'=\frac{g'f-f'g}{f^2}$
[商の微分公式の証明と例題](https://manabitimes.jp/math/1108)​

**【例題】**$y=f(x)=\frac{\cos (x)}{\sin (x)}$の微分

```math
\frac{dy}{dx}\cos(x)=-\sin(x)\\
\frac{dy}{dx}\sin(x)=\cos(x)\\
\cos^2(x)+\sin^2(x)=1
```

よって商の微分公式$\frac{g'f-f'g}{f^2}$に従って

```math
y'=\frac{\cos(x)'sin(x)-\sin(x)'\cos(x)}{\sin^2(x)}\\  
= \frac{-\sin^2(x)-\cos^2(x)}{\sin^2(x)}\\
= -\frac{1}{\sin^2(x)}
```
pythonによる実装結果

```python
import sympy as sp
x,y = sp.symbols('x,y')

#原始関数
f=sp.cos(x)/sp.sin(x)

#微分
g=sp.diff(f,x,1)

#積分
h=sp.integrate(g,x)

#画面表示
sp.init_printing()
display(f) 
print(sp.latex(f))
display(g) 
print(sp.latex(g))
display(h) 
print(sp.latex(h))
```
原始関数

```math
\frac{\cos{\left(x \right)}}{\sin{\left(x \right)}}
```
微分結果

```math
-1 - \frac{\cos^{2}{\left(x \right)}}{\sin^{2}{\left(x \right)}}
```
積分による復元

```math
\frac{\cos{\left(x \right)}}{\sin{\left(x \right)}}
```
<font color=green>まぁSympyの癖も分かってきたので、別にこの程度では驚きません。逆をいえば「**ここでピタゴラスの定理の導入を確実に思いつく**」のが人工知能とも?</font>

* 特にg=1の時，$(\frac{1}{f})'=-\frac{f'}{f^2}$(<font color=blue>逆数の微分</font>)

**【例題】**$\frac{1}{log(x)}$の微分

```math
\frac{dy}{dx}\log(x)=\frac{1}{x}
```
よって逆数の微分公式$-\frac{f'}{f^2}$に従って

```math
y'=-\frac{\log(x)'}{\log(x)^2}=-\frac{1}{x(\log(x)^2}
```
pythonによる実装結果

```python
import sympy as sp
x,y = sp.symbols('x,y')

#原始関数
f=1/sp.log(x)

#微分
g=sp.diff(f,x,1)

#積分
h=sp.integrate(g,x)

#画面表示
sp.init_printing()
display(f) 
print(sp.latex(f))
display(g) 
print(sp.latex(g))
display(h) 
print(sp.latex(h))
```
原始関数

```math
\frac{1}{\log{\left(x \right)}}
```
微分結果

```math
- \frac{1}{x \log{\left(x \right)}^{2}}
```
積分による復元

```math
\frac{1}{\log{\left(x \right)}}
```
###④合成関数の微分：
$(f(g(x))'=f'(g(x))g'(x)$
[合成関数の微分公式と例題7問](https://manabitimes.jp/math/936)
別の記法ではこうなる。

>y=f(u)がuについて微分可能であり、かつu=g(x)がxについて微分可能であれば、**合成関数**y=f(g(x))はxについて微分可能であり$\frac{dy}{dx}=\frac{dy}{du}\frac{du}{dx}$と表される。

**【例題】**関数$y=(5x+3)^2$の微分

```math
y=f(u)=u^2…①\\
u=g(x)=5x+3…②
```
この時、yは以下の形でxの関数となる。

```math
y=f(u)=u^2=(5x+3)^2…③
```
逆をいえば$y=(5x+3)^2$の式について$u=5x+3$と置いたのが①②ともいえる。③を展開して項別にxで微分すれば導関数が得られるが、③は①②の合成関数なので以下となる。

```math
\frac{dy}{dx}=\frac{dy}{du}\frac{du}{dx}\\
\frac{dy}{du}u^2=2u\\
\frac{du}{dx}5x+3=5\\
\frac{dy}{du}\frac{du}{dx}=2u*5=10u=10(5x+3)=50x+30
```
pythonによる実装結果

```python
import sympy as sp
x,y = sp.symbols('x,y')

#原始関数
f=(5*x+3)**2

#微分
g=sp.diff(f,x,1)

#積分
h=sp.integrate(g,x)

#画面表示
sp.init_printing()
display(f) 
print(sp.latex(f))
display(g) 
print(sp.latex(g))
display(h) 
print(sp.latex(h))
```
原始関数

```math
\left(5 x + 3\right)^{2}
```
微分結果

```math
50 x + 30
```
積分による復元

```math
25 x^{2} + 30 x
```
<font color=green>合成微分には積分による復元性に難がある?</font>

**【例題】**$\cos(ax+b)$の微分
$y=\cos(u),u=ax+b$と置く。

```math
\frac{dy}{dx}\cos(x)=-\sin(x)\\
```
よって合成関数の微分公式に従い

```math
\frac{dy}{dx}=\frac{dy}{du}\frac{du}{dx}\\
\frac{dy}{du}cos(x)=-sin(u)\\
\frac{du}{dx}ax+b=a\\
\frac{dy}{du}\frac{du}{dx}=a(-\sin(u))=-a(\sin(ax+b))
```
pythonによる実装結果

```python
import sympy as sp
x,a,b = sp.symbols('x,a,b')

#原始関数
f=sp.cos(a*x+b)

#微分
g=sp.diff(f,x,1)

#積分
h=sp.integrate(g,x)

#画面表示
sp.init_printing()
display(f) 
print(sp.latex(f))
display(g) 
print(sp.latex(g))
display(h) 
print(sp.latex(h))
```
原始関数

```math
\cos{\left(a x + b \right)}
```
微分結果

```math
- a \sin{\left(a x + b \right)}
```
積分による復元

```math
- a \left(\begin{cases} - \frac{\cos{\left(a x + b \right)}}{a} & \text{for}\: a \neq 0 \\x \sin{\left(b \right)} & \text{otherwise} \end{cases}\right)
```
<font color=green>やはり合成微分には積分による復元性に難がある?</font>

**【例題】**$e^{ax}$の微分
$y=e^u,u=ax$と置く。

```math
\frac{dy}{dx}e^x=e^x\\
```
よって合成関数の微分公式に従い

```math
\frac{dy}{dx}=\frac{dy}{du}\frac{du}{dx}\\
\frac{dy}{du}e^u=-e^u\\
\frac{du}{dx}ax=a\\
\frac{dy}{du}\frac{du}{dx}=a(e^u)=-ae^{ax}
```
pythonによる実装結果

```python
import sympy as sp
x,a= sp.symbols('x,a')

#原始関数
f=sp.exp(a*x)

#微分
g=sp.diff(f,x,1)

#積分
h=sp.integrate(g,x)

#画面表示
sp.init_printing()
display(f) 
print(sp.latex(f))
display(g) 
print(sp.latex(g))
display(h) 
print(sp.latex(h))
```
原始関数

```math
e^{a x}
```
微分結果

```math
a e^{a x}
```
積分による復元

```math
e^{a x}
```

###⑤（公式ではないが）対数微分法

* 「**両辺>0**」を確認して対数を取る（正でないときは絶対値を取る）
* 両辺を x で微分する
* y'について解く

[対数微分法のやり方と例題](https://manabitimes.jp/math/923)
**【例題】**$y=x^x$の微分。

>$y=x^x$の指数の底は正なので，そもそも関数の定義域がx>0xであることに注意する。よって，両辺は正なのでそのまま対数を取れる：

```math
\log(y)=xlog(x)
```
>両辺をxで微分する。左辺は前提知識で述べた公式，右辺は積の微分の公式を使う： 

```math
\frac{y'}{y}=log(x)+1
```

>これを y'について解く：

```math
y'=y(\log(x)+1)
```

>yをもとに戻す： 

```math
y'=x^x(\log(x)+1)
```
pythonによる実装結果

```python
import sympy as sp
x= sp.symbols('x')

#原始関数
f=x**x

#微分
g=sp.diff(f,x,1)

#積分
h=sp.integrate(g,x)

#画面表示
sp.init_printing()
display(f) 
print(sp.latex(f))
display(g) 
print(sp.latex(g))
display(h) 
print(sp.latex(h))
```
原始関数

```math
x^{x}
```
微分結果

```math
x^{x} \left(\log{\left(x \right)} + 1\right)
```
積分による復元

```math
x^{x}
```
<font color=green>Sympy、こういうのはケロっとこなすのでかえってその考え方の重要性を見逃す契機になってしまう?</font>

##冪関数の微分

```math
(x^{α})'=αx^{α-1}（αは任意の実数）
```
[べき関数（y=x^n）の微分公式の3通りの証明](https://manabitimes.jp/math/1132)

* 例えば$(x^2)'=2x,(x^{10})'=9x^9$

[【数理考古学】冪乗算の微積分](https://qiita.com/ochimusha01/items/0bd18370cfa96e83347a)

>対応する以下の積分公式がα-1の時に使えない。

```math
\int x^α dx=\frac{1}{α+1}x^{α+1}+C(α \in \mathbb{R} \land α \neq -1)
```

>それで以下の式で代用する事が多い。

```math
(α^x)'=α^x \log(α)\\
\int α^x dx = \frac{α^x}{\log(α)}
```

* α=−1とすると$(\frac{1}{x})'=-\frac{1}{x^2}$

pythonによる実装結果

```python
import sympy as sp
x= sp.symbols('x')

#原始関数
f=1/x

#微分
g=sp.diff(f,x,1)
g2=sp.diff(f,x,2)

#積分
h=sp.integrate(g,x)

#画面表示
sp.init_printing()
display(f) 
print(sp.latex(f))
display(g) 
print(sp.latex(g))
display(h) 
print(sp.latex(h))
display(g2) 
print(sp.latex(g2))
```
原始関数

```math
\frac{1}{x}
```
微分結果

```math
- \frac{1}{x^{2}}
```
積分による復元

```math
\frac{1}{x}
```
二階微分だと負号が消える

```math
\frac{2}{x^{3}}
```

* α=$\frac{1}{2}$とすると$(\sqrt{x})'=\frac{1}{2\sqrt{x}}$

pythonによる実装結果

```python
import sympy as sp
x= sp.symbols('x')

#原始関数
f=sp.sqrt(x)

#微分
g=sp.diff(f,x,1)

#積分
h=sp.integrate(g,x)

#画面表示
sp.init_printing()
display(f) 
print(sp.latex(f))
display(g) 
print(sp.latex(g))
display(h) 
print(sp.latex(h))
```
原始関数

```math
\sqrt{x}
```
微分結果

```math
\frac{1}{2 \sqrt{x}}
```
積分による復元

```math
\sqrt{x}
```

[【Pyrhon演算処理】同心集合①乗法的同心集合とは？](https://qiita.com/ochimusha01/items/9ed0fa5871b3f913306a)

>この演算は「**同心円集合**$r^{\frac{1}{2}}$**に内接/外接する正方形集合**」と縁深い事で知られている。

![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fcdn-ak.f.st-hatena.com%2Fimages%2Ffotolife%2Fo%2Fochimusha01%2F20210619%2F20210619222241.gif?ixlib=rb-4.0.0&auto=format&gif-q=60&q=75&w=1400&fit=max&s=7b472277a78ea87609240b8dabc8dbab)

##​指数対数関数の微分

微積分演算に対応しない形で以下の観測結果集合を構成する。
[【初心者向け】指数・対数関数の発見とそれ以降の発展について。](https://qiita.com/ochimusha01/items/3e28f4f06abd8cef0af1)
![imge.png](https://qiita-user-contents.imgix.net/https%3A%2F%2Fqiita-image-store.s3.ap-northeast-1.amazonaws.com%2F0%2F597772%2F216d94f6-b445-ae38-64d8-bd6f8db5a2c9.png?ixlib=rb-4.0.0&auto=format&gif-q=60&q=75&w=1400&fit=max&s=ca3acfaa779fa75fe0550daf03cdcf57)

<table border="1" class="dataframe">  <thead>    <tr style="text-align: right;">      <th></th>      <th>自然指数関数</th>      <th>自然対数関数</th>    </tr>  </thead>  <tbody>    <tr>      <th>0</th>      <td>y=exp(x)</td>      <td>y=log(x)</td>    </tr>    <tr>      <th>1</th>      <td>y=exp(-x)</td>      <td>y=log(-x)</td>    </tr>    <tr>      <th>2</th>      <td>y=-exp(x)</td>      <td>y=-log(x)</td>    </tr>    <tr>      <th>3</th>      <td>y=-exp(-x)</td>      <td>y=-log(-x)</td>    </tr>  </tbody></table>

* $(e^x)'=e^x$…積分の結果も$\int e^x dx = e^x$となる。
* $(e^{ax})'=ae^{ax}$…積分の結果は$\int e^{ax} dx = \frac{1}{a}e^{ax}$となる。
* $(\log(x))'=\frac{1}{x}$…これに対応する積分の結果は$\int \frac{1}{x} dx = \log(x)$となる。
* $(\log_{a}(x))'=\frac{1}{x \log(a)}$…一般の対数関数の場合。
* $(α^x)'=α^x \log(α)$…用例は上掲。<br>[【数理考古学】冪乗算の微積分](https://qiita.com/ochimusha01/items/0bd18370cfa96e83347a)

これに関連して

```math
\int \log(x) dx = x(log(x)-1)+C
```

pythonによる実装結果

```python
import sympy as sp
x,a= sp.symbols('x,a')

#原始関数
f=sp.log(x)

#微分
g=sp.integrate(f,x)

#積分
h=sp.diff(g,x,1)

#画面表示
sp.init_printing()
display(f) 
print(sp.latex(f))
display(g) 
print(sp.latex(g))
display(h) 
print(sp.latex(h))
```
原始関数

```math
\log{\left(x \right)}
```
積分結果

```math
x \log{\left(x \right)} - x
```
微分による復元

```math
\log{\left(x \right)}
```
そして虚冪演算の微積分は回転を通じて三角関数に関連付けられます。
[【Python演算処理】冪算と乗除算の関係について。](https://qiita.com/ochimusha01/items/8a66fcab8285a146265f)
![image.png](https://qiita-user-contents.imgix.net/https%3A%2F%2Fcdn-ak.f.st-hatena.com%2Fimages%2Ffotolife%2Fo%2Fochimusha01%2F20210504%2F20210504162020.png?ixlib=rb-4.0.0&auto=format&gif-q=60&q=75&w=1400&fit=max&s=a6955a0a4d3736fc57aed6753517f458)
![image.png](https://qiita-user-contents.imgix.net/https%3A%2F%2Fcdn-ak.f.st-hatena.com%2Fimages%2Ffotolife%2Fo%2Fochimusha01%2F20210504%2F20210504162153.png?ixlib=rb-4.0.0&auto=format&gif-q=60&q=75&w=1400&fit=max&s=2335f454fb0ffe93e583434dcb9060e4)

* $e^{ix}\frac{d^n}{dθ^n}=(i e^{i x}(-\log ix),-e^{ix},-i e^{i x}(\log ix),e^{ix},…)$
* $\int \int \int … \int e^{ix}(dθ)=(- i e^{i x}(\log ix),- e^{i x},i e^{i x}(-\log ix),e^{ix},…)$

さてドン尻に控えしは…

* $\(\log(x+\sqrt{x^2+a})\)'=\frac{1}{\sqrt{x^2+a}}$…何に用いられる式か不明だが、大学受験の出題範囲に入っているっぽくてネットに解答例が頻出。

pythonによる実装結果

```python
import sympy as sp
x,a= sp.symbols('x,a')

#原始関数
f=sp.log(x+sp.sqrt(x**2+a))

#微分
g=sp.diff(f,x,1)

#積分
h=sp.integrate(g,x)

#画面表示
sp.init_printing()
display(f) 
print(sp.latex(f))
display(g) 
print(sp.latex(g))
display(h) 
print(sp.latex(h))
```
原始関数

```math
\log{\left(x + \sqrt{a + x^{2}} \right)}
```
微分結果

```math
\frac{\frac{x}{\sqrt{a + x^{2}}} + 1}{x + \sqrt{a + x^{2}}}
```
積分による復元

```math
\log{\left(x + \sqrt{a + x^{2}} \right)}
```
<font color=green>Sympy完全に猫跨ぎ…微塵のモチベーションも感じられません。この猫っぽい仕草が本当に愛らしいのですよ。</font>

##三角関数の微分
[【Python演算処理】冪算と乗除算の関係について。](https://qiita.com/ochimusha01/items/8a66fcab8285a146265f)

>**三角関数**の微分（<font color=blue>それぞれが+90度の変遷に対応</font>）

```math
\cos(θ)\frac{d^n}{dθ^n}=(-\sin(θ),-\cos(θ),\sin(θ),\cos(θ),…)\\
\sin(θ)\frac{d^n}{dθ^n}=(\cos(θ),-\sin(θ),-\cos(θ),\sin(θ),…)
```

>**三角関数**の積分（<font color=blue>それぞれが-90度の変遷に対応</font>）

```math
\int \int \int … \int \cos(θ)(dθ)=(\sin(θ),-\cos(θ),-\sin(θ),\cos(θ),…)\\
\int \int \int … \int \sin(θ)(dθ)=(-\cos(θ),-\sin(θ),\cos(θ),\sin(θ),…)
```

[【Python演算処理】等速円運動についての物理学と数学の立場の違い?](https://qiita.com/ochimusha01/items/d5526764c7d1199671df)
![image.gif](https://camo.qiitausercontent.com/fa349dd64b85d0da84f63ee475e79148eb983e70/68747470733a2f2f73746167696e672d71696974612d757365722d636f6e74656e74732e696d6769782e6e65742f687474707325334125324625324671696974612d757365722d636f6e74656e74732e696d6769782e6e6574253246687474707325323533412532353246253235324671696974612d696d6167652d73746f72652e73332e61702d6e6f727468656173742d312e616d617a6f6e6177732e636f6d2532353246302532353246353937373732253235324631616332643937622d323937382d346230362d623131382d3930363732613536663538392e67696625334669786c696225334472622d312e322e322532366175746f253344666f726d61742532366769662d71253344363025323671253344373525323677253344313430302532366669742533446d61782532367325334430663266393364333236343937366634393336343733393831346436336537383f69786c69623d72622d342e302e30266175746f3d666f726d6174266769662d713d363026713d373526773d31343030266669743d6d617826733d3834616130626365363061653733623835616234313639343162323261653264)

* $\cos(θ)=\frac{e^{θi}+e^{-θi}}{2}$
* $\sin(θ)=\frac{e^{θi}-e^{-θi}}{2i}$

[tanxと1/tanxの微分公式のいろいろな証明](https://manabitimes.jp/math/1191)

* $(\tan x)'=\frac{1}{\cos^2(x)}$…タンジェントの微分結果

pythonによる実装結果

```python
import sympy as sp
x,a= sp.symbols('x,a')

#原始関数
f=sp.tan(x)

#微分
g=sp.diff(f,x,1)

#積分
h=sp.integrate(g,x)

#画面表示
sp.init_printing()
display(f) 
print(sp.latex(f))
display(g) 
print(sp.latex(g))
display(h) 
print(sp.latex(h))
```
原始関数

```math
\tan{\left(x \right)}
```
微分結果

```math
\tan^{2}{\left(x \right)} + 1
```
積分による復元

```math
\frac{\sin{\left(x \right)}}{\cos{\left(x \right)}}
```
<font color=green>Sympy「にゃー(タンジェントなんてオワコン)」圧を感じます。</font>
​
* $(\cot x)'=(\frac{1}{\tan(x)})'=-\frac{1}{\sin^2(x)}$…コタンジェント(タンジェントの逆数)の微分結果

pythonによる実装結果

```python
import sympy as sp
x,a= sp.symbols('x,a')

#原始関数
f=sp.cot(x)

#微分
g=sp.diff(f,x,1)

#積分
h=sp.integrate(g,x)

#画面表示
sp.init_printing()
display(f) 
print(sp.latex(f))
display(g) 
print(sp.latex(g))
display(h) 
print(sp.latex(h))
```
原始関数

```math
\cot{\left(x \right)}
```
微分結果

```math
- \cot^{2}{\left(x \right)} - 1
```
積分による復元

```math
\frac{\cos{\left(x \right)}}{\sin{\left(x \right)}}
```
<font color=green>Sympy「にゃー(コタンジェントなんてオワコン)」圧を感じます。</font>

[三角関数 - Wikipedia](https://ja.wikipedia.org/wiki/%E4%B8%89%E8%A7%92%E9%96%A2%E6%95%B0)
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/e5b88116-dcd7-6704-0c78-5456c9a84e0a.png)
三角関数のグラフ: **Sine**（<font color=blue>青実線</font>）、 **Cosine**（<font color=green>緑実線</font>）、 **Tangent**（<font color=red>赤実線</font>）、 **Cosecant**（<font color=blue>青点線</font>）、 **Secant**（<font color=green>緑点線</font>）、 **Cotangent**（<font color=red>赤点線</font>）

[逆三角関数の重要な性質まとめ](https://manabitimes.jp/math/787)

* $(\mathrm{Arcsin}(x))′=\frac{1}{\sqrt{1-x^2}}$

pythonによる実装結果

```python
import sympy as sp
x,a= sp.symbols('x,a')

#原始関数
f=sp.asin(x)

#微分
g=sp.diff(f,x,1)

#積分
h=sp.integrate(g,x)

#画面表示
sp.init_printing()
display(f) 
print(sp.latex(f))
display(g) 
print(sp.latex(g))
display(h) 
print(sp.latex(h))
```
原始関数

```math
\operatorname{asin}{\left(x \right)}
```
微分結果

```math
\frac{1}{\sqrt{1 - x^{2}}}
```
積分による復元

```math
\operatorname{asin}{\left(x \right)}
```

* $(\mathrm{Arccos}(x))'=-\frac{1}{\sqrt{1-x^2}}$

pythonによる実装結果

```python
import sympy as sp
x,a= sp.symbols('x,a')

#原始関数
f=sp.acos(x)

#微分
g=sp.diff(f,x,1)

#積分
h=sp.integrate(g,x)

#画面表示
sp.init_printing()
display(f) 
print(sp.latex(f))
display(g) 
print(sp.latex(g))
display(h) 
print(sp.latex(h))
```
原始関数

```math
\operatorname{acos}{\left(x \right)}
```
微分結果

```math
- \frac{1}{\sqrt{1 - x^{2}}}
```
積分による復元

```math
\operatorname{acos}{\left(x \right)}
```

* $(\mathrm{Arctan}(x))'=\frac{1}{1+x^2}$

pythonによる実装結果

```python
import sympy as sp
x,a= sp.symbols('x,a')

#原始関数
f=sp.atan(x)

#微分
g=sp.diff(f,x,1)

#積分
h=sp.integrate(g,x)

#画面表示
sp.init_printing()
display(f) 
print(sp.latex(f))
display(g) 
print(sp.latex(g))
display(h) 
print(sp.latex(h))
```
原始関数

```math
\operatorname{atan}{\left(x \right)}
```
微分結果

```math
\frac{1}{x^{2} + 1}
```
積分による復元

```math
\operatorname{atan}{\left(x \right)}
```
* $(\mathrm{Arccot}(x))'=-\frac{1}{1+x^2}$

pythonによる実装結果

```python
import sympy as sp
x,a= sp.symbols('x,a')

#原始関数
f=sp.acot(x)

#微分
g=sp.diff(f,x,1)

#積分
h=sp.integrate(g,x)

#画面表示
sp.init_printing()
display(f) 
print(sp.latex(f))
display(g) 
print(sp.latex(g))
display(h) 
print(sp.latex(h))
```
原始関数

```math
\operatorname{acot}{\left(x \right)}
```
微分結果

```math
-\frac{1}{x^{2} + 1}
```
積分による復元

```math
\operatorname{acot}{\left(x \right)}
```

[逆三角関数 - Wikipedia](https://ja.wikipedia.org/wiki/%E9%80%86%E4%B8%89%E8%A7%92%E9%96%A2%E6%95%B0)
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/c1fb71d5-7036-65e0-e132-3126a5e5b688.png)
平面上の直交座標系で図示された**arcsin(x)**（<font color=red>赤</font>）と **arccos(x)**（<font color=blue>青</font>）の通常の定義における主値。
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/c2075173-f043-200c-b46d-20fde9f4ead7.png)
平面上の直交座標系で図示された**arctan(x)**（<font color=red>赤</font>）と **arccot(x)** （<font color=blue>青</font>）の通常の定義における主値。

[双曲線関数にまつわる重要な公式まとめ](https://manabitimes.jp/math/617)

* $(\sinh(x))'=\cosh(x)$

pythonによる実装結果

```python
import sympy as sp
x,a= sp.symbols('x,a')

#原始関数
f=sp.sinh(x)

#微分
g=sp.diff(f,x,1)

#積分
h=sp.integrate(g,x)

#画面表示
sp.init_printing()
display(f) 
print(sp.latex(f))
display(g) 
print(sp.latex(g))
display(h) 
print(sp.latex(h))
```
原始関数

```math
\sinh{\left(x \right)}
```
微分結果

```math
\cosh{\left(x \right)}
```
積分による復元

```math
\sinh{\left(x \right)}
```

* $(\cosh x)'=\sinh(x)$

pythonによる実装結果

```python
import sympy as sp
x,a= sp.symbols('x,a')

#原始関数
f=sp.cosh(x)

#微分
g=sp.diff(f,x,1)

#積分
h=sp.integrate(g,x)

#画面表示
sp.init_printing()
display(f) 
print(sp.latex(f))
display(g) 
print(sp.latex(g))
display(h) 
print(sp.latex(h))
```
原始関数

```math
\cosh{\left(x \right)}
```
微分結果

```math
\sinh{\left(x \right)}
```
積分による復元

```math
\cosh{\left(x \right)}
```

* $(\tanh x)'=1-\tanh^2(x)$

pythonによる実装結果

```python
import sympy as sp
x,a= sp.symbols('x,a')

#原始関数
f=sp.tanh(x)

#微分
g=sp.diff(f,x,1)

#積分
h=sp.integrate(g,x)

#画面表示
sp.init_printing()
display(f) 
print(sp.latex(f))
display(g) 
print(sp.latex(g))
display(h) 
print(sp.latex(h))
```
原始関数

```math
\tanh{\left(x \right)}
```
微分結果

```math
1 - \tanh^{2}{\left(x \right)}
```
積分による復元

```math
\tanh{\left(x \right)}
```
[双曲線関数 - Wikipedia](https://ja.wikipedia.org/wiki/%E5%8F%8C%E6%9B%B2%E7%B7%9A%E9%96%A2%E6%95%B0)
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/7fff5f1a-eb25-b2f3-47af-c6bf357e80b4.png)

そして受験数学ではほとんど黙殺されているに等しいこの双曲線関数こそが機械学習概念の発端となったシグモイド曲線へとつながっていくのです。
[【Python演算処理】冪算と乗除算の関係について。](https://qiita.com/ochimusha01/items/8a66fcab8285a146265f)

```math
{{\displaystyle \varsigma _{a}(x)={\frac {1}{1+e^{-ax}}}={\frac {\tanh(ax/2)+1}{2}}}
}
```
![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fcdn-ak.f.st-hatena.com%2Fimages%2Ffotolife%2Fo%2Fochimusha01%2F20210726%2F20210726130329.gif?ixlib=rb-4.0.0&auto=format&gif-q=60&q=75&w=1400&fit=max&s=eb571ea3371cc5a19d1bef35354350c0)

##円と球を巡る微積分
[【Python演算処理】半径・直径・円周長・円の面積・球の表面積・球の体積の計算上の往復](https://qiita.com/ochimusha01/items/e180451a210227bb19a5)
[【Python演算処理】「N次球概念導入による円/球関係の数理の統合」？](https://qiita.com/ochimusha01/items/f7c9077eccfdbd0dd47a)

* **円の面積**$πr^2$$を**半径**rで微分すると**円周の長さ**2πrとなる。逆に**円周の長さ**2πrを**半径**rで積分すると**円の面積**$πr^2$になる。
* **球の体積**$4/3πr^3$を**半径**rで微分すると**球の表面積**$4πr^2$、逆に**球の表面積**$4πr^2$を**半径**rで積分すると**球の体積**$4/3πr^3$となる。
* **球の表面積**$4πr^2$を**半径**rで微分すると**円周長**2πr、逆に**円周長**2πrを半径rで積分すると**球の表面積**$4πr^2$となる。

>どうして微積分可能範囲がこんな形で「**離散的**」に現れるのか。それを説明するには階乗の概念を複素数全体に拡張した**Γ関数**の導入を必要不可欠とする。要するにこれらはそれぞれ構造全体を説明する巨大な単一方程式の部分抽出に過ぎないのである。

これで概ね大学受験数学までに習う微分の範囲は網羅した？　とはいえ、これらの知識を総動員したとしても「**微分とは何か**」自体はまるで見えてきません。それを理解するにはまた別のアプローチが必要となる様なのです。

#微分とは何か？(<font color=blue>理論編</font>)

[涌井良幸「高校生からわかるフーリエ解析(2019年)」](https://www.amazon.co.jp/dp/B081RVK1D1/ref=dp-kindle-redirect?_encoding=UTF8&btkr=1)

>関数y=f(x)の**定義域**(x)内の点αにおいて$Δx \to 0$の時、$\frac{f(α+Δx)-f(α)}{Δx}$が一定の値に収束する時「**関数**y=f(x)**は**x=α**で微分可能である**」とし、この一定の値をx=aにおける**微分係数**と呼びf'(α)で表す。

```math
f'(a)=\lim_{Δx \to 0}\frac{f(α+Δx)-f(α)}{Δx}
```

* Δx→0はΔx>0から0に近づけても、Δx<0から0に近づけても、また正負を交互に変化させつつ0に近づけても良い。いずれにせよα+Δxの数直線上の絶対値がαに漸近する事が重要なのである。

>関数f(x)がx=aの十分近くで定義されている時、Δy=f(a+Δx)-f(a)だとすると$\frac{Δy}{Δx}=\frac{f(a+Δx)-f(a)}{Δx}$は関数f(x)上の2点a,bを通る直線lの傾きを表す。すなわち「**x=αで微分可能である**」という事は、点a,bを限りなく近づけた時直線lの傾きが一定に近づく一点が存在する事を意味する。すなわち点A(a,f(a))(あるいは点B(b,f(b)))付近のグラフが滑らかで連続しておりドンドン拡大していくと直線とみなせる状態に到達する、すなわち「<font color=purple>拡大して直線と見做せた時に**微分可能**と見做され、その時の直線の傾きが**微分係数**となる</font>」のである。

```math
\lim_{Δx \to 0}\frac{Δy}{Δx}=\lim_{Δx \to 0}\frac{f(α+Δx)-f(α)}{Δx}
```

* ここで導入された「**十分に近い**」概念はおそらく抽象数学における「**近傍**(Neighbourhood)」概念と密接に関係している。<br>[近傍とは - コトバンク](https://kotobank.jp/word/%E8%BF%91%E5%82%8D-54489)


>ここで集合x=(a,b,c,…)にその微分係数(f'(a),f'(b),f'(c),…)を対応させる関数f'(x)を考え、**原始関数**(Primitive Function)f(x)の**導関数**(Derivative Function)と呼ぶ事にする。

```math
f'(x)=\lim_{Δx \to 0}\frac{Δy}{Δx}=\lim_{Δx \to 0}\frac{f(x+Δx)-f(x)}{Δx}
```

>**導関数**の記号は$y',\frac{dy}{dx},\frac{dy}{dx}f(x)$などでも良い。いずれにせよ関数f(x)の**導関数**を求める事を「**微分する**」と呼ぶ。

何だか分かるような分からない様な抽象的な定義ですね。しかしまぁ、ある意味それこそが厳密さを極めた現代数学の特徴とも?　そして…
[高校数学からヤコビアンに至るまで](https://qiita.com/bellbind/items/bd0411957e2363456ea1)

* 「**合成関数の微分**」概念は最小単位Δを共有する「**チェーンメール変換**」イメージを経て偏微分方程式の各変数を関数とする多重積分に拡張される。
* さらにここでいう共通最小単位Δの定義がイプシロンデルタ論法によって拡張される。
* こうして強化された共通最小単位Δの概念を直交座標と極座標の変換に導入した結果がヤコビアン行列となる。<br>![image.png](https://qiita-user-contents.imgix.net/https%3A%2F%2Fqiita-image-store.s3.amazonaws.com%2F0%2F79115%2F4c884000-71de-a6c4-9780-9d5cf1dae547.png?ixlib=rb-4.0.0&auto=format&gif-q=60&q=75&w=1400&fit=max&s=9a303a10bce8717703a5663cee1b9ff6)

大雑把にいうとこうした考え方が大学で習う「**専門数学**」の領域に入ってくる様です。まだまだ全体像は掴めないままですが、何となく道筋は掴めてきた感じがします。

* これまでの投稿で述べてきた 「**同心集合**(<font color=blue>Concentric Set</font>)=<font color=purple>中心0の1点でのみ固定された数直線上に置かれたそれぞれの値が描く同心円/同心球面t</font>」概念自体がここに登場しないのは、「**空間全体を構成する偏微分方程式を各変数ごとの多重積分に変換する過程**」においては、まず半径r=スカラー量を抽出し、しかるのちに一次トーラス=半径1の単位円(円周長2π)の直積によって全体を把握していくのが**直交座標**(r(スカラー量),φ(-π~+π),θ(-π~+π),…)の流儀だから。
* もちろんこの考え方だけでは「**メルカトル座標系**(φ=-π~+π,θ=0~π)と**二次元トーラス座標系**(φ=-π~+π,θ=-π~+π)」「**ヤコビアン座標系**(r=0~1,φ=-π~+π,θ=0~π)と**三次元トーラス**(四元数)**座標系**(i,j,kそれぞれが-1~+1の稼働範囲)」それぞれのメリットとデメリットが上手く説明出来ず、これからの克服課題として急浮上してくる。

現時点での自認は大体こんな感じ?

#微分とは何か？(<font color=blue>実践編</font>)

[【初心者向け】自然指数関数と自然対数関数の交点について。](https://qiita.com/ochimusha01/items/77b41156340b6f694373)

>導関数がＸ軸とY軸の間に間に作る面積は原始関数の値と一致する。

改めて具体的に「**微分とは何か?**」再掌握するには、この辺りが鍵となってくる様です。

###一次方程式<font color=purple>y=x</font>

```math
y=x\\
y'=1
```
導関数の導出過程

```math
f'(x)=\lim_{Δx \to 0}\frac{Δy}{Δx}=\lim_{Δx \to 0}\frac{f(x+Δx)-f(x)}{Δx}\\
=\lim_{Δx \to 0}\frac{(x+Δx)-x}{Δx}=\lim_{Δx \to 0}\frac{Δx}{Δx}=1
```
x=4の時、原始関数y=x→y=4、導関数y=1→y=1×4=4
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/b505fa28-c089-fcb1-edce-e23c4638e45b.png)

###二次方程式<font color=purple>y=x^2</font>

```math
y=\frac{1}{2}x^2 \\
y'=x
```
導関数の導出過程

```math
f'(x)=\lim_{Δx \to 0}\frac{Δy}{Δx}=\lim_{Δx \to 0}\frac{f(x+Δx)-f(x)}{Δx}\\
=\lim_{Δx \to 0}\frac{(x+Δx)^2-x^2}{Δx}\\
=\lim_{Δx \to 0}\frac{(x^2+2xΔx+Δx^2)-x^2}{Δx}\\
=\lim_{Δx \to 0}\frac{2xΔx+Δx^2}{Δx}\\
=\lim_{Δx \to 0}2x+Δx=2x
```
【参考】**二次形式行列**における「**マスキング**」

```python
import sympy as sp
a,b = sp.symbols('a,b')
#Rows
x=sp.Matrix([[a,1]])
#Columns
y=sp.Matrix([[a],[1]])
a00=sp.Matrix([[1,1],[1,1]])
a11=sp.Matrix([[1,0],[0,0]])
a12=sp.Matrix([[0,1],[0,0]])
a21=sp.Matrix([[0,0],[1,0]])
a22=sp.Matrix([[0,0],[0,1]])
adp=sp.Matrix([[1,0],[0,1]])
adm=sp.Matrix([[0,1],[1,0]])
sp.init_printing()
display(x) 
display(a00) 
display(a11) 
display(a12) 
display(a21) 
display(a22) 
display(adp)
display(adm)
display(y) 
display(sp.simplify(x*a00*y)) 
print(sp.latex(x)+sp.latex(a00)+sp.latex(y)+"="+sp.latex(sp.simplify(x*a00*y)))
display(sp.simplify(x*a11*y)) 
print(sp.latex(x)+sp.latex(a11)+sp.latex(y)+"="+sp.latex(sp.simplify(x*a11*y)))
display(sp.simplify(x*a12*y)) 
print(sp.latex(x)+sp.latex(a12)+sp.latex(y)+"="+sp.latex(sp.simplify(x*a12*y)))
display(sp.simplify(x*a21*y)) 
print(sp.latex(x)+sp.latex(a21)+sp.latex(y)+"="+sp.latex(sp.simplify(x*a21*y)))
display(sp.simplify(x*a22*y)) 
print(sp.latex(x)+sp.latex(a22)+sp.latex(y)+"="+sp.latex(sp.simplify(x*a22*y)))
display(sp.simplify(x*adp*y)) 
print(sp.latex(x)+sp.latex(adp)+sp.latex(y)+"="+sp.latex(sp.simplify(x*adp*y)))
display(sp.simplify(x*adm*y)) 
print(sp.latex(x)+sp.latex(adm)+sp.latex(y)+"="+sp.latex(sp.simplify(x*adm*y)))
```

```math
\left[\begin{matrix}a & 1\end{matrix}\right]\left[\begin{matrix}1 & 1\\1 & 1\end{matrix}\right]\left[\begin{matrix}a\\1\end{matrix}\right]=\left[\begin{matrix}a \left(a + 1\right) + a + 1\end{matrix}\right]\\=[a^2+2a+1]=[(a+1)^2]…①\\

\left[\begin{matrix}a & 1\end{matrix}\right]\left[\begin{matrix}1 & 0\\0 & 0\end{matrix}\right]\left[\begin{matrix}a\\1\end{matrix}\right]=\left[\begin{matrix}a^{2}\end{matrix}\right]…②\\

\left[\begin{matrix}a & 1\end{matrix}\right]\left[\begin{matrix}0 & 1\\0 & 0\end{matrix}\right]\left[\begin{matrix}a\\1\end{matrix}\right]=\left[\begin{matrix}a\end{matrix}\right]…③\\

\left[\begin{matrix}a & 1\end{matrix}\right]\left[\begin{matrix}0 & 0\\1 & 0\end{matrix}\right]\left[\begin{matrix}a\\1\end{matrix}\right]=\left[\begin{matrix}a\end{matrix}\right]…④\\

\left[\begin{matrix}a & 1\end{matrix}\right]\left[\begin{matrix}0 & 0\\0 & 1\end{matrix}\right]\left[\begin{matrix}a\\1\end{matrix}\right]=\left[\begin{matrix}1\end{matrix}\right]…⑤\\

\left[\begin{matrix}a & 1\end{matrix}\right]\left[\begin{matrix}1 & 0\\0 & 1\end{matrix}\right]\left[\begin{matrix}a\\1\end{matrix}\right]=\left[\begin{matrix}a^{2} + 1\end{matrix}\right]…⑥\\

\left[\begin{matrix}a & 1\end{matrix}\right]\left[\begin{matrix}0 & 1\\1 & 0\end{matrix}\right]\left[\begin{matrix}a\\1\end{matrix}\right]=\left[\begin{matrix}2 a\end{matrix}\right]…⑦
```
この考え方でいうと⑦のタイプのマスキングに該当する様です？

原関数がx=4(y=8)の時、導関数は$\frac{4×4}{2}=8$
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/05a64110-eeb7-b98f-236f-c5fd56e79687.png)

これは物理学の基礎に対応しています。

[落下運動](http://www.daiichi-g.co.jp/rika/subtextbook/data/46658/46658page14-19.pdf)
[運動エネルギーの基本式](https://www.structure.jp/column8/topic802.html)

>**自由落下運動**(<font color=blue>Free Fall Movement</font>)とは以下の様な運動を指す。
>①**重力加速度**<font color=green>0.98</font>(<font color=blue>=:1,単位m/s^2</font>)の力のかかる向きが進行方向と完全に一致する。
>②**初速**<font color=green>0</font>(<font color=blue>m/s</font>)で**時間**<font color=green>t</font>=<font color=green>0</font>(<font color=blue>s</font>)の時の**位置**<font color=green>y</font>=<font color=green>0</font>(<font color=blue>m</font>)とする。
>この時、
>①**時間**<font color=green>t</font>=<font color=green>x</font>(<font color=blue>s</font>)時点の**速度**<font color=purple>v</font>=<font color=purple>gt</font>(<font color=blue>m/s</font>)となる。<br>![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/06a2147b-2101-4044-469e-27a294c5ae53.png)
>②**時間**<font color=green>t</font>=<font color=green>x</font>(s</font>)時点の**位置**<font color=purple>y</font>=<font color=purple>1/2×gt^2</font>(<font color=blue>m</font>)となる。<br>![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/09eb6fb0-f1fa-80c4-40ea-63aa3ac75b3f.png)
>③<font color=purple>1/2</font>なる数字は(<font color=blue>三角形の面積(h(高さ)+v(幅))/2を求める要領で</font>)**時間**<font color=green>t</font>=<font color=green>0</font>時点の**速度**<font color=green>v</font>=<font color=green>0</font>(<font color=blue>m/s</font>)から**時間**<font color=green>x</font>(<font color=blue>s</font>)時点までの**速度**<font color=green>v</font>=<font color=green>gt</font>(<font color=blue>m/s</font>)までの**平均速度**(<font color=blue>m/s</font>)を求める過程で現れる。
従って<font color=purple>v^2</font>=<font color=purple>2gy</font>が成立。

ここで重力係数0.98を1に単純化すると見慣れた微分式$y=x^2→y=2x$の式が現れる訳ですね。そこでは単位m/sが面積に対応します。

* **速度**s(m/s)=**時間**t(s)×**重力係数**g(m)
* **初速**$v_0=0(m/s)$から**時間t時点の速度**$v_t=gt(m/s)$の**平均速度**vは$\frac{gt}{2}(m/s)$となる。さらにg=1と置くと$\frac{t}{2}(m/s)$となる。
* x(時間t)時点の位置y(m)は**平均速度**v(m/s)×**時間**t(s)=$\frac{gt^2}{2}(m)$となる。 さらにg=1と置くと$\frac{t^2}{2}(m)$となる。
* 従って重力加速度=1と置くと$x^2(t)=2y(m)$となる。

さらにいうと、ここで重力係数$g(0.98) \fallingdotseq 1$と置く事は、それを**実数環**(Real Ring)すなわち**実数列**(Real Sequence)を**垂直軸**(Vertical Axis)z、半径1の単位円を**水平軸**(Horizontal Axis)xyに取った**円筒座標系**(Cylindrical Coordinate System)へと射影(Projection)する事を意味します。そうすると嬉しい事にz軸(実数列)の値がそのまま勾配($\frac{z}{1}$)に対応する展開を迎える訳ですね。

###三次方程式<font color=purple>y=x^3</font>

```math
y=\frac{1}{3}x^3 \\
y'=x^2
```
導関数の導出過程

```math
f'(x)=\lim_{Δx \to 0}\frac{Δy}{Δx}=\lim_{Δx \to 0}\frac{f(x+Δx)-f(x)}{Δx}\\
=\lim_{Δx \to 0}\frac{(x+Δx)^3-x^3}{Δx}\\
=\lim_{Δx \to 0}\frac{(x^3+3x^2Δx+3xΔx^2+Δx^3)-x^3}{Δx}\\
=\lim_{Δx \to 0}\frac{3x^2Δx+3xΔx^2+Δx^3}{Δx}\\
=\lim_{Δx \to 0}3x^2+3xΔx+Δx^2\\
=3x^2+3*0+0^2=3x^2
```
これへも上掲の「**マスキング**」という概念を適用しようとすると、突如として「**行列のN乗**」なる概念が襲いかかってきます。

* 要するに**二次方程式**同様、Sympyで以下の様な行列を準備し…

```math
(a+1)^3=a^3+3a^2+3a+1\\
\left[\begin{matrix}\left[ a^{3}, \  a^{2}\right] & \left[ a^{2}, \  a^{2}\right]\\\left[ a, \  a\right] & \left[ a, \  1\right]\end{matrix}\right]
```
* 以下の様なマスキング演算"・"によって②が微分結果に一致する事を示したかったのですが、その演算"・"に該当する操作がどうしても見つからないのです(<font color=blue>ちなみに4次方程式以降もこの関係は続くが、二項定理と微分の関係について触れた文章をこれまで一度も見た事がなく、誰の関心も集めてこなかった袋小路の匂いがプンプンと…</font>)。<br>[【Python演算処理】パスカルの三角形と二項定理または二項展開](https://qiita.com/ochimusha01/items/d668b19d5389b52e0dd2)

```math
\left[\begin{matrix}\left[ a^{3}, \  a^{2}\right] & \left[ a^{2}, \  a^{2}\right]\\\left[ a, \  a\right] & \left[ a, \  1\right]\end{matrix}\right]・\left[\begin{matrix}\left[ 1, \  1\right] & \left[ 1, \  1\right]\\\left[ 1, \  1\right] & \left[ 1, \  1\right]\end{matrix}\right]=a^3+3a^2+3a+1\\
\left[\begin{matrix}\left[ a^{3}, \  a^{2}\right] & \left[ a^{2}, \  a^{2}\right]\\\left[ a, \  a\right] & \left[ a, \  1\right]\end{matrix}\right]・\left[\begin{matrix}\left[ 1, \  0\right] & \left[ 0, \  0\right]\\\left[ 0, \  0\right] & \left[ 0, \  0\right]\end{matrix}\right]=a^3…①\\
\left[\begin{matrix}\left[ a^{3}, \  a^{2}\right] & \left[ a^{2}, \  a^{2}\right]\\\left[ a, \  a\right] & \left[ a, \  1\right]\end{matrix}\right]・\left[\begin{matrix}\left[ 0, \  1\right] & \left[ 1, \  1\right]\\\left[ 0, \  0\right] & \left[ 0, \  0\right]\end{matrix}\right]=3a^2…②\\
\left[\begin{matrix}\left[ a^{3}, \  a^{2}\right] & \left[ a^{2}, \  a^{2}\right]\\\left[ a, \  a\right] & \left[ a, \  1\right]\end{matrix}\right]・\left[\begin{matrix}\left[ 0, \  0\right] & \left[ 0, \  0\right]\\\left[ 1, \  1\right] & \left[ 1, \  0\right]\end{matrix}\right]=3a…③\\
\left[\begin{matrix}\left[ a^{3}, \  a^{2}\right] & \left[ a^{2}, \  a^{2}\right]\\\left[ a, \  a\right] & \left[ a, \  1\right]\end{matrix}\right]・\left[\begin{matrix}\left[ 0, \  0\right] & \left[ 0, \  0\right]\\\left[ 0, \  0\right] & \left[ 0, \  1\right]\end{matrix}\right]=1…④
```

```python
import sympy as sp
a = sp.symbols('a')
a3=a**3
a2=a**2
a1=a
a0000=sp.Matrix([[[a3,a2],[a2,a2]],[[a1,a1],[a1,1]]])
a1111=sp.Matrix([[[1,1],[1,1]],[[1,1],[1,1]]])
a1000=sp.Matrix([[[1,0],[0,0]],[[0,0],[0,0]]])
a0100=sp.Matrix([[[0,1],[1,1]],[[0,0],[0,0]]])
a0010=sp.Matrix([[[0,0],[0,0]],[[1,1],[1,0]]])
a0001=sp.Matrix([[[0,0],[0,0]],[[0,0],[0,1]]])
sp.init_printing()
display(a0000) 
print(sp.latex(a0000))
display(a1111) 
print(sp.latex(a0000)+"・"+sp.latex(a1111)+"="+"a^3+3a^2+3a+1")
display(a1000) 
print(sp.latex(a0000)+"・"+sp.latex(a1000)+"="+"a^3")
display(a0100) 
print(sp.latex(a0000)+"・"+sp.latex(a0100)+"="+"3a^2")
display(a0010) 
print(sp.latex(a0000)+"・"+sp.latex(a0010)+"="+"3a")
display(a0001)
print(sp.latex(a0000)+"・"+sp.latex(a0001)+"="+"1")
```
そもそも行列演算では$(a+1)^3$を$(a+1)^2(a+1)$の形で扱う事しか考えてない様にも見えます。

[行列のn乗(3)](https://www.geisya.or.jp/~mwm48961/kou3/matrix_power.htm)

>高校数学の行列は旧教育課程の数学Cに含まれていたが，平成21年（2009年）告示の教育課程では数学Cはなくなっており数学Ⅲなどの他の科目にも行列は含まれなかったため，高校では行列計算は原則として習わない．
>
>高校の教育課程からなくなったということは大学入試問題では忠実に反映されるので，行列そのものを出題することはないが，高卒向けや大卒向けの就職試験となるとこの制限はあいまいになる．
>
>この頁では「高校の旧教育課程にあった数学Cのレベル」で「2×2行列に限定して」行列のｎ乗を扱う．

まさしく高校数学と大学数学の中間領域…

* それはそれとして、円錐の面積は円筒の体積の$\frac{1}{3}$となる。ここでも微積分との関係は連綿と継続している訳である。<br>![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fcdn-ak.f.st-hatena.com%2Fimages%2Ffotolife%2Fo%2Fochimusha01%2F20210412%2F20210412062849.gif?ixlib=rb-4.0.0&auto=format&gif-q=60&q=75&w=1400&fit=max&s=9d40f6ad1e1c4373a2244b3b32dac65c)<br>![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fcdn-ak.f.st-hatena.com%2Fimages%2Ffotolife%2Fo%2Fochimusha01%2F20210412%2F20210412062919.gif?ixlib=rb-4.0.0&auto=format&gif-q=60&q=75&w=1400&fit=max&s=b9344280100346dc40d9972fd3cd22b7)

```python
%matplotlib nbagg
import math as m
import cmath as c
import numpy as num
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
import matplotlib.animation as animation

#円柱データ作成
c0=num.linspace(0,m.pi*120,1201,endpoint = True)
s0=[]
for nm in range(len(c0)):
    s0.append(complex(m.cos(c0[nm]),m.sin(c0[nm])))
s1=num.array(s0)
z0=num.linspace(-1,1,1201,endpoint = True)

#「曲率」を計算
cv1=num.linspace(-1,1,1201,endpoint = True)
#cv1=num.sqrt(1-cv0**2)

#断面線(z)
cutz0=num.linspace(-1,1,61,endpoint = True)
cutz=cutz0[::-1]
cutx=num.abs(cutz)
cuty=num.repeat(0,61)

#単位円データ作成
u0=num.linspace(0,m.pi*2,61,endpoint = True)
u1=[]
for nm in range(len(u0)):
    u1.append(complex(m.cos(u0[nm]),m.sin(u0[nm])))
uc=num.array(u1)
uz0=num.repeat(-1,61)
uz1=num.repeat(-0,61)
uz2=num.repeat(1,61)

#グラフ表示
plt.style.use('default')
fig = plt.figure()
ax = Axes3D(fig)

#関数定義
def unit_cylinder(n):
    plt.cla()
    #円柱描画
    ax.plot(s1.real*cv1,s1.imag*cv1,z0,color="gray",lw=0.5)
    #スポーク描画
    #for num in range(len(uc)):
    #    ax.plot([0,uc[num].real],[0,uc[num].imag],[0,0],color="gray",lw=0.5)
    #for num in range(len(uc)):
    #    ax.plot([0,uc[num].real],[0,uc[num].imag],[1,1],color="gray",lw=0.5)
    for num in range(len(uc)):
        ax.plot([0,uc[num].real],[0,uc[num].imag],[0,0],color="olivedrab",lw=1)
    #単位円描画
    ax.plot(uc.real,uc.imag,uz0,color="red",lw=1)
    ax.plot(uc.real,uc.imag,uz1,color="green",lw=1)
    ax.plot(uc.real,uc.imag,uz2,color="blue",lw=1)
    #実数線追加
    ax.plot([0,0],[0,0],[-1,1],color="black",lw=1)
    ax.plot([0,1],[0,0],[-1,-1],color="red",lw=1)
    ax.plot([0,1],[0,0],[0,0],color="black",lw=1)
    ax.plot([0,1],[0,0],[1,1],color="blue",lw=1)
    ax.plot([1,1],[0,0],[-1,0],color="red",lw=1)
    ax.plot([1,1],[0,0],[0,1],color="blue",lw=1)
    #断面線描画
    ax.plot(cutx,cuty,cutz,color="black",lw=1)
    #諸元追加
    ax.set_ylim([-1.1,1.1])
    ax.set_xlim([-1.1,1.1])
    ax.set_zlim([-1.1,1.1])
    ax.set_title("Unit Cylinder")
    ax.set_xlabel("Real")
    ax.set_ylabel("Imaginal")
    ax.set_zlabel("Cycle")
    # グラフを回転(elev=0にすると水平表示に)
    ax.view_init(elev=25, azim=Time_code[n])

Time_code0=num.arange(0,360,6)  
Time_code=Time_code0[::-1] 
#unit_cylinder(len(s1))
#plt.show()

ani = animation.FuncAnimation(fig, unit_cylinder, interval=50,frames=len(Time_code))
ani.save("output512.gif", writer="pillow")
```
もしかしたら、この関係は4次方程式経由で四次元の世界へも継承されていく？ そんな可能性が示唆されたあたりで以下続報…