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