1 MatrixXd A
= MatrixXd::Random(6,6);
2 cout
<< "Here is a random 6x6 matrix, A:" << endl
<< A
<< endl
<< endl
;
4 EigenSolver
<MatrixXd
> es(A
);
5 MatrixXd D
= es
.pseudoEigenvalueMatrix();
6 MatrixXd V
= es
.pseudoEigenvectors();
7 cout
<< "The pseudo-eigenvalue matrix D is:" << endl
<< D
<< endl
;
8 cout
<< "The pseudo-eigenvector matrix V is:" << endl
<< V
<< endl
;
9 cout
<< "Finally, V * D * V^(-1) = " << endl
<< V
* D
* V
.inverse() << endl
;