Dynamics of an Eco-Epidemic Predator–Prey Model Involving Fractional Derivatives with Power-Law and Mittag–Leffler Kernel

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.