.求π的近似值

请利用“正多边形逼近”的方法求出π的近似值


*问题分析与算法设计
利用“正多边形逼近”的方法求出π值在很早以前就存在，我们的先人祖冲之就是用这种方法在世界上第一个得到精确度达小数点后第6位的π值的。
利用圆内接正六边形边长等于半径的特点将边数翻番，作出正十二边形，求出边长，重复这一过程，就可获得所需精度的π的近似值。
假设单位圆内接多边形的边长为2b，边数为i，则边数加倍后新的正多边形的边长为：
x=√──────
2-2*√───
1-b*b
──────
2
周长为： 
y=2 * i * x i:为加倍前的正多边形的边数

#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;
}

pai=3.141592653589794
The number of edges of required polygon:100663296