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