7-15 计算圆周率 (15 分)
根据下面关系式,求圆周率的值,直到最后一项的值小于给定阈值。
2π=1+31+3×52!+3×5×73!+⋯+3×5×7×⋯×(2n+1)n!+⋯
输入格式:
输入在一行中给出小于1的阈值。
输出格式:
在一行中输出满足阈值条件的近似圆周率,输出到小数点后6位。
输入样例:
0.01
输出样例:
3.132157
/*计算圆周率*/
#include <iostream>
#include <cmath>
using namespace std;
int main (void)
{
double m, n;
double end, pi;
int i = 1;
m = 1.0f;
n = 1.0f;
pi = 1.0f;
scanf("%lf", &end);
do{
m *= i;
n *= (1.0 + 2.0 * i);
i++;
pi += (m/n);
}while(m/n >= end);
printf("%.6lf", pi*2.0);
return 0;
}