蒙地卡罗法求PI

 解法蒙地卡罗的解法适用于与面积有关的题目，例如求PI值或椭圆面积，这边介绍如何求PI
值；假设有一个圆半径为1，所以圆面积就为PI，而包括此四分之一圆的正方形面积就
为1，如下图所示：
如果随意的在正方形中投射飞标（点）好了，则这些飞标（点）有些会落于四分之一圆内，假
设所投射的飞标（点）有n点，在圆内的飞标（点）有c点，则依比例来算，就会得到上图中最
后的公式。
至于如何判断所产生的点落于圆内，很简单，令乱数产生X与Y两个数值，如果X^2+Y^2不大于1
就是落在圆内。

#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#define N 9999999
int main(void)
{
 int i, sum = 0;
 double x, y;
 srand(time(NULL)+rand());
 for(i = 0; i < N; i++)
 {
  srand(time(NULL)+rand()+i);
  x = (double) rand() / RAND_MAX;
  y = (double) rand() / RAND_MAX;
  if((x * x + y * y) <= 1)
   sum++;
 }
 printf("PI = %f\n", (double) 4 * sum / N);
 getchar();
 return 0;
}