Fractional Order Liu System

using FractionalDiffEq, Plots

h=0.005
alpha = [0.95, 0.95, 0.95]
x0 = [0.2, 0, 0.5]
tspan = (0, 100)
function Liu!(du, u, p, t)
    a, b, c, e, n, m = 1, 2.5, 5, 1, 4, 4
    du[1] = -a*u[1]-e*u[2]^2
    du[2] = b*u[2]-n*u[1]*u[3]
    du[3] = -c*u[3]+m*u[1]*u[2]
end
prob = FODESystem(Liu!, alpha, x0, tspan)
sol = solve(prob, h, GL())

plot3d(sol, vars=(1,2,3), title="Fractional Order Liu System")

This page was generated using DemoCards.jl.