1 Matrix3f m
= Matrix3f::Random();
2 Matrix3f y
= Matrix3f::Random();
3 cout
<< "Here is the matrix m:" << endl
<< m
<< endl
;
4 cout
<< "Here is the matrix y:" << endl
<< y
<< endl
;
6 x
= m
.colPivHouseholderQr().solve(y
);
7 assert(y
.isApprox(m
*x
));
8 cout
<< "Here is a solution x to the equation mx=y:" << endl
<< x
<< endl
;