/* モンテカルロ法で円の面積を求める */
#include <stdio.h>
#include<stdlib.h>					// randを使うときのヘッダ
#define N 400						// 繰り返し回数
#define D 10.0						// 正方形の辺の長さ

main()
{
	unsigned long i;
	unsigned long in=0,out=0;
	float ratio,x,y,r,square;

	square=D*D;						// 正方形の面積
	r=D/2.0;						// 正方形に内接する円の半径

	for(i=0;i<N;++i) {
		x= (double)rand()/(double) RAND_MAX *r;		// x座標を一様乱数で求める
		y= (double)rand()/(double) RAND_MAX *r;		// y座標を一様乱数で求める
		if( (x*x+y*y)<=r*r ) ++in;		// (x,y)が円の内側か外側かを判定 ←判定条件
		else ++out;
		printf("回数=%d 内側の回数=%d 外側の回数=%d x=%5.2f y=%5.2f\n",i,in,out,x,y);
	}
	ratio=(float)in/(float)(in+out);				// (x,y）が円の内側のときと外側のときの比

	printf("範囲内の点数= %d   範囲外の点数= %d\n",in,out);
	printf("正方形の面積＝%f\n",square);
	printf("円の面積（＝正方形の面積×ｒの内外の比）＝%f\n",square*ratio);
	printf("円周率を使って求めた円の面積（3.14*r*r）＝%f\n",3.14*r*r);
}