00001 MatrixXf m = MatrixXf::Random(3,5); 00002 cout << "Here is the matrix m:" << endl << m << endl; 00003 MatrixXf ker = m.fullPivLu().kernel(); 00004 cout << "Here is a matrix whose columns form a basis of the kernel of m:" 00005 << endl << ker << endl; 00006 cout << "By definition of the kernel, m*ker is zero:" 00007 << endl << m*ker << endl;