00001 MatrixXd ones = MatrixXd::Ones(3,3); 00002 SelfAdjointEigenSolver<MatrixXd> es(ones); 00003 cout << "The first eigenvector of the 3x3 matrix of ones is:" 00004 << endl << es.eigenvectors().col(1) << endl;