« 任意精度実数計算用にPARI/GPを導入してみる。 | トップページ | ルンゲクッタ法4次、8次、14次の比較 »

2010年5月24日 (月)

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桁で計算してますよ。

「Lorenz.gp」をダウンロード

使い方は、PARI/GPを立ち上げて、

\r Lorenz

でロード。

rk14("ファイル名")

で実行されて、ファイル名のところに保存される。散布図がどうも描きにくそうなので、ファイルに落としてExcelで描いてみた。それがこれ。

Lorenz14thparigp

すばらしく安定感のある計算ができるなあ。これからはこれでしばらくいろいろ計算していってみよう(続く)。

« 任意精度実数計算用にPARI/GPを導入してみる。 | トップページ | ルンゲクッタ法4次、8次、14次の比較 »

学問・資格」カテゴリの記事

コメント

コメントを書く

(ウェブ上には掲載しません)

トラックバック


この記事へのトラックバック一覧です: 35段14次のルンゲクッタ法をPARI/GPで計算(Lorenz方程式):

« 任意精度実数計算用にPARI/GPを導入してみる。 | トップページ | ルンゲクッタ法4次、8次、14次の比較 »

最近の記事

2021年3月
  1 2 3 4 5 6
7 8 9 10 11 12 13
14 15 16 17 18 19 20
21 22 23 24 25 26 27
28 29 30 31      
フォト
無料ブログはココログ