#include <iostream>
#include <synaps/linalg.h>
#include <synaps/upol.h>       // univariate polynomials

typedef UPolDse<double>    pol_t;
typedef MatrDse<double,lapack::rep2d<double> > matr_t;


matr_t companion(const pol_t & p)
{
  ...
}
int main(int argc, char** argv) 
{
  using std::cout;
  using std::endl;
  pol_t p("2*x^3-34*x+23");                     cout<<p<<endl;

  matr_t C=companion(p);
  cout <<  C <<endl;
  cout <<Eigenval(C)<<endl;
  
}