精确度问题

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

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

1. double s=0.0,s1=1.0;
   
2. int n=1;
   
3. while(s1>=eps)
   
4. {
   
5. s=s+s1;
   
6. s1=s1*n/(2*n+1);
   
7. n++;
   
8. }
   
9. return 2*s;

复制代码

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

1. #include <stdio.h>
   
2. #include <math.h>
   
3. double  fun ( double  eps)
   
4. {
   
5.         double pi,a=1.0,b=1.0,i=1.0,t;
   
6.         pi=1.0;
   
7.         t=1.0;
   
8.         while(t>=eps)
   
9.         {
   
10.                 a*=i;
    
11.                 b*=(2*i+1);
    
12.                 t=a/b;
    
13.                 pi+=t;
    
14.                 i++;
    
15.         }
    
16.         return (2*pi);
    
17. 
    
18. 
    
19. }
    
20. 
    
21. void main( )
    
22. { double  x;void NONO ();
    
23.   printf("Input eps:") ;
    
24.   scanf("%lf",&x); printf("\neps = %lf, PI=%lf\n", x, fun(x));
    
25.   NONO();
    
26. }
    
27. 
    
28. void NONO ()
    
29. {/* 本函数用于打开文件，输入数据，调用函数，输出数据，关闭文件。 */
    
30.   FILE *fp, *wf ;
    
31.   int i ;
    
32.   double x ;
    
33. 
    
34.   fp = fopen("in.dat","r") ;
    
35.   wf = fopen("out.dat","w") ;
    
36.   for(i = 0 ; i < 10 ; i++) {
    
37.     fscanf(fp, "%lf", &x) ;
    
38.     fprintf(wf, "%lf\n", fun(x)) ;
    
39.   }
    
40.   fclose(fp) ;
    
41.   fclose(wf) ;
    
42. }

复制代码

人造人

?
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);
}