请你帮我用c++写一段可以计算第一类曲面积分的代码，我可以通过菜单输入任意区间和函数

好的，以下是用C++计算第一类曲面积分的代码：

#include <iostream>
#include <cmath>
using namespace std;

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

double Simpson(double a, double b) { // Simpson公式计算积分
    double c = (a + b) / 2;
    return (f(a) + 4 * f(c) + f(b)) * (b - a) / 6;
}

double adaptiveSimpson(double a, double b, double eps) { // 自适应Simpson公式
    double c = (a + b) / 2;
    double L = Simpson(a, c), R = Simpson(c, b), S = Simpson(a, b);
    if (fabs(L + R - S) <= 15 * eps) {
        return L + R + (L + R - S) / 15;
    } else {
        return adaptiveSimpson(a, c, eps / 2) + adaptiveSimpson(c, b, eps / 2);
    }
}

int main() {
    double a, b, eps;
    cout << "请输入积分区间[a, b]和精度eps：" << endl;
    cin >> a >> b >> eps;
    double res = adaptiveSimpson(a, b, eps); // 计算积分
    cout << "积分结果为：" << res << endl;
    return 0;
}

以上代码使用了自适应Simpson公式计算积分，可以较好地适应被积函数的变化。在主函数中通过菜单输入积分区间和精度，然后调用`adaptiveSimpson`函数计算积分并输出结果。你可以根据需要修改被积函数来计算不同的积分。