key: cord-1003439-kjjljbl5
authors: Abdo, Mohammed S.; Panchal, Satish K.; Shah, Kamal; Abdeljawad, Thabet
title: Existence theory and numerical analysis of three species prey–predator model under Mittag-Leffler power law
date: 2020-05-27
journal: Adv Differ Equ
DOI: 10.1186/s13662-020-02709-7
sha: ba34d1c92b3b8a9e464949ffa971612ea11e438a
doc_id: 1003439
cord_uid: kjjljbl5

In this manuscript, the fractional Atangana–Baleanu–Caputo model of prey and predator is studied theoretically and numerically. The existence and Ulam–Hyers stability results are obtained by applying fixed point theory and nonlinear analysis. The approximation solutions for the considered model are discussed via the fractional Adams Bashforth method. Moreover, the behavior of the solution to the given model is explained by graphical representations through the numerical simulations. The obtained results play an important role in developing the theory of fractional analytical dynamic of many biological systems.

A predator-prey model is a two-component system, where one of them lives at the expense of the other. A diversity of mathematical techniques is applied at modeling a predator-prey system due to numerous factors that may affect its evolution. In this regard there have been introduced some models in [1] [2] [3] [4] [5] [6] [7] in which the first model, which regards in a specified way only substantial phenomena (gluttony and fertility), is of the type

where p(t) and q(t) are the number of prey, the number of predators, respectively, and a 1 , a 2 and a 3 are the average of death of predators, the measurement of the tendency of prey to predation, and the predatory capability, respectively. The model (1) has a unique solution. However, the solutions of (1) are not structurally stable w.r.t. perturbation of the initial conditions. Inside the restricted scope of quadratic differential equations, those which cover competition and also predation, must be slightly more realistic. A second model with competition within preys is formulated as ⎧ ⎨ ⎩ p (t) = p(t)[a 1a 2 q(t) + a 5 p(t)],

where a 1 a 3 > a 4 a 5 , a 5 > 0 describes the competition of the prey. According to biologically sensible hypotheses, there exists a unique positive solution of model (2) which is asymptotically stable.

Dai, and Zhao investigated the dynamic complexities of a predator-prey model with state dependent on impulsive influences as

The authors used the analogue of the Poincaré norm to obtain the existence and stability of the model (3) . For details, see [8] . The following dynamic model, which addresses the case of predatory prey with disease, was analyzed by Das et al. [9] :

dy dt = (1 -y+x k )r 2 y + a 2 xyα 2 x(z + w)my,

where x(0), y(0), z(0), w(0) > 0. The authors showed that the model is globally stable on every side of the internal equilibrium point according to certain standard conditions. So, their analysis shows that the force of infection and predation average are the main parameters on the dynamics of the model. Fractional calculus deals with differentiation and integration involving fractional order, which is advantageous over the ordinary integer order in the explanation of real-world problems, as also in the modeling of real phenomena due to characterization of memory and hereditary properties [10, 11] . Further, the integer-order derivative does not describe the dynamics between two various points. Various types of fractional-order or nonlocal derivatives were proposed in the present literature to deal with the reduction of a traditional derivative. For instance, based on a power-law, Riemann-Liouville introduced the idea of a fractional derivative. Afterwards Caputo-Fabrizio in [12] have proposed a new fractional derivative utilizing the exponential kernel. This derivative has a few problems related to the locality of the kernel. Newly, to overcome Caputo-Fabrizio's problem, Atangana and Baleanu (AB) in [13] have proposed a new modified version of a fractional derivative with the aid of a generalized Mittag-Leffler function (MLF) as a nonsingular kernel and being nonlocal. Since the generalized MLF is used as the kernel it is guaranteed to have no singularity. Furthermore, the AB fractional derivative supplies a description of memory as discussed in [14] [15] [16] [17] [18] [19] [20] .

Most of the published work describes the mathematical system of predators and prey as a problem of Cauchy type of a system of classical differential equations [21] [22] [23] [24] [25] . However, recently, there has been great interest in studying the behavior of the solution for some biological systems using fractional differential equations involving the Atangana-Baleanu operator by several authors for the purpose of investigating several real-world systems and modeling infectious diseases; see [26] [27] [28] [29] [30] [31] [32] [33] [34] [35] [36] . Some fractional-order models have been investigated via the new operators recently. For instance its use has been suggested for the dynamics of smoking in [32] . Along the same line, the transference model for the Ebola virus together with AB operator was studied in [31] . A fractional-order model of leptospirosis infection was considered in [26] . The dynamical behavior of coronavirus (COVID-19) epidemic infection model through the ABC derivative has been studied in [33] . Also, the existence results and analytic solutions of fractional-order dynamics of COVID-19 with ABC derivative has been obtained in [34] . There is no literature available on prey-predator fractional models with three species under the aforesaid derivative. Just some fractional models have been found in the previous years; however, they have been confined to a standard fractional derivative. Furthermore, in the presence of the mentioned derivatives, recently some fruitful results have been published in [37] [38] [39] .

Due to the success of this operator in modeling the biological systems and infectious diseases, we have studied the dynamical behavior of the mathematical model which describes three prey-predator species by a nonlocal Atangana-Baleanu-Caputo (ABC) derivative operator with 0 < α ≤ 1 as

with the initial conditions

where ABC D α 0 + (·) is the ABC fractional derivative of order α, P 0 is the initial population density of prey, S 0 is the initial population density of susceptible predator, and I 0 initial population density infected predator. Here a denotes the saturation constant whereas susceptible predators threaten the prey, b is a search rate of the prey across a susceptible predator, c is the conversion rate of the susceptible predator due to prey, and d is the disease transmission coefficient. The symbol k represents the carrying capacity of the prey population, the proportionality constant is denoted by b 1 , the growth rate of the prey population is represented as r 1 . In the proposed model, m and n define the death rate of sensitive predator and death rate of the infected predator, respectively. Further, we ramark that the right hand sides of our considered model (4) under ABC fractional derivtive are assumed to vanish at zero, (for details, see Theorem 3.1 in [28] ).

The main aim of the paper is to demonstrate the existence, uniqueness and Ulam stability of the solution for the model (4)-(5) by using the Picard and fixed point techniques. Moreover, the numerical simulations via the fractional version of the Adams Bashforth technique to approximate the ABC fractional operator are performed. Graphical presentations are also given of the numerical results. This paper is organized as follows: Sect. 1 presents an introduction which contains a survey of the literature. Section 2 consists of some foundational preliminaries related to fractional calculus and nonlinear analysis. The existence and Ulam stability results on a proposed model are obtained in Sects. 3, 4. The numerical solution and numerical simulations of the model at hand are presented in Sect. 5.

For the next analysis, let

Definition 1 ([13] ) Let α ∈ (0, 1] and σ ∈ H 1 (0, T). Then the left-sided ABC fractional derivative with the lower limit zero of order α for a function σ is defined by

where ABC[α] is known as the normalization function which is defined as ABC[α] = α 2-α , 0 < α ≤ 1 and satisfies the result ABC(0) = ABC(1) = 1, and E α is called the Mittag-Leffler function defined by the series

here Re(α) > 0 and Γ (·) is a gamma function.

Let α ∈ (0, 1] and σ ∈ L 1 (0, T). Then the left-sided AB fractional integral with the lower limit zero of order α for a function σ is defined by

L ABC D α 0 + σ (t) = ABC[α] s α (1 -α) + α s α L σ (t) -s α-1 σ (0) .

The solution of the proposed problem for α ∈ (0, 1]

is given by

here κ is the Lipschitz constant for Π . If κ < 1 we say that Π is a contraction. 

Now, we address the existence and uniqueness results of the model (4)-(5) by utilizing the fixed point technique. Let us reformulate model (4) in the appropriate form

where

Utilizing Lemma 1, the model (8) can be turned to the fractional integral equation in the sense of AB fractional integral as follows:

Theorem 2 The kernels W ( = 1, 2, 3) agree with the contraction and Lipchitz conditions if there exists a constant L such that 0 ≤ L < 1, = 1, 2, 3.

For W 1 , let P and P * be two functions, then we have

where L 1 := (b 1 + b 1 k (A 1 + A * 1 ) + abC 1 + r 1 D 1 ), and P , P * , S , I are functions bounded by the constants A 1 , A * 1 , C 1 , D 1 , respectively. Consequently

Obviously, the Lipschitz condition is verified for W 1 . Besides, W 1 leads to a contraction due to 0 ≤ L 1 < 1. Likewise, we can show that W 2 and W 3 admit the contraction and Lipschitz condition, i.e.,

where L 2 := ( cb a+A 1 A 1 + dD 1 + m(C 1 + C * 1 )) and L 3 := (n + d + ckA 1 ). (11)-(13) hold. If

Then the solution of the fractional model given in (4)-(5) exists and is unique.

The initial conditions and the recurrence form of the model (10) are, respectively,

The successive difference between the terms is defined as

Clearly

Taking the norm of Eqs. (15) , it follows from the conditions (11)-(13) that

Let us consider P, S and I as bounded functions that comply with the Lipschitz condition. It follows from Eqs. (16) and (17) that

This shows the existence for the solutions. Moreover, to prove that Eqs. (18) are solutions for the model (4)-(5), we consider

On using recessive techniques, we get

As n → ∞, M 1n (t) → 0. In a similar way, we conclude that M 2n (t) and M 3n (t) tends to 0. Next, we address the uniqueness of the solution to the proposed mode (4)- (5) . To this end, let P * (t), S * (t) and I * (t) be other solutions. Then

It means that

From our hypothesis

It follows that P(t)-P * (t) = 0. Likewise, we conclude that S(t)-S * (t) = 0 and I(t)-I * (t) = 0.

For the notion of Ulam stability, see [42, 43] . The aforesaid stability has been scrutinized for classical fractional derivatives in many of the research articles; we refer to some of them like [44] [45] [46] [47] . Additionally, since stability is a prerequisite in respect of approximate solution, we endeavor on Ulam type stability for the model (4) via using nonlinear functional analysis.

Ulam-Hyers stable if there exists λ = max(λ 1 , λ 2 , λ 3 ) > 0 and = max( 1 , 1 , 1 ) > 0, for each P, S, I ∈ E × E × E, with the following inequalities:

then there exists (P, S, I) ∈ E × E × E satisfying the coupled system (4) with the following initial conditions:

such that ( P, S, I) -(P, S, I) Ω ≤ λ .

Remark 1 Consider a small perturbation g 1 ∈ C[0, T] that depends only on the solution such that g 1 (0) = 0 with the following properties: 1 |g 1 (t)| ≤ 1 , for t ∈ [0, T] and 1 > 0.

2 Furthermore, one has

Note that we will only discuss the first equation from the proposed system and the rest of the equations are similar in technique, i.e. P -P E ≤ λ 1 1 .

satisfies the relation

where P g 1 (t) is a solution of (21), P(t) satisfies (19-a) and κ := ( Γ (α)-Γ (α+1)+T α ABC[α]Γ (α) ).

Proof 3 Thanks to Remark 1, and Lemma 1, the solution of (21) is given by

Also, we have

It follows from Remark 1 that (11), the system (4)-(5) will be Ulam-Hyers stable in Ω.

Let P ∈ E be the solution of the inequality (19-a) and the function P ∈ E be a unique solution of Eq. (4-a) with the condition

That is,

Due to (22) , P 0 = P 0 . Hence Eq. (23) becomes

Thus by condition (11) and Lemma 2, we obtain P(t) -P(t) ≤ P(t) -P g 1 (t) + P g 1 (t) -P(t)

Γ (α) )L 1 < 1. For λ 1 = 2κ 1-Λ 1 , we get P -P E ≤ λ 1 1 . Similarly, we conclude that S -S E ≤ λ 2 2 , and I -I E ≤ λ 3 3 , where λ = 2κ 1-Λ ( = 2, 3). For some , λ > 0, ( P, S, I) -(P, S, I) Ω ≤ λ .

Hence the model (4)-(5) is Ulam-Hyers stable.

In this part, we give approximation solutions of the ABC fractional model (4)- (5) . Then the numerical simulations are acquired via the suggested scheme. To this aim, we employ the modified fractional version for AMB [48] to approximate the fractional integral in the AB sense. To procure an iterative scheme, we go ahead with the first equation of the model (10) as follows:

Set t = t n+1 , for n = 0, 1, 2, . . . , it follows that

Now, we approximate the function W 1 (θ , P) on the interval [t , t +1 ] through the interpolation polynomial as follows:

which implies

Now, we compute the integrals I -1,α and I ,α as follows:

and

Put t = , we get

α(α + 1) -2(α + 1)(n -) α + (α + 1)(n + 1 -) α -(n -) α+1 + (n + 1 -) α+1 = α+1 α(α + 1) (n -) α -2(α + 1) -(n -) + (n + 1 -) α (α + 1 + n + 1 -) = α+1 α(α + 1) (n + 1 -) α (n -+ 2 + α) -(n -) α (n -+ 2 + 2α) (26) and

Substituting (26) and (27) into (25), we get

Similarly

Γ (α + 2) α (n + 1 -) α (n -+ 2 + α) -(n -) α (n -+ 2 + 2α) -W 2 (t -1 , S(t -1 )) Γ (α + 2) α (n + 1 -) α+1 -(n -) α (n -+ 1 + α) (29) and

Now, to present the numerical simulations of the ABC fractional model (4)-(5), we apply the iterative solution contained in (28)- (30) . Take the time range up to 100 units. The numerical values of the parameters applied in the simulations are specified in Table 1 .

The graphical representations of numerical solution for species P, S, I at various fractional orders, α = 0.4, 0.6, 0.8, 1.0, of the considered model (4) are given in Figs. 1-3 , respectively. From Figs. 1-3, we observe that species I depends on species P and S. Therefore the papulation density of specie P and S gradually go on decreasing with different rate due to the fractional order in the first 50 days. The lower the fractional order, the faster the decay rate and hence the more rapidly the system becomes stable and vice versa. On the other hand, the species I is going on increasing with different rate, the lower the order the slower is the growth rate until it becomes stable and vice versa. The fractional order greatly Figure 3 Graphical representation of numerical solution for specie I at various fractional orders of the considered model (4) affects the stability of the system and also provides the global nature of the dynamics of the considered model.

Here we claim that the established numerical technique is powerful and converges for the ABC fractional derivative. Meanwhile the iterative techniques like perturbation and decomposition methods do not show the perfect behavior for the said derivatives for approximate solutions in many cases.

In this paper, the population density model of prey and sensitive predatory and infected predatory has been studied theoretically and numerically. Theoretically, the existence and stability results in the sense of Ulam-Hyers have been obtained through the help of fixed point theory and nonlinear analysis. Numerically, the approximation solution of the ABC fractional model (4) has been discussed via the use of a fractional Adam Bashforth method. Moreover, the behavior of the solutions of the model (4) has also been explained through graphs using some numerical values for the parameter. The obtained results play an important role in developing the theory of fractional analytical dynamics of various phenomena of real-world problems.

Théorie mathématique de la lutte pour la vie

Elements of Physical Biology

Sulla theoria di Volterra della lotta per l'esistenza

Mathematical Biology. Harrap

Models in Ecology

Mathematical Biology

Nonlinearities in mathematical ecology: phenomena and models, would we live in Volterra's world

Mathematical and dynamic analysis of a prey-predator model in the presence of alternative prey with impulsive state feedback control

A predator-prey mathematical model with both the populations affected by diseases

Theory and Applications of Fractional Differential Equations

Fractional Differential Equations

A new definition of fractional derivative without singular kernel

New fractional derivatives with non-local and non-singular kernel: theory and application to heat transfer model

Non validity of index law in fractional calculus: a fractional differential operator with Markovian and non-Markovian properties

Fractional derivatives with no-index law property: application to chaos and statistics

Decolonisation of fractional calculus rules: breaking commutativity and associativity to capture more natural phenomena

Dynamical study of fractional order mutualism parasitism food web module

Stability and numerical simulation of a fractional order plant-nectar-pollinator model

Green function's properties and existence theorems for nonlinear singular-delay-fractional differential equations

A fractional order HIV-TB coinfection model with nonsingular Mittag-Leffler law

A predator-prey mathematical model with both the populations affected by disease

The dynamical analysis of a prey-predator model with a refuge-stage structure prey population

A fractional calculus approach to Rosenzweig-MacArthur predator-prey model and its solution

An impulsively controlled three-species prey-predator model with stage structure and birth pulse for predator

A hybrid predator-prey model with general functional responses under seasonal succession alternating between Gompertz and logistic growth

Fractional Derivatives with Mittag-Leffler Kernel

A study of behaviour for immune and tumor cells in immunogenetic tumour model with non-singular fractional derivative

On a class of ordinary differential equations in the frame of Atangana-Baleanu fractional derivative

Solution for fractional generalized Zakharov equations with Mittag-Leffler function

Analysis of lakes pollution model with Mittag-Leffler kernel

Modelling the spread of Ebola virus with Atangana-Baleanu fractional operators

On a nonlinear fractional order model of dengue fever disease under Caputo-Fabrizio derivative

Semi-analytical study of pine wilt disease model with convex rate under Caputo-Fabrizio fractional order derivative

Study on Krasnoselskii's fixed point theorem for Caputo-Fabrizio fractional differential equations

Hyers-Ulam stability and existence criteria for coupled fractional differential equations involving p-Laplacian operator

On the existence of positive solutions for a non-autonomous fractional differential equation with integral boundary conditions

Analysis of some generalized ABC-fractional logistic models

Fractional logistic models in the frame of fractional operators generated by conformable derivatives

A fractional-order epidemic model with time-delay and nonlinear incidence rate

Discrete fractional differences with nonsingular discrete Mittag-Leffler kernels

Basic Theory of Fractional Differential Equations

Problems in Modern Mathematics

A Collection of Mathematical Problems. Interscience

Investigation of Ulam stability results of a coupled system of nonlinear implicit fractional differential equations

On Ulam's stability for a coupled systems of nonlinear implicit fractional differential equations

Ulam stability to a toppled systems of nonlinear implicit fractional order boundary value problem

Existence and Ulam-Hyers stability for Caputo conformable differential equations with four-point integral conditions

New numerical approximation of fractional derivative with non-local and non-singular kernel: application to chaotic models

The authors are very thankful to the reviewers for their useful suggestions. The fourth author would like to thank Prince Sultan University for funding this research work.

Data sharing is not applicable to this article as no data sets were generated or analyzed during the current study.

The authors declare that they have no competing interests.

All authors equally contributed to this manuscript and approved the final version. 

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.Received: 15 April 2020 Accepted: 19 May 2020