Python＋Scipyでルンゲクッタ8次のDOP853(Dormand＆Prince)を使う（その12） 誤差設定を1e-12にしてローレンツ方程式でちょっとだけ初期値が違うものを並べてGIFアニメにする。
DOP853での設定でatol, rtol（絶対誤差、相対誤差）があるがどちらも1E-12くらいにしておくと、ピタゴラスの三体問題なんかが簡単に計算できる。今回はローレンツ方程式で初期値をx0=0.1, 0.1001, 0.100001にしてみた（y0=z0=0.1は共通）。
計算結果はこちら。クリックでGIFアニメが始まります。
« 新型コロナウイルス、中国、日本、韓国、アメリカでの感染者数を指数関数とロジスティック関数でフィッティングした。（10/18更新）アメリカと日本は再び増加傾向。アメリカはトランプがあれだから、ですが日本も油断している人多すぎでは、、、 | トップページ | ココイチでスパイスカレーTHEチキベジ（3辛大盛）を食す。特製スパイスと赤玉ねぎのピクルスがついてくる。 »
「学問・資格」カテゴリの記事
- TensorFlow(Keras)でモデルをsaveで保存してload_modelで読み込むときに突然エラー(utf8で読めないとかディレクトリが存在しないとか）が出始めた。なんで？といろいろやると、単に日本語がフォルダ名に使われているときにエラーになるだけだった…（Windowsネイティブの場合）(2023.05.28)
- 時代に逆行してCOBOL(GnuCOBOL)を学んでみる(4) 10進31桁まで計算できることを生かして4段4次のルンゲクッタ法でローレンツ方程式を計算して図示してみる。テキストベース（アスキーアート）で！(2023.05.23)
- 時代に逆行してCOBOL(GnuCOBOL)を学んでみる(3) ロジスティック写像の分岐図を31桁まで10進計算ができることを生かして描いてみる。テキストベースで！今回は2次元配列の練習。(2023.05.19)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース(5/18) IEEE Microwave Magazineはグラフェンのマイクロ波応用、Microwave Journalはテラヘルツ波の発生、Next G allianceの新6Gレポート、Hexa-Xの新動画、ESAが軌道上に6Gラボ、De-embeddingって何？など。(2023.05.18)
- カシオの高精度計算サイトkeisan.casio.jpが復旧してログインや自作式の作成ができるようになった。ということで作った自作式のアクセスベスト10とお勧めを紹介。(2023.05.17)
« 新型コロナウイルス、中国、日本、韓国、アメリカでの感染者数を指数関数とロジスティック関数でフィッティングした。（10/18更新）アメリカと日本は再び増加傾向。アメリカはトランプがあれだから、ですが日本も油断している人多すぎでは、、、 | トップページ | ココイチでスパイスカレーTHEチキベジ（3辛大盛）を食す。特製スパイスと赤玉ねぎのピクルスがついてくる。 »
コメント