[已解决]两段程序的比较

fairylee83

题目：用π/4 ≈ 1-1/3+1/5-1/7+…公式求π的近似值，直到某一项的绝对值小于10^6为止。
小甲鱼：
#include <stdio.h>
#include<math.h>
void main()
{
    int s;    float n，t，pi;    t=1；pi=0;   n=1.0;    s=1;   
    while( fabs(t) > 1e-6 )   {
        pi = pi+t;   n = n+2;    s = -s;    t = s/n;
   }   
pi = pi*4;    printf(″pi=%10.6f＼n″，pi);}
运行结果：pi=3.141594
我：
#include <stdio.h>
#include <math.h>
void main()
{
        double s,k,pi;
        int i,j;
        k=-1.0;       
        i=1;
        j=1;
        s=1.0;
        while(fabs(k/j)>10e-6)
        {       
                j=2*i+1;
                s+=k/j;
                k*=-1;
                i++;
        }
        pi=4.0*s;
        printf("%lf\n",pi);
}
运行结果：3.141613
二者的差别好大，哪位大神能帮忙解释下产生这么大误差的原因出在哪里？为什么小甲鱼的程序最优最接近正确答案？

最佳答案

月排行榜 / 总排行榜

U201010009

2018-1-9 17:12:24

①fabs(k/j)>10e-6 判断应该是 1e-6 不是 10e-6
②小甲鱼老师和你的操作有些不一样。小甲鱼老师是先计算 +(k/j)之后再判断  k/(j+2) > 1e-6,而你的是 +(k/j)之后再判断k/j > 1e-6 算法都不一样的

跳转到最佳答案楼层

U201010009

①fabs(k/j)>10e-6 判断应该是 1e-6 不是 10e-6
②小甲鱼老师和你的操作有些不一样。小甲鱼老师是先计算 +(k/j)之后再判断  k/(j+2) > 1e-6,而你的是 +(k/j)之后再判断k/j > 1e-6 算法都不一样的

爱晴海

float 与double型数据运算时它们精度相差的问题，精度高系统储存是要求越高，达到时就会舍去。具体我也不太清楚。望有用

fairylee83

嗯。谢谢您。把10e-6改成了1e-6就得到了3.141595。看来是我实数型数据的表达方法没有学明白。又返工去复习了一下。
另外，好好奇。虽然您的答案最佳，但怎么不像是我评的呢？我失忆了？