時代に逆行してCOBOL(GnuCOBOL)を学んでみる(5) 中二病のような名前のブルースカイカタストロフィ(Blue-sky catastrophe)を生じるGavrilov Shilnikov modelをドルマン・プリンス法（ode45)で計算する。テキストベース（アスキーアート）で！
さて前回はCOBOLで４段４次のルンゲクッタ法でローレンツ方程式を計算してみた。
今回はさらに高次のドルマンプリンス法を使ってみよう。COBOLでやる意味があるのかどうなのかというのはあるが、まあ31桁で10進計算ができるということで。
前にPythonでやったのを参考に：
Python＋Scipyでルンゲクッタ8次のDOP853(Dormand Prince)を使う（その11） 中二病用語、ブルースカイ・カタストロフィを生じるGavrilov Shilnikov modelを計算してGIFアニメに。
ソースコードはこんな感じで（長い…配列にすればよかった）。
結果：ちょっと薄い…がちゃんと望む形が計算出来た。
これまでのシリーズ：
時代に逆行してCOBOL(GnuCOBOL)を学んでみる(3) ロジスティック写像の分岐図を31桁まで10進計算ができることを生かして描いてみる。テキストベースで！今回は2次元配列の練習。
« 河童ラーメン本舗で河童ラーメン煮卵付きをいただく。私はここのラーメン一番好き。替え玉も無料。 | トップページ | 平均が0でなく相関がある正規分布の比(商)Z=X/Yの分布 Fieller-Hinkleyの式（平均が0ならコーシー分布）をカシオの高精度計算サイトkeisan.casio.jpにあげていたが式にミスがあったのでやり直し。Pythonでモンテカルロして全然違う値になったのでもっと新しい論文探した。 »
「パソコン・インターネット」カテゴリの記事
「学問・資格」カテゴリの記事
- 高周波・RFニュース 2025年6月12日 QorvoがKバンド衛星通信向けパワーアンプ発表、QuectelがWi-SUNモジュール発表、Kymetaがマルチバンドアンテナ発表、iFixitがAppleがiPadの修理は悪夢と認めたと語る、中国の6G戦略解説、iFixitがNintendo Switch2のチップ同定(2025.06.12)
- 高周波・RFニュース 2025年6月11日 Siversの28GHzモジュールがaiRadarに採用、LG Innotekが車載5G衛星通信モジュール発表、Nuvotronicsがフィルタ設計のためのStrataWorksプラットフォーム発表、TechInsigtsのXiaomi 15S Pro分解、QualcommがAlphawave Semi買収(2025.06.11)
「日記・コラム・つぶやき」カテゴリの記事
- 高周波・RFニュース 2025年6月12日 QorvoがKバンド衛星通信向けパワーアンプ発表、QuectelがWi-SUNモジュール発表、Kymetaがマルチバンドアンテナ発表、iFixitがAppleがiPadの修理は悪夢と認めたと語る、中国の6G戦略解説、iFixitがNintendo Switch2のチップ同定(2025.06.12)
- 高周波・RFニュース 2025年6月11日 Siversの28GHzモジュールがaiRadarに採用、LG Innotekが車載5G衛星通信モジュール発表、Nuvotronicsがフィルタ設計のためのStrataWorksプラットフォーム発表、TechInsigtsのXiaomi 15S Pro分解、QualcommがAlphawave Semi買収(2025.06.11)
« 河童ラーメン本舗で河童ラーメン煮卵付きをいただく。私はここのラーメン一番好き。替え玉も無料。 | トップページ | 平均が0でなく相関がある正規分布の比(商)Z=X/Yの分布 Fieller-Hinkleyの式（平均が0ならコーシー分布）をカシオの高精度計算サイトkeisan.casio.jpにあげていたが式にミスがあったのでやり直し。Pythonでモンテカルロして全然違う値になったのでもっと新しい論文探した。 »
コメント