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

void main(void)
{
	float radius = 1.0;
	float pi = 3.1415926;
	int count=10000,inside_count = 0,inside_probability;
	int i, x, y;
	for(i=0;i<count;i++){
		x = (rand()%32768)/32768;
		y = (rand()%32768)/32768;
		if(sqrt(x*x+y*y) < radius)
			inside_count++;
	}		
	inside_probability = inside_count/count;
	printf("pi difference = %f\t<error rate>\n",
abs(inside_probability-pi);
}