00001 Matrix3f m = Matrix3f::Random(); 00002 Matrix3f y = Matrix3f::Random(); 00003 cout << "Here is the matrix m:" << endl << m << endl; 00004 cout << "Here is the matrix y:" << endl << y << endl; 00005 Matrix3f x; 00006 x = m.colPivHouseholderQr().solve(y); 00007 assert(y.isApprox(m*x)); 00008 cout << "Here is a solution x to the equation mx=y:" << endl << x << endl;