Hi,
You are right, something is wrong with the predictor from the transcitical bifurcation. In case you need it, this works:
parameters2 = @set parameters.g=br1.specialpoint[2].param + 0.01
solution2 = br1.specialpoint[2].x
br3, = continuation(f, solution2, parameters2, (@lens _.g), options,
plot = true,
recordFromSolution = (x,p) -> (C = x[2], M = x[1]),
plotSolution = (x, p; k...) -> begin
plot!(br1, subplot = 1)
end,
)