Abstract
This paper investigates a fractional-order mathematical model of predator-prey interaction in the ecology considering the fear of the prey, which is generated in addition by competition of two prey species, to the predator that is in cooperation with its species to hunt the preys. At first, we show that the system has non-negative solutions. The existence and uniqueness of the established fractional-order differential equation system were proven using the Lipschitz Criteria. In applying the theory of Routh-Hurwitz Criteria, we determine the stability of the equilibria based on specific conditions. The discretization of the fractional-order system provides us information to show that the system undergoes Neimark-Sacker Bifurcation. In the end, a series of numerical simulations are conducted to verify the theoretical part of the study and authenticate the effect of fear and fractional order on our model's behavior.