复化Simpson公式求定积分的C语言代码

下面是使用复化Simpson公式求定积分的C语言代码：

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

double f(double x) {
    return sin(x); // 定义被积函数
}

double simpson(double a, double b) {
    int n = 1000; // 分割区间的个数
    double h = (b - a) / n;
    double sum = 0;
    for (int i = 0; i < n; i++) {
        double x1 = a + i * h;
        double x2 = a + (i + 1) * h;
        sum += f(x1) + 4 * f((x1 + x2) / 2) + f(x2);
    }
    return sum * h / 6;
}

int main() {
    double a = 0, b = M_PI; // 积分区间 [a, b]
    double result = simpson(a, b);
    printf("The result is: %lf\n", result);
    return 0;
}

其中，`f(x)` 函数用于定义被积函数，`simpson(a, b)` 函数用于计算积分值。在 `main()` 函数中，调用 `simpson(a, b)` 函数计算积分值，并输出结果。