拡散方程式をDOP853（ドルマンプリンス-ルンゲクッタ8次）で計算（Excel VBA）
ルンゲクッタ8次のルーチン、DOP853は偏微分方程式にも適用できる。まずは拡散方程式でもやってみよう。
∂u(x,t)/∂t = D ∆u(x,t)
を
du_i/dt = D (u_i+1 + u_i-1 - 2*u_i) / h^2
のように空間部分を差分化すると、これは普通の常微分方程式になった。でそのままDOP853で計算すると、
と一瞬で計算できた。ではこれから偏微分方程式シリーズを少しやってみよう（続く）。
« WiiのビリーズブートキャンプのCMがセイバンの「天使のはね」だというのは度肝を抜かれた！ | トップページ | KdV方程式をDOP853（ルンゲクッタ8次）でExcel VBAで計算 »
「学問・資格」カテゴリの記事
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2019年11月11日 Huawei Mate 30 Pro 5G分解、AppleのUWBチップU1解析、Qualcommの5G RFモジュール解説など。(2019.11.11)
- 又吉直樹のヘウレーカ 「世界最速！？量子コンピューターって何？」を見てました。リアルタイムツイートのまとめ。(2019.11.07)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2019年10月28日 AppleのUWBモジュール（U1チップ内蔵）解析、GSAのミリ波レポート、Infineonの60GHzレーダセンサの配線とか。(2019.10.28)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2019年10月24日 Qualcommが5G関連企業に200Mドル出資＆5Gの利点を プレゼン、放熱機構がすごいMicrosoft Surface Laptopや60GHzセンサが搭載のGoogle Pixel 4 分解(2019.10.24)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2019年10月12日 Qualcommのセルラー通信解説、Skyworks/BroadcomのWi-Fi6製品、台風19号（ハギビス）の118GHzでの観測(2019.10.12)
« WiiのビリーズブートキャンプのCMがセイバンの「天使のはね」だというのは度肝を抜かれた！ | トップページ | KdV方程式をDOP853（ルンゲクッタ8次）でExcel VBAで計算 »
コメント