LambertのW関数 (第k解）をカシオの高精度計算サイト keisan.casio.jpにUP!
z=w*exp(w)の逆関数、w=W(z)であるLambert（ランベルト）のW関数については、何回か自作式を作ったけど、結構間違っていたりした、、、
今回、MapleやMatlabでの実装を参考に、WolframAlphaでも検証したものに作り替えたので再UP!
ちゃんと第k分岐の解まで計算できるようになった。
リンクはこちら。
説明文：
ランベルトのW関数は z=w*exp(w)の逆関数 w=W(z)です。主要解の他に第k解まで計算できます。
|
最終的に画面イメージはこんな感じに。
参考文献は、大本はこれ。
”On the Lambert W function”
Matlabのも参考にしました。
そしてWolframAlphaについては
ProductLog[k,z]
で計算できる。
例えばこんな感じ。
« キャプテンサンダーボルト（上下、文庫版）（阿部和重＆伊坂幸太郎）を読んだ。期待通りめちゃくちゃ面白かった。 | トップページ | 橿原神宮へ行ってきた。神武天皇の絵巻物が公開されていた。 »
「パソコン・インターネット」カテゴリの記事
- 高周波回路シミュレータQucsStudioを使ってみる（その３）Mixed Mode S parameterを計算(2018.10.12)
- 高周波回路シミュレータQucsStudioを使ってみる（その２）SパラメータのTouchStoneフォーマットで出力するには？(2018.10.11)
- 高周波回路シミュレータQucsStudioを使ってみる（その１）まずは何をさておきμの文字化けだけには注意。(2018.10.10)
- 円の弧長,弦長,矢高,半径のどれか2つを与えて残りを計算（カシオの高精度計算サイト自作式）で180°以上、複数解に対応。(2018.10.09)
- macbook proをmacOS Mojaveにアップデート。せっかくなんでダークモードにしてみる。(2018.09.26)
「学問・資格」カテゴリの記事
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2021年2月16日 IEEE Microwave Magazineの特集はオールデジタルのRFID、Microwave JournalはＥバンド ミリ波通信に衛星や気球を使う話、アメリカの半導体企業がバイデンに投資を迫る、(2021.02.17)
- カオスを生じる電気回路、Chua’s circuitをLTspiceで回路シミュレーションしてみる。(2021.02.19)
- Labyrinth Chaos(迷宮カオス）を生むThomas-Rössler方程式のパラメータbを色々変えて、Python＋Scipyでルンゲクッタ8次のDOP853(Dormand＆Prince)を使って計算してGIFアニメ(2021.02.16)
- フィッツヒュー・南雲 (FitzHugh-Nagumo) 方程式をPython+Scipyでルンゲクッタ8次のDOP853(Dormand Prince)で計算。(2021.02.23)
- 「水晶振動子の等価回路計算」をカシオの高精度計算サイトkeisan.casio.jpの自作式としてUP！ インピーダンスの大きさと位相がグラフ化できる。(2021.02.12)
「日記・コラム・つぶやき」カテゴリの記事
- ExcelのLET関数＋SEQUENCE関数で数値計算シリーズ（その5）ワンライナーでラマヌジャンの円周率公式を使って計算する。n=2で収束した。(2021.03.01)
- mRNAワクチンやCRISPR cas9のことが全然わからんので昔買ったブルーバックスの「アメリカ版 大学生物学の教科書」を3巻まで読む、、、と思ったら11年ぶりに完全改訂されたものが出てた！(2021.03.02)
- 新型コロナウイルス、中国、日本、韓国、アメリカ、ドイツ、フランス、イギリスでの感染者数を指数関数＆ロジスティック関数＆Log-Logプロットでべき関数フィッティングした（2/28更新）さすがにどの国も伸びは鈍ったか。(2021.03.01)
コメント
« キャプテンサンダーボルト（上下、文庫版）（阿部和重＆伊坂幸太郎）を読んだ。期待通りめちゃくちゃ面白かった。 | トップページ | 橿原神宮へ行ってきた。神武天皇の絵巻物が公開されていた。 »
民視眼鏡 - 媒體認證台北市士林最便宜眼鏡配到好 http://www.people-eye.com.tw/
投稿: Arnoldkasia | 2017年11月29日 (水) 23時10分