7-15 计算圆周率 (15分)
根据下面关系式,求圆周率的值,直到最后一项的值小于给定阈值。
输入格式:
输入在一行中给出小于1的阈值。
输出格式:
在一行中输出满足阈值条件的近似圆周率,输出到小数点后6位。
输入样例:
0.01
输出样例:
3.132157
#include <stdio.h>
int main()
{
double frac = 1, tip = 1, tips = 1,sum = 0, n = 1, a = 1;
int cnt = 0;
double x;
scanf("%lf", &x);
while(frac >= x) //frac表示分式
{
if(cnt)
tip *= n++; //n的阶乘
else //cnt为0,即次数为0时
cnt++;
tips *= a; //等差数列2n+1的和
a += 2;
frac = tip / tips;
sum += frac;
}
printf("%.6lf\n", sum * 2.0);
return 0;
}
n的阶乘和等差数列的变化,条件控制。