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辛大盛）を食す。特製スパイスと赤玉ねぎのピクルスがついてくる。 »
「学問・資格」カテゴリの記事
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2021年2月16日 IEEE Microwave Magazineの特集はオールデジタルのRFID、Microwave JournalはＥバンド ミリ波通信に衛星や気球を使う話、アメリカの半導体企業がバイデンに投資を迫る、(2021.02.17)
- カオスを生じる電気回路、Chua’s circuitをLTspiceで回路シミュレーションしてみる。(2021.02.19)
- Labyrinth Chaos(迷宮カオス）を生むThomas-Rössler方程式のパラメータbを色々変えて、Python＋Scipyでルンゲクッタ8次のDOP853(Dormand＆Prince)を使って計算してGIFアニメ(2021.02.16)
- フィッツヒュー・南雲 (FitzHugh-Nagumo) 方程式をPython+Scipyでルンゲクッタ8次のDOP853(Dormand Prince)で計算。(2021.02.23)
- 「水晶振動子の等価回路計算」をカシオの高精度計算サイトkeisan.casio.jpの自作式としてUP！ インピーダンスの大きさと位相がグラフ化できる。(2021.02.12)
« 新型コロナウイルス、中国、日本、韓国、アメリカでの感染者数を指数関数とロジスティック関数でフィッティングした。（10/18更新）アメリカと日本は再び増加傾向。アメリカはトランプがあれだから、ですが日本も油断している人多すぎでは、、、 | トップページ | ココイチでスパイスカレーTHEチキベジ（3辛大盛）を食す。特製スパイスと赤玉ねぎのピクルスがついてくる。 »
コメント