Fractional Order Duffing System

using FractionalDiffEq, Plots

h=0.005
alpha = [0.9, 1]
x0 = [0.21, 0.31]
tapan = (0, 100)
function Duffing!(du, u, p, t)
    α, δ, ω = 0.15, 0.3, 1
    du[1] = u[2]
    du[2] = u[1]-α*u[2]-u[1]^3+δ*cos(ω*t)
end
prob = FODESystem(Duffing!, alpha, x0, tspan)
sol = solve(prob, h, GL())

plot(sol, vars=(1,2), title="Fractional Order Duffing System")

This page was generated using DemoCards.jl.