sympy - symbolically solve equations with float power
<p>Using sympy, I define the symbols,</p>

<pre><code>a, b, c = sympy.symbols(['a', 'b', 'c'])
</code></pre>

<p>Then, when I try to solve the following system of equations,</p>

<pre><code>sympy.solve([sympy.Eq(b - a**2.552 - c), sympy.Eq(a, 2)])
</code></pre>

<p>I get the solution,</p>

<pre><code>[{b: c + 5.86446702875684, a: 2.00000000000000}]
</code></pre>

<p>But, when I try to solve,</p>

<pre><code>sympy.solve([sympy.Eq(b - a**2.552 - c), sympy.Eq(b, 2)])
</code></pre>

<p>It just seems to keep running (for ~4hrs), with no solution. Any help would be appreciated!</p>

<p>I don't know why but <code>rational=False</code> helps</p>

<pre><code>sympy.solve([sympy.Eq(b - a**2.552 - c), sympy.Eq(b, 2)], rational=False)
</code></pre>

<p>see: <a href="http://stackoverflow.com/questions/17087629/sympy-hangs-when-trying-to-solve-a-simple-algebraic-equation">sympy hangs when trying to solve a simple algebraic equation</a></p>

python sympy