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