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