陰的ルンゲクッタ(ラダウIIA)でOregonatorを計算(Excel VBA)
さて、スティフな(硬い)方程式の例としてまずはオレゴネータかな。8次のルンゲクッタで計算したときに、すさまじく計算しにくかったです。
http://www.scholarpedia.org/article/Oregonator
εdx/dt = q*y - x*y + x*(1-x)
ε'dx/dt = -q*y - x * y + f *z
dz/dt = x - z
で計算結果はこちら。何も考えなくてもふつーに計算できる!!!!すごい!!!!
これもソースをつけておきます。
「implicit-Oregonator1.txt」をダウンロード
« フィリピンの刑務所で受刑者が踊るThey don't care about usは感動的。 | トップページ | くじら雲 »
「学問・資格」カテゴリの記事
- 高周波・RFニュース 2025年1月13日 IEEE Microwave Magazineの特集はニューラルネットワークとマイクロ波、Siversがミリ波ビームフォーマー開発を受注、バイデン・ハリス政権が ワイヤレス革命に1億 1,700 万ドル、HoneywellとNXPが航空機技術で提携(2025.01.13)
- UnityでVisual C#用の数値計算ライブラリMath.NET numericsを使う(1) まずはNuGetForUnityを使ってインストール。2Dゲーム画面に連立方程式を解いた結果を表示。(2025.01.14)
- 高周波・RFニュース 2025年1月9日 CES2025に合わせて各社プレスリリース、特にQualcomm、NVIDIA、INTELが大量。SEMCOのC0G MLCC (1210 inch, 22nF, 1000V)解説、TIのAI搭載60GHz車内レーダ、MarvellのCPO、Qorvoの車載UWB SoC、TDKのセンサがAI白杖に採用(2025.01.09)
- NHK パンサー尾形さんの笑わない数学 微分・積分 スペシャルがもうすぐ始まる。これから見てリアルタイムでポストしたのでそのスレッドを残す。(2024.12.29)
« フィリピンの刑務所で受刑者が踊るThey don't care about usは感動的。 | トップページ | くじら雲 »
コメント