eigen-2.1.1: eigen3/doc/snippets/Matrix_setIdentity_int_int.cpp

MatrixXf m;
m.setIdentity(3, 3);
cout << m << endl;