eigen/doc/snippets/MatrixBase_identity_int_int.cpp

cout << MatrixXd::identity(4, 3) << endl;