# 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.*