1MatrixXf A = MatrixXf::Random(3, 2);
2VectorXf b = VectorXf::Random(3);
3cout << "The solution using the QR decomposition is:\n"
4     << A.colPivHouseholderQr().solve(b) << endl;
5