title: sympy solve() unable to find solution
body: <p>I have a fairly basic code right now that's supposed to solve the equation shown below but it's not converging to a solution. It just hangs on a pulsing "_"</p>
  
  <pre><code>from sympy import *
  
  gamma = 1.4
  M_a = 1.0
  
  y = Symbol('y', real=True) 
  eqn = Eq((1.0/((gamma/2.0)*(y**2.0))) * ((((1.0 + ((gamma - 1.0)/2.0)*(y**2.0))/(1.0 + ((gamma - 1.0)/2.0)*(M_a**2.0)))**(gamma/(gamma - 1.0))) - 1.0) ,(-0.5704/((1.0 - (y**2.0))**(1.0/2.0))))
  print solve(eqn, y)
  </code></pre>
  
  <p>simply printing the equation yields</p>
  
  <pre><code>1.42857142857143*y**(-2.0)*((0.166666666666667*y**2.0 + 0.833333333333333)**3.5- 1.0) == -0.5704*(-y**2.0 + 1.0)**(-0.5)
  </code></pre>
  
  <p>which when plugged into wolfram or maple yields the correct solution. ~= 0.696256</p>
  
  <p>So I'm trying to figure out why sympy is unable to solve the equation.</p>
  
  <p>The equation should look like this <a href="http://i.imgur.com/fk7wPQM.png" rel="nofollow">Picture</a></p>
  
  <p>If sympy can not be used to solve this equation what can I use instead?</p>
  
  <p>Thanks</p>
  
  <p>Phil</p>