時代に逆行して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ニュース 2026年7月22日 IEEE Journal of Microwaves7月号発行、メクテックが薄型フレキシブル電波吸収体を開発、ミニサーキットが33GHz信号発生器発表、Redmi Note 17 Pro分解動画など(2026.07.22)
- Qwen 3.8 Max PreviewでPythonのscikit-rfを使ってTouchstoneフォーマットのSパラメータファイルを読んでdB, 位相, スミスチャート, TDRを表示するGUIアプリを作成。何度か修正お願いするとものすごく長い思考とものすごい数のウェブ検索をして完成した。(2026.07.22)
- RF Weekly Digest (Codexで作ったOllama経由のgemma4-26B-A4B QAT、DuckDuckGo Searchを使った高周波・RF情報の週刊まとめアプリ）2026/7/12-7/19(2026.07.19)
「日記・コラム・つぶやき」カテゴリの記事
- 高周波・RFニュース 2026年7月22日 IEEE Journal of Microwaves7月号発行、メクテックが薄型フレキシブル電波吸収体を開発、ミニサーキットが33GHz信号発生器発表、Redmi Note 17 Pro分解動画など(2026.07.22)
- Qwen 3.8 Max PreviewでPythonのscikit-rfを使ってTouchstoneフォーマットのSパラメータファイルを読んでdB, 位相, スミスチャート, TDRを表示するGUIアプリを作成。何度か修正お願いするとものすごく長い思考とものすごい数のウェブ検索をして完成した。(2026.07.22)
- RF Weekly Digest (Codexで作ったOllama経由のgemma4-26B-A4B QAT、DuckDuckGo Searchを使った高周波・RF情報の週刊まとめアプリ）2026/7/12-7/19(2026.07.19)
« 河童ラーメン本舗で河童ラーメン煮卵付きをいただく。私はここのラーメン一番好き。替え玉も無料。 | トップページ | 平均が0でなく相関がある正規分布の比(商)Z=X/Yの分布 Fieller-Hinkleyの式（平均が0ならコーシー分布）をカシオの高精度計算サイトkeisan.casio.jpにあげていたが式にミスがあったのでやり直し。Pythonでモンテカルロして全然違う値になったのでもっと新しい論文探した。 »
コメント