7-1 用格里高利公式求给定精度的PI值 (15分)
本题要求编写程序,计算序列部分和 4∗(1−1/3+1/5−1/7+…) ,直到最后一项的绝对值小于给定精度eps。
输入格式:
输入在一行中给出一个正实数eps。
输出格式:
在一行中按照“Pi = pi”的格式输出部分和的值pi,精确到小数点后四位。题目保证输入数据和计算结果均不超过双精度范围。
输入样例1:
0.0001
输出样例1:
Pi = 3.1418
输入样例2:
0.2
输出样例2:
Pi = 2.8952
#include<stdio.h>
#include<math.h>
int main(void)
{
int flag = 1,fenmu = 1;
double eps,item = 1.0,sum = 0;
scanf("%lf",&eps);
if(eps>1){ //精度>1,输出首项
printf("Pi = 4.0000");
}else{
while(fabs(item)>=eps) //调运绝对值函数fabs()
{
item = flag*1.0/fenmu;
sum += item;
flag = -flag;
fenmu += 2;
}
printf("Pi = %.4lf\n",sum*4);
}
return 0;
}
//ps1: eps=1时,item(首项)>=eps=1,满足循环条件进行一次循环 -> item = 1/3,循环结束
//ps2: 最后一项的绝对值<给定精度 <==> 最后一项的绝对值>=给定精度eps