实验6-1 近似求PI_c语言_PT

本题要求编写程序，根据下式求π的近似值，直到最后一项小于给定精度eps。

输入格式：

输入在一行中给出精度eps，可以使用以下语句来读输入：

scanf("%le", &eps);

输出格式：

在一行内，按照以下格式输出π的近似值（保留小数点后5位）：

PI = 近似值

输入样例：

1E-5

输出样例：

PI = 3.14158

思路：这题真的有点折磨，最大的坑在于要求最后一项是要小于eps的，不论倒数第二项是大于还是等于eps。发现了这个点就不难了。 

#include <stdio.h>

double fact(double n){//阶乘函数
	
	if(n==0||n==1)
		return 1;
	else
		return n*fact(n-1);
}


int main(int argc, char *argv[]) {
	
	double eps;
	scanf("%le", &eps);
	double sum=0;
	double fenzi=1,fenmu=3;
	int i=1;
	
	while(fenzi/fenmu>=eps)
	{
		sum=sum+fenzi/fenmu;
		i++;
		fenzi=fact(i);
		fenmu=1.0*fenmu*(2*i+1);
		
	}
	if(eps<=1)
	sum=sum+fenzi/fenmu;//因为要求最后一项小于eps，但是由于不符合while循环，所以这里要补加上。
//但是如果eps非常大，那么只保留第一项1就可以，与分子分母无关，while循环不会执行，所以不需要修正sum的大小。
	
	sum=sum+1;
	printf("PI = %.5lf\n",2*sum);

}