# ## Brusselator

# ## Model

# The Brusselator is a two ODEs theoretical model for autocatalytic reactions,
# where $x$ and $y$ are chemical concentrations:


module brusselator      #jl

using ReachabilityAnalysis

const A = 1.0
const B = 1.5
const B1 = B + 1

@taylorize function brusselator!(du, u, p, t)
    x, y = u
    x² = x * x
    aux = x² * y
    du[1] = A + aux - B1*x
    du[2] = B*x - aux
    return du
end

function model(X0)
    S = @system(x' = brusselator!(x), dim:2)
    return IVP(S, X0)
end

end #jl

# ## Reachability settings

# ## Results



# ## References
