1MatrixXf A = MatrixXf::Random(4,4); 2RealSchur<MatrixXf> schur(4); 3schur.compute(A, /* computeU = */ false); 4cout << "The matrix T in the decomposition of A is:" << endl << schur.matrixT() << endl; 5schur.compute(A.inverse(), /* computeU = */ false); 6cout << "The matrix T in the decomposition of A^(-1) is:" << endl << schur.matrixT() << endl; 7