All in all:
DO = DeflationOperator( 2, # Algorithm parameter required when using deflated continuation, see BifurcationKit documentation.
dot, # Algorithm parameter required when using deflated continuation, see BifurcationKit documentation.
1., # Algorithm parameter required when using deflated continuation, see BifurcationKit documentation.
[fill(0.,length(rn.states))]); # Guess(es) of the fixed point for the initial parameter set. Do not need to be exact.
params_input = setindex!(copy(params),p_span[1],p_idx) # The input parameter values have to start at the first index of our parameter span.
branches, = continuation(F, J, params_input, (@lens _[p_idx]), opts , DO, # Gives our input.
verbosity = 0, plot=false, # We do not want to display, or plot, intermediary results.
recordFromSolution = (x, p) -> x[plot_var_idx], # How we wish to print the output in the diagram. Here we simply want the value of the target varriable.
perturbSolution = (x,p,id) -> (x .+ 0.8 .* rand(length(x))), # Parameter for the continuation method, see BifurcationKit documentation.
callbackN = BifurcationKit.cbMaxNorm(1e7), # Parameter for the continuation method, see BifurcationKit documentation.
)
plot(branches...)