Fallen Physicist, Rising Engineer まとめリンク：

<ドラマに出てきた数式シリーズ>

新ガリレオ（ドラマ）に出てきた数式まとめ　+　α

物理学者が変身する仮面ライダービルドに出てきた数式解説まとめ

＜カシオの高精度計算サイトの自作式＞

カシオの高精度計算サイトkeisan.casio.jpに投稿した自作式まとめ：

　　ランキング2位：「円の弧長,弦長,矢高,半径のどれか2つを与えて残りを計算」をkeisan.casio.jpにUP!

　　ランキング1位：運命数の計算（やりすぎコージー都市伝説より）をカシオの高精度計算サイト、keisan.casio.jpの自作式として作ってみたよ。

＜Excel VBAで数値計算＞

Excelマクロ（VBA)で数値計算-複素TDGL方程式、蔵本シバシンスキー方程式、ピタゴラスの三体問題、シュレーディンガー方程式、FPU問題などなど

Dormand-Prince(ルンゲ・クッタ8次）の有名なルーチンDOP853をExcel VBAに移植

Excel VBAで複素数/FFTが使えるライブラリ

＜超高次ルンゲクッタ＞

35段14次のルンゲクッタ法をPARI/GPに実装、ローレンツ方程式を計算し、通常の4次、そして8次、オイラー法と比較してみる。

＜Scratch＞

Scratchまとめ：Numerical Calculation

＜もろもろ所感＞

電子レンジは、水分子の固有振動数（共振周波数）を利用しているのではないです。

世界で最も恐ろしい試験は「マクローリン展開の剰余項を求めよ」→GeoGebraでマクローリン展開してみる。

Powerpointで数式入力をLaTeX形式で行う。WordはデフォルトでできるけどPowerpointはちょっと初期設定が必要。」

京都ｖｓ滋賀の決め文句、琵琶湖の水止めてやる！は、、、京都の疏水事務所で管理してるから滋賀県民には止められない

長谷工コーポレーションのCM、一体何に”気づいた？”なのかというと、、、

41個以上の正方形を見つけられるか？？？　-　わたしも思いついた！　42個目も思いついたぞ（Squareは正方形だけじゃない）。

出雲大社大阪分祠へ行ってきた。やはり大きなしめ縄があった。

南海の初芝駅から歩いてここまで。

分祠とは言ってもとても立派です。

もちろん大きなしめ縄もある。

南海初芝駅付近の前ヶ池に神々しく翼を広げる鳥が石の上に止まっていた。Google画像検索するとミミヒメウと出た。

最初、オブジェかな？と思っていたら動いてびっくりした。こんなかっこよく止まるのはすごいな。

Google画像検索するとミミヒメウとでた。

ここには本当のオブジェもある。パンダが怖い、、、

南海の萩原天神駅（はぎはらてんじん）で降りて萩原天神（萩原神社）でお参り。

久しぶりに南海電車に乗った。で降りたのは萩原天神駅。はぎわらでもおぎわらでもなくてはぎはらなんですね。

まだ桜が残っていた。

さてここから（続く）

高校事変Xを読んだ。もう魁男塾とかドラゴンボールとかそんな感じになってるな。。。そして次巻ではおそらく松岡作品の全ヒロインが集合すると思う。

前に原爆まで出てもう「強い奴のインフレ」（by　さるマン）に行きついたこのシリーズですが、ラスボスを出したのと、敵を傭兵にして舞台を海外に移して結衣が勝てない感を出してリセットされてます。が、2km走るだけで対応してたりするので、ドラゴンボールの精神と時の部屋に入ったようだ、、、それと、本作で「探偵の探偵」との繋がりがはっきりしました。

あらすじは

「慧修学院高校3年の生徒たちが国際交流のためにホンジュラスを訪問。そのさなか、メキシコの過激派組織ゼッディウムに襲撃される。背後には結衣の兄、架禕斗（この漢字出せないよ、、、）の姿があった。そのころ日本では、凛香の母親の市川凛（これが探偵の探偵で出てきた悪人）が怪しい動きを、、、」

というもの。ラスボスが兄弟姉妹というのはたまにありますが、確かに行きついた結衣の強さに対抗しようと思うと兄弟姉妹を出すしかないという納得の展開。そして探偵の探偵と結びついたことで、これも探偵の鑑定のような感じで次のお話、松岡作品の全ヒロイン（万能鑑定士Q,水鏡推理などなど）でそうな雰囲気。この本も最後に過去作品の最新作の宣伝も出てます。千里眼の復活も。

さすがを「小説家になって億を稼ごう」を書いた松岡さんで売り上げが上がる方法を知ってますな。

松屋でさばの味噌煮御前（定食にとろろ追加と豚汁変更したもの）、ご飯特盛でいただく。生姜がとてもきいた味噌煮でこれはかなり美味しい。

松屋で鯖味噌？とちょっと訝しく思ったが珍しいので注文。

もろみ味噌がかかっているのもありますが、それよりしょうがの利き方が非常によく、おいしかった。

さらに松屋さん豚汁はいつも熱々で熱々でいいです。

とろろとの相性も良く、ご飯特盛で正解でした。最後はご飯に余った味噌を乗せて全部いただいた。

円と接線を使った1枚の絵でsin/cos/tan/sec/csc/cotを表すもの、複素数を使って接線のパラメータ表示(z=exp(iθ)+i exp(iθ) t )で計算してみる。

昨日、円と接線を使った1枚の絵でsin/cos/tan/sec/csc/cotを表すのを見た。

All the basic trigonometric identities in one figure pic.twitter.com/5tFml7VFZY

— Fermat's Library (@fermatslibrary) April 3, 2021

Geogebraでも描いてみた。

もちろん普通に相似とか利用して計算すればいいですが、今回は複素数を使った計算をやってみよう。和、積の公式とか全部複素数にすると簡単に求まるのでその延長で。

まず円上の点は

z=exp(iθ)とかけて、

接線は、傾きはexp(i(θ+π/2))=i exp(iθ)になるので（大きさは1になっている）、

円上の点を通ってこの傾きの直線をパラメータ表示（tを使う）すると、

z=exp(iθ)+i exp(iθ) t = (cosθ - t sinθ) + i (sinθ + t cosθ)

となる。

y=0になる点はt=-sinθ/cosθ = - tanθ

なので、接点と、接線を伸ばしてx軸とぶつかる間の距離はtanθ。

原点から接線を伸ばしてx軸とぶつかる点の間の距離は

cosθ - (-tanθ)*sinθ = cosθ + (sinθ)^2/cosθ =((cosθ)^2+ (sinθ)^2)/cosθ = 1/cosθ = secθ

x=0になる点はt=cosθ/sinθ = cotθ

なので、接点と、接線を伸ばしてy軸とぶつかる間の距離はcotθ。

原点から接線を伸ばしてy軸とぶつかる点の間の距離は

sinθ + (cotθ)*cosθ = sinθ + (cosθ)^2/sinθ =((sinθ)^2+ (cosθ)^2)/sinθ = 1/sinθ = cscθ

とかなり簡単に求められる。

新型コロナウイルス、中国、日本、韓国、アメリカ、ドイツ、フランス、イギリスでの感染者数を指数関数＆ロジスティック関数＆Log-Logプロットでべき関数フィッティングした（4/4更新） フランス・ドイツの増え方がエグイが日本・韓国もやはり増加中。イギリスの減速がすごいな。

まずはリニアスケール。イギリスが目に見えて減速しているのは続いている。しかしフランス、ドイツがひどい。

日本も報道されているように増えている。韓国も実は何げに全然減速してない。

各国のログログプロット。なんやかんやで各国の傾きこのスケールだと似てきてる。

日本の詳細ログログプロット。このままいくと第四波がすぐやってくるな。まん延防止等重点措置がどこまで効果がでるのか。

長岡天神（八条ヶ池）の夜桜ライトアップをiPhone 12 Pro Maxで撮影(2021/4/3)

桜は散り始めていますが、十分綺麗なライトアップ。

さすがはiPhone 12 Pro Maxのカメラ。何の設定もしなくても普通に綺麗に撮影できた（リサイズしてますが）。

ココイチで牛カレー（ぎゅうかれー）＋やさい、400gの２辛を食す。かなりコクがあって美味しい。

一見していつものポークカレーと違う、濃い色で味も予想通りコクがあって美味しかった。

肉も結構はいってる。やさいを入れてさらに美味しく（やさいいれないとちょっと具が少ないかも）。

分数階微分(Fractional Calculus)の計算をPythonのdifferint（Riemann-Liouvilleの定義で計算するもの）でやってみる。

たまたま非整数階微分のことを調べていて、

整数階でない微積分法について

を見ていた。昔、非整数階の拡散方程式（∂u/∂t　=　D∂^α u / ∂u^α ）を計算したことがありますが、

その時とは違うRiemann-Liouvilleの定義で計算しようと思った。

こういうのは絶対Pythonでやっている人がいるだろう、と思ったら、やっぱりあった。こちら。differint。

https://pypi.org/project/differint/

では早速やってみる。

f(x)=x

f(x)=sin(x)

f(x)=exp(-x^2)

なるほど直感的にもそんな感じだ。なぜこれを調べていたかというと（続く）。

京都の猿ラインへ行く（その３） 猿ヶ辻（京都御所） 何度も見ていた光景なのに、なんでこんな形か全然疑問に思ったことがなかった。

最後は猿ヶ辻。御所には何度も言っているのに、そしてこれよく考えるとおかしい形なのに疑問にも思わず、理由を知らなかった。
こういう理由です。

https://ja.kyoto.travel/glossary/single.php?glossary_id=387

そう知ってみると見方が違うな。

猿ライン３か所回ってかなり歩き回って疲れたがすっきりした。まあこれが御利益と言えばそうか。

滋賀にもいつかいきたい。

京都の猿ラインへ行く（その２）幸神社（さいのかみやしろ） 出町柳の近く、、、大学から近かったのに全然しらなかった。

赤山禅院の次は幸神社へ。さいのかみやしろと読むそうです。

大学から近かったのに全然知らなかった。Googleマップを頼りに歩きますが、見落としそうな位置にありました。かなりこじんまりした神社です。

石神さんがいた。

続く。

ピタゴラミングスイッチ（ピタゴラスイッチ＋プログラミング）を録画して見た。笹の葉のミシン穴に驚く、、、

『ピタゴラミングスイッチ』3月30日朝9時15分放送　プログラミング的思考をはぐくむ

というのを知って録画してさっき見ました。

百科おじさんが登場してプログラミングは

・じゅんびがだいじ

・よそうがだいじ

・じゅんばんがだいじ

だという。

まずはじゅんびがだいじ、ということでもんくたれぞうが登場。

ズボンと長靴が脱ぎっぱなし、、、ではなくて消防士さんがすぐに着替えられるように。

http://www.chunou-119.jp/diary/3084

なるほどー、納得。

でこれを応用して、白い枠の中にピタゴラ装置を準備するプログラミング装置その１が紹介され（文字で説明しにくい、、、）

次はぼてじん。

コンベアに乗ってもらう。

直角に3つあって、継ぎ目のところでころりと回転する。ゴールで正面を向くにはどの向きで乗ったらいい？

って絵がないとわからないのでパワーポイントでかいてみた。

下を向いた形では、、、残念。

逆から考えるとわかる。最後に正面を向いていたいなら、、、とか考える。

左を向いてかつ下を向くと。正面を向いた。

次はプログラミング装置その２，

４３２１２１装置。サイコロに書かれた数字がこの順番で並ぶ。

１＝ら

２＝か

３＝ど

４＝の

にすると？また

１＝え

２＝め

３＝ぎ

４＝や

だったら？というもの。

という問題でした。

次はフローチャート。

左側のは何？

→猿の虱取り。

右側は？

→馬飛び。

次は想像してみよう。

びっくりする笹の葉のミシン穴の理由。。。

こういうことらしいです。リンク参照。

https://sambuca.jp/bamboo-leaf/

次は２枚の扉でのキリンの通路（言葉で説明しにくい、、、）を挟んで、、、

順番が大事、ということで鉄琴階段。こういうやつで、入れ替えると違う曲になる。

https://www.momotoys.jp/2019/12/02/creativeselection-02/

そして蜘蛛の糸のすごい張り方でおしまい。すごいなこれ何度見ても。

https://www.nacsj.or.jp/archive/2012/08/1180/

京都の猿ラインへ行く（その１）赤山禅院で猿を屋根のてっぺんに見る。呪術廻戦の禪院真希先輩の名前も入ってる。

先日、鳥居の密室　世界にただひとりのサンタクロース（島田荘司さん）を読んだ。

そこに猿ラインというのが出ていた。なんでも平安京を守るために猿を使うが、滋賀県の日吉大社→京都の赤山禅院→幸神社→猿ヶ辻が一直線上に並んでいるところが防衛ラインだというような話。

面白そうなので行ってみよう。滋賀は遠いのでちょっと後回しでまずは京都、赤山禅院だ。ちょうど呪術廻戦の禪院真希先輩の名前も入ってるし。

あれ？猿はどこだ？と思ったら、、、屋根の上でした！

続く。

iPhoneで使える数学ソフトMaple Calculatorで手書き入力でMathematicaの常微分方程式の例を試す（その１）

MapleというとMaxima、Mathematica、reduceなどと同じくらい有名な数式処理ソフトウェアだと思いますが、いつの間にかスマホで（一部機能を）無償で使えるようになっていた。

https://jp.maplesoft.com/products/maplecalculator/index.aspx?L=J

どの程度使える？というのを見るために、まずは常微分方程式を解かせてみよう。

例題はWolfram言語（Mathematica）のやつを使う。

https://reference.wolfram.com/language/tutorial/DSolveOrdinaryDifferentialEquations.html.ja?source=footer

大量にあるので、まずは順番に上から4つ。

全部手書きで書いて、カメラで撮って解かせた。認識度も高い。

結果はこちら：

おお、ちゃんと解いてくれてるじゃないの！

ぼちぼち続きをやっていこう。