Monte Carlo Method

mc.c

#include <stdio.h>
#include <stdlib.h>
#include <time.h>

int main(void)
{
	double area, pi, x, y, n, N, i;
	
	n = 0;
	N = 100000000;
	
	srand((unsigned)time(NULL));
	for(i=0;i<N;i++)
	{
		x = rand() / (double)RAND_MAX;
		y = rand() / (double)RAND_MAX;
		
		if((x*x + y*y) <= 1)
		{
			n++;
		}
	}
	
	area = 1 * n / N;
	pi = area * 4;
	
	printf("%f\n",pi);
	
}