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

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