(14) \(\left\{\begin{matrix}\dot{x}(t)=\sigma-\delta x(t)+\rho x(t)\frac{y(t)}{y(t)+n}-\mu x(t)y(t)-K_{x}M(t)x(t)\\ \dot{y}\left(t\right)=\alpha y\left(t\right)\left(1-\beta y\left(t\right)\right)-x\left(t\right)y\left(t\right)-K_{y}M\left(t\right)y\left(t\right)\text{\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ }\\ \dot{M}\left(t\right)=-\xi M\left(t\right)+v_{M}\left(t\right)\text{\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ }\\ \end{matrix}\right.\ \)
(15)
(16)