# Fractional Order Lotka Volterra System

```
using FractionalDiffEq, Plots
h=0.005
alpha = [0.95, 0.95, 0.95]
x0 = [1, 1.4, 1]
tspan = (0, 50)
function LotkaVolterra!(du, u, p, t)
a, b, c, d, e, p, s = 1, 1, 1, 1, 2, 3, 2.7
du[1] = a*u[1]+e*u[1]^2-b*u[1]*u[2]-s*u[3]*u[1]^2
du[2] = -c*u[2]+d*u[1]*u[2]
du[3] = -p*u[3]+s*u[3]*u[1]^2
end
prob = FODESystem(LotkaVolterra!, alpha, x0, tspan)
sol = solve(prob, h, GL())
plot(sol, vars=(1,2,3), title="Fractional Order Lotka Volterra System")
```

*This page was generated using DemoCards.jl.*