割圆术计算圆周率

include

include

using namespace std;
int main()
{
int n;
double h,x=1,s=3*sqrt(3)/2;//从六边形开始计算
for(n=6;n<=5999;n+=n)//运用n边形与2n边形之间的关系求解
{
h=sqrt(1-x*x/4);
s+=n*x*(1-h)/2;
x=sqrt((x*x)/4+(1-h)*(1-h));
}
cout<<”pi=”<