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