Returns polynomial roots as a Vector.
Unlike solve(), which returns a SolutionSet with solver metadata, this function preserves
the roots API's Vector return type. Numeric constants retain the historical square-root
interpretation. Univariate polynomials delegate to PolynomialSolver and therefore return
numerical roots.
