I have the set of nonlinear differential equations describing a system I modeled for my research (spread of epidemics or information for instance):

$$\begin{array}{rl} \dot{p}(t) &= \gamma r(t)-u(t)p(t)\\ \dot{q}(t) &= u(t)p(t)-\delta q(t)\\ \dot{r}(t) &= \delta q(t)-\gamma r(t)\end{array}$$

Since $\dot{p}+\dot{q}+\dot{r}=0$, I only consider the last two equations:

$$\begin{array}{rl} \dot{q}(t) &= u(t)(1-q(t)-r(t))-\delta q(t)\\ \dot{r}(t) &= \delta q(t)-\gamma r(t)\end{array}$$

with input $$u(t)\in \{{u}:u\textrm{ is Lebesgue integrable, } 0\leq u(t)\leq u_{\max}\}$$ and the set $$\Omega=\{(q,r)\vert 0 \leq {q},\,{r}\leq 1, {q}+{r}\leq 1\}$$

My goal is to show that $\Omega$ is positively invariant with respect to the nonlinear dynamical system expressed above, i.e., a solution starting from any initial condition $(q,r)\in\Omega$ remains confined to $\Omega$.

Edit: $p(t),\,q(t),$ and $r(t)$ represent probabilites that sum up to $1$: ${p}(t)+{q}(t)+{r}(t)=1$.