Consider the system without prey harvesting, that is \begin{equation} \begin{split} x' &= x(a - cx -dy) \\ y' & = -y(b - ex) \end{split}\label{eq1} \end{equation} where all of $a, b, c, d, e$ are positive, and $a/c>b/e.$
The following simulation shows different solutions to the equations (\ref{eq1}) for different initial conditions and fixed values $d = 0.02$ and $e=0.02$.
Things to try: