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