ex_compan.m

p =  [-1 -1 2 1];   %polynme -x^3-x^2 + 2x + 1
a =  compan(p)      %matrice compagnon
poly(a)             %polynme caractristique

r0 = [1 2 3 4];     %vecteur de racines
p =  poly(r0)       %polynme ayant ces racines
r1 =  roots(p)      %racines du polynme