title: Scipy - Non-linear Equations System with linear constraints (beginner)
body: <p>I have seen this <a href="https://stackoverflow.com/questions/8739227/how-to-solve-a-pair-of-nonlinear-equations-using-python">amazing example</a>.
  But I need to solve system with boundaries on X and F, for example:</p>
  
  <pre><code>    f1 = x+y^2 = 0
      f2 = e^x+ xy = 0
  -5.5&lt; x &lt;0.18
  2.1&lt; y &lt; 10.6
     # 0.15&lt; f1 &lt;20.5 - not useful for this example
     # -10.5&lt; f2 &lt; -0.16 - not useful for this example
  </code></pre>
  
  <p>How could I set this boundary constrains to fsolve() of scipy? Or may be there is some other method?
  Would You give me a Simple code example?</p>