东软学习小组成员:时雾
用格里高利公式求给定精度的PI值
本题要求编写程序,计算序列部分和 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
分析:循环计算一个绝对值,直到在某个位置其分子分母小于eps,退出算式,fabs作为双浮点函数的绝对值计算。另外abs是作为int函数的绝对值计算,前后只差个f。
#include<stdio.h>
#include<math.h>
int main(){
double eps=0,Pi=0;
scanf("%lf",&eps);
int fm=1,zf=1;//fm为分母缩写,zf为正负缩写
Pi=Pi+zf*1.0/fm;
while(fabs(zf*1.0/fm)>=eps){
fm+=2;
zf=-zf;
Pi=Pi+zf*1.0/fm;
}
printf("Pi = %.4lf",4*Pi);
return 0;
}