35段14次のルンゲクッタ法をPARI/GPで計算(Lorenz方程式)
先日導入したPARI/GPに35段14次のルンゲクッタ法を実装してみた。
係数はこちら↓60桁分あります。
http://sce.uhcl.edu/rungekutta/
まずはいつものようにローレンツ方程式をやってみよう。式はこちら。
dx/dt = -σ*(x-y)
dy/dt=-y-x*z+r*x
dz/dt=x*y-b*z
でプログラムのソースコードはこちら。60桁で計算してますよ。
使い方は、PARI/GPを立ち上げて、
\r Lorenz
でロード。
rk14("ファイル名")
で実行されて、ファイル名のところに保存される。散布図がどうも描きにくそうなので、ファイルに落としてExcelで描いてみた。それがこれ。
すばらしく安定感のある計算ができるなあ。これからはこれでしばらくいろいろ計算していってみよう(続く)。
« 任意精度実数計算用にPARI/GPを導入してみる。 | トップページ | ルンゲクッタ法4次、8次、14次の比較 »
「学問・資格」カテゴリの記事
- Wordleを100回やってみて正解率の分布を見た。1回目で当たる確率をpとして、2回目からはk倍ずつ確率が上がっていく(kp,k²p,k³p…)とすると、私の結果はp=0.0089, k=2.4になった。(2022.05.27)
- 高周波回路シミュレータQucsStudioを使ってみる(その5)マルチポート(4ポートだとs4pみたいな)のSパラメータデータをTouchstoneファイルに出力しようとするとExport to SnPでは2ポートになってしまう。Project→Import Dataからポート数を手で入れないと。(2022.05.24)
- 高周波(RF・マイクロ波・ミリ波・5G)関連ニュース2022年5月19日 IEEE Microwave Magazineで磁性体を使わない非可逆デバイス、Microwave Journalは5Gスモールセル、THz、Movanoがミリ波血圧&血糖値計発表、三菱の3Dプリンタ衛星アンテナ、など。(2022.05.19)
- シン・ウルトラマンをIMAXで観てきた。面白かった!物理学者の有岡君がラグランジアン,AdS,プランクブレーンとかしゃべってる!Randall–Sundrumも。マグカップは標準理論+重力で、物理監修は予想通り橋本幸士さんでした!ブラックホールの構造もタイムリー!(2022.05.14)
- 個人の正解確率pとして多数決をとったときの正解確率が、pが0.5に近づくとどのくらいの人数で0.9を超えるか計算してみた。p=0.50005だと日本の人口くらいの多数決。(神とさざなみの密室読んで正当性確率が気になったので)(2022.05.10)
« 任意精度実数計算用にPARI/GPを導入してみる。 | トップページ | ルンゲクッタ法4次、8次、14次の比較 »
コメント