DifferentialEquations.jl, Error: Exponentiation yielding a complex result requires a complex argument

General Usage
6

Thanks mauro3 for your relpy.

Is i correct that isoutofdomain=(u,p,t) -> any(x -> x < 0, u) ensure that all variables to be within the positive domain?

If so, is there a way to specify that u[2] > 0 and is free u[2]?