Numerical bifurcation of predator-prey fractional differential equations with a constant rate harvesting

In this article saddle and Hopf bifurcation points of predator-prey fractional differential equations system with a constant rate harvesting are investigated. The numerical results based on Grunwald-Letnikov discretization for fractional differential equations together with the Mickens' non-standard discretization method agree with those found by the corresponding ordinary differential equation system.