Discrete Dynamical Systems and Beautiful Graphs

Sharing my SciML homework

using DifferentialEquations
using Flux
using Plots

function calc_attractor!(f, x0, r, warmup, num_attract)
    x = Vector{typeof(x0)}(undef,warmup+num_attract)
    x[1] = x0   
    for i in 1:warmup+num_attract-1
        x[i+1] = f(x[i], r)
        display(i)
    end
    x[warmup:warmup+num_attract]
end

f(x,r) = r.*x.*(1.0.-x)
r = 2.9:0.001:4

sol = calc_attractor!(f, [0.25], r, 400, 150)
plot(r, sol, label=["x" "r"])