title: Solve polynomial equation with negative power or non-integer power with python
body: <p>I tried to solve the equation</p>
  
  <pre><code>1600 = 0.41 + 6.31*d**-1.54 + 2.42*d**-3 
  </code></pre>
  
  <p>Mathematica can give me <code>d=6.4673</code> in less than 1 second. But I could not get the answer with python symbolic calculation. </p>
  
  <p>Using "solve" of sympy take forever. Is there any way to solve this equation by using python symbolic calculation? It seems that the issues comes mainly from the non-integer negative power. </p>