Fractional Order Van der Pol Oscillator

using FractionalDiffEq, Plots

h=0.005
alpha = [1.2, 0.8]
x0 = [0.2, -0.2]
tspan = (0, 60)
function VanderPol!(du, u, p, t)
    ϵ = 1
    du[1] = u[2]
    du[2] = -u[1]-ϵ*(u[1]^2-1)*u[2]
end
prob = FODESystem(VanderPol!, alpha, x0, tspan)
sol = solve(prob, h, GL())

plot(sol, vars=(1,2), title="Fractional Order Van der Pol Oscillator")

This page was generated using DemoCards.jl.