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

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