1c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan KamathMatrixXf m = MatrixXf::Random(3,5);
2c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathcout << "Here is the matrix m:" << endl << m << endl;
3c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan KamathMatrixXf ker = m.fullPivLu().kernel();
4c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathcout << "Here is a matrix whose columns form a basis of the kernel of m:"
5c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath     << endl << ker << endl;
6c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathcout << "By definition of the kernel, m*ker is zero:"
7c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath     << endl << m*ker << endl;
8