A fractional model for predator-prey with omnivore

We consider the model of interaction of predator and prey with omnivore using three different waiting time distributions. The first waiting time is induced by the power law distribution which is the generator of Pareto statistics. The second waiting time is induced by exponential decay law with a particular property of Delta Dirac distribution when the fractional order tends to 1, this distribution is link to the Poison distribution. While the last waiting distribution, induced by the Mittag-Leffler distribution, presents a crossover from exponential to power law. For each model, we presented the conditions under which the existence of unique set of exact solutions is reached using the fixed-point Picard's method. Making use of a recent suggested numerical scheme, we solved each model numerically and some numerical simulations were generated for different values of fractional orders. We noticed a new attractor which can be considered as a combination of the Brownian motion and power law distribution in the model with the Atangana-Baleanu fractional derivative. With the aim to capture more attractors, we modified the model and presented also some numerical simulations. Our new model provides more attractors than the existing one even for fractional differential cases. We presented finally the Maximal Lyapunov exponent and the bifurcation diagrams. The comparative study shows that modeling with non-local and non-singular kernel fractional derivative leads to more attractors as this kernel is able to capture more physical problems.