doc/snippets/FullPivLU_kernel.cpp - eigen - Git at Google

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