“…Therefore, if the density of the prey population falls below the threshold ET , that is, H(Z) < 0, then ε = 0, which indicates that the prey may avoid the predator via a habitat shift by moving to the refuge and the density of the predator will decrease [2]; if the density of the prey population increases and exceeds the threshold ET , that is, H(Z) > 0, then ε = 1, which indicates that the prey population may re-appear and once again become accessible to predators [2,16,17,24]. According to the above definition, the extended classical Lotka-Volterra model with a piecewise saturating function can be defined as the following Filippov system dx(t) dt = rx(t) − εbx(t)y(t) 1 + bhx(t) , dy(t) dt = εkbx(t)y(t) 1 + bhx(t) − δy(t),…”