Stability Analysis and Numerical Computation of the Fractional Predator–Prey Model with the Harvesting Rate

Abstract

In this work, a fractional predator-prey model with the harvesting rate is considered. Besides the existence and uniqueness of the solution to the model, local stability and global stability are experienced. A novel discretization depending on the numerical discretization of the Riemann–Liouville integral was introduced and the corresponding numerical discretization of the predator–prey fractional model was obtained. The net reproduction number ${\mathcal{R}}_0$ was obtained for the prediction and persistence of the disease. The dynamical behavior of the equilibria was examined by using the stability criteria. Furthermore, numerical simulations of the model were performed and their graphical representations are shown to support the numerical discretizations, to visualize the effectiveness of our theoretical results and to monitor the effect of arbitrary order derivative. In our investigations, the fractional operator is understood in the Caputo sense.