Yes, NumPy does this too. See https://github.com/numpy/numpy/blob/v1.8.1/numpy/lib/polynomial.py#L126 The first problem still works because computing the polynomial's coefficients and then computing the roots from that isn't stable either.