求π的近似值(2)

利用随机数法求π的近似值


*问题分析与算法设计
随机数法求π的近似值的思路：在一个单位边长的正方形中，以边长为半径，以一个顶点为圆心，在政权方形上作四分之一圆。随机的向正方形内扔点，若落入四分之一圆内则计数。重复向正方形内扔足够多的点，将落入四分之一圆内的计数除以总的点数，其值就是π值四分之一的近似值。
按此方法可直接进行编程，注意：本方法求出的π值只有统计次数足够多时才可能准确。

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

int main()
{
    double e = 0.1,b = 0.5, c, d;
    long int i; 
    for (i = 6; ; i *= 2) 
    {
        d = 1.0 - sqrt(1.0 - b*b); 
        b = 0.5*sqrt(b*b + d*d);
        if (2*i*b - i*e < 1e-15) 
        {
            break;
        } 
        e = b; 
    }
    printf("pai = %.15lf\n", 2*i*b); 
    printf("The number of edges of required polygon:%ld\n",i);

    return 0;
}

3.122267
3.139733
3.133733