35段14次ルンゲクッタ法（適応刻み幅付き）でLorenz方程式を計算(PARI/GP)
おなじみのローレンツ方程式
dx/dt = -σ*(x-y)
dy/dt=-y-x*z+r*x
dz/dt=x*y-b*z
をちゃんと精度見ながら計算しようと思う。PARI/GPで60桁精度で、35段14次のルンゲクッタ法で、tol = 1.E-20で計算してみた。
図示は久々にgnuplotを使用。
5回以上このブログではこのローレンツ方程式を計算してきたような気がするが、これが一番精度がいいはず。
残念ながら4段4次のルンゲクッタと形が大幅にずれるわけでなく、あまり面白くはないが。。。
« 「鹿男あをによし」を読んだ。 | トップページ | 新田辺駅をぶらぶら。 »
「学問・資格」カテゴリの記事
- 高周波・RFニュース 2025年1月23日 5G Americasの新ホワイトペーパー「AI時代のセルラーネットワークの信頼性とセキュリティ」、KyoceraAVXの新薄膜フィルタ、TDKの車載／一般用C0G特性1,250V 3225サイズMLCC、Semtechの5G LPWAモジュール(2025.01.23)
- 高周波・RFニュース 2025年1月22日 everythingRFマガジンにMarkiの宇宙向けミリ波部品の記事、NordicのRF52810を使った太陽電池で動き暗闇でも3週間持つアセットトトラッカー、KnowlessのMRIの技術解説記事、Broadcomの3.5Dパッケージング解説(2025.01.22)
- UnityでVisual C#用の数値計算ライブラリMath.NET numericsを使う(3) 3D画面に補間（Interpolate) を行って表示する。リニア、３次スプライン、有理関数などいろいろ使える。(2025.01.23)
コメント