1 MatrixXf matA(2,2);
2 matA << 2, 0,  0, 2;
3 matA = matA * matA;
4 cout << matA;
5