1 #include <unsupported/Eigen/Polynomials>
2 #include <iostream>
3 
4 using namespace Eigen;
5 using namespace std;
6 
main()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