Insight into Caputo Fractional-Order Extension of Lotka–Volterra Model with Emphasis on Immigration Effect

For over a hundred years, the model of prey–predator has been handled extensively. It has been proved that Lotka–Volterra (LV) models are not stable from all previous studies. In these models, there is divergent extinction of one species or cyclic oscillation. There is no solution for asymptotic stability for the prey–predator models. In this paper, a Caputo fractional-order one prey-two predators’ model with an immigration effect is investigated. Then, the effect of adding a small immigration term to the prey or predators’ populations is studied. The extension is both the inclusion of a second predator and the use of Caputo derivatives. Although LV systems have been extensively investigated, there are no previous studies on our proposed nonlinear modifications. The Grünwald–Letnikov method was applied to find numerical solutions for the proposed model to validate our findings. Several different cases were studied to reach confirmed conclusions. The effect of the variation of fractional-order derivative is handled in this study. From all the results, the effect of adding the small immigration terms is positive for stability. Hence, it can be concluded that the proposed small immigration terms can stabilize the natural populations of the one prey-two predators’ model. An optimal control strategy for the proposed model is shown at the end of the paper.