用c++写Simpson算法

以下是用 C 语言实现 Simpson 算法的代码：

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

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

double simpson(double a, double b, int n) {
    double h = (b - a) / n;
    double x, sum = 0;
    int i;

    for (i = 1; i < n; i += 2) {
        x = a + i * h;
        sum += 4 * f(x);
    }

    for (i = 2; i < n - 1; i += 2) {
        x = a + i * h;
        sum += 2 * f(x);
    }

    return (h / 3) * (f(a) + f(b) + sum);
}

int main() {
    double a = 0, b = M_PI; // 积分区间 [a, b]
    int n = 100; // 分割数

    printf("Simpson's rule: %lf\n", simpson(a, b, n));

    return 0;
}

这个程序实现了 Simpson 算法，用于计算被积函数在给定区间上的定积分。其中，`f(x)` 函数是被积函数，`simpson(a, b, n)` 函数是 Simpson 算法的实现，`main()` 函数则是程序的入口。