22 #include "sinfo_solve_poly_root.h"
40 sinfo_set_companion_matrix(
const double *a,
size_t nc,
double *m)
44 for (i = 0; i < nc; i++)
45 for (j = 0; j < nc; j++)
46 MAT (m, i, j, nc)= 0.0;
48 for (i = 1; i < nc; i++)
49 MAT (m, i, i - 1, nc)= 1.0;
51 for (i = 0; i < nc; i++)
52 MAT (m, i, nc - 1, nc)= -a[i] / a[nc];