title: Creating a bifurcation diagram in python
body: <p>I am trying to come up with a code that will allow me to plot a diagram for period doubling bifurcation. </p>
  
  <p>I am using the equation <code>x = rx − 1(1 − x)</code>, and am trying to model it with r values from 0.5 to 4.  Here is code that I am working with</p>
  
  <pre><code>startr = 0.5
  finalr = 4
  max_time = 200
  x = [0.1]
  r= np.linspace(.5,4,200)
  
  for n in range(0,200):
  
      x = np.append(r * x[n] * (1-x[n]))
  
  
  plt.plot(x, label='x');
  plt.xlabel('t');
  </code></pre>
  
  <p>This keeps getting kicked out</p>
  
  <pre><code>TypeError: append() missing 1 required positional argument: 'values'
  </code></pre>