時代に逆行して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でモンテカルロして全然違う値になったのでもっと新しい論文探した。 »
「パソコン・インターネット」カテゴリの記事
「学問・資格」カテゴリの記事
- Interface2025年1月号はMATLABで1ニューロンから手作り 数学＆図解でディープ・ラーニング。初歩からAlexNetの転移学習、CNNまで話題が豊富で、なんとMatlab Onlineの半年ライセンスがついてくる。Simulinkや各種toolboxも使える。早速MATLAB入門オンラインコース修了した。(2024.12.03)
- 高周波・RFニュース 2024年11月29日 NuvotronicsのPolystrataがQorvoのGaNパワーアンプに採用、WürthがBalunのラインアップ拡張、MinewSemiのモジュールにNordicのnRF54Lが採用、AIチップにPCIe 7.0が必要な理由、SIJの特集はThe Road from 1 Gbps-NRZ to 224 Gbps-PAM4(2024.11.29)
- 高周波・RFニュース 2024年11月27日 Ericsson Mobility Report Nov.2024発行、NordicのnRF54L15がRaytacモジュールに搭載、SIMComの5G RedCapモジュール、北京の清華大学でLevel4の自律ネットワークが稼働、AmpleonのDoherty RF Power Transistor、Rohmのシリコンキャパシタ(2024.11.27)
「日記・コラム・つぶやき」カテゴリの記事
- Interface2025年1月号はMATLABで1ニューロンから手作り 数学＆図解でディープ・ラーニング。初歩からAlexNetの転移学習、CNNまで話題が豊富で、なんとMatlab Onlineの半年ライセンスがついてくる。Simulinkや各種toolboxも使える。早速MATLAB入門オンラインコース修了した。(2024.12.03)
- 高周波・RFニュース 2024年11月29日 NuvotronicsのPolystrataがQorvoのGaNパワーアンプに採用、WürthがBalunのラインアップ拡張、MinewSemiのモジュールにNordicのnRF54Lが採用、AIチップにPCIe 7.0が必要な理由、SIJの特集はThe Road from 1 Gbps-NRZ to 224 Gbps-PAM4(2024.11.29)
- 高周波・RFニュース 2024年11月27日 Ericsson Mobility Report Nov.2024発行、NordicのnRF54L15がRaytacモジュールに搭載、SIMComの5G RedCapモジュール、北京の清華大学でLevel4の自律ネットワークが稼働、AmpleonのDoherty RF Power Transistor、Rohmのシリコンキャパシタ(2024.11.27)
« 河童ラーメン本舗で河童ラーメン煮卵付きをいただく。私はここのラーメン一番好き。替え玉も無料。 | トップページ | 平均が0でなく相関がある正規分布の比(商)Z=X/Yの分布 Fieller-Hinkleyの式（平均が0ならコーシー分布）をカシオの高精度計算サイトkeisan.casio.jpにあげていたが式にミスがあったのでやり直し。Pythonでモンテカルロして全然違う値になったのでもっと新しい論文探した。 »
コメント