21 lines
635 B
C++
21 lines
635 B
C++
|
#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();
|
||
|
|
}
|