HomotopyContinuation: Duffing Revisited: Plotting Solution from Multiple Paths

[Domenico_Lahaye]

Using

Avec = Vector{Float64}[]
Bvec = Vector{Float64}[]
Cvec = Vector{Float64}[]
for (i,ri) in enumerate(myrealsols)
    Aveci = []
    Bveci = [] 
    Cveci = [] 
    for (j,rj) in enumerate(ri)
        Avecj = rj[1]
        Bvecj = rj[2] 
        Cvecj = sqrt(Avecj^2+Bvecj^2)
        push!(Aveci,Avecj) 
        push!(Bveci,Bvecj) 
        push!(Cveci,Cvecj) 
    end 
    push!(Avec,Aveci)
    push!(Bvec,Bveci)
    push!(Cvec,Cveci)
end 

X = [ i*ones(length(res)) for (i,res) in pairs(myrealsols)];

scatter(X, Cvec, ms=5, legend =:false)

I now obtain

matching well with the previous solution.