title: Scipy's fsolve alternatives?
body: <p>I'm trying to solve numerically a system of non-linear equations:</p>
  
  <pre><code>def func(p):
      x, f = p
      return (math.exp(-x/O)-f,
              L - L*((1 - math.exp(-x/O))**W) - x*math.exp(-x/O))
  </code></pre>
  
  <p>and I'm currently using scipy.fsolve for it in the following way:</p>
  
  <pre><code>x, f = fsolve(func, (10, 0.2))
  </code></pre>
  
  <p>I'm sure my way of using fsolve is correct: it works perfectly for certain range of parameters. However, it totally fails for another (e.g., O=8, L=1.67, W=8) with the following error:</p>
  
  <pre><code>RuntimeWarning: The number of calls to function has reached maxfev = 600.
  </code></pre>
  
  <p>I'm sure it is well solvable -- there are (at least) matlab tools, which do it. Is there anything I'm doing wrong or is there any other solver I could try?</p>
  
  <p>Thanks in advance for any hints!</p>