compan(p)
length(p)-1
rows and columnsp
with -p[2:length(p)]/p[1]
as first row.The eigenvalues of this matrix are the roots of the polynomial.
roots
p <- c(1, 0, -7, 6)
compan(p)
# 0 7 -6
# 1 0 0
# 0 1 0
Run the code above in your browser using DataLab