B. 实验6-1 近似求PI

B. 实验6-1 近似求PI
题目描述

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

​​​π/2=1+1/3+2!/(3*5)+3!/(3*5*7)+4!/(3*5*7*9)+...​​


输入

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

scanf("%le", &eps);

输出

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

PI = 近似值

样例查看模式 
正常显示
查看格式
输入样例1 
1E-5
输出样例1
PI = 3.14158

double factorial(double i)
{
	int sum = 1, j;
	for (j = 1; j <= i; j++)
	{
		sum *= j;
	}
	return sum;

}

#include<stdio.h>
int main()
{
    double eps, z, sum = 0, i = 0, one=1, a = 1, b = 1;
    scanf("%le", &eps);
    while (1)
    {
        b = 1;
        one = factorial(i) / b;
        for (z = 1; z <= a; z++)
        {
            one = one / b;
            b += 2;
        }

        i++; a++;
        if (one < eps)
            break;
        sum += one;
    }
    printf("PI = %.5lf", sum * 2);
    return 0;
}

上面这是我写的代码，这道题我写了好久一直最后的输出是3.14148  后来感觉心态都炸了

这里只能说自己学的还不够，学的还不扎实 （哭泣）

下面给出两种解决办法

首先是使用递推方法解决

#include<stdio.h>
int main()
{
	double sum=0,re=1,eps,i=1;
	scanf("%le",&eps);
	while(re>=eps)
	{
		sum+=re;
		re=re*i/(2*i+1);
		i+=1;
	}
	sum+=re;
	printf("PI = %.5lf",sum*2);
	
	return 0;
}

很简短吧！！！但是要注意最后要加一个re   因为循环里面没有加上！！！！

然后是递归方法解决

double one(double n)
{
    if(n==1)
    return 1;
    else
    {
        return (n-1)/(2*n-1)*one(n-1);
     } 
}

#include<stdio.h>
int main()
{
    double eps,i=1,re=1,sum=0;
    scanf("%le",&eps);
    while(eps<=re)
    {
        re=one(i);
        sum+=re;
        i++;
    }
    
    printf("PI = %lf",sum*2);
    return 0;
}

两个方法都很类似，值得反复揣摩