I'm trying to create a bifurcation plot for a driven damped pendulum. In particular, I'm trying to recreate the plot found in Taylor's 'Classical Mechanics' (page 484) for a driving strength $\gamma$ in the range $1.060 \leq \gamma \leq 1.087$:
I believe I have the code down and I can reproduce the main properties of the plot, however for $\gamma \geq 1.081$ the points shift down a couple of units which causes a large discontinuity. I'd appreciate any help on the matter. I've attached an image of my result below.