乱数を使用した積分

/****************************/
/* 乱数を使用した積分       */
/*      coded by Y.Suganuma */
/****************************/
#include <stdio.h>
#include <stdlib.h>

main()
{
	double x, y;
	int i1, count = 0, n = 1000000;

	srand(123);   // 乱数の初期設定

	for (i1 = 0; i1 < n; i1++) {
		x = (double)rand() / RAND_MAX;
		y = (double)rand() / RAND_MAX;
		if (y <= x*x)
			count++;
	}

	printf("面積： %f\n", (double)count/n);

	return 0;
}