モンテカルロ

http://ogata.plala.jp/main/hatena_pictures/pi_monte_1000s.gif
 明日のレポート、モンテカルロ法で円周率を求める? らしい。
なんか数年前にも同じようなことをやった気がするんだが(笑
とりあえず3分ハッキングで完成。
あほらしいので、gnuplotにpipeしてグラフも表示するようにしてみた。
20分ハッキングですね。(gnuplotの使い方がいまいち分かってなかったので)
あとはせこせこと出力を見やすく調整してみたり とか。
まぁ相変わらず適当なソースだけど晒しておきます。
好きに使ってください。(たぶん自分で一から書いた方が早いと思うけど…)
http://ogata.plala.jp/main/public_cs-kobe/pi_monte_open.c