title: Python: Finding multiple roots of nonlinear equation
body: <p>Assume the following function:</p>
  
  <p><code>f(x) = x * cos(x-4)</code></p>
  
  <p>With <code>x = [-2.5, 2.5]</code> this function crosses <code>0</code> at <code>f(0) = 0</code> and <code>f(-0.71238898) = 0</code>.</p>
  
  <p>This was determined with the following code:</p>
  
  <pre><code>import math
  from scipy.optimize import fsolve
  def func(x):
      return x*math.cos(x-4)
  x0 = fsolve(func, 0.0)
  # returns [0.]
  x0 = fsolve(func, -0.75)
  # returns [-0.71238898]
  </code></pre>
  
  <p>What is the proper way to use <code>fzero</code> (or any other Python root finder) to find both roots in one call? Is there a different <code>scipy</code> function that does this?</p>
  
  <p><a href="http://docs.scipy.org/doc/scipy/reference/generated/scipy.optimize.fsolve.html#scipy.optimize.fsolve" rel="noreferrer"><code>fzero</code> reference</a></p>