求π的近似值

int main()
{
    int i,n,sum=0; 
    double x,y;
    printf("输入点的数量:");
    scanf("%d",&n); 
    srand(time(NULL));    
    //点落在阴影里面的次数除了总次数，得出概率，就为阴影部分的面积。
    for(i=1;i<n;i++)
    {
        x=(double)rand()/RAND_MAX; //产生0~1之间的一个随机数 
        y=(double)rand()/RAND_MAX;//产生0~1之间的一个随机数 
        if((x*x+y*y)<=1) //若在阴影区域 
            sum++; //计数 
    }
    printf("PI=%f\n",(double)4*sum/n); //输出结果 
    getch();
    return 0; 
}