1 #include <unsupported/Eigen/Polynomials> 2 #include <iostream> 3 4 using namespace Eigen; 5 using namespace std; 6 7 int main() 8 { 9 Vector4d roots = Vector4d::Random(); 10 cout << "Roots: " << roots.transpose() << endl; 11 Eigen::Matrix<double,5,1> polynomial; 12 roots_to_monicPolynomial( roots, polynomial ); 13 cout << "Polynomial: "; 14 for( int i=0; i<4; ++i ){ cout << polynomial[i] << ".x^" << i << "+ "; } 15 cout << polynomial[4] << ".x^4" << endl; 16 Vector4d evaluation; 17 for( int i=0; i<4; ++i ){ 18 evaluation[i] = poly_eval( polynomial, roots[i] ); } 19 cout << "Evaluation of the polynomial at the roots: " << evaluation.transpose(); 20 } 21
