7-50 近似求PI (15分)

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

​π / 2​​=1+​1! / 3​​+​2! / 3×5​​​​+3! / ​3×5×7​​​​+⋯+​​​i! / 3×5×⋯×(2×i+1)​​+⋯

输入格式：

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

scanf("%le", &eps);

 

输出格式：

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

PI = 近似值

 

输入样例：

1E-5

 

输出样例：

PI = 3.14158

 

作者

张高燕

单位

浙江大学城市学院

代码长度限制

16 KB

时间限制

400 ms

内存限制

64 MB

---

编译器 (1)

C (gcc)

#include <stdio.h>
double fact(int n);
int main(void)	//如果不需要从命令行中获取参数，请用int main(void)；否则请用int main(int argc, char *argv[])
{
	double eps;
	scanf("%le", &eps);

	int i;
	double sum=0.0,sub=1.0;
	for(i=0;;i++){
		sub*=2*i+1;
		sum+=fact(i)/sub;
		if(fact(i)/sub<eps){
			printf("PI = %.5f",sum*2);	// \n
			break;
		}
	}
	
	return 0;
}

double fact(int n)
{
	int i;
	double factorial=1.0;
	if(n==0){	//为什么此处的if不能省略? 
		return 1.0;
	}
	for(i=2;i<=n;i++){
		factorial*=i;
	}
	return factorial;
}