| #include <unsupported/Eigen/Polynomials> |
| #include <iostream> |
| |
| using namespace Eigen; |
| using namespace std; |
| |
| int main() |
| { |
| Vector4d roots = Vector4d::Random(); |
| cout << "Roots: " << roots.transpose() << endl; |
| Eigen::Matrix<double,5,1> polynomial; |
| roots_to_monicPolynomial( roots, polynomial ); |
| cout << "Polynomial: "; |
| for( int i=0; i<4; ++i ){ cout << polynomial[i] << ".x^" << i << "+ "; } |
| cout << polynomial[4] << ".x^4" << endl; |
| Vector4d evaluation; |
| for( int i=0; i<4; ++i ){ |
| evaluation[i] = poly_eval( polynomial, roots[i] ); } |
| cout << "Evaluation of the polynomial at the roots: " << evaluation.transpose(); |
| } |