fortran66のブログ

fortran について書きます。

2011-11-22から1日間の記事一覧

Monte Carlo で円周率

半径1の1/4円中に入るか入らないかで円周率を求める古典的な例です。一千万回の試行をしているけれども、組み込み乱数がちゃんとしているかは不明w 試行回数Nに対して、おおむね 1/sqrt(N) に比例して誤差が減るはずなので目安の数字も出力しています。…