SymPy and square roots of complex numbers
<p>When using <code>solve</code> to compute the roots of a quadratic equation, SymPy returns expressions which could be simplified but I can't get it to simplify them.  A minimal example looks like so:</p>

<pre><code>from sympy import *
sqrt(-24-70*I)
</code></pre>

<p>Here, SymPy just returns <code>sqrt(-24-70*I)</code> while Mathematica or Maple will answer with the equivalent of <code>5-7*I</code>.</p>

<p>I'm aware that there are two square roots, but this behavior entails that SymPy will, for example, return pretty complicated solutions from</p>

<pre><code>z = symbols("z")
solve(z ** 2 + (1 + I) * z + (6 + 18 * I), z)
</code></pre>

<p>while, again, Maple and Mathematica will both happily give me the two Gaussian integers that solve this equation.</p>

<p>Is there an option or something that I'm missing?</p>