title: How to perform non-linear optimization with scipy/numpy or sympy?
body: <p>I am trying to find the optimal solution to the follow system of equations in Python:</p>
  
  <pre><code>(x-x1)^2 + (y-y1)^2 - r1^2 = 0
  (x-x2)^2 + (y-y2)^2 - r2^2 = 0
  (x-x3)^2 + (y-y3)^2 - r3^2 = 0
  </code></pre>
  
  <p>Given the values a point(x,y) and a radius (r):</p>
  
  <pre><code>x1, y1, r1 = (0, 0, 0.88)
  x2, y2, r2 = (2, 0, 1)
  x3, y3, r3 = (0, 2, 0.75)
  </code></pre>
  
  <p>What is the best way to find the optimal solution for the point (x,y)
  Using the above example it would be:<br>
   ~ (1, 1)</p>