1*bf2c3715SXin Li #include <iostream>
2*bf2c3715SXin Li #include <Eigen/Dense>
3*bf2c3715SXin Li
4*bf2c3715SXin Li using namespace std;
5*bf2c3715SXin Li using namespace Eigen;
6*bf2c3715SXin Li
main()7*bf2c3715SXin Li int main()
8*bf2c3715SXin Li {
9*bf2c3715SXin Li Matrix3f A;
10*bf2c3715SXin Li Vector3f b;
11*bf2c3715SXin Li A << 1,2,3, 4,5,6, 7,8,10;
12*bf2c3715SXin Li b << 3, 3, 4;
13*bf2c3715SXin Li cout << "Here is the matrix A:\n" << A << endl;
14*bf2c3715SXin Li cout << "Here is the vector b:\n" << b << endl;
15*bf2c3715SXin Li Vector3f x = A.colPivHouseholderQr().solve(b);
16*bf2c3715SXin Li cout << "The solution is:\n" << x << endl;
17*bf2c3715SXin Li }
18