1 MatrixXd ones = MatrixXd::Ones(3,3); 2 SelfAdjointEigenSolver<MatrixXd> es(ones); 3 cout << "The eigenvalues of the 3x3 matrix of ones are:" 4 << endl << es.eigenvalues() << endl; 5