« 2020年12月 | トップページ | 2021年2月 »

2021年1月

2021年1月31日 (日)

特性インピーダンスが高くて線路長さが短い伝送線路はインダクタンスに見える(プローブのキャリプレーションキットのロードの設定)

フォームファクター社(カスケードマイクロテックを買収)の高周波プローブを使うとき、キーサイトのネットアナが伝統的に校正キットのロードにインダクタンス成分が入れられなくしているので、代わりにディレイの項目で、Delay=Lterm/500, Zo=500にしてくださいと書いてある。・・・が若手とか、がなぜかわからんとか、何も考えずに500Ωを50Ωにしていたりする。

今年に入ってからも4回この話を聞いた、、、ので説明してみる。

※それもあってカシオの高精度計算サイトにこれをUPしていた。

「伝送線路のついた負荷の入力インピーダンス」をカシオの高精度計算サイトkeisan.casio.jpに自作式としてUP!

Zl_20210130181601

ではやってみる。

伝送線路の特性インピーダンスをZo, 位相定数をβ(今はロスなしとするのでγ=j*β), 線路内の速度をv=c/√εeff、線路長さをL、負荷インピーダンスをZLとする。

ロスなしの場合の入力インピーダンスはこうかける。

               ZL+j*Zo*tan(β*L)
Zin=Zo* ───────────
       Zo+j*ZL*tan(β*L)

で、通常はZoが50Ω。

でもここでZoが負荷インピーダンスよりずっと大きい(|Zo|>>|ZL|)とする。(プローブの場合は500Ωに設定で、負荷インピーダンスの方が50Ω)

すると入力インピーダンスはZL/Zoの一次のオーダーで

               ZL/Zo+j*tan(β*L)
Zin=Zo* ─────────────
       1+j*(ZL/Zo)*tan(β*L)

    ≒ZL+j*Zo*tan(β*L)

と近似できる。さらにLが非常に短いときはtan(x)≒x の近似を使って

Zin≒ZL+j*Zo*β*L=ZL+j*Zo*(ω/v)*L

ここで、L/vというのが線路内を進む遅延時間(ディレイ)になるのでこれをτと置くと

Zin≒ZL+j*ω*Zo*τ

これとインダクタンスL(線路のLと重なるが我慢して)がついた負荷のインピーダンス

Zin≒ZL+j*ω*L

と比べると

L=Zo*τ

となっている。つまり、τ=L/Zoをオフセットディレイに入れて、Zoを特性インピーダンスにするとインダクタンスが再現できるということになる。

あくまで線路長が短くて、特性インピーダンスが負荷インピーダンスよりずっと大きいときのみに成り立つことに注意。

 

2021年1月30日 (土)

「積層セラミックコンデンサ(MLCC)の加速試験による寿命予測」をカシオの高精度計算サイトkeisan.casio.jpの自作式としてUP!

久しぶりに最先端のMLCCの話を聞く機会があった。それで何か自分でも計算してみようということで、積層セラミックコンデンサ各社が出している寿命予測の経験式をやってみた。

有名どころ三社のリンク:

https://product.tdk.com/info/ja/contact/faq/faq_detail_D/1432655794672.html

https://www.murata.com/ja-jp/support/faqs/products/capacitor/ceramiccapacitor/qlty/0010

https://www.yuden.co.jp/jp/product/support/faq/q020.html

これらを参考に作ったのがこれ:

 積層セラミックコンデンサ(MLCC)の加速試験による寿命予測

 

Mlcc1

積層セラミックコンデンサ(MLCC)の市場での寿命予測に使われる経験式です。
Ln:通常の使用条件での寿命(時間)
La:加速条件での寿命(時間)
Vn:通常の使用電圧(V)
Va:加速条件での電圧(V)
Tn:通常の使用温度(℃)
Ta:加速条件での温度(℃)
n:電圧加速係数
θ: 温度加速係数

Mlcc2

2021年1月29日 (金)

「伝送線路のついた負荷の入力インピーダンス」をカシオの高精度計算サイトkeisan.casio.jpに自作式としてUP!

伝送線路のついた負荷の入力インピーダンスの計算、回路シミュレータでやればいいのだがちょっとしたものを計算するのにお気軽なものはないかな、と思って自分で作ってみた。

こんな感じの状況で、

Zl

 

負荷インピーダンスZL(ωの関数で入力)に対し、特性インピーダンスZo, 長さL、有効比誘電率εeff, 減衰定数αの伝送線路がついたとき、どういう入力インピーダンスになるかを計算します。

というもの。ZLとαがωの関数として入力できるのはあんまりないんじゃないか。

リンクはこちら:

伝送線路のついた負荷の入力インピーダンス

画面とグラフ表示はこんな感じ。

Zl3

Zl2

なか卯でとろろ付き漬けまぐろ鉄火丼(普通盛)+サラダ+味噌汁をいただく。

豪快盛もあるけど、ちょっとそこまでは、、、ということで普通盛にして代わりにとろろを追加。

20210111-114321 20210111-114322

やはり鉄火丼ととろろ(+わさび)はとても合うので美味しい。

20210111-114445

2021年1月28日 (木)

10000000000000000÷99999999999999989999999999999999でフィボナッチ数列76項までわかる。PARI-GPで2000桁で見てみる。

というツイート見た。
おお、これはフィボナッチ数列を生み出す母関数 f(x)=x/(1-x-x^2)のx=1/100の場合ですね。

ではもっと進めて、x=1/10000000000000000の場合はどうなる?というと、、、(普通には計算できないのでPARI-GPを使って桁数を2000桁まで上げる)

Fibo20210127

これで76項まで計算できてる!

元データはこちら:

(22:38) gp > \p 2000
realprecision = 2003 significant digits (2000 digits displayed)
(22:41) gp > 10000000000000000/99999999999999989999999999999999.
%14 = 1.0000000000000001000000000000000200000000000000030000000000000005000000000000000800000000000000130000000000000021000000000000003400000000000000550000000000000089000000000000014400000000000002330000000000000377000000000000061000000000000009870000000000001597000000000000258400000000000041810000000000006765000000000001094600000000000177110000000000028657000000000004636800000000000750250000000000121393000000000019641800000000003178110000000000514229000000000083204000000000013462690000000002178309000000000352457800000000057028870000000009227465000000001493035200000000241578170000000039088169000000006324598600000001023341550000000165580141000000026791429600000004334944370000000701408733000000113490317000000018363119030000002971215073000000480752697600000077787420490000012586269025000002036501107400000329512800990000053316291173000008626757127200001395838624450000225851433717000036543529616200005912867298790000956722026041000154800875592000025047307819610004052739537881000655747031984200106102098577230017167680177565002777789003528800449455702128530072723460248141011766903046099401903924907091350308061521170129049845401187926408065155330493931304969544928657211148507797805034164546229067075527939700884757894439432379146544723340246762233416728348467688788906237314391213057907216116009194853094755513050064381636711396954969111226270196140727489740989163763861236800877783661021089979416004714477006719437081658600466103755310630113804746347649016041512187871202742198682263610434634990105073070885485833143411434898484365071852344706767941299583455520444848481792619723897844013817176838269219307914922805362068963260663228399975475294376460687180136069930068472766550757613719078015775062056635467085082342854326866258854851787357476708913722004410259439890074015793033126227445681897711523484726120102414622929430987356697140204299759815943314739849548565733516982553016006499096750787257985079500608885863498917568761166200686762964975255057851984109187512652828060671301843802647159 E-16

ココイチでベジタブルスープカレー(3辛)を食べた。豆腐とコーンのナゲットが予想以上に美味しかった。

ローストチキンスープカレーは以前も食べたし、ここは野菜のみのベジタブルスープカレーにした。
じゃがいも、キャベツ、人参、オクラ、なす、ごぼう、かぼちゃ、れんこん、カリフラワーのたっぷりな野菜とコーン入りの豆腐ナゲット

が入っているとのこと。特に、オクラ、レンコンが美味しいのと予想以上にコーン入りの豆腐ナゲットが美味しかった。

スパイスもついてくるのでかなり多めにかけていただいた。

20210104-123708 20210104-123710 20210104-123711 20210104-123712

2021年1月27日 (水)

「レーダ最大検知距離(レーダ方程式より)」をカシオの高精度計算サイトkeisan.casio.jpの自作式としてUP!

ちょっと計算したいことがあって、一番簡単なレーダー方程式

Pr=Pt*G^2*λ^2*σ/[(4π)^3*R^4]

を使って最大検知距離を計算する自作式を作った。

これです。リンクをクリック!

 レーダ最大検知距離(レーダ方程式より)

画面はこんな感じで。いろいろな単位を選べるようにはしています。

ちなみに、入っているデフォルトの値はあえてデタラメなものにしています。いろいろやばいからね。民生でも。

Radar_eq Radar_eq2

2021年1月26日 (火)

松のやで「シュクメルリチーズBigメンチハンバーグ定食」を食す。ものすごくニンニクが効いていて、かつハンバーグがでかい!

松屋ではシュクメルリ鍋定食が始まってますが、こちらの松のやでもシュクメルリチーズBigメンチハンバーグ定食が始まった。

メンチカツ(松のやさんではハンバーグとつけるこだわりが、、、)が巨大で、かつチーズとニンニクが効いていてこれはボリュームがあっていいな。このお店はご飯・みそ汁おかわり無料な店なので1回ご飯をおかわり。ものすごくお腹いっぱいになった。

もちろん最後はご飯を鍋にINで。

20210121-193135 20210121-193137 20210121-193136

 

 

2021年1月25日 (月)

三宮神社(京都西京区。桂駅近く)でお参り。ムクノキがあった。

阪急の桂駅から歩いてぶらぶらと。

20210117-134101 20210117-134113 20210117-134121 20210117-134140 20210117-134144

御神木としてムクノキがあった。

20210117-134149 20210117-134152

2021年1月24日 (日)

#NHKスペシャル 「わたしたちの“目”が危ない 超近視時代サバイバル」を見てTwitterでリアルタイムメモしてました。そのリンク。

私も目がどんどん悪くなっているので(コロナでテレワークが増えてからさらにひどく)、興味深く見てました。

眼軸という言葉は初めて知った。で結論は

・合ったメガネに変えろ

・光を浴びなさい

・3つの20(近業をする(30㎝以下の作業)場合は20分に一回20秒、20フィート(6m)離れたところを見る)

でした。

 

新型コロナウイルス、中国、日本、韓国、アメリカ、ドイツ、フランス、イギリスでの感染者数を指数関数&ロジスティック関数&Log-Logプロットでべき関数フィッティングした(1/24更新)日本の増加率、アメリカ超えて韓国に近づいている。今回から各波ごとのべきを描いてみた。

緊急事態宣言が出ても大きくは減少してない。アメリカはもともと多いからあれですが、増加率だけ見るとアメリカ以上になってる。イギリスやドイツも超えてる。中国も、、、さすがにごまかしきれないのかちょっとだけ増加させてる。

リニアスケール:

Coronaww202101241

Log-Logプロット:

Coronaww202101242

今回から、日本の各時期でのべき乗フィッティングしてみた。

現在の第三波の増え方は第二波を超えて、第一波の後期くらいの急激な増加になってる。

Coronaww202101243

2021年1月23日 (土)

「ヘモグロビンの吸光スペクトル(パルスオキシメーターの原理」をカシオの高精度計算サイトkeisan.casio.jpに自作式としてUP!

新型コロナウイルスの影響で、一般の人までパルスオキシメーター(指に挟むもの)を入手して、実際に必要な人に届かないということが起きていると聞いた。Apple Watch Series6とかにも機能がついているのは知っていたが、原理はむかーし聞いたことがあるけど完全に忘れていた。

コニカミノルタさんのサイトが参考になる。

パルスオキシメーター知恵袋 基礎編

なるほど酸素を持っていたり離したりするヘモグロビンの吸光スペクトルを、波長の違う2つのLED&センサで受けて測定しているということなのか。

これをカシオの高精度計算サイトに作ってみよう。

元データは、こちら:

Optical Absorption of Hemoglobin

ただ関数形が複雑、、、なので、区分的に分けてそれぞれで6次関数でフィッティングすることにした。

できた自作式のリンクがこちら:

 ヘモグロビンの吸光スペクトル(パルスオキシメーターの原理)

画面はこんな感じで、

0

こんな感じで250nmから1000nmまで図示できる。

1

実際に使われている範囲がこちら。

2

うーん、区分的にフィッティングしているので、つなぎ目でずれてるなあ、、、

まあ実用のためというよりは雰囲気をつかむためのものなので勘弁してもらうとするか、、、

2021年1月22日 (金)

三輪神社(大阪の高槻市富田)でお参り。

やはり奈良の三輪神社と縁がある神社でした。この辺は事故が多い(富田)ので、その犠牲者へのご冥福の祈りも込めてお参り。

20201107-1442232 20201107-144223 20201107-144140 20201107-144102 20201107-144032 20201107-144027

2021年1月21日 (木)

向日市の体育館近くの巨大鬼瓦をiPhone 12 Pro Maxで撮影。

ここ、大昔(12年前)に記事にしていたのを思い出した。

その時に思ったのが仮面ライダーアマゾンの十面鬼?でした。今も同じことを思った。ポケモンGOのスキャン対象になっていたのでスキャンしました。

20201220-135021

2021年1月20日 (水)

1/20の又吉直樹のヘウレーカ 今日のテーマは「友達の友達はみな友達って本当?」をTwitterでリアルタイムメモしてました。東京大学生産技術研究所教授の羽田野直道さんがゲスト。そのリンク。

以下のリンクからたどれます。

 

関連:

Zachary's karate club

「友達の友達はみな友達だ-複雑ネットワークの科学」

<講義ノート>複雑ネットワーク:統計物理学の視点

Communicability and Communities in Complex Socio-Economic Networks

Social_network_model_of_relationships_in Images

木嶋神社(蚕ノ社)でお参り2021。徒歩で。

西院春日神社⇒西院野々宮神社⇒山王神社⇒猿田彦神社の次(最後)はこちら。

標識と鳥居がだいぶ前にある。

20210111-125955 20210111-130000

ちょっと歩くとまた鳥居。

20210111-130027 20210111-130334 20210111-130336 20210111-130416

ここでお参りです。

20210111-130447 20210111-130535

 

20210111-130538 20210111-130710

2021年1月19日 (火)

”ストライドとピッチ、タイム(2021年共通テスト数学IA第二問)”をカシオの高精度計算サイトkeisan.casio.jpに自作式をUP!

今年の共通テストで面白かったのが、

https://www.toshin.com/kyotsutest/suugaku-1a_question_2.html

2021

ピッチとストライド、タイムを計算させる問題。イラストも面白い。

で、カシオの高精度計算サイトにこれを解いてくれる自作式をUPした。

100mには限らないのでマラソンでも、、、

これです。

ストライドとピッチ、タイム(2021年共通テスト数学) 

こんな感じで計算してくれる。

2021_20210119224901

西院野々宮神社⇒山王神社⇒猿田彦神社(京都)でお参り2021。徒歩で。

西院春日神社の後はこの流れで。

まずは野々宮神社。

 

20210111-121800 20210111-121742 20210111-121835 20210111-121836

 

そして山王神社。

20210111-123909 20210111-123915

夫婦岩と、

20210111-124019 20210111-123929

大きな楠がある。

20210111-12410120210111-124103

そして猿田彦神社。伊勢神宮の近くにある猿田彦神社には毎年参っているのでこちらも。

20210111-125044 20210111-125224

ここにも区民の誇りの木(庚申楠)がある。

20210111-125207 20210111-125258 20210111-125229 20210111-125259

猿がかわいい。

20210111-125317

そしてまだ続く、、、

2021年1月18日 (月)

高周波(RF・マイクロ波・ミリ波・5G)関連ニュース2021年1月18日 IEEE Microwave Magazineの特集はパッシブ部品。負の群遅延を持つ線路など。Microwave Journalの特集はレーダとアンテナ。MIMOレーダとか。Cバンドのオークションは$80.9B,5G americasの3GPP Releases 16, 17 & Beyond

今月号のMicrowave Magazineの特集はPassive Component Technologies and Applications。

Negative Group Delay Circuits and Applications

Waveguide Components Based on Multiple-Mode Resonators

Tunable Balanced Power Dividers

https://ieeexplore.ieee.org/xpl/mostRecentIssue.jsp?punumber=6668

帯域は狭いけれど負の群遅延を持つ線路とか。

Microwavemagazine_202101

Microwave Journalの特集はRadar and antennas。

https://www.microwavejournal.com/publications/1

MIMOレーダなど。

Microwavejournal_202101

CoilcraftのRFインダクタの資料。

https://www.microwavejournal.com/articles/35291-key-parameters-for-selecting-rf-inductors

 

Cバンドのアメリカでのオークションが$80.9B。。。って一兆円くらい?すごいな。
https://www.eetimes.com/c-band-auction-brings-home-80-9b/

 

5Gアメリカズの3GPP rel.16, 17のホワイトペーパー。

https://www.5gamericas.org/lte-and-5g-wireless-standards-continue-to-evolve-with-each-new-3gpp-release/

 

西院春日神社でお参り2021。

神社参りには2020年のうちに沢山行った(密にならないように)のですが、

そろそろ2021年も行っとこうと訪問。まずは西院の春日神社。

20210111-11565820210111-115709

20210111-115738 20210111-115824 20210111-115810 20210111-115913

ここは一願蛙がいます。

20210111-115948

2021年1月17日 (日)

新型コロナウイルス、中国、日本、韓国、アメリカ、ドイツ、フランス、イギリスでの感染者数を指数関数&ロジスティック関数&Log-Logプロットでべき関数フィッティングした(1/17更新) 日本の増加率がイギリスやドイツレベルに達している。緊急事態宣言の効果はこれからか、、、

日本の増え方ひどいな。もうグラフが連続して見えないくらいに増えている。中国もそろそろ偽情報でごまかしきれないくらいには増加し始めた。

Coronaww202101172

 

べき乗で見ると、日本はもうイギリスやドイツと同じ増加率。緊急事態宣言が出てもこれじゃあ、、、

Coronaww202101171

三宮の東急ハンズ、やっぱり閉店していた、、、

ニュースで見たから当然なんですが、

https://sannomiya.tokyu-hands.co.jp/item/post-1455525674.html

去年末で三宮の東急ハンズが閉店。ここにはかなりお世話になったので残念。。。

ありがとう、神戸の垂れ幕が出ていた。こちらこそありがとう。

20210109-154715

そのすぐそばの生田神社は相変わらずすごい人出でした。

20210109-154823 20210109-154911 20210109-154946

 

2021年1月16日 (土)

いつの間にか京都駅前に精密な羅城門のミニチュアがあった。2016年から移設されていたというのを全く知らなかった、、、がよみがえる羅城門という計画が最近発表されたそう。

JR京都駅前に

平安京羅城門模型デジタルサイネージ

というのがあった。あれ?2016年から移設したということだが、全然気が付かなかった、、、(こちらの方はあまり行かない方向で、、、)

20210110-141009

20210110-141012 20210110-141035

でもここにあるQRコード読むと、、、

「よみがえる羅城門」を呼びかける記者会見のご報告

おお、模型じゃなくて本物をよみがえらせる試みがあるのか。それは実現したら面白いな。

 

 

2021年1月15日 (金)

無限に続く円周率Πで、1桁ずつ違う方向に動かすとランダムウォークのようになるというRandom Walk of Pi をScratch(プログラム言語)で作った。

このTweetみて、

 

これは面白いな、自分でもやってみようと思った。ただオリジナルのVennさんのはなぜか0-7までのみ。(絵が描きにくかったのかな?)
0-9までScratchでやってみた。しかも、10000桁以上で。
これです。

RandomWalk of PI

Randomwalk_of_pi

動かすとこんな感じ(GIFアニメです。クリックすると始まる)。

Randomwalk_scratch

こう見るとなんか左右に広がっているのは気のせいだろうか。

 

2021年1月14日 (木)

ノーマンズランド(誉田哲也さん、姫川玲子シリーズ文庫最新作)を読んだ。まさかあの国家犯罪と結びつくとは、、、犯人に対して複雑な思いを抱く作品。そして最後まで正体が明かされない人物とは?

これは最初は何の変哲もない殺人事件(殺人が変哲もないというのもおかしいですが)から始まり、政治家が絡んできてあれ?と思ったらあの国家犯罪(詳細はネタバレになるのでかきません)と結びつくというすごい作品。

あらすじは「東京葛飾区のマンションで女子大生が殺害された。特捜本部入りした姫川玲子班だが、容疑者として浮上した男はすでに別の男を殺した件で逮捕されていた。しかしその情報は不自然なほどに遮断され、捜査は行き詰まる。一体事件の背後に何があるのか?そして二十年前の少女失踪事件との関係はあるのか?」

というもの。その少女失踪事件との結びつきと言うのが驚く。それを小説に持ってくるとはさすが誉田さん。なので犯人がその立場に自分が立ったらどうするかと考えると複雑な気持ちに。シリーズの恒例?のグロシーンもあるけど、、、

しかし死体写真を見ながらハンバーガーを普通に食べる姫川さんにも笑う。

そして最後まで正体が明かされない(ずっと出ているのに)もまあその人しかいないと途中で思うけど驚く。

今後のシリーズでも重要人物になるのだろう。楽しみ。

20210111-154903

2021年1月13日 (水)

高周波(RF・マイクロ波・ミリ波・5G)関連ニュース2021年1月13日 IEEE Journal of Microwaves創刊、CES2021でQorvoやSkyworksがFEM関連の出展、ミニサーキットの60GHz LTCCパッケージ、など。

いよいよ新しいマイクロ波関係の雑誌が創刊された。

IEEE Journal of Microwavesだ。通常のMTTでは扱う範囲がかなり狭いので、それを垣根を取っ払って広げた感じかな。

https://ieeexplore.ieee.org/xpl/mostRecentIssue.jsp?punumber=9171629

Ieee_journal_of_microwave

内容は、CMB(宇宙マイクロ波背景放射 )から化学応用、5G/6G、レーダ、センシング、イメージング、SIW、フィルタ、量子コンピュータなどマイクロ波の基礎、応用に関わること全般。結構面白いです。

そして今やっているCES2021。日本では家電見本市のように扱われますが、RFの製品もたくさん出ている。

Qorvoや

https://www.qorvo.com/newsroom/trade-shows/ces-2021

Skyworksは

https://www.skyworksinc.com/en/Thought-Leadership/CES-2021

WiFi 6E用のFEMが新しいかな。

次はミニサーキットの60GHzまで対応したLTCCパッケージの特許。

Mini-Circuits Awarded Patent for Low-Cost LTCC Packaging Technology up to 60 GHz

1ポートキャリブレーションの基礎。こういうのはみんな知っておいた方がいい。

One-Port VNA Calibration: A Look Under the Hood

かつやで胡麻担々チキンカツ定食(ご飯大盛)を食す。麻辣感はあまりないけど胡麻坦々はいい感じ。

気になっていた本格麻辣ということで食べてみた。丼を食べている人が多いので

豚汁がついている定食のほうにした。

うーん、あまり痺れはないな。辛さもないが、胡麻坦々そのものは美味しいです。ネギと青梗菜もいい感じで入っている。

20210110-125158 20210110-125156

2021年1月12日 (火)

帰宅途中でJRが止まって90分遅れになったので、電車の中でiPhoneのメモでTwitterで聞かれた微分方程式を解く、、、と間違って削除したので、シェイクや3本指タッチで取り消すことを初めて知った、、、

今日は朝は雪で電車が遅れ、午後は人身事故で電車が90分遅れた、、、

 

 

あまりにも暇なので、たまたま私のTwitterに返信があったこれ、

https://www.wolframalpha.com/input/?i=%28y%27%27%29%5E%282%29-y%27y%27%27%27%3D%28%28y%27%29%2F%28x%29%29%5E%282%29

を計算しようと。(この方、WolframAlphaProと書いている人不特定多数に同じ文面で出しているような人なので普段ならほっとくがあまりにもひまなので、、、)

ただ、大混雑で筆記具とメモとかで計算するわけにもいかず、ただ片手でiPhoneなら入力できる、ということでメモを使って計算しようとして。。。で途中まで計算してコピペしようとして間違って消した!

うわー、、、と思ったが取り消す方法あるはず、と調べた。

昔からあるのは降る(シェイク)。今は3本指タップすればいいらしい。

【iOS 13】「シェイクで取り消し」はもう古い! 3本指でタップすればiPhoneを落とす心配なし

お、戻った。

20210112-214214

ということで計算したのがこれ。

 

 

Due to quotient rule,
y’’^2-y’y’’’=(y’’y’’-y’y’’’)/y’^2 *y’^2
=- (y’’/y’)’ *y’^2
Then
-(y’’/y’)’*y’^2=(y'/x)^2
(y’’/y’)’=-1/x^2
y’’/y’=c1+1/x
(log(y’))’=c1+1/x
log(y’)=c1*x+log(x)+log(c2)
y’=c2*x*exp(c1*x)
Due to integration by parts,
y=c2*exp(c1*x)*(c1*x-1)/c1^2+c3

 

片手で入力しつつ計算するのが死ぬほどめんどくさかった、、、紙と鉛筆の3倍くらい時間がかかったが暇が潰せた。

騙し絵の牙(塩田武士さん)を読んだ。大泉洋さんを主人公にあてがきし、映画版も今年公開。確かに大泉さんにしか思えない主人公で、いったい何に騙されているのかが最後までわからない。

前から気になっていったんですが、映画を観に行ったときに予告編が流れて、先に原作読んでおこうと買ってきた。

確かに大泉洋さんとしか顔が浮かばない主人公だ。

あらすじは「出版大手、薫風社でカルチャー誌 トリニティの編集長を務める速水(これが大泉さん)。彼は笑顔とユーモア、ウィットにとんだ会話(&ものまね)で周囲を魅了する。しかしこの本が売れない時代。ある夜、上司から廃刊の可能性をにおわされたことを機に組織に翻弄されていく。社内抗争、大物作家の大型連載、企業(パチンコ業界)とのタイアップ。速水は廃刊の危機を救えるのか?」

というもの。

だが、もうこのあらすじの時点で騙されていることがエピローグになって初めてわかる。この騙しは、ミステリ小説のトリック(叙述トリックとかいろいろ)と全く違う騙され方で、これは新鮮。騙されたというか、そうみんなが思い込んでいるだけで。

誰しもこういうところはあるだろう。面白かったです。

映画は登場人物の名前や、関係性を結構変えているみたいですが、

https://movies.shochiku.co.jp/damashienokiba/

また別の何かにだまされるのだろうか、、、見に行こう。

20210104-150225

 

2021年1月11日 (月)

文化財修理の最先端@京都国立博物館 平成知新館を観てきた。あの騎馬武者像(足利尊氏じゃなかった)の裏で弓の下絵があったり、絹を古くするのに電子ビームを使ったり、と知らなかったことがたくさん。国宝の五智如来坐像はすごい迫力。丑づくしも。

ここは予約は不要でした。あまり派手な宣伝などはしないで、今あるもので何とかいい展示をしようという

心意気が感じられてよかったです。ものすごく展示の数は多いけれど、人はものすごく少なくてゆっくり見られる。

ただ、ゆっくり見すぎると結構疲れるので、体調整えていくのがいいかも。

20210110-154853 20210110-142910 20210110-142928 20210110-143242

騎馬武者像が、昔は足利尊氏だと学校で習っていたのがひっくり返った、その絵がありました。

裏を見ると、弓の下絵があったり。

ただ、掛け軸については用語の予備知識があったほうがいいかも。私はなかったのですぐ検索した。

https://kakejikuya-nagoya.jp/yougo/

あと、絹を古くして修復するのに電子ビームつかうとかも面白かった。

膠が50年くらいしか持たないとか、いろんな知識も仕入れられた。そしてやはり1Fの仏像のコーナーは迫力がすごい。

丑づくしのコーナーも入って700円の入場料とお得です。ゆっくり見られるし。

2021年1月10日 (日)

新型コロナウイルス、中国、日本、韓国、アメリカ、ドイツ、フランス、イギリスでの感染者数を指数関数&ロジスティック関数&Log-Logプロットでべき関数フィッティングした(1/10更新) ここ数日の伸びでいうと日本が一番ひどいか。

まずはリニアスケール。ここ数日の陽性者数の伸びでいうと日本が一番ひどい(増えすぎて点がつながって見えない)。

中国も(なぜかデータが一週間遅れのものしかないが)ようやく伸びてきたことを認めたのか。

Coronaww202101101

ログログプロットすると、やはり韓国の係数が高い。イギリスとドイツが同じくらいに見える。

フランスが伸びてないのが意外と言えば以外。

Coronaww202101102

つなぐ TSUNAGU―THE POWER OF KOBE CITY MUSEUM@神戸市立博物館を観てきた。栗花落(鬼滅の刃のカナヲの苗字)が神戸に実際にあって、かつその由来に驚く。また伊能忠敬の北海道地図の正確さにも。予約制でものすごくすいていました。

予約制でメールアドレスを登録する必要があります。そのおかげか、ものすごく人が少なくて全く密にならないでゆっくり見て回れました。

栗花落っていう苗字が実際に神戸にあったのを知って驚く。ただ、読み方はカナヲのつゆり、じゃなくてつゆ、でした。

20210109-145449 20210109-140954

伊能忠敬の北海道地図もあったが、めっちゃ正確。昔の世界地図では日本が三角形でZimpaguになっていたが、、、

20210109-141014 20210109-143206 20210109-140724

こんなにゆっくり見て回れたのは久しぶりでよかった。ただ、もう緊急事態宣言が出そうなのでぜひ行ってとは言いにくいですが、、、

2021年1月 9日 (土)

NHKスペシャル 2030 未来への分岐点 (1)「暴走する温暖化 “脱炭素”への挑戦」を観てました。Twitterで実況していたリンク貼っときます。

こちらからスレッドたどれます。26ツイートしていた。

 

諏訪神社(京都)でお参り。

のぼりが立っていたのでお参りしようとしたが、結構細い路地を行くのでちょっと迷ったり。

七不思議ってなんだろう、と調べるとこういうのらしい。

https://ja.wikipedia.org/wiki/%E8%AB%8F%E8%A8%AA%E5%A4%A7%E7%A4%BE%E4%B8%83%E4%B8%8D%E6%80%9D%E8%AD%B0

20210103-125807 20210103-125905 20210103-125914

2021年1月 8日 (金)

「友達以上探偵未満」(麻耶雄嵩さん)をたまたま三重県に行く電車の中で読んだ。伊賀ももと上野あおの2人の女子高生コンビの推理だが、麻耶さんらしく一筋縄ではいかないトリック。読者への挑戦状も。

読者への挑戦状があるのが最近では珍しい、と思ったらもともとは第一話目はNHKのBSプレミアムの犯人当ての番組「謎解きLIVE」から来ているそうでした。女子高生コンビということでノリは軽く見えても謎解きは驚きの展開の全3話。

伊賀の里殺人事件、は色の違う忍者装束と芭蕉の格好をしたミステリーツアーの参加者が殺される。

よーく色と登場人物の数を数えておかないと真相にはたどり着けない、、、というもの(いや数えていてもたどり着けないか)

夢うつつ殺人事件は、美術部の女の子が夢うつつで聴いた殺人事件が本当になる、、、というものですが、

これはありなのか?と思ったりするが、よーく最初から読み直すと、だったら仕方ないか、という(笑)。

そして最後の夏の合宿殺人事件は、探偵役とワトソン役の立場についての話が出ていてそれも面白い。最初の2作で絶妙のコンビと見えたももとあおが昔は、、、というもの。

どれも麻耶さんらしい一筋縄ではいかないトリックで面白かったです。

20201228-203924

すき家でゆず香る チーズ豆乳牛鍋定食~やみつき塩麹ダレ付き~(肉二倍、ご飯大盛)を食す。うーん、チーズは美味しいけどチーズなしの方が好みかな。

この前はチーズなしのを食べてとても美味しかったのでチーズありはどうか?ということで注文。

うーん、確かにチーズはきいてるけどちょっとゆずの微妙な味とかが消されてる感があり、チーズなしの方が好きかも。

でも最後に雑炊風にするときはチーズ入りの方がコクがあってよかったり。

20201219-144820 20201219-144822 20201219-145346

2021年1月 7日 (木)

かつやでコクうま味噌カツ定食(ご飯大盛)をいただく。味噌がものすごく甘い。

ちょっと卵が残念なことになっている(黄身が沈んでいる)けれどなかなかおいしそうなビジュアル。

20201220-141417 20201220-141420

で食べると、、、うわ、めっちゃ味噌が甘い。ちょっと甘すぎでは、、、という感じ。なので卓上の七味と組み合わせるといいかと。

ご飯と豚汁は安定のおいしさですが。

20201220-141429 20201220-141427

2021年1月 6日 (水)

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はちょっと初期設定が必要。」

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

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

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

Fallenphysicistrigingengineer

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

こちらから計算式を検索してください:
https://keisan.casio.jp/

ストリップラインのSパラメータ(ロス込み)の計算をカシオの高精度計算サイトhttp://keisan.casio.jp
に自作式としてUP!Morganモデルで表面粗さを考慮し、複素特性インピーダンスを計算。
マイクロストリップラインのSパラメータ(ロス込み)の計算をカシオの高精度計算サイトkeisan.casio.jpに自作式としてUP!HammerstadとJensenモデルでインピーダンスを、Morganモデルで表面粗さを考慮し、複素特性インピーダンスを計算。
π型回路の2ポートSパラメータ計算&図示をカシオの高精度計算サイトkeisan.casio.jpの自作式としてUP!
T型回路の2ポートSパラメータの計算&グラフ化をカシオの高精度計算サイトkesan.casio.jpの自作式としてUP!
「感電の危険性(人体電流と通電時間による症状)」の計算の自作式をカシオの高精度計算サイトkeisan.casio.jpにUP!IEC60479-1に基づくもので、低電圧(<600V)の研修を受けてきたので。WebPlotDigitizerを利用した。
雨のときは走ると歩くのどちらが濡れる?(雨に走れば公式)をカシオの高精度計算サイトkeisan.casio.jpに自作式としてUP!
コロナウイルス、中国、日本、韓国、アメリカでの感染者数を指数関数とロジスティック関数でフィッティングした。(7/5更新)そりゃ連日東京で100人越えすればそうなるよな、、、という増加。アメリカは本格的に指数関数に近づいている。韓国も全然終息しない(一時終息とかいってなのはなんだったのか)
SEIRモデル(コロナウイルス感染の微分方程式その2)をカシオの高精度計算サイトkeisan.casio.jpの自作式としてUP!ルンゲクッタ法を使用。
SIRモデル(感染症の微分方程式、コロナウイルスなども)の計算をカシオの高精度計算サイトkeisan.casio.jpにUP! 4段4次のルンゲクッタ法を使用。 
”楕円の円周と一方の径からもう一方の径を求める。”をカシオの高精度計算サイトkeisan.casio.jpに自作式としてUP.
映画「アルキメデスの大戦」観てきた。面白かった!菅田将暉くん演じる天才数学者が数学で戦争を止めようとするお話。そのクライマックスのあの数式の自作式をカシオの高精度計算サイトkeisan.casio.jpにUP!
御香宮でお参り2019.算額も観てきた。そしてその問題を解いてみる。
円の弧長,弦長,矢高,半径のどれか2つを与えて残りを計算(カシオの高精度計算サイト自作式)で180°以上、複数解に対応。
東京医大の入試点数シミュレータ(男女で点数が変わる)をカシオの高精度計算サイト
keisan.casio.jp にUP!
合同数の判定(タネルの定理 Tunnell's
theorem)をする自作式をカシオの高精度計算サイトにUP!
導体の表面粗さによる電気伝導率モデル3種(Hammerstad-Morgan,Groisse
,Huray)を計算する自作式をカシオの高精度計算サイトにUP!
樹脂に粉を混ぜたときなど混合物の誘電率の計算に使われるMaxwell-Garnettモデルをカシオの高精度計算サイトkeisan.casio.jpにUP!
積を和から計算するProsthaphaeresisと対数表の動作を実際に計算してみる自作式をカシオの高精度計算サイトにUP!
矩形マイクロストリップパッチアンテナの設計式をカシオの高精度計算サイトkeisan.casio.jpにUP!
ある集まりに同じ誕生日の人がいる確率が高い、と言うのは知られていますが、ではN人中k人が一致する確率は?(今日は私の誕生日なので再掲)
6/7のGoogleトップページはヴァージニア・アプガーさん記念、ということでカシオの高精度計算サイトkeisan.casio.jpに新生児の健康を示すアプガースコアの計算をUP!
円柱導体の高周波インピーダンス(表皮効果含む)の計算をカシオの高精度計算サイトkeisan.casio.jpに自作式として作った。
お好きな正の整数の各桁を2乗して足すのを繰り返すとどうなるか?をカシオの高精度計算サイトkeisan.casio.jpの自作式としてUP!あるサイクルになる。
二重指数関数型積分公式(DE公式)は端点の特異点に強い、∫_0^1
dx/√(1-x^2)を計算すると?(カシオの高精度計算サイトで台形・シンプソン・ガウスクロンロッドなどと比較)
信頼性(寿命)評価に使われる”ワイブル分布のパラメータの算出”をカシオの高精度計算サイトに自作式として作った。
可解カオスの例として、楕円関数(ヤコビのSn)が出てくるシュレーダー写像(Schröder
map)を厳密解と精度保証付き102桁までの計算で比較(カシオの高精度計算サイト利用)
カオスを生じるロジスティック写像を厳密解と精度保証付き桁数の多い計算(カシオの高精度計算サイトkeisan.casio.jp)で桁数102桁まで変えて比較してみる。
#仮面ライダービルド アッカーマン関数の計算をカシオの高精度計算サイトkeisan.casio.jpの自作式として作ってみた。
完璧なゆで卵の作り方方程式(Peter
Barhamによる)をカシオの高精度計算サイトkeisan.casio.jpにUP!
”フーリエ級数を計算”をカシオの高精度計算サイトkeisan.casio.jpにUP!
メリークリスマスになる確率の公式をカシオの高精度計算サイトにUP!(Daily
Mailの記事より)
LambertのW関数 (第k解)をカシオの高精度計算サイト
keisan.casio.jpにUP!
オイラーのφ関数(トーシェント関数)の計算をカシオの高精度計算サイトkeisan.casio.jpにUPした。
#仮面ライダービルド
番組で出てきたロジスティック写像を50桁・精度保証で3つの異なる初期値で計算する自作式をカシオの高精度計算サイトにUP!
三角形の五心(重心・内心・外心・垂心・傍心)を計算する自作式をカシオの高精度計算サイトkeisan.casio.jpにUP!
インピーダンス⇒反射係数の変換をPower wavesとTraveling
wavesに対応したものをカシオの高精度計算サイトkeisan.casio.jpにUP!
#仮面ライダービルド 第10話の話数を表す数式は、ラマヌジャンの円周率の式
π^2=10-Σ1/(n^3*(n+1)^3)か!
余りの計算ができるカシオの電卓が大人気、と聞いて「余りの計算(調剤薬局や物流倉庫)」ができる自作式をカシオの高精度計算サイト
keisan.casio.jpにUP!
2変数のニュートン・ラフソン法をカシオの高精度計算サイトの自作式としてUPした。
4変数の常微分方程式をルンゲクッタ法(4段4次)で計算する自作式をカシオの高精度計算サイトkeisan.casio.jpにUP!
算数マジック(6個の数字を書いてもらって、隣り合うものを足して1桁目とり、、、を繰り返したものを最初から当てる)をカシオの高精度計算サイトkeisan.casio.jpにUP!
”円錐台の体積を半分にする高さを計算”をカシオの高精度計算サイトにUP! 人参の2分割のような問題。
ニンジンを同じ体積で2分割するには太いほうから1:4に分ける、というのを数学的に、そして円錐台に拡張
#すまたん の大人の数学特集見て
円の弧長,弦長,矢高,半径のどれか2つを与えて残りを計算するExcelシート(VBA利用)を作った。ダウンロードできます。カシオの高精度計算サイトのExcel移植版。
ハートを描く関数をカシオの高精度計算サイト keisan.casio.jpにUP!
ピアソンのχ^2検定(適合度)を計算する自作式をカシオの高精度計算サイトkeisan.casio.jpにUP!
ピアソンのχ^2検定(独立性の検定)の自作式をカシオの計算サイトkeisan.casio.jpにUP!
「全ての爆弾の母」と呼ばれるMOAB(Massive Ordnance Air
Blast)という大規模爆風爆弾兵器をアメリカがイスラム国に使用。サンシャイン池崎か!どんな爆弾か?
「単回帰分析(t検定と信頼区間)」を計算する自作式をカシオの計算サイト
keisan.casio.jpにUP!
2つの比率の差の検定、をカシオの計算サイトkeisan.casio.jpにUP!
2つの平均の差の検定のサンプルサイズ決定・検出力曲線をカシオの高精度計算サイトkeisan.casio.jpにUP!
サンプルサイズの決定(1つの母平均の検定)をカシオの高精度計算サイトkeisan.casio.jpの自作式として作った。
「検出力曲線(1つの平均の場合)」を描画する自作式をカシオの高精度計算サイトkeisan.casio.jpにUP!
「1つの平均のt検定」を計算する自作式をカシオの高精度計算サイトkeisan.casio.jpにUP!
「母比率の信頼区間の推定」計算式をカシオの高精度計算サイトにUP!出口調査で何人投票してたら当選確実になるか?のような。
「2群の平均の差のt検定」をカシオの高精度計算サイトkeisan.casio.jpに自作式として作ってみた。
ABテストに必要なサンプル数の計算式をカシオの高精度計算サイトkeisan.casio.jpに自作式として作ってみた。
ジェイムズ・ジョイスのユリシーズに出てくる数、ジョイス数列(9^(9^9)の桁数みたいな)をカシオの高精度計算サイトkeisan.casio.jpにUP!
誕生日のパラドクスはよく知られてますが、ではN人中K人の誕生日の一致する確率の計算は? #大人のピタゴラスイッチ
#ピタゴラスイッチ
シュリクマー関数(象の身長と前足の太さから表面積を求める)をカシオの高精度計算サイトkeisan.casio.jpの自作式としてUP!
cos(40°)^(1/3) + cos(80°)^(1/3) -
cos(20°)^(1/3) = ((3/2)*(9^(1/3)-2))^(1/3) を証明できる?(ラマヌジャンより)
ラマヌジャンのようなどうやって思いついたかわからないような円周率Πの計算公式を次々だすJesús
Guilleraさんの公式をカシオの高精度計算サイトkeisan.casio.jpで試してみる。
黄金比(フィボナッチ数列から)、白銀比(ペル数列から)、そしてプラスチック数(パドヴァン数列orペラン数列から)をカシオの高精度計算サイトkeisan.casio.jpに作ってみた。
分割数とラマヌジャン・ハーディの漸近公式の計算をカシオの高精度計算サイトkeisan.casio.jpに自作式として作った。奇蹟がくれた数式を観てきたので。
”奇蹟がくれた数式”(数学者ラマヌジャンとハーディの映画)を観てきた。よかった!出てくる数式を紹介。
ペレリマン数列(2乗するともとに戻る数)を10000桁までPari/GPで計算(天に向かって続く数、加藤文元さん著)
ペレリマン数列(2乗するともとに戻る数)の計算をカシオの高精度計算サイトkeisan.casio.jpにUP!(加藤文元さんの天に向かって続く数より)
カプレカ数(カプレカルーチン)を検証する自作式をカシオの高精度計算サイトkeisan.casio.jpにUp
御香宮(京都・伏見)の算額の問題を数値計算(ニュートン・ラフソン法)してみる(そしてカシオの高精度計算サイトkeisan.casio.jpにUP!)
#数楽
国立がん研究センターによる心筋梗塞・脳梗塞のリスクの計算式
ガチャを何回もやっても全部外れる確率、をカシオの高精度計算サイトkeisan.casio.jpにUP!
Scratchでキヨシゲーム(キヨシ関数のゲーム版)を作ってみた。ズン・ズン・ズン・ズンドコ キ・ヨ・シ!
キヨシ関数(ズン・ズン・ズン・ズンドコ)をカシオの高精度計算サイトkeisan.casio.jpに作ってみた。
[円周率の日なので再掲]円の弧長,弦長,矢高,半径のどれか2つを与えて残りを計算&円周率記事
平方剰余の計算を、カシオの高精度計算サイトkeisan.casio.jpにUP!
「あなたの飼っている犬や猫の年齢を人間に換算すると?」という計算式をカシオの高精度計算サイトkeisan.casio.jpにUP!
原子番号137以上の元素はないとファインマンが言っていた(原子番号113が理研が命名権獲得!で思い出した)
あけましておめでとうございます。2016、今年の年賀状は東山動植物園のイケメンゴリラのシャバーニにした。
浅草寺のおみくじシミュレータをカシオの高精度計算サイトkeisan.casio.jpにUp!2016年の運勢を占おう。
海辺の美女の問題(結婚相手の最初の4割捨てろ理論?)をカシオの高精度計算サイトに作った(林先生が驚く初耳学?)
#初耳学
水没したドアにかかる力の計算をカシオの高精度計算サイトにUP!(朝日新聞の記事・訂正記事見て)
標高と水の沸点の関係の計算(カシオの高精度計算サイト自作式)を再紹介-カップ麺の3分は注ぎ始めか終わりかの記事を見て。。。
2つの正規分布の比(商)の分布はコーシー分布とも限らない(平均が0のときのみ)。Fieller-Hikley分布になる。カシオの高精度計算サイトにUPしてみた。
平成27年センター試験 数学IIBの第6問、BASICプログラム(連分数)をカシオの高精度計算サイトにUP!だが、、、旧課程だけしかプログラムの問題ない。。。
友愛数の計算をカシオの高精度計算サイトkeisan.casio.jpにUP!
生年月日を1桁ずつ区切って足すと(運命数を出す一歩手前)、ものすごくきれいな正規分布になる(中心極限定理の応用)
運命数の分布の計算(1900~2015年生まれ)をやってみた(やりすぎコージーの運命数の記事のアクセスがめっちゃのびてるので)
運命数の計算(やりすぎコージー都市伝説より)をカシオの高精度計算サイト、keisan.casio.jpの自作式として作ってみたよ。
サーミスタの温度特性(Steinhart-Hart式)&そのパラメータを計算する自作式をカシオの高精度計算サイトkeisan.casio.jpにUP!【再掲】
n次方程式(複素数、n<=20)を計算する自作式をカシオの高精度計算サイトkeisan.casio.jpにUP.
健康診断で腹囲測定(メタボ検診)があって、ふと超楕円(ラメ曲線)の円周の長さはどうなるか計算した。でカシオの高精度計算サイトkeisan.casio.jpに自作式をUP.
球の抗力係数CDとレイノルズ数Reの関係
をカシオの高精度計算サイトkeisan.casio.jpの自作式にUP.
5次方程式と6次方程式をDKA法で解く自作式を再度、カシオの高精度計算サイトkeisan.casio.jpにアップしてみた。
ノイマン・コレクション 数理物理学の方法を買った。ホルツマーク分布の計算をカシオの高精度計算サイトkeisan.casio.jpにUP!
楕円の円弧の話 → カシオの高精度計算サイトkeisan.casio.jpに”楕円の円弧・弦長・矢高を長軸半径・短軸半径・中心角から求める”をUP!
Lychrel
numberは196が最小?という記事を見て、カシオの高精度計算サイトにそのあたりを計算する自作式を作ってみた。
N人中k人の誕生日が一致する確率は?(R言語のpbirthday,
qbirthday)【再掲】
電磁波の波長,反射率,透過率,電力半減深度,表皮厚みを計算するツールをカシオの高精度計算サイトに作ってみた。
脳卒中のリスクの計算式をExcelで作った。アップロードしときます。
LambertのW関数をkeisan.casio.jpの自作式として作ってみた。
平成25年センター試験 数学II・数学B 第6問(BASICの問題)をkeisan.casio.jpにUP!
Fourier係数を積分で求める自作式を作ってkeisan.casio.jpにUP!
ルンゲクッタ法(4段4次)で3変数の常微分方程式を解く自作式をkeisan.casio.jpにUP!
超楕円(ラメ曲線)を描くという自作式をkeisan.casio.jpにUP!
[再掲]”円の弧長,弦長,矢高,半径のどれか2つを与えて残りを計算”をkeisan.casio.jpの自作式で作成
[再掲]ドレイク方程式をkeisan.casio.jpにUP!
標高と水の沸点の関係をkeisan.casio.jpにUP!
Eulerの定理(数論)をkeisan.casio.jpにUP!
Eulerのφ関数をkeisan.casio.jpにUP.
ロトカ・ヴォルテラ方程式をkeisan.casio.jpにUP!
「円の弧長,弦長,矢高,半径のどれか2つを与えて残りを計算」をkeisan.casio.jpにUP!
”平成24年度センター試験 数学II・数学B 第6問”をkeisan.casio.jpにUP!
誕生日の一致する確率(N人中k人)の計算をkeisan.casio.jpにUP!
「富士山の図示」をkeisan.casio.jpにUP!
ダイポールアンテナの入力インピーダンス(ダイポール長さ変化)をkeisan.casio.jpに。
マグニチュードの差で地震規模がどう変わるかをkeisan.casio.jpに。
「NTCサーミスタのB定数の算出」をkeisan.casio.jpにUP!
ツィオルコフスキーの公式(ロケットの方程式)をkeisan.casio.jpにUP!
「双子素数」をkeisan.casio.jpにUP!
「素数階段」をkeisan.casio.jpにUP!
相対論的水素原子のエネルギー準位(Dirac方程式の解)をkeisan.casio.jpにUP!
Emden方程式の計算をkeisan.casio.jpにUP!
ポアソン乱数をkeisan.casio.jpにUP!
地震の発生頻度(グーテンベルグ・リヒター則)をkeisan.casio.jpにUP!
ゼルニケ多項式(Zernike
Polynomials)をkeisan.casio.jpにUP!
「3点を通る円」「4点を通る球」をkeisan.casio.jpにUP!
正規乱数をkeisan.casio.jpにUP!
平成22年センター試験 数学II・数学B 第6問をkeisan.casio.jpにUP

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

Kamen_rider_build_physics02

6話以降の記事のリンクはこちらのTwitterのモーメントからたどれます。

仮面ライダービルド 数式解説まとめ

#仮面ライダービルド 最終話、第49話の話数の数式667^2=444889,6667^2=44448889,66667^2=4444488889など48が間に挟まれることの証明。
#仮面ライダービルド 最終話、第49話の話数を表す数式は6667^2=44448889からだんだん6がなくなっていって最後に7^2=49になる。48話から49話への移り変わりかな。
#仮面ライダービルド 第48話の話数を表す数式はn>=48ならnと9/8nの間に必ず素数がある、というもの。nと2nの間に素数があるのはエルデシュが高校生で証明した。(ベルトラン・チェビシェフの定理)
#仮面ライダービルド
第47話の話数を表す数式はarctan(1/73)+arctan(1/132)=arctan(1/47)
#仮面ライダービルド 第46話の話数を表す数式は3^3+36^3+37^3=46^3
#仮面ライダービルド 第46話の黒板の数式はペンローズ・ホーキングの特異点定理?話の中でワームホールが出てきたのに合わせて。
#仮面ライダービルド 第45話の黒板の数式は量子情報関係かな。Popescu-Rohrlich
boxなどのワードが見えた。
#仮面ライダービルド 第45話の話数を表す数式は、オイラーの定数γを使って3^(2/γ)=44.9968...≒45
#仮面ライダービルド 第44話の話数を表す数式はモンモール数か!でも!5って記号なのかな。
#仮面ライダービルド 第44話の黒板の数式はKPZ方程式か!
#仮面ライダービルド 第43話の黒板の数式はクォークを記述するQCD(量子色力学)!閉じ込めと漸近的自由、QGP(クォーク・グルーオンプラズマ)などなど。
#仮面ライダービルド 第43話の話数を表す数式はexp(π√43)≒884736744でした。有名なexp(π√163)がほとんど整数なのと同じ理由でほとんど整数。
#仮面ライダービルド 第42話の黒板の数式はKAM理論(Kolmogorov-Arnold-Moser)に関するもの?
#仮面ライダービルド 第42話の話数を表す数式は分割数p(n)でp(10)=42。ラマヌジャンの映画でもおなじみ。
#仮面ライダービルド 第41話の黒板の数式はスピングラスかな?RSB(レプリカ対称性の破れ)が見えたり。
#仮面ライダービルド 第41話の話数を表す数式はオイラーの素数生成多項式ですね!q^2+q+41がq=0から39で全部素数になる。
#仮面ライダービルド 第40話の黒板の数式は、、、一瞬しか映らずわからなかった、、、ハミルトニアンと、L,Rが見えたので量子もつれ関連の何か?→全然違いました。密度行列繰り込み群か!
#仮面ライダービルド 第40話の話数を表す数式はペル方程式x^2-77y^2=1の最小解がx=351,y=40ということ。
#仮面ライダービルド 第39話の話数を表す数式、知らなかった!でも確かに4π^2-39-Σ16/(n^2(n+1)^2(n+2)^2)は成り立ってた!
#仮面ライダービルド 第39話の黒板の数式はワインバーグ・サラム理論(&ヒッグス機構)?
#仮面ライダービルド 第38話でクリアボードに書かれていた式はワイトマン公理系とCPT定理だった。
#仮面ライダービルド 第38話の話数を表す数式は魔六角陣(魔方陣の6角形版で縦横斜めの和が等しい)の和が38ということ!
#仮面ライダービルド 第37話の話数を表す数式は、非正則素数の最小が37ということ!
#仮面ライダービルド 第36話の黒板の数式はホーキング輻射!ボゴリューボフ変換が使われている!またベリーの位相も!
#仮面ライダービルド 第36話の話数を表す数式は立方和が三角数の二乗に等しいこと!Σk^3=(Σk)^2
#仮面ライダービルド 第35話の話数を表す数式は結び目理論!交点数8以下結び目数が35個!
#仮面ライダービルド 第34話の黒板の数式は、宇宙マイクロ波背景放射(CMB)とインフレーション!
#仮面ライダービルド 第34話の話数を表す数式はまたも巨大数、スキューズ数(第一スキューズ数SK1<10^(10^(10^34)))!
#仮面ライダービルド 第33話の黒板の数式は計算複雑性理論(PCP)に関するもの?CSP(制約充足問題)とかNP困難のワードが見えた。
#仮面ライダービルド 第33話の話数を表す数式は、ラマヌジャンの円周率公式を変形したもの!スーパーコンピュータでの円周率計算にも使われている。
#仮面ライダービルド 第32話の話数を表す数式は、結晶族点群が32種類あるということ!
#仮面ライダービルド 第31話の話数を表す数式はメルセンヌ素数!M5=2^5-1=31
#仮面ライダービルド に続き、数式がいっぱいでて科学者が主人公のアニメ、#重神機パンドーラ を観た。またこれも解説しようかな、、、
#仮面ライダービルド
第30話の黒板の数式は、量子テレポーテーション(&量子暗号?)?ちょっとしか映らなかったけどEPRペアが見えた。先週と合わせて。
#仮面ライダービルド 第30話の話数を表す数式は、自分以下の数で互いに素なものが全部素数、で最大なものが30ということ(totative
of nの最大が30)
#仮面ライダービルド アッカーマン関数の計算をカシオの高精度計算サイトkeisan.casio.jpの自作式として作ってみた。
#仮面ライダービルド #数式 新オープニングの最後に流れる数式解説まとめ
#仮面ライダービルド 第29話の黒板の数式は量子コンピュータ関係の計算と思いますがちょっとしか映らなかった、、、q-bitとかパウリ行列(Z,X)は見えた。
#仮面ライダービルド
第29話の話数を表す数式はアッカーマン関数A(3,2)=A(2,A(3,1))=29。A(m,n)はmが4より大きくなると巨大数に。
#仮面ライダービルド 第28話の話数を表す数式は、28が完全数であるということ!天才がタンクでやってくる、は馬鹿が戦車でやってくるのもじり。
#仮面ライダービルド 第27話の話数を表す数式は、Milnorのエキゾチック7次元球面は28個、、、あれ?通常の球面引いて27個か。
#仮面ライダービルド 第26話の黒板の数式はボース・アインシュタイン凝縮か!
#仮面ライダービルド 第26話の話数を表す数式は、有限単純群の分類で散在型単純群(Sporadic
group)が26個ということ。
#仮面ライダービルド 第25話の黒板の数式は、AKLTモデル&ハルデン予想?
#仮面ライダービルド 第25話の話数を表す数式はブロカールの問題(n!+1=m^2)の変形? (n-1)!+1=n^kになるのはn=5のときだけ。
#仮面ライダービルド 第24話の黒板の数式は超伝導のBCS理論!ギャップΔと3.53という転移温度の係数が見えた!
#仮面ライダービルド 第24話の話数を表す数式は平方数の和が平方数になる?というLUCAS’
SQUARE PYRAMID PROBLEM。Σn^2=m^2となるのは(1,1), (24,70)しかない。
#仮面ライダービルド 第23話の話数を表す数式はまたもやラマヌジャン!√(23-2√(23+2√(・・)=1+4√3*sin20°!
#仮面ライダービルド 第22話の話数を表す数式はまたもやラマヌジャン!2143/π^4≃22、あるいはπ≃(9^2+19^2/22)^(1/4)
#仮面ライダービルド 第21話の話数を表す数式は、連続したハーシャッド数は21個ない、ということ。
#仮面ライダービルド 第20話で戦兎が持っていた数式マグカップ、売っているのを見つけた!頭がメッチャよくなる?
マグカップ 数学
#仮面ライダービルド 第20話の黒板の数式は、量子コンピュータの量子誤り訂正(トポロジカル表面符号)か!
#仮面ライダービルド 第20話の話数を表す数式は、初めて近似式。”ほとんど整数”としてe^π√163と並んで有名なe^π-π≃20。
#仮面ライダービルド 第19話の黒板の数式は、、、端っこしか映らなかったのでわからなかった、、、量子力学の計算(|a>が見えた。何かの準位?)とは思いますが、、、⇒量子コンピュータのGrover探索アルゴリズムとのこと!
#仮面ライダービルド 第19話の話数を表す数式は、ラマヌジャンの見つけた(7(20)^1/3-19)^1/6=(5/3)^1/3-(2/3)^1/3。ガロア理論とも関連。
#仮面ライダービルド 第18話の黒板の数式は、量子ゼノン効果かな。
#仮面ライダービルド 第18話の話数を表す数式は、有限群の族の数、、、って18個だっけ?数え方によるのか、、、
#仮面ライダービルド 第17話の黒板の数式は、暗号の公開鍵で、、、Little
Dragonとか?Dragon 型多変数公開鍵暗号
#仮面ライダービルド 第17話の話数を表す数式は、フェルマー数Fn=2^2^n+1でn=2としたもの!
#仮面ライダービルド と千葉先生のおかげで蔵本由紀先生がクローズアップ。蔵本転移以外にも蔵本・シバシンスキーや複素TDGL方程式などもうちのブログで何回も取り上げていたので再掲。
#仮面ライダービルド 日経サイエンス2018年2月号買ってきた。仮面ライダーの方程式、という特集で監修の白石さんがなんと11ページ!(ゆらぎ、蔵本転移と暗号)の記事を書かれている。
#仮面ライダービルド 第16話の話数を表す数式は、クヌースの矢印記法!グラハム数と寿司 虚空編を思い出す。。。
#仮面ライダービルド 第15話の黒板の数式は、、、黒板をはみ出して部屋中に。流体力学極限で粘性をミクロから計算とかそういうやつ?⇒佐々さんの「ハミルトン粒子系から流体方程式へ」でした!
#仮面ライダービルド 第15話の話数を表す数式は、ラマヌジャン・スコーレムの定理!2^n-7=m^2を満たす自然数は?
#仮面ライダービルド 第14話の必殺技の時に出てきたのはワームホールの計量か!
#仮面ライダービルド 第14話の黒板の数式はエントロピー生成と揺らぎの定理など。
#仮面ライダービルド 第14話の話数を表す数式は、オイラーのトーシェント関数?でも14が何かわからない、、、⇒ノントーシェントでした。
#仮面ライダービルド 第13話の黒板の数式はほとんど映らず(位相空間の積分が見えた)。その代わりに必殺技がMaxwell方程式!
#仮面ライダービルド 第13話の話数を表す数式はアルキメデスの立体の数が13個、ということか!
#仮面ライダービルド
番組で出てきたロジスティック写像を50桁・精度保証で3つの異なる初期値で計算する自作式をカシオの高精度計算サイトにUP!
#仮面ライダービルド 第12話の黒板の数式は、、、ほとんど映らなかったけどロジスティック写像のリアプノフ指数の計算のように見えた。
#仮面ライダービルド 第12話の話数を表す数式はゼータ関数ζ(s)=∑1/m^sでζ(-1)=-1/12、つまり1+2+3+…=-1/12という衝撃の式。
#仮面ライダービルド 第11話の話数を表す数式、ミルの公式でn=12までPARI/GPで計算してみた。
#仮面ライダービルド 第11話の話数を表す数式は、素数を生み出すミルの定数(ミルの公式)か!
floor(A^(3^n)
#仮面ライダービルド 第11話の黒板の数式は、蔵本モデルの同期現象の”相転移”!
#仮面ライダービルド 第10話の話数を表す数式は、ラマヌジャンの円周率の式
π^2=10-Σ1/(n^3*(n+1)^3)か!
#仮面ライダービルド 第9話の黒板の数式はドラゴン曲線ということで
#Scratch で縮小写像系を使って描いてみる。
#仮面ライダービルド
第9話の黒板の数式は、、、一瞬だけ、ドラゴンカーブみたいなのが映ったけど本当に一瞬でよくわからない。
#仮面ライダービルド 第9話の話数を表す数式はカタラン予想か!x-y=1(x,y∈m^n)ということで、べき乗の差が1になるのは3^2-2^3=1でx=9だけということ。
#仮面ライダービルド 第八話の黒板の数式は、楕円曲線暗号と思いきや、RSA暗号か?先週のディッフィーヘルマン鍵交換に続き。
#仮面ライダービルド 第八話の話数を表す数式は、Hurwitzの定理(1,2,4,8の定理)?なんてマニアック。
#仮面ライダービルド 第7話の黒板の数式はRSA暗号かな?⇒ディフィー・ヘルマン鍵共有でした。劇中でアナグラム暗号が出てきたのに合わせて。
#仮面ライダービルド 第7話の話数を表す数式は、トーラスの塗り分けは7色いるというものか!(平面は4色で四色問題になる、というのが第四話でした)
#仮面ライダービルド 第六話の黒板の数式は、、、わからなかった。磁気遮蔽とかそういうやつ?→メタマテリアルの透明マントでした。
#仮面ライダービルド 第六話の話数を表す数式はバーゼル問題(Σ1/n^2=π^2/6)でした。
#仮面ライダービルド 第五話の黒板に書かれていた式は重力波の一般相対性理論(アインシュタイン方程式)からの摂動での導出かな。
#仮面ライダービルド
第五話の話数を表す数式は、五次方程式以上の代数方程式は代数的に解けない、ということ?アーベル・ルフィニの定理やガロア理論。
仮面ライダービルド、第四話の話数を表す式は、、、四色問題?
仮面ライダービルド、第四話の黒板に書かれている数式は特殊相対性理論のローレンツ変換。
仮面ライダービルド、第三話の黒板の数式はLeggett-Garg
Inequalityに関わるもの?
仮面ライダービルド、第二話の黒板の式はArthurs-Kelly-Goodmanの不等式とのこと。知らなかった、、、
仮面ライダービルド、オープニング曲で数式(エントロピー、ハイゼンベルク方程式、ローレンツ変換)が壁を突き破ってくる斬新な演出
物理学者、仮面ライダービルドのライダーキックのグラフの位置を取得して関数形を推定。サイクロイド(最速降下曲線)じゃなかった。
物理学者が主人公の仮面ライダービルド、第一話に出てきた式は、、、不確定性原理の一般化(ロバートソンの不等式)やガウス積分などなど。
物理学者が仮面ライダーになる、という仮面ライダービルド。黒板に書いている数式を解説。ファインマン経路積分、フォッカープランク方程式、Lindblad方程式、分配関数などなど(追記しました)

新年マグロの初競り価格を1999年からグラフにプロット(リニア&ログ)。2021年は去年の1/10近い2084万円に。こう見ると去年、一昨年、2013年は異常。

毎年マグロの初競り価格をプロットしている。

今年の値段は2084万円。ああ、めちゃくちゃ下がったな。コロナの影響か。

https://news.yahoo.co.jp/articles/427616d7bf85c0a5829fdbc93407d66bf87073f0

ということで1999年からの価格の推移をプロット。

リニア:

Maguro20211

ログ:

Maguro20212_20210106080701

こう見ると今年の落ちがすごいけど、逆に去年、一昨年、2013年以前がおかしくて

元に戻ったように見える。ただ、マグロ漁師さんはがっかりだろう、、、

十輪寺(京都大原野、なりひら寺)でお参り。

ここも善峯寺へ向かう道すがらでお参り。在原業平が晩年いたそうで、なりひら寺とも呼ばれているそうです。

 

20210102-135451 20210102-135515 20210102-135546 20210102-135552 20210102-135606じゅ20210102-14114120210102-135748

2021年1月 5日 (火)

ジャンプの編集者(呪術廻戦の担当)で情報幾何の修士持っているというのは早稲田卒の高野健さんみたい。研究会賞までもらっている本格的な研究されてたようです。

さっき

呪術廻戦 【呪術廻戦】五条悟の無下限呪術、情報幾何修士号のジャンプ編集に間違っていると指摘される 【呪術廻戦】五条悟の無下限呪術、情報幾何修士号のジャンプ編集に間違っていると指摘される

と言うのを見た。ほぉ、ジャンプの編集者でそんな方がいるのか、と調べてみる。

『鬼滅の刃』担当編集者座談会

にも出てきた高野健さんのようだ。以下の写真とたぶん同じ人みたい。

研究会賞 | 情報論的学習理論と機械学習 (IBISML) 研究会


マーク付き点過程の距離計算手法と判別への応用

という発表で、2014年の研究会賞までもらっている。

それだけでなくて

非負値行列分解の情報幾何

など7件も発表。

https://www.ieice.org/ken/search/index.php?search_mode=form&year=39&psort=1&pskey=author%3A%22%E9%AB%98%E9%87%8E+%E5%81%A5%22&ps2=1&layout=&lang=&term=AUTHOR

すごいな、そういう人が集英社に就職してジャンプの編集者やっているとは!

 

 

大歳神社(京都大原野)でお参り。

善峯寺へ向かう道すがらで見つけた神社。初めてお参りです。

 

20210102-133130 20210102-133140 20210102-133148

2021年1月 4日 (月)

全部が素数である等差数列の見つかっている最長の記録は27個でa(n)=224584605939537911 + 81292139·23#·nだそうなのでPARI/GPで見てみる。

素数の列が任意の長さの等差数列を含んでいるということが証明されているそうだ。

グリーン・タオの定理

では実際に数値的に確認されている記録は?というと最長で27個で、019年に見つかった。

Primes in Arithmetic Progression Records

数列はこんな形で、ここで23#は素数階乗(23までの素数の積)

224584605939537911 + 81292139·23#·n

疑い深いので、PARI/GPで確かめてみる。23#=223092870で最初に計算している。
結果はこれ。数字の後ろの1は素数(isprime(a)を使っている)。0は合成数。

n , 数列   , 素数なら1
0 , 224584605939537911 , 1
1 , 242720302537486841 , 1
2 , 260855999135435771 , 1
3 , 278991695733384701 , 1
4 , 297127392331333631 , 1
5 , 315263088929282561 , 1
6 , 333398785527231491 , 1
7 , 351534482125180421 , 1
8 , 369670178723129351 , 1
9 , 387805875321078281 , 1
10 , 405941571919027211 , 1
11 , 424077268516976141 , 1
12 , 442212965114925071 , 1
13 , 460348661712874001 , 1
14 , 478484358310822931 , 1
15 , 496620054908771861 , 1
16 , 514755751506720791 , 1
17 , 532891448104669721 , 1
18 , 551027144702618651 , 1
19 , 569162841300567581 , 1
20 , 587298537898516511 , 1
21 , 605434234496465441 , 1
22 , 623569931094414371 , 1
23 , 641705627692363301 , 1
24 , 659841324290312231 , 1
25 , 677977020888261161 , 1
26 , 696112717486210091 , 1
27 , 714248414084159021 , 0
28 , 732384110682107951 , 0
29 , 750519807280056881 , 0
30 , 768655503878005811 , 0

おお、本当だ。n=0~26までは素数で27から合成数になった。

20201231/20210101は両方とも素数なのに年は2020=2*2*5*101,2021=43*47とも合成数。こういうことが起きるのは次は25331231/25340101(PARI/GPで確認)

6年前に

20201231/20210101は素数、20291231/20300101も素数、そして21611231/21620101も素数

というのを書いて、まだまだ先かと思ったらもう来てしまった、、、年月の経つのは早い。

一方、このTweetみて面白いなと思った。

 

なるほど、年は素数っぽく見えるけど合成数なのか。

そんなことが次に起きるのはいつなのか、またPARI/GPで計算してみよう。3000年末まで見てみた。

参考にしたのはこちら:

pari-gpで初等数論

計算結果:

20201231 素数  2020=2*2*5*101
20210101 素数 2021=43*47

20291231 素数  2029  素数
20300101 素数 2030=2*5*7*29

21611231 素数  2161  素数
21620101 素数  2162=2*23*47

23831231 素数  2383  素数
23840101 素数  2384=2*2*2*2*149

24371231 素数  2437  素数
24380101 素数  2438=2*23*53

25331231 素数  2533=17*149
25340101 素数  2534=2*7*181

25811231 素数  2581=29*89
25820101 素数  2582*2*1291

26021231 素数  2602=2*1301
26030101 素数  2603=19*137

26501231 素数  2650=2*5*5*53
26510101 素数  2651=11*241

26891231 素数  2689 素数
26900101 素数  2690=2*5*269

29321231 素数  2932=2*2*753
29330101 素数  2933=7*419

ちょっと見やすくした表:

Prime_date

次に年月日が連続して素数なのは

20291231と20300101だけれど、その時の年は片方が素数。次もそう。

次に年が両方とも合成数になるのは、なんと25331231/25340101。

さすがにここまでは生きてない、、、か?(脳だけになって生きる時代が数十年で来れば、、、)

善峯寺でお参り2021。景色がすごい。しかし下からずっと徒歩で登っていったのでさすがに疲れた、、、

そういや善峯寺には1回しかいったことないな、と思い立って訪問。

下から、というかはるかかなたの自宅から歩いて行った。かなり登るので最後はへろへろになった。。。

20210102-142754 20210102-142931

日本一の松がある。

20210102-143318 20210102-143404 20210102-143517 20210102-143544

ああ、松の日本一って横に伸びるほうか。

20210102-143553 20210102-143602

幸福を招くお地蔵さんがある。

20210102-143937 20210102-143942

登るだけあって、景色はどんどん壮大になっていきます。

20210102-144054 20210102-144122 20210102-144313 20210102-144328 20210102-144812 20210102-144820

ここからまた徒歩で帰る。疲れた、、、

2021年1月 3日 (日)

すき家でわさび山かけ牛丼 塩麹とん汁セット(特盛)を食す。山かけと塩麴がなかなか合って美味しい。

塩麴の鍋が美味しかったのでこちらも注文。

半分を豚汁に、半分を牛丼にかけていただく。

20201227-124110 20201227-124112 20201227-124114

わさび山かけ初めて食べたが塩麹とも合って美味しかった。

2021年1月 2日 (土)

京都に続いて2021年長岡京初詣も、密にならないように2020年12月末に徒歩で6社お参り。

昨日は京都市内の初詣で15社回った。今日は長岡京市。

まずは、やはり長岡天神から。

20201231-090723 20201231-090027 20201231-090038 20201231-090059

次は走田神社。何年か前に木がなくなってから眺めがよくなった。

20201231-094620 20201231-094256 20201231-094506 20201231-094120

次は長方稲荷神社。道が濡れていて急な坂でこけそうに、、、

20201231-095838 20201231-095744

次は長法寺へ。ここも結構坂を登る。

20201231-101318 20201231-101343 20201231-101348 20201231-101349

次は子守勝手神社。さらに登る。

20201231-102234 20201231-102243 20201231-102313 20201231-102345 20201231-102341

そして最後は光明寺。やっぱり除夜の鐘は一般の人はつけなくて(コロナの影響で)僧侶のみとなっている。

20201231-103744 20201231-104127 20201231-104148 20201231-104202 20201231-104334 20201231-104431

 

 

 

2021年1月 1日 (金)

2021年京都初詣は、密にならないように2020年12月末に15社徒歩で回ってお参り。Googleマップで神社を検索して一番近いところから回る。

初詣ではないのですが、まあ仕方ない。超メジャーどころも避けて、密にならないようにお参りした。

出発は四条新町から(たまたまそうなった)。ルールは、Googleマップで神社を検索して、一番近いところを回っていく。

では最初から。

龍田稲荷神社

20201230-121519

越後神社

20201230-122646

幾世稲荷大明神

20201230-123413

玉藏稲荷大明神

20201230-124304

高松神明神社 鳥居が白い!

20201230-124934

御金神社 鳥居が金ぴか!

20201230-125553

源正稲荷大明神

20201230-130031

福永弁財天 これは探すのに難儀する、、、

20201230-130755

蛭子大神宮・末永大明神 これも探すのに難儀。

20201230-132021

菅原院天満宮神社、梅丸大明神 菅原道真の産湯の井戸がある。

20201230-132947 20201230-133047 20201230-133104 20201230-133106

護王神社 ここは足腰の神社と共に、いのししの神社でもある。なので、鬼滅の刃の伊之助のフレームがあったり。

20201230-133541 20201230-133548 20201230-133723 20201230-133605

福長神社

20201230-134952

霊光殿天満宮

20201230-135350

白峯神宮 ここはサッカーで有名。色んなボールが奉納してある。

 20201230-140340

20201230-140111 20201230-140316

晴明神社 最後にここの桃を触って終了。

20201230-141106 20201230-141129 20201230-141139

これで2時間以上かかっているので結構疲れた、、、

目標の20社には届かなかった。

グーグルマップのタイムラインはこちら:

Googlemaps20201230

2021 あけましておめでとうございます。Happy New Yearの方程式(w + log |N| = log|ap^2H + e^aR| + log|y|) を解いてみよう。あと Σ1/(5^n*(1/3*n^3+3*n^2-1/3*n+8))を計算すると20210101に。

あけましておめでとうございます。去年と同じですがこれからスタート。
w + log |N| = log|ap^2H + e^aR| + log|y|
を計算してみましょう。GIFアニメにしてみた。
Happynewyear

 

元ネタはこのあたりから:
で、Scratchでもこれをアップロードしておきました。
Scratch_new_year
Happynewyearscratch

 

で終わるのもあれなんで、ISC(Inverse Symbolic Caluculator)で20210101を探すと、、、

20210101

になった。これの10^9乗を見てみると?

f1 f1
1 18181818.181818181818181818181818181818181818181818
2 20000000
3 20186046.511627906976744186046511627906976744186047
4 20207099.143206854345165238678090575275397796817625
5 20209700.769223114507766864694350737877023813077788
6 20210044.855244619884110950715856114221109834583164
7 20210092.795319526251152149217728773397139797129981
8 20210099.751841265381586931826424425571052840608242
9 20210100.794611122815395485797911187281847137960584
10 20210100.955112690213514608054964478817897294700396
11 20210100.980365464566165656143743763651436135637511
12 20210100.98441289539620518183544336839451913959008
13 20210100.985071946080034626726835162441180442888551
14 20210100.985180737580698637351005149160967959488817
15 20210100.985198911734886102686557007175388369915883
16 20210100.9852019798996801101771937112952385571818
17 20210100.985202502723772651062714253776291608637724
18 20210100.985202592560647222687114527936675433175078
19 20210100.985202608113539539507049265403302546791216
20 20210100.985202610824439229269200247823157769128341

 

うーん、あともう少しだが今年はこれで。

« 2020年12月 | トップページ | 2021年2月 »

最近の記事

最近のコメント

2021年7月
        1 2 3
4 5 6 7 8 9 10
11 12 13 14 15 16 17
18 19 20 21 22 23 24
25 26 27 28 29 30 31
フォト
無料ブログはココログ