Romain Veltz

I suggest that you first use Natural() continuation to get à feeling. I think the problem is from the scaling of k. I would write k = 1e6 * k2...

You can pass the jacobian to a `BifurcationProblem` with `; J = ...)`

I am surprised Natural did not work, it is a basic Newton method

Oh are using complex numbers? You can only use reals

You can look at `br.branch` or do `plot(br)` to check Beta(k). I am surprised natural does not work. You should check the function newton in verbose mode. Once this works...

> suggestions for overcoming this difficulty? Yes I think you should decrease the option `maxIter`, somehow Newton works too well. Indeed it found a solution far from the guess as...

> Regarding Natural, here is the output with verbosity = 3, but I'm not sure how to modify the code based on this result: I am wondering what happens if...

That is in NewtonPar

can you share the function F please so I can help more?

can we progress on this?