| #include <iostream> | 
 | #include <Eigen/Dense> | 
 |  | 
 | int main() | 
 | { | 
 |    Eigen::Matrix2f A; | 
 |    A << 1, 2, 2, 3; | 
 |    std::cout << "Here is the matrix A:\n" << A << std::endl; | 
 |    Eigen::SelfAdjointEigenSolver<Eigen::Matrix2f> eigensolver(A); | 
 |    if (eigensolver.info() != Eigen::Success) abort(); | 
 |    std::cout << "The eigenvalues of A are:\n" << eigensolver.eigenvalues() << std::endl; | 
 |    std::cout << "Here's a matrix whose columns are eigenvectors of A \n" | 
 |         << "corresponding to these eigenvalues:\n" | 
 |         << eigensolver.eigenvectors() << std::endl; | 
 | } |