Consider the following predator-prey system:
\begin{cases} \dfrac{dN}{dT} &= rN\left(1-\dfrac{N}{K(T)}\right) -aN^\alpha P -\dfrac{qEN}{mE+nN},\\ \dfrac{dP}{dT} &=-d_1P+aeN^\alpha P. ...Proving positivity the solution of a Predator-Prey model – math.stackexchange.com
user1661695