Found this unfortunate behavior in polyroots:


v := [-1 1 0 0 0]


polyroots(v) =
index out of bounds


the 4th order polynomial degenerates to a first order with single real root x=1. Polyroots couldn't figure that out. Should it not analyze the coefficient vector to determine correct polynomial order and then proceed to find the roots?