生成一段封装好的Jacobi矩阵的C++代码

好的，以下是一段计算Jacobi矩阵的示例C++代码：

#include <iostream>
#include <vector>

using namespace std;

// 计算Jacobi矩阵
vector<vector<double>> jacobi_matrix(vector<double> (*f)(vector<double>), vector<double> x, double h=1e-6) {
    int n = x.size();
    vector<vector<double>> J(n, vector<double>(n));
    vector<double> fx = f(x);
    for (int i = 0; i < n; i++) {
        vector<double> xh = x;
        xh[i] += h;
        vector<double> fxh = f(xh);
        for (int j = 0; j < n; j++) {
            J[j][i] = (fxh[j] - fx[j]) / h;
        }
    }
    return J;
}

// 定义目标函数
vector<double> f(vector<double> x) {
    vector<double> fx(2);
    fx[0] = x[0] * x[0] - x[1];
    fx[1] = x[0] + x[1] * x[1] - 7;
    return fx;
}

int main() {
    vector<double> x = {1.0, 2.0};  // 初始点
    vector<vector<double>> J = jacobi_matrix(f, x);  // 计算Jacobi矩阵
    for (int i = 0; i < J.size(); i++) {
        for (int j = 0; j < J[i].size(); j++) {
            cout << J[i][j] << " ";
        }
        cout << endl;
    }
    return 0;
}

其中，`jacobi_matrix` 函数接受目标函数的指针、初始点和步长 `h` 等参数，返回计算得到的Jacobi矩阵。在 `main` 函数中，我们定义了目标函数 `f`，并指定了初始点 `x` 和步长 `h` 来计算Jacobi矩阵。你可以根据具体的问题，将目标函数传递给函数，并指定适当的初始点和步长来计算Jacobi矩阵。