プログラミング

・円周率を求めるための当たり外れの試行を、10000回行うことによって、円周率の値の1回分を求める。
・円周率を100回求めて、その平均値を答えとして表示

これをモンテカルロ法によって求めるプログラムを教えてください。
説明も加えてください。お願い致します！

No.1 回答者： yhr2 回答日時： 2021/07/18 19:50

＞これをモンテカルロ法によって求めるプログラムを教えてください。

そのアルゴリズムを考えるところからがプログラミングなんじゃない？

No.2 回答者： satoumasaru 回答日時： 2021/07/18 19:52

ググってみましょう。

モンテカルロ法で円周率を求める
http://www.cmpt.phys.tohoku.ac.jp/open-campus/20 …

モンテカルロ法を用いて円周率を求めるプログラム
https://www.rs.kagu.tus.ac.jp/yamalab/2010/haman …

No.3 回答者： ありものがたり 回答日時： 2021/07/19 12:47

コードをここに書けというのなら、

使用言語を指定しなきゃ話になるまい。

No.4 回答者： kamiyasiro 回答日時： 2021/07/22 10:24

企業で統計を推進する立場の者です。

既に回答は#2さんによって示されていますが、オーソドックスなアルゴリズムを説明し、Rで書いてみました。それをご使用中の言語に直して下されば幸甚です。

円周率をπとする。
x軸y軸とも±ｒとなる正方形を考える。その面積は4r^2
半径ｒの円を考える。その面積はπr^2
πは面積比で与えられる。すなわち、π＝４・(円の面積)／(正方形の面積)
この面積比を、±ｒの範囲に打点された一様乱数の座標数で代用する。
（簡単のためにｒ＝１とする）

sum((x1^2 + x2^2) <= 1) は、本来であれば、x座標^2＋y座標^2の平方根が半径ですから√が必要ですが、半径を１とすることによって計算負荷を下げています。学校の課題であれば、
sum(sqrt(x1^2 + x2^2) <= 1)
として下さいね。

～～～～～～～～～～～～～～～～～～～～～

pai <- NULL
i <- 0

while(i < 100){

x1 <- runif(10000, -1, 1)
x2 <- runif(10000, -1, 1)

p <- 4 * sum((x1^2 + x2^2) <= 1) / 10000

pai <- append(pai, p)

i <- i + 1

}

mean(pai)