1Tridiagonalization<MatrixXf> tri; 2MatrixXf X = MatrixXf::Random(4,4); 3MatrixXf A = X + X.transpose(); 4tri.compute(A); 5cout << "The matrix T in the tridiagonal decomposition of A is: " << endl; 6cout << tri.matrixT() << endl; 7tri.compute(2*A); // re-use tri to compute eigenvalues of 2A 8cout << "The matrix T in the tridiagonal decomposition of 2A is: " << endl; 9cout << tri.matrixT() << endl; 10