Conflict with states (DifferentialEquations and HybridSystems)

[schillic]

Using the model from https://github.com/JuliaReach/ReachabilityAnalysis.jl/blob/5279c0195c2954b96292ba34c853c12c3eaa5d3a/test/models/hybrid/thermostat.jl and initial_condition = X0:

using ReachabilityAnalysis
using HybridSystems, MathematicalSystems, LazySets, LinearAlgebra
using LazySets: HalfSpace  # resolve name-space conflicts with Polyhedra

function thermostat()
    c_a = 0.1

    # automaton structure
    automaton = LightAutomaton(2)

    # mode on
    A = hcat(c_a)
    inv = HalfSpace([1.0], 22.0)  # x ≤ 22
    # @vars x v
    # @set x ≤ 22
    m_on = ConstrainedLinearContinuousSystem(A, inv)

    # mode off
    A = hcat(-c_a)
    inv = HalfSpace([-1.0], -18.0)  # x ≥ 18
    m_off = ConstrainedLinearContinuousSystem(A, inv)

    # modes
    modes = [m_on, m_off]

    # transition from on to off
    add_transition!(automaton, 1, 2, 1)
    guard = HalfSpace([-1.0], -21.0)  # x ≥ 21
    t_on2off = ConstrainedIdentityMap(2, guard)

    # transition from off to on
    add_transition!(automaton, 2, 1, 2)
    guard = HalfSpace([1.0], 19.0)  # x ≤ 19
    t_off2on = ConstrainedIdentityMap(2, guard)

    # transition annotations
    resetmaps = [t_on2off, t_off2on]

    # switching
    switchings = [AutonomousSwitching()]

    ℋ = HybridSystem(automaton, modes, resetmaps, switchings)

    # initial condition
    X0 = Singleton([18.0])
    initial_condition = X0

    return InitialValueProblem(ℋ, initial_condition)
end;

S = thermostat();
T = 5.0;

sol = solve(S, T=T);  # works fine

If we instead first load DifferentialEquations, there is a conflict that leads to a crash:

julia> import DifferentialEquations
julia> using ReachabilityAnalysis

[...]

julia> sol = solve(S, T=T);
WARNING: both DifferentialEquations and HybridSystems export "states"; uses of it in module ReachabilityAnalysis must be qualified
ERROR: UndefVarError: states not defined