用c++写复合Simpson算法，精度为0.0000001

以下是用 C 语言编写的复合 Simpson 算法，精度为 0.0000001 的代码：

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

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

double simpson(double a, double b) {
    double c = (a + b) / 2;
    return (f(a) + 4 * f(c) + f(b)) * (b - a) / 6;
}

double adaptive_simpson(double a, double b, double eps, double S) {
    double c = (a + b) / 2;
    double S1 = simpson(a, b);
    double S2 = simpson(a, c) + simpson(c, b);
    if (fabs(S2 - S1) <= 15 * eps) {
        return S2 + (S2 - S1) / 15;
    } else {
        return adaptive_simpson(a, c, eps / 2, S1) + adaptive_simpson(c, b, eps / 2, S2);
    }
}

int main() {
    double a = 0, b = 1, eps = 0.0000001;
    double S = simpson(a, b);
    double result = adaptive_simpson(a, b, eps, S);
    printf("The result is: %lf\n", result);
    return 0;
}

这个程序使用递归的方式实现了自适应 Simpson 积分算法，可以达到精度为 0.0000001。其中 `f` 函数是被积函数，可以根据需要进行修改。