用lapacke_dgesv求解cpp中的线性方程组

#include <iostream>
#include <lapacke.h>

int main() {
    // 定义线性方程组的系数矩阵A和右侧向量b
    double A[3][3] = {{1, 2, 3}, {4, 5, 6}, {7, 8, 10}};
    double b[3] = {3, 6, 9};

    // 定义必要的变量
    int n = 3; // 线性方程组的维度
    int lda = n; // 矩阵A的列数
    int ipiv[n]; // 主元素索引
    int info; // 返回值

    // 调用lapacke_dgesv函数求解线性方程组
    info = LAPACKE_dgesv(LAPACK_ROW_MAJOR, n, 1, *A, lda, ipiv, b, 1);

    // 检查求解结果
    if (info == 0) {
        std::cout << "线性方程组求解成功！" << std::endl;
        std::cout << "解向量x：" << b[0] << ", " << b[1] << ", " << b[2] << std::endl;
    } else {
        std::cout << "线性方程组求解失败！" << std::endl;
    }

    return 0;
}