1MatrixXd X = MatrixXd::Random(5,5); 2MatrixXd A = X * X.transpose(); 3X = MatrixXd::Random(5,5); 4MatrixXd B = X * X.transpose(); 5 6GeneralizedSelfAdjointEigenSolver<MatrixXd> es(A,B,EigenvaluesOnly); 7cout << "The eigenvalues of the pencil (A,B) are:" << endl << es.eigenvalues() << endl; 8es.compute(B,A,false); 9cout << "The eigenvalues of the pencil (B,A) are:" << endl << es.eigenvalues() << endl; 10