We consider a model of predator–prey interaction at fractional-order where the predation obeys the ratio-dependent functional response and the prey is linearly harvested. For the proposed model, we show the existence, uniqueness, non-negativity and boundedness of the solutions. Conditions for the existence of all possible equilibrium points and their stability criteria, both locally and globally, are also investigated. The local stability conditions are derived using the Magtinon’s theorem, while the global stability is proven by formulating an appropriate Lyapunov function. The occurrence of Hopf bifurcation around the interior point is also shown analytically. At the end, we implemented the Predictor–Corrector scheme to perform some numerical simulations.

In this paper, we consider a fractional-order eco-epidemic model based on the Rosenzweig–MacArthur predator–prey model. The model is derived by assuming that the prey may be infected by a disease. In order to take the memory effect into account, we apply two fractional differential operators, namely the Caputo fractional derivative (operator with power-law kernel) and the Atangana–Baleanu fractional derivative in the Caputo (ABC) sense (operator with Mittag–Leffler kernel). We take the same order of the fractional derivative in all equations for both senses to maintain the symmetry aspect. The existence and uniqueness of solutions of both eco-epidemic models (i.e., in the Caputo sense and in ABC sense) are established. Both models have the same equilibrium points, namely the trivial (origin) equilibrium point, the extinction of infected prey and predator point, the infected prey free point, the predator-free point and the co-existence point. For a model in the Caputo sense, we also show the non-negativity and boundedness of solution, perform the local and global stability analysis and establish the conditions for the existence of Hopf bifurcation. It is found that the trivial equilibrium point is a saddle point while other equilibrium points are conditionally asymptotically stable. The numerical simulations show that the solutions of the model in the Caputo sense strongly agree with analytical results. Furthermore, it is indicated numerically that the model in the ABC sense has quite similar dynamics as the model in the Caputo sense. The essential difference between the two models is the convergence rate to reach the stable equilibrium point. When a Hopf bifurcation occurs, the bifurcation points and the diameter of the limit cycles of both models are different. Moreover, we also observe a bistability phenomenon which disappears via Hopf bifurcation.

In this paper, a modified Leslie–Gower predator-prey model with Beddington–DeAngelis functional response and double Allee effect in the growth rate of a predator population is proposed. In order to consider memory effect on the proposed model, we employ the Caputo fractional-order derivative. We investigate the dynamic behaviors of the proposed model for both strong and weak Allee effect cases. The existence, uniqueness, non-negativity, and boundedness of the solution are discussed. Then, we determine the existing condition and local stability analysis of all possible equilibrium points. Necessary conditions for the existence of the Hopf bifurcation driven by the order of the fractional derivative are also determined analytically. Furthermore, by choosing a suitable Lyapunov function, we derive the sufficient conditions to ensure the global asymptotic stability for the predator extinction point for the strong Allee effect case as well as for the prey extinction point and the interior point for the weak Allee effect case. Finally, numerical simulations are shown to confirm the theoretical results and can explore more dynamical behaviors of the system, such as the bi-stability and forward bifurcation.

We consider a mathematical model to describe the interaction between predator and prey that includes predator cannibalism and refuge. We aim to study the dynamics and its long-term behavior of the proposed model, as well as to discuss the effects of crucial parameters associated with the model. We first show the boundedness and positivity of the solution of the model. Then, we study the existence and stability of all possible equilibrium points. The local stability of the model around each equilibrium point is studied via the linearized system, while the global stability is performed by defining a Lyapunov function. The model has four equilibrium points. It is found that the equilibrium point representing the extinction of both prey and predator populations is always unstable, while the other equilibrium points are conditionally stable. In addition, there is forward bifurcation phenomena that occur under certain condition. To support our analytical findings, we perform some numerical simulations.