« 御影近くのNTTのアンテナ | トップページ | 「アトポス」を読んだ。 »

2010年6月21日 (月)

超高精度でロジスティック写像を計算(PARI/GP)

カオスを生じる写像として非常によく知られているロジスティック写像

xn+1 = a xn ( 1- xn)

ですが、計算の桁数が少ないとすぐ信用されない結果しか出さないのでも有名。

で、PARI/GPで100000桁!!!精度で計算して

Excelと50桁計算できるサイトkeisan.casio.jpと比較してみた。

a=3.9, 初期値は0.1としておく。

Log

誤差の結果がこれ。

Excelはもはや60くらいの繰り返しでアウト(倍精度実数のCでも同じことっす)。keisan.casio.jpは善戦するも。200超えたくらいでアウト。

ということで、こういうのの計算にPARI/GPはぴったりです。

« 御影近くのNTTのアンテナ | トップページ | 「アトポス」を読んだ。 »

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

コメント

コメントを書く

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

トラックバック


この記事へのトラックバック一覧です: 超高精度でロジスティック写像を計算(PARI/GP):

« 御影近くのNTTのアンテナ | トップページ | 「アトポス」を読んだ。 »

最近の記事

最近のコメント

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