Weekend Mathematics／情報／モデル化とシミュレーション／円周率を豆まきで



円周率を豆まきで



・・・モンテカルロ法を用いてπの値を求める・・・




区間[0,1)の乱数により、ｘ座標とｙ座標を与える。 ｘ２＋ｙ２≦１ならば、円内にあると見なす。 試行回数ｍに対して、円内に入った回数ｐをカウントすることにより、 π＝ｐ×４／ｍで近似する。 （ブラウザの更新ボタンを押してください。）





Java Script によるプログラム

<Script Language="JavaScript">
<!--
m=10000
p=0
for (i=1;i<=m;i++){
	x=Math.random();
	y=Math.random();
	if(x*x+y*y<=1 ){
		p++;
        }
}
pi=p*4/m;
document.write("π=",pi);
//-->
</Script>




戻る