乱数

質問

乱数ｘ、ｙ(0〈ｘ〈1)、(0〈ｙ〈1)をｎ回発生させて、これらを座標とする点をｘ‐ｙ座標平面上にプロットしてこの点と原点との距離が1以下となるときの回数をカウントしｍ回としてn:m＝（長さ1の正方形の面積）：（円の4分の1の面積）となる。これを利用して円周率を求めるプログラムなんてどう作るんですか？

回答No.2

imogasi

モンテカルロ法というやつでしょう。
１．第１系列の乱数発生させる。
そのままで０．Ｘ・・の大きさの数のはず。
２．第2系列の乱数発生させる。
大きさは０．ｘ・・。
以上でｐ＝（ｘ、ｙ）という一点が決まる。これは
第１象限の辺１の１ｘ１の正方形の中の１点です。
３．てんｐと原点との距離を計算する
　　ｄ＝ＳＱＲ（Ｘ＾２＋ｙ＾２）
またはｄ＝ＳＱＲ（Ｘ＊Ｘ＋Ｙ＊Ｙ）
４．もしｄ＜＝１であれば円内の点としてカウント
５．円内の点の数／試行回数が３や３．１、３．１４付近になれば、この試行は納得というところ。

回答No.1

BLUEPIXY

質問文に、ほぼ書かれているとおりです。
距離を求めるには
sqr(x*x+y*y)
とします。