[已解决]C语言求Π的近似值的问题

肖-肖

我的问题是为什么输出的是:5.135300
我感觉我写的没错啊，代码如下

1. /**编写程序,根据以下公式求Π的近似值,要求累加到某项小于1e-6次方时为止
   
2. Π/2=1 + 1/3 + 1*2/3*5 + 1*2*3*4/3*5*7*9 + ......+ 1*2*...*n/3*5*...*(2n+1)
   
3. */
   
4. #include <stdio.h>
   
5. int main()
   
6. {
   
7.     float i,j,n,m,pi;//i:分子 j:分母
   
8.     n=1.0;m=1.0;i=1.0;j=1.0;pi=0.0;
   
9.     while(n/m>=exp(-6))
   
10.     {
    
11.         pi+=n/m;//求和
    
12.         i++;
    
13.         j=j+2;
    
14.         n=n*i;
    
15.         m=m*j;
    
16.     }
    
17.     pi=pi*2;
    
18.     printf("%f",pi);
    
19.     return 0;
    
20. }
    

复制代码

请帮我分析一波,谢谢，尽可能详细点~~~

最佳答案

月排行榜 / 总排行榜

yuxijian2020

2021-4-17 16:54:47

题目的意思是 到某项小于 1e-6 为止，但是第一个小于1e-6的项还是要的

1. int main()
   
2. {
   
3.     double pi = 1.0;        //第一项 为1
   
4.     double cur = 0.0;       //当前项
   
5.     int i = 1;              //分子
   
6.     int j = 2 * i + 1;;     //分母
   
7.     double up = 1;          //计算后的分子
   
8.     double down = 1;        //计算后的分母
   
9. 
   
10.     do
    
11.     {
    
12.         up *= i;
    
13.         down *= j;
    
14.         cur = up / down;
    
15. 
    
16.         pi += cur;
    
17. 
    
18.         i++;
    
19.         j = 2 * i + 1;
    
20. 
    
21.     } while (cur >= 1e-6);
    
22.    
    
23. 
    
24.     printf("%lf\n", 2 * pi);
    
25. 
    
26.     return(0);
    
27. }

复制代码

结果

跳转到最佳答案楼层

wp231957

第二次循环貌似不对吧
期望1/3   你的貌似2/3

CodeCloud

第二次循环，i变成2，m变成3，加成了2/3而不是1/3

肖-肖

wp231957 发表于 2021-4-16 21:54
第二次循环貌似不对吧
期望1/3   你的貌似2/3

嗯 我把i的初值改为了0.0
结果是: 3.137130
      和: 3.1415926
差不多了,这么改行吧,应该没什么么问题了吧~

yuxijian2020

题目的意思是 到某项小于 1e-6 为止，但是第一个小于1e-6的项还是要的

1. int main()
   
2. {
   
3.     double pi = 1.0;        //第一项 为1
   
4.     double cur = 0.0;       //当前项
   
5.     int i = 1;              //分子
   
6.     int j = 2 * i + 1;;     //分母
   
7.     double up = 1;          //计算后的分子
   
8.     double down = 1;        //计算后的分母
   
9. 
   
10.     do
    
11.     {
    
12.         up *= i;
    
13.         down *= j;
    
14.         cur = up / down;
    
15. 
    
16.         pi += cur;
    
17. 
    
18.         i++;
    
19.         j = 2 * i + 1;
    
20. 
    
21.     } while (cur >= 1e-6);
    
22.    
    
23. 
    
24.     printf("%lf\n", 2 * pi);
    
25. 
    
26.     return(0);
    
27. }

复制代码

结果