title: Orbit Diagrams with DynamicalSystems Initial Condition not Working
body: <p>I am trying to using the DynamicalSystems julia package to plot an orbit diagram of a discrete dynamical system governed by the following equation:</p>
  
  <p><code>y_{n+1} = 1- by_n^2</code></p>
  
  <p>I'm able to get it to plot however I can't seem to get the initial condition to work, I am trying to use the initial condition of <code>y_0 = -0.68</code>. </p>
  
  <p>Here is my code:</p>
  
  <pre><code>using DynamicalSystems
  using PyPlot
  
  
  function eom(dx, x, p, n)
      dx[1] = 1 - p[1] * x[1]^2
  end
  
  ds = DiscreteDynamicalSystem(eom, [-0.68], [2.0, 0.0])
  
  i = 1
  pvalues = 0.0:0.001:2.0
  n = 2000
  Ttr = 2000
  p_index = 1
  output = orbitdiagram(ds, i, p_index, pvalues; n = n,
                        Ttr = Ttr, u0 = get_state(ds))
  
  L = length(pvalues)
  x = Vector{Float64}(undef, n*L)
  y = copy(x)
  for j in 1:L
      x[(1 + (j-1)*n):j*n] .= pvalues[j]
      y[(1 + (j-1)*n):j*n] .= output[j]
  end
  
  figure()
  PyPlot.title("total points: $(L*n)")
  plot(x, y, ls = "None", ms = 0.5, color = "black", marker = "o", alpha = 0.05)
  xlim(pvalues[1], pvalues[end]);
  xlabel("\$b\$"); ylabel("\$y\$")
  tight_layout()
  </code></pre>
  
  <p>And here is the figure it generates:
  <a href="https://i.stack.imgur.com/jkxTA.png" rel="nofollow noreferrer"><img src="https://i.stack.imgur.com/jkxTA.png" alt="enter image description here"></a></p>
  
  <p>If possible, could someone also tell me if there is a way to get the bifurcation points using DynamicalSystems.</p>