辛普森规则错误

该代码使用辛普森法则来计算X * SIN（X）与（1,2）边界的积分。我遇到的问题是，它的实际价值非常接近。即使有999次迭代，它仍然没有达到要求。虽然我有一个单独的程序，它使用梯形法则来处理同样的事情，并且它在1000次迭代后恰好达到了这一点。它应该打的一点是“1.440422”辛普森规则错误

是什么应该发生的事情对辛普森的规则？或者我的代码有问题吗？

#include <stdlib.h> 
#include <stdio.h> 
#include <math.h> 
double f(double x); 
int main() 
{ 
    double x,result,y,z,h,s1,s2; 
    s1 = 0; 
    s2 = 0; 
    int i,n; 
    printf("\nHow many points to you want it evaluated at (odd number)? And what are the bounds? lower bound,upper bound >\n"); 
    scanf("%d %lf,%lf",&n,&y,&z); 
    h = (z-y)/n; 
    result = 0; 
    if(n%2!=0) 
    { 
     for(i=0;i<n;i++) 
    { 
     if(i%2==0) 
     { 
      s1 = s1+f(y+i*h); 
     } 
     else 
     { 
      s2 = s2+f(y+i*h); 
     } 
    } 
     result = (h/3)*(f(y)+f(z)+4*s2+2*s1); 
     printf("\nThe value is %lf with %d interations\n",result,i); 
    } 
    else 
    { 
     printf("\n The number of points has to be odd, try again\n"); 
    } 
} 


double f(double x) 
{ 
    return(x*sin(x)); 
}

来源

2015-10-14 user3908631

欢迎近似的世界和浮点运算。如果您使用近似公式，则会得到近似结果。这可能是一种侥幸，这个梯形法则最终会以'正确'的正确答案结束。 –

是的，那令我困惑的部分是我认为辛普森的近似值比梯形误差更小。但我想这可能只是我正在使用的功能而已。谢谢 – user3908631

在99999次迭代中，我得到了你需要的值'1.440422'。一些函数比其他函数更快地收敛，例如有几种方法可以计算'pi'，但并非全部都有用。 –

您看到的问题可能是因为用于读取数字的格式字符串。

scanf("%d %lf,%lf",&n,&y,&z); 
     // ^^^ Is the , there on purpose?

尝试移除,从格式字符串，看看问题是否会消失。

它不能强调 - 总是检查返回值scanf。

if (scanf("%d %lf %lf", &n, &y, &z) != 3) 
{ 
    // Deal with error. 
}

为了确保数字阅读是准确的，添加回声输入回stdout线。

printf("n: %d, y: %lf, z: %lf\n", n, y, z);

我注意到几个错误，在你的代码：

间隔h是不正确的。由于您使用的是n分，因此有n-1间隔。因此，h必须是：
```
h = (z-y)/(n-1); 
```
既然你要添加在最后陈述f(y)和f(z)，循环必须是：
```
// Not good to use for(i=0;i<n;i++) 
for(i=1;i<n-1;i++) 
{ 
```

随着这些补丁，我得到输出1.440422使用n = 1001。

来源

2015-10-14 15:47:50

scanf中的逗号表示当您输入用逗号分隔它们的值时。但我仍然尝试过，并没有奏效。并用printf检查出它的值是正确的。 – user3908631

谢谢你，这些修复是一个更好的近似。 – user3908631

@ user3908631，我很高兴能够提供帮助。 –

辛普森规则错误

回答

相关问题