用格里高利公式求给定精度的PI值

本题要求编写程序,计算序列部分和 4∗(1−1/3+1/5−1/7+…) ,直到最后一项的绝对值小于给定精度eps。

输入格式:

输入在一行中给出一个正实数eps。

输出格式:

在一行中按照“Pi = pi”的格式输出部分和的值pi,精确到小数点后四位。题目保证输入数据和计算结果均不超过双精度范围。

输入样例1:

1
0.0001

输出样例1:

1
Pi = 3.1418

输入样例2:

1
0.2

输出样例2:

1
Pi = 2.8952

分析:循环计算一个绝对值,直到在某个位置其分子分母小于eps,退出算式,fabs作为双浮点函数的绝对值计算。另外abs是作为int函数的绝对值计算,前后只差个f。

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
#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;
}