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ニュース 2025年7月17日 IEEE Microwave Magazineで地上監視レーダ解説など、Microwave Journalでミリ波AiPのOTA測定解説など、5G Americasが5G Advancedホワイトペーパー発行、TDKのAI 微細欠陥検出edgeRX Vision、IEEE Trans. on THz の今月号発行(2025.07.17)
- Xのタイムラインが何かランダムウォークのガウス分布じゃなくて逆正弦法則のように見えてきている…仕方ないのでJulia言語でモンテカルロシミュレーションで試行回数を変えて逆正弦法則とガウス分布を見る。(2025.07.16)
- 高周波・RFニュース 2025年7月13日 Pythonの高周波ライブラリscikit-rfがv1.8.0に、SamsungがGalaxy Z Fold7など発表→QualcommがSnapdragon 8 Eliteが使われていると発表、NGMNが基地局アンテナの推奨事項をまとめる、STMicroとMetalenzがメタサーフェス光学のライセンス締結(2025.07.14)
- Interface2025年8月号Pythonで体験！はじめての暗号を買った。上杉暗号からRSA、AES、DHなど、特に楕円曲線暗号についてはコードも実際に動かすところまで詳しくかかれていた。耐量子暗号や聞いたことなかったY-00暗号や関数型暗号も記載。(2025.07.10)
- 高周波・RFニュース 2025年7月8日 NordicとSercommのセルラーIoTモジュール、iFixitがFairphone 6を分解、スコアは10／10、RCR wireless newsのウェビナー2件（6GとIndustry4.0)、SEMCOが高耐圧C0G MLCCを車載急速充電に提案(2025.07.09)
« 新型コロナウイルス、中国、日本、韓国、アメリカでの感染者数を指数関数とロジスティック関数でフィッティングした。（10/18更新）アメリカと日本は再び増加傾向。アメリカはトランプがあれだから、ですが日本も油断している人多すぎでは、、、 | トップページ | ココイチでスパイスカレーTHEチキベジ（3辛大盛）を食す。特製スパイスと赤玉ねぎのピクルスがついてくる。 »
コメント