精确度问题

h-b-z-d-j-s-m

题目是根据公式求Π，要求精确度为0.0005，即某项小于0.0005时停止
为什么我的运行出来结果跟答案给的结果不一样，哪里错了？
fun函数里面为自己写的
答案写的是

double s=0.0,s1=1.0;
int n=1;
while(s1>=eps)
{
s=s+s1;
s1=s1*n/(2*n+1);
n++;
}
return 2*s;

fun函数里面是我自己写的，麻烦帮我看下

#include <stdio.h>
#include <math.h>
double  fun ( double  eps)
{
        double pi,a=1.0,b=1.0,i=1.0,t;
        pi=1.0;
        t=1.0;
        while(t>=eps)
        {
                a*=i;
                b*=(2*i+1);
                t=a/b;
                pi+=t;
                i++;
        }
        return (2*pi);


}

void main( )
{ double  x;void NONO ();
  printf("Input eps:") ;
  scanf("%lf",&x); printf("\neps = %lf, PI=%lf\n", x, fun(x));
  NONO();
}

void NONO ()
{/* 本函数用于打开文件，输入数据，调用函数，输出数据，关闭文件。 */
  FILE *fp, *wf ;
  int i ;
  double x ;

  fp = fopen("in.dat","r") ;
  wf = fopen("out.dat","w") ;
  for(i = 0 ; i < 10 ; i++) {
    fscanf(fp, "%lf", &x) ;
    fprintf(wf, "%lf\n", fun(x)) ;
  }
  fclose(fp) ;
  fclose(wf) ;
}

人造人

?
https://fishc.com.cn/forum.php?m ... 582&pid=5842010

h-b-z-d-j-s-m

人造人 发表于 2022-5-20 08:45
?
https://fishc.com.cn/forum.php?mod=redirect&goto=findpost&ptid=213582&pid=5842010

那个没有解决我的问题，我想看看我的方法哪里错了

人造人

h-b-z-d-j-s-m 发表于 2022-5-20 10:42
那个没有解决我的问题，我想看看我的方法哪里错了

???
https://fishc.com.cn/forum.php?m ... 582&pid=5842010

#include <stdio.h>
#include <math.h>

double fun(double eps) {
    double pi, a = 1.0, b = 1.0, i = 1.0, t;
    pi = 1.0;
    t = 1.0;
    //while(t >= eps) {
    while(1) {
        a *= i;
        b *= (2 * i + 1);
        t = a / b;
        if(t < eps) break;
        pi += t;
        i++;
    }
    return (2 * pi);
}

// void main( )
int main(void) {
    double x;
    void NONO();
    //printf("Input eps:");
    //scanf("%lf", &x);
    x = 0.0005;
    printf("\neps = %lf, PI=%lf\n", x, fun(x));
    //NONO();
}

void NONO() { /* 本函数用于打开文件，输入数据，调用函数，输出数据，关闭文件。 */
    FILE *fp, *wf;
    int i;
    double x;

    fp = fopen("in.dat", "r");
    wf = fopen("out.dat", "w");
    for(i = 0; i < 10; i++) {
        fscanf(fp, "%lf", &x);
        fprintf(wf, "%lf\n", fun(x));
    }
    fclose(fp);
    fclose(wf);
}