小编roc*_*man的帖子

我在C++中使用Eigen库,我试图找到矩阵的行列式.根据我如何初始化矩阵,我会得到不同的结果.

方法一:

MatrixXd a(3, 3);
for (int n = 0; n < 3; n++)
    for (int m = 0; m < 3; m++)
        a(n,m) = (double) (n + m*m + 2.5)/3;

cout << "Matrix a: " << endl;
cout << a << endl;
cout << "Determinat of matrix a is: " << a.determinant() << endl;

这部分代码打印出来

Matrix a:
0.8333333  1.166667  2.166667
1.166667       1.5       2.5
1.5  1.833333  2.833333
Determinat of matrix a is: -7.401487e-17

方法二:

MatrixXd b(3, 3);
b << 0.8333333, …