This paper is concerned with a Holling-II stage-structured predator-prey system in which predators are divided into an immature and mature predators. The aim is to explore the impact of the prey's fear caused by the dread of mature predators in a prey-predator model including intraspecific competitions and prey shelters. The theoretical study includes the local and global stability analysis for the three equilibrium points of the system and shows the prey's fear may lead to improving the stability at the positive equilibrium point. A numerical analysis is given to ensure the accuracy of the theoretical outcomes and to testify the conditions of stability of the system near the non-trivial equilibrium points.

    The relationship between prey and predator populations is hypothesized and examined using a mathematical model. Predation fear, cannibalism among the prey population, and a refuge reliant on predators are predicted to occur. This study set out to look at the long-term behavior of the proposed model and the effects of its key elements. The solution properties of the model were investigated. All potential equilibrium points' existence and stability were looked at. The system's persistence requirements were established. What circumstances could lead to local bifurcation near equilibrium points was uncovered. Suitable Lyapunov functions are used to study the system's overall dynamics. Numerical simulations were conducted to verify the

An ecological model consisting of prey-predator system involving the prey’s fear is proposed and studied. It is assumed that the predator species consumed the prey according to prey square root type of functional response. The existence, uniqueness and boundedness of the solution are examined. All the possible equilibrium points are determined. The stability analysis of these points is investigated along with the persistence of the system. The local bifurcation analysis is carried out. Finally, this paper is ended with a numerical simulation to understand the global dynamics of the system.

A mathematical eco-epidemiological model consisting of harvested prey–predator system involving fear and disease in the prey population is formulated and studied. The prey population is supposed to be separated into two groups: susceptible and infected. The susceptible prey grows logistically, whereas the infected prey cannot reproduce and instead competes for the environment’s carrying capacity. Furthermore, the disease is transferred through contact from infected to susceptible individuals, and there is no inherited transmission. The existence, positivity, and boundedness of the model’s solution are discussed. The local stability analysis is carried out. The persistence requirements are established. The global behavior of th

Using a mathematical model to simulate the interaction between prey and predator was suggested and researched. It was believed that the model would entail predator cannibalism and constant refuge in the predator population, while the prey population would experience predation fear and need for a predator-dependent refuge. This study aimed to examine the proposed model's long-term behavior and explore the effects of the model's key parameters. The model's solution was demonstrated to be limited and positive. All potential equilibrium points' existence and stability were tested. When possible, the appropriate Lyapunov function was utilized to demonstrate the equilibrium points' overall stability. The system's persistence requirements were spe

        In this paper, the effects of prey’s fear on the dynamics of the prey, predator, and scavenger system incorporating a prey refuge with the linear type of functional response were studied theoretically as well as numerically approach. The local and global stabilities of all possible equilibrium points are investigated. The persistence conditions of the model are established. the local bifurcation analysis around the equilibrium points, as well as the Hopf bifurcation near the positive equilibrium point, are discussed and analyzed. Finally, numerical simulations are carried out, and the obtained trajectories are drowned using the application of Matlab version (6) to explain our found analytical