Local Stability Dynamics of Equilibrium Points in Predator-Prey Models with Anti-Predator Behavior

This article describes the dynamics of local stability equilibrium point models of interaction between prey populations and their predators. The model involves response functions in the form of Holling type III and anti-predator behavior. The existence and stability of the equilibrium point of the model can be obtained by reviewing several cases. One of the factors that affect the existence and local stability of the model equilibrium point is the carrying capacity (k) parameter. If x3∗, y3∗  > 0 is a constant solution of the model and ∈ (0,x3∗), then there is a unique boundary equilibrium point Ek (k , 0). Whereas, if k ∈ (x4∗, y4∗], then Ek (k, 0) is unstable and E3 (x3∗, y3∗) is stable. Furthermore, if k ∈ ( x4∗, ∞), then Ek ( k, 0) remains stable and E4 (x4∗, y4∗) is unstable, but the stability of the equilibrium point E3 (x3∗, y3∗) is branching. The equilibrium point E3 (x3∗, y3∗) can be stable or unstable depending on all parameters involved in the model. Variations of k parameter values are given in numerical simulation to verify the results of the analysis. Numerical simulation indicates that if k = 0,92 then nontrivial equilibrium point Ek (0,92 ; 0) stable. If k = 0,93 then Ek (0,93 ; 0) unstable and E3∗(0,929; 0,00003) stable. If k = 23,94, then Ek (23,94 ; 0) and E3∗(0,929; 0,143) stable, but E4∗(23,93 ; 0,0005) unstable. If k = 38 then Ek(38,0) stable, but E3∗(0,929; 0,145) and E4∗(23,93 ; 0,739) unstable.Keywords: anti-predator behavior, carrying capacity, and holling type III.