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