【数理考古学】ピタゴラスの定理あるいは三平方の定理からの出発

* 2020年03月06日 初投稿
* 2020年06月13日 Tex対応 

**ピタゴラスの定理**(<font color="blue">Pythagorean theorem</font>)あるいは**三平方の定理**(<font color="blue">Three-square theorem</font>)の証明から出発します。
![20190918061420.gif](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/c7cb8211-2ac7-8acb-c35a-573af944e931.gif)

```R
Three_square_theorem04<-function(x){
c0<-c(0,pi/2,-1*pi,-1*pi/2,0)
c0_cos<-cos(c0)
c0_sin<-sin(c0)
plot(c0_cos,c0_sin,type="l",main="Regular polygon rotation",xlab="Cos(θ)",ylab="Sin(θ)")
text(c0_cos,c0_sin, labels=c("a","b","c","d",""),col=c(rgb(1,0,0),rgb(1,0,0),rgb(1,0,0),rgb(1,0,0),rgb(1,0,0)),cex=c(2,2,2,2,2))
line_scale_01_cos<-seq(cos(c0[1]),cos(c0[2]),length=15)
line_scale_01_sin<-seq(sin(c0[1]),sin(c0[2]),length=15)
line_scale_02_cos<-seq(cos(c0[2]),cos(c0[3]),length=15)
line_scale_02_sin<-seq(sin(c0[2]),sin(c0[3]),length=15)
line_scale_03_cos<-seq(cos(c0[3]),cos(c0[4]),length=15)
line_scale_03_sin<-seq(sin(c0[3]),sin(c0[4]),length=15)
line_scale_04_cos<-seq(cos(c0[4]),cos(c0[1]),length=15)
line_scale_04_sin<-seq(sin(c0[4]),sin(c0[1]),length=15)
text(c(line_scale_01_cos[x],line_scale_02_cos[x],line_scale_03_cos[x],line_scale_04_cos[x]),c(line_scale_01_sin[x],line_scale_02_sin[x],line_scale_03_sin[x],line_scale_04_sin[x]), labels=c("ab","bc","cd","da",""),col=c(rgb(0,0,1),rgb(0,0,1),rgb(0,0,1),rgb(0,0,1),rgb(0,0,1)),cex=c(2,2,2,2,2))
#塗りつぶし
polygon(c(c0_cos[1],line_scale_04_cos[x],line_scale_01_cos[x],c0_cos[1]), #x
c(c0_sin[1],line_scale_04_sin[x],line_scale_01_sin[x],c0_sin[1]), #y
density=c(30), #塗りつぶす濃度
angle=c(45),　 　　　#塗りつぶす斜線の角度
col=rgb(0,1,0)) 　 

#塗りつぶし
polygon(c(c0_cos[2],line_scale_01_cos[x],line_scale_02_cos[x],c0_cos[2]),#x
c(c0_sin[2],line_scale_01_sin[x],line_scale_02_sin[x],c0_sin[2]),#y
density=c(30),#塗りつぶす濃度
angle=c(45),#塗りつぶす斜線の角度
col=rgb(0,1,0)) 　 

#塗りつぶし
polygon(c(c0_cos[3],line_scale_02_cos[x],line_scale_03_cos[x],c0_cos[3]),#x
c(c0_sin[3],line_scale_02_sin[x],line_scale_03_sin[x],c0_sin[3]), #y
density=c(30),#塗りつぶす濃度
angle=c(45),#塗りつぶす斜線の角度
col=rgb(0,1,0)) 　

#塗りつぶし
polygon(c(c0_cos[4],line_scale_03_cos[x],line_scale_04_cos[x],c0_cos[4]),#x
c(c0_sin[4],line_scale_03_sin[x],line_scale_04_sin[x],c0_sin[4]), #y
density=c(30), #塗りつぶす濃度
angle=c(45),#塗りつぶす斜線の角度
col=rgb(0,1,0)) 
}

#アニメーション
library("animation")
Time_Code=c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15)
saveGIF({
for (i in Time_Code){
　Three_square_theorem04(i)
}
}, interval = 0.1, movie.name = "TEST04.gif") 
```
#ピタゴラスの定理(三平方の定理)を証明

![20190919170846.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/c993c7f7-2812-10e4-5705-402141a017db.png)

1. **直角三角形**を４個集めると**正方形**が出来る。
2. 大きな正方形の**面積**<font color="green">S</font>は一辺が<font color="green">{a,ab}={b,bc}={c,cd}={d,da}</font>（<font color="blue">Aとする</font>）。かつ<font color="green">{ab,b}+{bc,c}+{cd,c}={da,a}</font>（<font color="blue">Bとする</font>）にして<font color="green">{a,ab}+{ab,b}={b,bc}+{bc+b}={c,cd}+{cd,d}={d,da}+{da,a}</font>なので<font color="green">S=(A+B)^2</font>。
3. また<font color="green">S</font>は一辺が<font color="green">{ab,bc}={bc,cd}={cd,da}={da=ab}</font>（<font color="blue">Cとする</font>）の正方形と<font color="green">ABの直角三角形4つ</font>の和でもあるので<font color="green">S=C^2+4*1/2*A*B</font>。よって<font color="green">$A^2+2AB+B^2=C^2+2AB$</font>。
4. これを整理すると<font color="green">$A^2+B^2=C^2$</font>となる。

#<font color="green">a^2+b^2=c^2</font>から<font color="green">sinθ^2+cosθ^2=1</font>へ

ところで、ここに登場する<font color="green">直角三角形ABC</font>は（<font color="blue">プログラム中での扱いを見ても分かる様に</font>）、**単位円**(<font color="blue">Unit Circle=半径１の円弧</font>)上においては「**円弧上の任意の点の中心からの距離が全て１となる**」事を表す式と重なってくるのです。
[【初心者向け】「円そのもの」の近似から派生した角度と経度の概念の起源](https://qiita.com/ochimusha01/items/80b70db93cbf2612bd18)
![20191028123539.gif](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/6ecec94f-f795-b69a-df0e-7cdc9fa8dde7.gif)

1. <font color="green">a=cos(θ)</font>
2. <font color="green">b=sin(θ)</font>
3. <font color="green">c=1</font>

従って同時に**任意の角度**<font color="green">θ</font>に対して以下が成立。

1. <font color="green">$cos(θ)^2+sin(θ)^2=1$</font>

![20190916055039.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/fa5fb658-6452-326c-9f58-1ee2f3767629.png)

```R
f0<-function(x){cos(x)^2}
plot(f0,col=rgb(0,0,1),xlim=c(0,2*pi),ylim=ｃc(0,1),main="Three square theorem",xlab="radian(π)",ylab="Amplitude")
par(new=T) #上書き
f1<-function(x){sin(x)^2}
plot(f1,col=rgb(0,1,0),xlim=c(0,2*pi),ylim=c(0,1),main="",xlab="",ylab="")
par(new=T) #上書き
f2<-function(x){cos(x)^2+sin(x)^2}
plot(f2,,col=rgb(1,0,0),xlim=c(0,2*pi),ylim=c(0,1),main="",xlab="",ylab="")
legend("bottomright", legend=c("y=cos(θ)^2","y=sin(θ)^2","y=cos(θ)^2+sin(θ)^2"), lty=c(1,1,1), col=c(rgb(0,0,1),rgb(0,1,0),rgb(1,0,0)))
```

#外接/内接多角形が生成する同心円集合
ところで以下の図は正方形における内接円/外接円の半径の推移を示唆してもいるのです。
![20190918061420.gif](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/c7cb8211-2ac7-8acb-c35a-573af944e931.gif)
例えば半径1の単位円に内接する正方形の辺長は$\sqrt{2}$であり、さらにこれに内接する正方形の辺長が最小となるポイントに注目すると…

```math
\cos(\frac{π}{4})=\sin(\frac{π}{4})=\frac{\sqrt{2}}{2}
```
これはまさに外接/内接正方形の比率連鎖から生じる同心円集合の一部に他なりません。その半径rの推移は以下となり、この比率はそれぞれの外接/内接正方形の辺長/対角線長としても現れてくる形となります。
![img.png](https://camo.qiitausercontent.com/f0d3bf8ddc9a019b0d4e786666b1a2e0cade6f9d/68747470733a2f2f71696974612d696d6167652d73746f72652e73332e61702d6e6f727468656173742d312e616d617a6f6e6177732e636f6d2f302f3539373737322f38623639616562642d313661342d373962302d323565352d3030386234646661393730652e706e67)

ここで単位円の**半径**(Radius)1を$r_0$,単位円に内接する正方形の**対角線**(Diagonal)2を$d_0$、その**辺長**(Edge)$\sqrt{2}$を$e_0$と置くと…

```math
r_n,e_n,d_n(n \in \mathbb{Z})=\\
(…,\\
r_{-4},e_{-5},d_{-6}=2^{-2}=\frac{1}{4},\\
r_{-3},e_{-4},d_{-5}=2^{-1.5}=\frac{\sqrt{2}}{4},\\
r_{-2},e_{-3},d_{-4}=2^{^-1}=\frac{1}{2},\\
r_{-1},e_{-2},d_{-3}=2^{0.5}=\frac{\sqrt{2}}{2},\\
r_0,e_{-1},d_{-2}=2^{0}=1,\\
r_1,e_0,d_{-1}=2^{0.5}=\sqrt{2},\\
r_2,e_1,d_0=2^{1}=2,\\
r_3,e_2,d_1=2^{1.5}=2\sqrt{2},\\
r_4,e_3,d_2=2^{2}=4\\
…）
```
[【オイラーの多面体定理と正多面体】内接円/球面の半径と外接円/球面の半径の狭間](https://qiita.com/ochimusha01/items/1312cc54694cf1a27e17)

<!-- html table generated in R 3.5.2 by xtable 1.8-4 package -->
<!-- Sat May  9 17:05:36 2020 -->
<table border=1>
<tr> <th>  </th> <th> Target_size </th> <th> Target_names </th> <th> Target_values </th>  </tr>
  <tr> <td align="right"> 1 </td> <td> 2^-1 </td> <td> 2^-1(2^-0.5r) </td> <td> 0.5 </td> </tr>
  <tr> <td align="right"> 2 </td> <td> 2^-1 </td> <td> 2^-1d=2^-1*4 </td> <td> 1 </td> </tr>
  <tr> <td align="right"> 3 </td> <td> 2^-1 </td> <td> 2^-1a1=2^0a0 </td> <td> sqrt(2)=1.414214 </td> </tr>
  <tr> <td align="right"> 4 </td> <td> 2^-1 </td> <td> 2^-1a1(2^-0.5a0)*4 </td> <td> 4sqrt(2)=5.656854 </td> </tr>
  <tr> <td align="right"> 5 </td> <td> 2^-0.5 </td> <td> 2^-0.5(2^0r) </td> <td> sqrt(2)/2=0.7071068 </td> </tr>
  <tr> <td align="right"> 6 </td> <td> 2^-0.5 </td> <td> 2^-0.5d=2^-0.5*4 </td> <td> sqrt(2)=1.414214 </td> </tr>
  <tr> <td align="right"> 7 </td> <td> 2^-0.5 </td> <td> 2^-0.5a1=2^0a0 </td> <td> 1 </td> </tr>
  <tr> <td align="right"> 8 </td> <td> 2^-0.5 </td> <td> 2^-1a1(2^0.5a0)*4 </td> <td> 4 </td> </tr>
  <tr> <td align="right"> 9 </td> <td> 2^0 </td> <td> 2^0(2^-0.5R,2^0.5r) </td> <td> 1 </td> </tr>
  <tr> <td align="right"> 10 </td> <td> 2^0 </td> <td> 2^0d=2^0*4 </td> <td> 2 </td> </tr>
  <tr> <td align="right"> 11 </td> <td> 2^0 </td> <td> 2^0a1=2^0.5a0 </td> <td> sqrt(8)=2sqrt(2)=2.828427 </td> </tr>
  <tr> <td align="right"> 12 </td> <td> 2^0 </td> <td> 2^0a1(2^0.5a0)*4 </td> <td> sqrt(32)=4sqrt(2)=5.656854 </td> </tr>
  <tr> <td align="right"> 13 </td> <td> 2^0.5 </td> <td> 2^0.5(2^-0R,2^1r) </td> <td> sqrt(2)=1.414214 </td> </tr>
  <tr> <td align="right"> 14 </td> <td> 2^0.5 </td> <td> 2^0.5d=2^0.5*4 </td> <td> 2sqrt(2)=2.828427 </td> </tr>
  <tr> <td align="right"> 15 </td> <td> 2^0.5 </td> <td> 2^0.5a1=2^1a0 </td> <td> 4 </td> </tr>
  <tr> <td align="right"> 16 </td> <td> 2^0.5 </td> <td> 2^0.5a1(2^0a0)*4 </td> <td> 8 </td> </tr>
  <tr> <td align="right"> 17 </td> <td> 2^1 </td> <td> 2^1(2^1-0.5R,2^0.5r) </td> <td> 2 </td> </tr>
   </table>

同様に外接/内接多角形の比率連鎖が同心円空間を生成するのは他には正三角形の場合だけとなります。

![image.gif](https://camo.qiitausercontent.com/70d583e3b26928d2c5deff8201340994f34e18db/68747470733a2f2f71696974612d696d6167652d73746f72652e73332e61702d6e6f727468656173742d312e616d617a6f6e6177732e636f6d2f302f3539373737322f39666662633532312d656339642d376331302d356162662d6635626533303738303236372e676966)

正三角形に内接する正三角形の辺長比は1:1(最大=同サイズ)から1:2(最小)にかけて推移するので、それが生み出す同心円集合の半径比は綺麗に$2^n$乗($n \in \mathbb{Z}$)で連鎖していくのですが、ややこしい事に各同心円に内接する正三角形の一辺はその円の半径の$\sqrt{3}$倍となるのです。

<!-- Sat May  9 16:53:13 2020 -->
<table border=1>
<tr> <th>  </th> <th> Target_size </th> <th> Target_names </th> <th> Target_values </th>  </tr>
  <tr> <td align="right"> 1 </td> <td> 3^-1 </td> <td> 3^-1(3^0r) </td> <td> 1/2=0.5 </td> </tr>
  <tr> <td align="right"> 2 </td> <td> 3^-1 </td> <td> 3^-1h=3^-1*3 </td> <td> 3/2=1.5 </td> </tr>
  <tr> <td align="right"> 3 </td> <td> 3^-1 </td> <td> 3^-1a1=3^0a0 </td> <td> 2*sqrt(1.5)=2.44949 </td> </tr>
  <tr> <td align="right"> 4 </td> <td> 3^-1 </td> <td> 3^-1a1(3^0a0)*3 </td> <td> 6*sqrt(1.5)=7.348469 </td> </tr>
  <tr> <td align="right"> 5 </td> <td> 3^0 </td> <td> 3^0(3^-1R,3^1r) </td> <td> 1.0 </td> </tr>
  <tr> <td align="right"> 6 </td> <td> 3^0 </td> <td> 3^0h=3^0*3 </td> <td> 3.0 </td> </tr>
  <tr> <td align="right"> 7 </td> <td> 3^0 </td> <td> 3^0a1=3^1a0 </td> <td> 2*sqrt(3)=3.464102 </td> </tr>
  <tr> <td align="right"> 8 </td> <td> 3^0 </td> <td> 3^0a1(3^1a0)*3 </td> <td> 6*sqrt(3)=10.3923 </td> </tr>
  <tr> <td align="right"> 9 </td> <td> 3^1 </td> <td> 3^1(3^0R) </td> <td> 3.0 </td> </tr>
   </table>

ところで$\sqrt{3}$といえば1辺が1の単位立方体の立方対角線$\sqrt{1^2+1^2+1^2}=\sqrt{3}$の計算結果でもあったりします。全方向に対する立方体の無限連鎖が構成する同心円集合においてはここで挙げた「**$2^n$の数理**」と「**3の数理**」が統合されるとも考えられるかもしれません。

```math
\cos(π)=-1\\
\sin(2π)=sin(π)=\cos(\frac{π}{2})=0\\
\cos(\frac{π}{4})=2\sqrt{2}(0.7071068)\\
\cos(\frac{π}{3})=\sin(\frac{π}{6})=\frac{1}{2}(0.5)\\
\cos(\frac{π}{6})=\sin(\frac{π}{3})=\frac{\sqrt{3}}{2}(0.8660254)\\
\cos(2π)=\sin(\frac{π}{2})=1
```

```R
> cos(2*pi)
[1] 1
> round(sin(2*pi),digit=8)
[1] 0
> cos(pi)
[1] -1
> round(sin(pi),digit=8)
[1] 0
> round(cos(pi/2),digits=8)
[1] 0
> sin(pi/2)
[1] 1
> cos(pi/3)
[1] 0.5
> sin(pi/3)
[1] 0.8660254
> sqrt(3)/2
[1] 0.8660254
> cos(pi/4)
[1] 0.7071068
> sqrt(2)/2
[1] 0.7071068
> cos(pi/6)
[1] 0.8660254
> sin(pi/3)
[1] 0.8660254
> sqrt(3)/2
[1] 0.8660254
> sin(pi/6)
[1] 0.5
```
[【オイラーの多面体定理と正多面体】とある「球面幾何学(Spherical Geometry)」の出発点…](https://qiita.com/ochimusha01/items/5f8dbf33bac5b7227a11)
![img.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%252F20201121%252F20201121223022.gif%3Fixlib%3Drb-1.2.2%26auto%3Dformat%26gif-q%3D60%26q%3D75%26w%3D1400%26fit%3Dmax%26s%3D3f3105fa287092825b39b75fa7a36f21?ixlib=rb-4.0.0&auto=format&gif-q=60&q=75&w=1400&fit=max&s=54e1688127dec878f58461ec6d247c25)
![img.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-4.0.0&auto=format&gif-q=60&q=75&w=1400&fit=max&s=5f959037ec8d18f2a24f65e77e7cce11)
![img.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fcdn-ak.f.st-hatena.com%2Fimages%2Ffotolife%2Fo%2Fochimusha01%2F20201129%2F20201129055848.gif?ixlib=rb-4.0.0&auto=format&gif-q=60&q=75&w=1400&fit=max&s=61e634a0808fc8df078f8c1ff744232c)

* 半径$\frac{\sqrt{3}}{2}$(直径$\sqrt{3}$)の球面に内接する立方体は辺長1、平方対角線長$\sqrt{2}$,立方対角線長=$\sqrt{3}$となる。さらにこれに内接する球面は半径$\frac{1}{2}$(直径1)となる。
* これに内接する正四面体の辺長は$\sqrt{2}$となる。
* これに内接する正八面体の辺長は$\frac{\sqrt{2}}{2}$となり半径$\frac{1}{2}$(直径1)の球面に内接する。

この辺りについてはまだまだ全然調べ込みが足りていません。何しろこの段階ではまだいわゆる「**プラトン立体**」５個のうち3個が登場するに過ぎないのです。
[【オイラーの多面体定理と正多面体】とある「球面幾何学(Spherical Geometry)」の出発点…](https://qiita.com/ochimusha01/items/5f8dbf33bac5b7227a11)
![image,gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fcdn-ak.f.st-hatena.com%2Fimages%2Ffotolife%2Fo%2Fochimusha01%2F20210125%2F20210125160959.gif?ixlib=rb-4.0.0&auto=format&gif-q=60&q=75&w=1400&fit=max&s=63f3be7d9ee3e2411b0d1d2a2cbcf9c3)
![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%252F20210125%252F20210125162938.gif%3Fixlib%3Drb-4.0.0%26auto%3Dformat%26gif-q%3D60%26q%3D75%26s%3Dc537d6648d806c18a199c8628fff3ea3?ixlib=rb-4.0.0&auto=format&gif-q=60&q=75&w=1400&fit=max&s=11943feee24090ad422d1d638fc49e1d)
[【オイラーの多面体定理と正多面体】正六面体と正十二面体の間の往復(「屋根掛け法」とは何か？)](https://qiita.com/ochimusha01/items/37e73f243a0b2f2cba9d)
![image.png](https://camo.qiitausercontent.com/a6f177167a6099d331ad6b12e4b886200357da08/68747470733a2f2f71696974612d696d6167652d73746f72652e73332e61702d6e6f727468656173742d312e616d617a6f6e6177732e636f6d2f302f3539373737322f62626165623038612d333161372d303866302d663861332d3466306239666366653461332e706e67)
![image.gif](https://cdn-ak.f.st-hatena.com/images/fotolife/o/ochimusha01/20210613/20210613193801.gif)

#二項展開(<font color="blue">Binomial Expansion</font>)の限界
ところで**平方ピタゴラスの定理**<font color="green">$a^2=b^2+c^2$</font>は一般の**2項展開**(<font color="blue">Binomial Expansion</font>)<font color="green">(a±b)^2</font>では扱い切れません。
[【初心者向け】パスカルの三角形と二項定理または二項展開](https://qiita.com/ochimusha01/items/d668b19d5389b52e0dd2)

その証明には「**角度が**<font color="green">90度</font>(<font color="blue">π/2ラジアン</font>)**という特別条件**」あるいは「**複素平面**(<font color="blue">Complex Plane</font>)**における単位円上の操作、すなわち虚数**<font color="green">i^2=-1</font>**の概念の導入と共益複素数**(<font color="blue">Conjugate Complex Number</font>)**の相互打ち消し作用**」を必要とするのです。

```R
#余弦定理(cosine formula)で角度が90度の場合
a^2=(b-c)^2=b^2+c^2−2bc*cos(α)
#ところでα=90度(π/2ラジアン)の場合、cos(α)=0
a^2=b^2+c^2−2bc*0=b^2+c^2 

#オイラーの公式(Euler's Formula)あるいは
#共益複素数(Conjugate Complex Number)同士の打ち消し作用を用いた証明。
1=e^0=e^(iθ-iθ)=e^iθ-e^iθ)
=(cos(θ)+sin(θ)i)(cos(θ)-sin(θ)i)
=cos(θ)^2+sin(θ)i-sin(θ)i-sin(θ)i^2
=cos(θ)^2+sin(θ)^2 #i^2=-1
```
そういえば受験数学では不思議とこういう肝心の事を習いません。**自明の場合**(<font color="blue">Trivial Case</font>)だから？

#三角比概念の限界

ただし**三角比**(<font color="blue">Triangle Ratio</font>)すなわち**直角三角形**（<font color="blue">Right triangle</font>）の概念に拘っている限り、X座標もY座標も<font color="green">0</font>以上、すなわち面積が<font color="green">0</font>にならない範囲にしか想像が及ぶ事がなかったのでした。
![CFS01.gif](https://cdn-ak.f.st-hatena.com/images/fotolife/o/ochimusha01/20210129/20210129040756.gif)

```R
#CFS=Circle Function Single
CFS<-function(Radian){
cx<-seq(-1,1,length=60)
f0<-function(x) sqrt(1-x^2)
cy<-f0(cx)
plot(cx,cy,asp=1,type="l",xlim=c(0,1),ylim=c(0,1),main="Cos(θ)(ad)^2+Sin(θ)(bd)^2=1",xlab="X(1→0)",ylab="sqrt(1-X^2)(0→1)")

#塗りつぶし（背景）
polygon(cx, #x
cy, #y
density=c(30), #塗りつぶす濃度
angle=c(45),　 　　　#塗りつぶす斜線の角度
col=c(200,200,200)) 　#塗りつぶす色

#塗りつぶし（三角）
polygon(c(0,cx[Radian],1), #x
c(0,cy[Radian],0), #y
density=c(30), #塗りつぶす濃度
angle=c(45),　 　　　#塗りつぶす斜線の角度
col=rgb(0,0,1)) 　#塗りつぶす色
#線追加（三角）
segments(cx[Radian],cy[Radian],cx[Radian],0,col=c(100,100,100))
#文字追加
text(c(0,cx[Radian],1,cx[Radian]),c(0,cy[Radian],0,0), labels=c("a","b","c","d"),col=c(rgb(0,0,0),rgb(0,0,0),rgb(0,0,0),rgb(0,0,0)),cex=c(2,2,2,2))
}
#アニメーション
library("animation")
Time_Code=rev(seq(30,60, length=30))
saveGIF({
for (i in Time_Code){
　CFS(i)
}
}, interval = 0.1, movie.name = "CFS01.gif")
```
xの値が<font color="green">1</font>から<font color="green">0</font>にかけて遷移する**関数**<font color="green">$y=sqrt(1-x^2)$</font>の内容は、**三角関数**(<font color="blue">Trigonometric Function</font>)的に**表現**(<font color="blue">Expression</font>)すると、**角度**<font color="green">θ</font>の<font color="green">0ラジアン</font>(<font color="blue">0度</font>)から<font color="green">π/2ラジアン</font>(<font color="blue">90度</font>)に向かう変遷を扱った形となります。**頂点**<font color="green">b</font>および**頂点**<font color="green">d</font>のx値が<font color="green">Cos(θ)</font>,**頂点**<font color="green">b</font>のy値が<font color="green">Sin(θ)</font>に対応。
[【Rで球面幾何学】「半円しか描けなかった」世界の思い出？](https://qiita.com/ochimusha01/items/cf74b1fff650edb0a496)
以下続報…