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

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

コメント

コメントを書く

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

トラックバック

この記事のトラックバックURL:
http://app.cocolog-nifty.com/t/trackback/512682/48411280

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

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

最近のコメント

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