1 #include <unsupported/Eigen/Polynomials>
9 Vector4d roots = Vector4d::Random();
10 cout <<
"Roots: " << roots.transpose() << endl;
13 cout <<
"Polynomial: ";
14 for(
int i=0;
i<4; ++
i ){ cout << polynomial[
i] <<
".x^" <<
i <<
"+ "; }
15 cout << polynomial[4] <<
".x^4" << endl;
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();