根据以下公式求PI值。例如，给指定精度的变量eps输入0.0005时，应当输出Pi＝3.140578
          PI/2＝1+1/3+(1/3)*(2/5)+(1/3)*(2/5)*(3/7)+(1/3)*(2/5)*(3/7)*(4/9)+...

#include<stdio.h>
main()
{
    double r,eps,temp;
    int m=1;
    printf("\nPlease enter a precision: ");
    scanf("%lf",&eps);
    r=0.0;
    /************found*********/
    temp=0;
    /************found*********/
    while(temp<=eps)
    {
        r+=temp;
        temp=(temp*m)/(2*m+1);
        m++;
    }
    /************found*********/
    printf("\neps=%lf,Pi=%lf\n\n",eps,r);
}

这个程序的答案我知道，但是不明白为什么这么改，请教各位大侠，谢谢！

[此贴子已经被作者于2016-4-17 11:11编辑过]

#include <stdio.h>

int main()
{
    double r,eps,temp;
    int m=1;

    printf("\nPlease enter a precision: ");
    scanf("%lf",&eps);
    r=0.0;

    /************found*********/
    temp=1;  //=1,下面的while是从第二项1/3开始的。
    /************found*********/
    while(temp>=eps)  //>=,小于精度才停止迭加
    {
        r+=temp;
        temp=(temp*m)/(2*m+1);
        m++;
    }
    /************found*********/
    printf("\neps=%.4f,Pi=%.6f\n",eps,r*2);  //r*2,这个你看公式，PI/2=....上面的结果只是我们想要的一半

    return 0;
}

[此贴子已经被作者于2016-4-17 12:27编辑过]

temp=(temp*m)/(2*m+1);
这个语句执行后temp始终等于0

PI/2＝1+ 1*1/3+(1/3)*(2/5)+(1/3)*(2/5)*(3/7)+(1/3)*(2/5)*(3/7)*(4/9)+...
PI/2＝1+1/3+(1/3)*(2/5)+(1/3)*(2/5)*(3/7)+(1/3)*(2/5)*(3/7)*(4/9)+...
PI/2＝1+1/3+(1/3)*(2/5)+(1/3)*(2/5)*(3/7)+(1/3)*(2/5)*(3/7)*(4/9)+...
PI/2＝1+1/3+(1/3)*(2/5)+(1/3)*(2/5)*(3/7)+(1/3)*(2/5)*(3/7)*(4/9)+...

不是应该 temp = temp * (m / (2 * m +1));

[此贴子已经被作者于2016-4-17 12:46编辑过]

不错，我那样写是一样的，数学上的转换。

恩恩，我范傻了