愛の方程式、i^i^i^i^i^...(虚数iのi乗のi乗のi乗の、、、つまり愛の愛情の愛情の愛情の、、、)をPython(SciPy)のLambertのW関数で計算する。 そしてexp(π/2)^exp(π/2)^exp(π/2)^… = -i (虚数)になる。
ずっと前に書いた話ですが、、、
愛の方程式、i^i^i^i^i^...(iのi乗のi乗のi乗の、、、つまり愛の愛情の愛情の愛情の、、、)の答えは?
というのをPythonでやろうと思う。こういうやつ。
普通にプロットすると、、、(10万回繰り返す)
どこかに収束している。その値は、
(0.4382829367270323+0.3605924718713855j)
になった。これをLambertのW関数を使って求めよう。
c=zzzzzzz…
と置いて、両辺のlogをとる。
log(c) = c*log(z)
ここで
c = exp(log(c))
なので、
log(c) = exp(log(c)) * log(z)
log(z) = exp(-log(c))*log(c)
ということで
(-log(z)) = (-log(c)) * exp( -log(c))
となる。
ところでLambertのW関数というのは
y = w*exp(w) となるようなw=W(y)のことをいう。
-log(c) = W(-log(z))
c = 1/exp(W(-log(z))
これで終わってもいいですが、
W(-log(z))*exp(W(-log(z))) = -log(z)
でもあるので、
c = W(-log(z)) / (-log(z))
ともかける。
ではこれをSciPyで計算しよう。簡単に、
from scipy.special import lambertw
z = 1j
c = lambertw(-np.log(z))/(-np.log(z))
print(c)
で答えは、
となってさっきと一致!
ちなみに、z=√2 とすると、
from scipy.special import lambertw
z = np.sqrt(2)
c = lambertw(-np.log(z))/(-np.log(z))
print(c)
(2.0000000000000004-0j)
2ですね。これを参照。
√2の√2乗の√2乗の、、、、を無限に繰り返すと? (ちょっと追記)
そしてz = exp(π/2)とすると、、、なんと
from scipy.special import lambertw
z = np.exp(np.pi/2)
c = lambertw(-np.log(z))/(-np.log(z))
print(c)
(-1.0028537113604063e-17-1j)
-iになる!
« tan1°が有理数か無理数かという入試問題があったが、Python+Sympyで実際に計算してみる。 | トップページ | 京都の向日市のJR西日本の車両基地?で黄色い電車を見た。キヤ141系「ドクターWEST(ウェスト)」というやはり検査用の車両だそう。 »
「日記・コラム・つぶやき」カテゴリの記事
- 高周波・RFニュース 2025年2月7日 フジクラがMWCバルセロナでミリ波(28GHz,60GHz)モジュールデモ、Qorvo、NXP、AMDらのヘルスケアウェビナー、Samsung Galaxy S25もやはり5Gミリ波AiPは1個、u-bloxがBluetooth LEモジュール4種発表(2025.02.07)
- 高周波・RFニュース 2025年2月6日 3GPPがRel-20(5G Advancedと6G)の計画を公開、SEMCOが産業用MLCC(0805 inch, X6S, 22uF, 25V)と車載MLCC(1005,2.2uF,10V)を発表、KYOCERA AVXのスーパーキャパシタシミュレーションソフト、I向けフォトニックチップ解説(2025.02.06)
- 高周波・RFニュース 2025年2月5日 Mavenirの5G NTNホワイトペーパー、ミニサーキットの1.5mm角QFNパッケージMMIC、Keysigntがスペインに6G研究ラボ設立、アイルランドとスコットランドの大学がRIS(reconfigurable intelligent surfaces)で協力(2025.02.05)
- 高周波・RFニュース 2025年2月4日 Samsung Galaxy S25plus、S25ultra分解、やっぱり5Gミリ波AiPは1個だけ、スマートワイヤレス聴診器にNordicのICが採用、TDKのTMR角度センサ解説(2025.02.04)
« tan1°が有理数か無理数かという入試問題があったが、Python+Sympyで実際に計算してみる。 | トップページ | 京都の向日市のJR西日本の車両基地?で黄色い電車を見た。キヤ141系「ドクターWEST(ウェスト)」というやはり検査用の車両だそう。 »
コメント