请你帮我用c++写一个计算任意区间任意函数的第一类曲线积分的代码

好的，以下是一个计算任意区间任意函数的第一类曲线积分的C++代码示例：

#include <iostream>
#include <cmath>

using namespace std;

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

double integral(double a, double b, double eps) {
    // a 和 b 分别为积分区间的上下限，eps 为精度
    int n = 100; // 初始分割数
    double Tn = 0.5 * (func(a) + func(b)); // Tn 为梯形公式的初值
    double Sn = 0.0; // Sn 为 Simpson 公式的初值
    double x = a;
    double h = (b - a) / n;
    for (int i = 1; i < n; i++) {
        x += h;
        Tn += func(x);
        if (i % 2 == 0) {
            Sn += func(x);
        }
    }
    Tn *= h;
    Sn = (4 * Sn + 2 * Tn) * h / 3;
    double In = Sn;
    n *= 2;
    h /= 2;
    do {
        Tn = 0.5 * (Tn + In);
        Sn = 0.0;
        x = a + h;
        for (int i = 1; i < n; i++) {
            x += h;
            if (i % 2 == 1) {
                Sn += func(x);
            }
        }
        Sn *= h * 4 / 3;
        double temp = In;
        In = (4 * Sn - Tn) / 3;
        n *= 2;
        h /= 2;
    } while (fabs(In - temp) > eps);
    return In;
}

int main() {
    double a = 0.0, b = M_PI; // 积分区间
    double eps = 1e-6; // 精度
    double result = integral(a, b, eps);
    cout << "The result is: " << result << endl;
    return 0;
}

该代码使用了梯形公式和 Simpson 公式对被积函数进行数值积分，以达到计算曲线积分的目的。