トップ «前の日記(2008-10-04) 最新 次の日記(2010-01-01)» 編集

にっき

2005|07|08|09|10|11|12|
2006|01|02|05|06|07|10|11|12|
2007|01|06|07|12|
2008|01|03|04|10|
2010|01|03|05|10|11|12|
2011|06|07|09|10|
2012|03|04|05|
2013|03|05|08|
2014|10|11|
2015|07|
2016|04|06|07|10|12|
2017|05|06|08|09|

2008-10-07

_ ローレンツ方程式を数値計算してみた

手慰みにローレンツ方程式をC言語で計算してみました。割とキレイな画になったのでうpしてみます。

画像の説明

_ 数値計算には4次ルンゲクッタを用い10000回計算しています。尚、同様のことをschemeでもやってみました。

画像の説明

_ なぜここまで差がでるのだろう?

_ 両者共に、初期値もルンゲクッタのスケールhも同じだけど、C言語はdouble小数点以下6桁(double)なのに対して、schemeは少数点以下16桁まである数字があるため10000行×3個の数字の8行目からC言語での結果と差が生まれている(少数第7位の数字を四捨五入した場合)。その結果、図のような違いが生じたようだ。

[]


あなたは本日人目のお客様です!
総計: