Impact of Fear on a Crowley–Martin Eco-Epidemiological Model with Prey Harvesting †

: In this paper, we develop a three-species food web model that incorporates the use of interactions between diseased predator–prey models. The logistically growing prey populations are susceptible and diseased prey. Prey populations are assumed to grow logistically in the absence of predators. We investigate the effect of fear on susceptible prey through infected prey populations. In Crowley–Martin-type interactions, it is assumed that interdependence between predators happens regardless of whether an individual predator is searching for prey or handling prey at the time. Also, the prey harvesting of susceptible and infected prey has been considered. The existence of all possible equilibrium points for biological systems has been established. The criteria for the local and global stability of equilibrium points are examined. Additionally, we look at Hopf-bifurcation analysis for the suggested model in relation to the existence of harvesting rate ( h 1 ) . Numerical simulations are provided in order to explain the phenomenon and comprehend the complex interactions between predators and prey.


Introduction
The predator-prey models developed by Lotka [1] and Volterra [2] are regarded as the earliest developments in contemporary mathematical ecology in coupled systems of nonlinear differential equations.Since Kermack and Mckendrick's pioneering work on SIRS, [3], epidemiological models have attracted much interest from researchers.Mathematical modeling of predator-prey interactions, known as "functional response," is among the most important factors in predator-prey population modeling.Crowley-Martin functional responses take into account both prey and predators.In the recent era, some eminent authors, [4][5][6][7], have studied to understand the importance and interactions of prey.To make the model system more realistic and feasible in the ecosystem, they incorporated some functional responses, i.e., Crowley-Martin-type functional responses.Kadhim and Azhar [8] use a type II Holling function to represent two disease types in a predator population model with a linear functional response.In [9], a nonlinear analysis of a discrete effects predator-prey model is investigated.Prey refuge and prey harvest [10,11], with ratio-dependent and Holling type II functional responses.Several investigations have been conducted on the dynamic behavior of Crowley-Martin diseased predator-prey models.To our knowledge, only a few researchers have looked into three-species preypredator models that take into account species interactions, including Crowley-Martin disease in prey populations.This study examines how fear affects a Crowley-Martin ecoepidemiological model with prey harvesting.The rest of the paper is structured as follows: In Section 2, we describe how the study's model formation was created.In Sections 3 and 4, we talk about boundary equilibrium points and their stability.In Section 5, the Hopf-bifurcation is the positive equilibrium point E * (u * , v * , w * ).Numerical simulations of the suggested model are examined in Section 6.The paper is concluded in Section 7, which also discusses the biological consequences of our mathematical findings.

Mathematical Model Formation
The model explains the relationship between the structure of the infected prey and the following equations.The proposed framework was used to discover a non-linear prey-predator mathematical model.
and the positive conditions are described as S 0 ≥ 0, I 0 ≥ 0 and P 0 ≥ 0. The Table 1 displays the specific biological meanings of the parameters.The condition for the impact of fear is F(ρ, i) = 1 1+ρi .This refers to the infected's fear effect on susceptible prey.Here, f is the amount of fear.It is appropriate to modify the variables as follows in order to decrease the number of systems (1) variables s = S K , i = I K , p = P K , and to consider the dimension time t = λKT.Now, we apply the following transformations 1) can be represented in dimensionless form using the above transformations.

Existence of Equilibrium Points
The model (2) exhibits the following equilibrium points based on observation: 1.
. Here, î is the unique positive root of the equation

Stability Analysis
In order to determine local stability around various equilibrium points, we compute the Jacobian matrix.At each given point (s, i, p), the Jacobian matrix is given by where, Theorem 1.The trivial equilibrium point E 0 (0, 0, 0) is always unstable.Proof.
Here, the eigenvalues of Hence, E 0 is locally asymptotically stable only if r < h 1 and unstable otherwise. where where, In the above characteristic equation, we obtain one of the eigenvalue is N 7 , which is negative as δ > bc î 1+η î + βc ŝ 1+η ŝ and the other two eigenvalues should be negative if X 11 > 0 and X 12 > 0.
Theorem 5. E * is locally asymptotically stable if Z 1 > 0, Z 3 > 0, and Z 1 Z 2 − Z 3 > 0. Proof. where, The characteristic equation is where According to the Routh-Hurwitz criterion, Z 1 , Z 3 , and Z 1 Z 2 − Z 3 must all be positive, and the characteristic of all the roots must be negative.Hence, E * is locally asymptotically stable.

Global Stability Analysis
Theorem 6.If E * is the endemic equilibrium point is globally asymptotically stable in G = {(s, i, p) : s > s * , i > i * and p > p * ) or s < s * , i < i * and p < p * }.
Differentiating L 1 with respect to t, we obtain ].
Now, we see that dL 1 dt ≤ 0. Whenever G = {(s, i, p) : s > s * , i > i * and p > p * ) or s < s * , i < i * and p < p * } and Consequently, for all solutions in G, L is a Lyapunov function.

Hopf-Bifurcation Analysis
Theorem 7. If the critical value for the bifurcation parameter h 1 is exceeded, the model (2) will experience the Hopf-bifurcation.The following Hopf-bifurcation requirements are present for h where S is the zeros of the characteristic equation corresponding to the non-negative equilibrium point.

Conclusions
We investigated a three-species food-web model that involved the use of the interactions between diseased prey-predator model.The local and global stability of ( 2) is used for each set of biologically possible equilibrium points in the model.It is used to modify the harvesting rate (h 1 ) and the level of fear (ρ) as control parameters.In addition, we investigated the stability analysis of the model ( 2) and studied the Hopf-bifurcation phenomenon.As a result, we found that modifying the harvesting rate h 1 significantly affects the stability of the system (2).The analytical and numerical findings demonstrate that the harvesting rate has a significant impact on every population.A decrease in the population of susceptible prey and an increase in infected prey population density are the effects of increasing the harvesting rate.This study shows the complex behavior of the proposed model.

Table 1 .
Biological representation of the model.