拡散方程式を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ニュース 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)
« WiiのビリーズブートキャンプのCMがセイバンの「天使のはね」だというのは度肝を抜かれた！ | トップページ | KdV方程式をDOP853（ルンゲクッタ8次）でExcel VBAで計算 »
コメント