龙贝格算法C++代码，要求被积函数是可以输入的

以下是一个简单的龙贝格算法的 C 代码，其中被积函数可以在运行程序时输入：

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

// 定义被积函数
double func(double x) {
    // 这里可以根据需要修改被积函数
    return sin(x);
}

// 定义龙贝格算法函数
double romberg(double a, double b, double eps) {
    int i, j, k, n = 1, m = 20;
    double h, s, q, ep, p[m + 1][m + 1];
    h = b - a;
    p[0][0] = h * (func(a) + func(b)) / 2.0;
    do {
        s = 0.0;
        for (i = 0; i < n; i++) {
            q = a + (i + 0.5) * h;
            s += func(q);
        }
        p[1][0] = 0.5 * (p[0][0] + h * s / n);
        for (k = 1; k <= m; k++) {
            n = pow(2, k);
            h = (b - a) / n;
            s = 0.0;
            for (i = 0; i < n; i++) {
                q = a + (i + 0.5) * h;
                s += func(q);
            }
            p[k + 1][0] = 0.5 * (p[k][0] + h * s / n);
            for (j = 1; j <= k; j++) {
                ep = (p[k + 1][j - 1] - p[k][j - 1]) / (pow(4, j) - 1);
                p[k + 1][j] = p[k + 1][j - 1] + ep;
            }
            if (fabs(ep) < eps) {
                return p[k + 1][k + 1];
            }
        }
    } while (k <= m);
    return -1.0;
}

int main() {
    double a, b, eps, res;
    printf("请输入积分区间 [a, b] 和精度 eps：\n");
    scanf("%lf %lf %lf", &a, &b, &eps);
    res = romberg(a, b, eps);
    printf("积分结果为 %.6f\n", res);
    return 0;
}

在运行程序时，会提示用户输入积分区间和精度，用户可以根据需要输入。程序会输出计算得到的积分结果。