1MatrixXd ones = MatrixXd::Ones(3,3);
2EigenSolver<MatrixXd> es(ones);
3cout << "The first eigenvector of the 3x3 matrix of ones is:"
4     << endl << es.eigenvectors().col(1) << endl;
5