« 任意精度実数計算用に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次の比較 »

最近の記事

最近のコメント

2019年10月
    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    
フォト
無料ブログはココログ