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 Weekly Digest (Google AI Studio BuildによるAIで高周波・RF情報の週刊まとめアプリ） 2025/11/9-2025/11/16(2025.11.16)
- Visual Studio 2026がリリースされたので早速新しいPCにインストール。全面的にGitHub Copilotを使うようになっている。とりあえずC#でMath.NET numericsを使って連立方程式を計算するコードを書いてもらったら一発で動く。他の例として固有値や非線形計算もコードを出してくれた。(2025.11.14)
- 高周波・RFニュース 2025年11月13日 QorvoがTWTA置き換えの広帯域SSPA発表、iFixitがiPad Pro M5分解、KYOCERA AVXが0.9Vで動く超小型クロック発表、Mini-Circuitsが様々なBALUNやトランス解説、軍用5G解説、imecの110GHz可能な300mm GeSiウェハー(2025.11.13)
- 高周波・RFニュース 2025年11月12日 Qualcommが語る技術の標準化、STMicroelectronicsが語るシリコンフォトニクス、HuaweiがイノベーションとIPフォーラム開催、6G SummitでFCCのコミッショナーが米国が6Gをリードすべしと語る(2025.11.12)
« 新型コロナウイルス、中国、日本、韓国、アメリカでの感染者数を指数関数とロジスティック関数でフィッティングした。（10/18更新）アメリカと日本は再び増加傾向。アメリカはトランプがあれだから、ですが日本も油断している人多すぎでは、、、 | トップページ | ココイチでスパイスカレーTHEチキベジ（3辛大盛）を食す。特製スパイスと赤玉ねぎのピクルスがついてくる。 »
コメント