1 MatrixXd
X = MatrixXd::Random(5,5);
2 MatrixXd
A =
X +
X.transpose();
3 cout <<
"Here is a random symmetric 5x5 matrix:" << endl <<
A << endl << endl;
4 Tridiagonalization<MatrixXd>
triOfA(
A);
6 cout <<
"The orthogonal matrix Q is:" << endl << Q << endl;
8 cout <<
"The tridiagonal matrix T is:" << endl << T << endl << endl;
9 cout <<
"Q * T * Q^T = " << endl << Q * T * Q.transpose() << endl;
cout<< "Here is a random symmetric 5x5 matrix:"<< endl<< A<< endl<< endl;Tridiagonalization< MatrixXd > triOfA(A)