实验6-1 近似求PI (15分)

实验6-1 近似求PI (15分)

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

输入格式：

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

scanf("%le", &eps);

输出格式：

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

PI = 近似值

输入样例：

1E-5

输出样例：

PI = 3.14158

代码:

#include <stdio.h>

double fact(int n);
double oddfact(int n);
	
int main()
{
    int i;
	double sum=1,eps,pi,item=1;

    scanf("%le", &eps);
    
    for(i=1;item>eps;i++)
    {
    	item=fact(i)/oddfact(i);
    	//printf("PI = %.5f",item);	
		sum=sum+item;
		//printf("PI = %.5f",sum*2);	
	}
	
	printf("PI = %.5f",sum*2);
    return 0;
}

double fact(int n)
{
	int i;
	double sum=1;
	
	for(i=2;i<=n;i++)
		sum=sum*i;
	
	return sum;
}

double oddfact(int n)
{
	int i;
	double sum=1;
	
	for(i=1;i<=n;i++)
		sum=sum*(2*i+1);
	
	return sum;	
}