時代に逆行して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でモンテカルロして全然違う値になったのでもっと新しい論文探した。 »
「パソコン・インターネット」カテゴリの記事
- Julia言語でタッパーの自己言及式（不等式を計算して図示するとまた不等式になる）を描いてみる。543桁の数を含む計算が必要だが、デフォルトで任意精度演算が可能なので容易にできた。(2025.03.11)
- UnityでVisual C#用の常微分方程式ソルバーOpen Solving Library for ODEs(OSLO)を使う(3)三体問題の周期解としてまずは有名な8の字を描くものをやってみる。(2025.03.14)
- UnityでVisual C#用の常微分方程式ソルバーOpen Solving Library for ODEs(OSLO)を使う(2)ブルースカイ・カタストロフィを生じるGavrilov Shilnikov modelを計算してDormand＆PrinceのRK547Mで計算して玉を動かして軌跡を付ける。ぐるぐる回っていたと思ったら突然広がって戻る。(2025.03.12)
「学問・資格」カテゴリの記事
- 高周波・RFニュース 2025年3月14日 Microwave Journalでテラヘルツデバイス製造解説、IgnionのRF、アンテナ向けAIツールOxionが2.0に、TIが1.38mm²の超小型MCU発表、IntelのCEOはLip-Bu Tanに。(2025.03.14)
- 高周波・RFニュース 2025年3月13日 NordicとQorvoがAliroとMatterのリファレンスアプリケーション提供、TSMCとMediaTekがパワーアンプと電源管理ユニット統合、3GPPのTSG RANの議長がSamsungの人に、Silicon Labsが超小型Bluetooth Soc, QuectelがWi-Fi/BTモジュール発表(2025.03.13)
- 高周波・RFニュース 2025年3月12日 iFixitが任天堂Alarmoを分解、なんとSocionextの24GHzミリ波センサ’積んでる！Next G Allianceが6Gに向けたデジタルツインとFWAのホワイトペーパー発行、Qorvoが統合型UWB SoC発表、SemtechがLora、u-bloxとTelitがGNSSモジュール発表(2025.03.12)
「日記・コラム・つぶやき」カテゴリの記事
- 高周波・RFニュース 2025年3月14日 Microwave Journalでテラヘルツデバイス製造解説、IgnionのRF、アンテナ向けAIツールOxionが2.0に、TIが1.38mm²の超小型MCU発表、IntelのCEOはLip-Bu Tanに。(2025.03.14)
- 高周波・RFニュース 2025年3月13日 NordicとQorvoがAliroとMatterのリファレンスアプリケーション提供、TSMCとMediaTekがパワーアンプと電源管理ユニット統合、3GPPのTSG RANの議長がSamsungの人に、Silicon Labsが超小型Bluetooth Soc, QuectelがWi-Fi/BTモジュール発表(2025.03.13)
- 高周波・RFニュース 2025年3月12日 iFixitが任天堂Alarmoを分解、なんとSocionextの24GHzミリ波センサ’積んでる！Next G Allianceが6Gに向けたデジタルツインとFWAのホワイトペーパー発行、Qorvoが統合型UWB SoC発表、SemtechがLora、u-bloxとTelitがGNSSモジュール発表(2025.03.12)
« 河童ラーメン本舗で河童ラーメン煮卵付きをいただく。私はここのラーメン一番好き。替え玉も無料。 | トップページ | 平均が0でなく相関がある正規分布の比(商)Z=X/Yの分布 Fieller-Hinkleyの式（平均が0ならコーシー分布）をカシオの高精度計算サイトkeisan.casio.jpにあげていたが式にミスがあったのでやり直し。Pythonでモンテカルロして全然違う値になったのでもっと新しい論文探した。 »
コメント