doc/snippets/MatrixBase_identity.cpp

cout << Matrix<double, 3, 4>::identity() << endl;