몬테카를로 방법으로 원주율 구하기

몬테카를로 방법으로 원주율 구하기

몬테카를로 방법 원주율 계산

1 방법

• 예제(500회 시행): http://zetawiki.com/ex/php/Monte-Carlo-Pi.php
• 많이 시행할수록 실제 값에 근접할 확률이 높아짐

<?php
$inner_count = 0;
for($i=1; $i<=500; $i++) {
	$x = 200 * mt_rand() / mt_getrandmax() - 100;
	$y = 200 * mt_rand() / mt_getrandmax() - 100;
	$r = sqrt( $x*$x + $y*$y );
	if( $r <= 100 ) $inner_count++;
	$pi = 4 * $inner_count / $i;
	echo "#$i ($x, $y) r=$r, [$inner_count/$i] pi=$pi<br>";
}

실행예시

... (생략)
#99997 (47.072822855354, -8.7733809411402) r=47.883429959746, [78581/99997] pi=3.143334300029
#99998 (-10.641246713065, 0.96602174498422) r=10.685004895652, [78582/99998] pi=3.1433428668573
#99999 (-50.901550031687, -29.689785339725) r=58.927507576235, [78583/99999] pi=3.1433514335143
#100000 (-76.904958196406, 65.226926265856) r=100.84108540313, [78583/100000] pi=3.14332

2 같이 보기

• 몬테카를로 방법
• 원주율

3 참고

• http://navercast.naver.com/contents.nhn?rid=2871&contents_id=90216