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