1#include <iostream>
2#include <Eigen/Dense>
3
4using namespace std;
5using namespace Eigen;
6
7int 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