C 语言的复合辛普森规则

Question

已关闭。这个问题是 not reproducible or was caused by typos 。目前不接受答案。

---

这个问题是由拼写错误或无法再重现的问题引起的。虽然类似的问题可能是 on-topic在这里，这个问题的解决方式不太可能帮助 future 的读者。

已关闭 9 年前。

Improve this question

我正在尝试用 C 语言编写复合辛普森法则。公式为:

哪里x_j=a+jh对于 j=0, 1, ..., n-1 , n与 h=(b-a)/n ;特别是x_0=a和x_n=b .

由于某种原因，第一个和第二个循环具有相同的值。我检查了很多次，但似乎找不到我的错误。

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

float f(float);
float a;
float b;
float x;
float h;
int n;
int j;


 a=0;
 b=2;
 n=8;
 h = (n - j) / b; 

float first;
float second;
int main() {


  sum = (h / 3.0f) * (f(h) + f(n));

  printf("%f\n", sum);

  second =  (4.0f) * h * f(a);
  }
  printf("second sum: %f\n",second );

  sum = sum + first + second;
  printf("%f\n", sum); 
  return 0;
}

答案应该在 3.1 左右(最终的值 sum )

Answer 1

您的部门可能不会达到您的预期:

(2 / 3) == 0

将 int 除以 int 将得到 int。

使用浮点常量(2.0f/3.0f)

编辑:

您仍然遇到与其他 n/2 相同的问题。

打印 float 时应该使用%f:printf("first sum: %f\n",first);

The value of the integral is: 3.109337