#include <unsupported/Eigen/Polynomials> #include <cmath> #include <iostream> int main() { Eigen::PolynomialSolver<double, Eigen::Dynamic> solver; Eigen::VectorXd coeff(3); coeff[0] = 1.0; coeff[1] = 1.0; coeff[2] = 1.0; solver.compute(coeff); const Eigen::PolynomialSolver<double, Eigen::Dynamic>::RootsType & r = solver.roots(); for(int i =0;i<r.rows();++i) { std::cout << r[i] << std::endl; } return 0; }