title: Solving polynomial equations in Python
body: <p>Up to now I have always Mathematica for solving analytical equations. Now however I need to solve a few hundred equations of this type (characteristic polynomials)</p>
  
  <pre><code>a_20*x^20+a_19*x^19+...+a_1*x+a_0=0 (constant floats a_0,...a_20)
  </code></pre>
  
  <p>at once which yields awfully long calculation times in Mathematica.</p>
  
  <p>Is there like a ready to use command in numpy or any other package to solve an equation of this type? (up to now I have used Python only for simulations so I don't know much about analytical tools and I couldn't find anything useful in the numpy tutorials).</p>