Решение полиномиальных уравнений в Python

До сих пор у меня всегда была Mathematica для решения аналитических уравнений. Однако теперь мне нужно решить несколько сотен уравнений этого типа (характеристических полиномов)

a_20*x^20+a_19*x^19+...+a_1*x+a_0=0 (constant floats a_0,...a_20)

одновременно, что приводит к ужасно долгому времени вычислений в системе Mathematica.

Есть ли готовая к использованию команда в numpy или любом другом пакете для решения уравнения такого типа? (До сих пор я использовал Python только для моделирования, поэтому я мало разбираюсь в аналитических инструментах и ​​не смог найти ничего полезного в учебниках по numpy ).