请帮我用C++写一个计算二重积分的程序，适用于任意函数，我可以使用它输入某一个函数，在某一区间上得到他的二重积分。

好的，以下是一个简单的C++程序，可以计算给定函数在矩形区域上的二重积分。代码中使用了拉格朗日4点公式进行数值积分。

#include <iostream>
#include <cmath>

using namespace std;

double f(double x, double y) {
    // 这里输入你想要计算的函数
    return x*x + y*y;
}

double integrate(double a, double b, double c, double d, int n) {
    double h1 = (b-a)/(2*n), h2 = (d-c)/(2*n);
    double sum = 0;
    for (int i = 0; i < 2*n; i++) {
        double x1 = a + i*h1, x2 = a + (i+1)*h1;
        for (int j = 0; j < 2*n; j++) {
            double y1 = c + j*h2, y2 = c + (j+1)*h2;
            double fx1y1 = f(x1, y1), fx2y1 = f(x2, y1), fx1y2 = f(x1, y2), fx2y2 = f(x2, y2);
            double Lx1 = (fx1y1 + fx2y1)*h1/2 + (fx1y2 + fx2y2)*h1/2;
            double Lx2 = (fx2y1 + f(x2+h1, y1))*h1/2 + (fx2y2 + f(x2+h1, y2))*h1/2;
            sum += (Lx1 + Lx2)*h2/2;
        }
    }
    return sum;
}

int main() {
    double a, b, c, d;
    int n;
    cout << "Enter the limits of integration (a, b, c, d): ";
    cin >> a >> b >> c >> d;
    cout << "Enter the number of subintervals (must be even): ";
    cin >> n;
    double result = integrate(a, b, c, d, n);
    cout << "The result is: " << result << endl;
    return 0;
}

代码中的 `f` 函数用于输入要计算的函数，可以根据自己的需要修改。`integrate` 函数用于计算二重积分，输入矩形区域的边界以及将矩形分成的小矩形数量。代码中使用了拉格朗日4点公式进行数值积分。

请注意，这个程序只是一个简单的示例，对于某些函数，可能需要更高级的数值积分方法才能得到准确的结果。