Solving Fractional-Order Logistic Equation Using a New Iterative Method

A fractional version of logistic equation is solved using new iterative method proposed by Daftardar-Gejji and Jafari (2006). Convergence of the series solutions obtained is discussed. The solutions obtained are compared with Adomian decomposition method and homotopy perturbation method.


Introduction
The following model describing growth of population was first studied by Pierre Verhulst in 1938 1 where N t is population at time t, and r > 0 is Malthusian parameter describing growth rate and K is carrying capacity.Defining x N/K gives the following differential equation: which is called as logistic equation.Logistic equation of fractional order has been discussed in the literature 2, 3 .El-Sayed et al. 2 have investigated the equation D α x t rx t 1 − x t , where D α is Caputo International Journal of Differential Equations fractional derivative of order 0 < α ≤ 1. Momani and Qaralleh 3 have employed Adomian decomposition method ADM for solving fractional population growth model in a closed system.
In the present paper we use New Iterative Method NIM introduced by Daftardar-Gejji and Jafari 4 to solve fractional version of logistic equation.
NIM is useful for solving a general functional equation of the form where f is a given function, L and N linear and nonlinear operators, respectively.The NIM has fairly simple algorithm and does not require any knowledge of involved concepts such as Adomian polynomials, homotopy, or Lagrange multipliers.Rigorous convergence analysis of NIM has been worked out recently 5 .This method has been applied by present authors successfully for solving partial differential equations 6 , evolution equations 7 , and fractional diffusion-wave equations 8 .NIM has been further explored by many researchers.Several numerical methods with higher order convergence can be generated using NIM.M. A. Noor and K. I. Noor 9, 10 have developed a three-step predictor-corrector method for solving nonlinear equation f x 0. Further, they have shown that this method has fourth-order convergence 11 .Some new methods 12, 13 are proposed by these authors using NIM.Mohyud-Din et al. 14 solved Hirota-Satsuma coupled KdV system using NIM.These authors 15 also have applied NIM in solutions of some fifth order boundary value problems.Noor and Mohyud-Din 16 have used NIM to solve Helmholtz equations.NIM is applied to solve homogeneous and inhomogeneous advection problems 17 , diffusion equations 18 , Schrödinger equations 19 , time fractional partial differential equations 20 , and so on.Yaseen and Samraiz 21 proposed modified NIM and used it to solve Klein-Gordon equations.Srivastava and Rai 22 have proposed a new mathematical model for oxygen delivery through a capillary to tissues in terms of multiterm fractional diffusion equation.They have solved the multi-term fractional diffusion equation using NIM and ADM and have shown that the results are in perfect agreement.
Recently Usman et al. 23 have solved a wide range of physical problems using NIM and various other methods and have shown that NIM has better performance as compared to other methods.
In the present paper, we solve fractional-order logistic equation.We compare the results obtained by NIM, ADM 3 , and homotopy perturbation method HPM with exact solution.Further, we propose sufficient condition for the convergence of NIM solution of fractional order logistic equation.

Preliminaries and Notations
In this section, we set up notation and recall some basic definitions from fractional calculus 24 .

Fractional-Order Logistic Equation
Consider Operating I α on both sides of 3.1 and using 3.2 , we get x t x 0 rI α x t − rI α x 2 t .3.3

Adomian Decomposition Method
Adomian decomposition method 25 is one of the most powerful methods used to solve 1.3 .In this method, solution u is assumed to be of the form ∞ i 0 u i , and the nonlinear operator N is represented as where , n 0, 1, 2, . . .

3.5
are the Adomian polynomials.The terms of the series are given as

International Journal of Differential Equations
For convergence of this series, we refer the reader to 26 .Fractional order logistic equation is solved using ADM in 3 .

New Iterative Method
Recently Daftardar-Gejji and Jafari 4, 8 have used a different decomposition of the term N u .They have decomposed N u as follows: The terms of the series u ∞ i 0 u i are determined by the following recurrence relation:

3.8
Discussion regarding convergence of this series and comparison of ADM and NIM can be found in 5 .
In view of NIM, the solution of 3.3 is given by x t ∞ i 0 x i where , n 1, 2, . . . .

3.9
The k-term approximate solution is given by 3.10

Homotopy Perturbation Method
Homotopy perturbation method HPM proposed by He 27, 28 is a useful technique to solve nonlinear problems.In view of HPM, we construct the following homotopy for solving fractional order logistic equation 3.1 : where p ∈ 0, 1 is an embedding parameter.For p 0, the homotopy 3.11 gives a linear equation D α x 0 and for p 1, it gives the original equation 3.1 .The solution of homotopy equation 3.11 is assumed to be The solution of 3.1 is obtained from series 3.12 by substituting p 1. Substituting 3.12 in 3.11 and equating like powers of p, we get

3.13
and so on.The linear equations given in 3.13 give x 4 0, 3.14 and so forth.It can be seen from 3.11 that the HPM solutions are the same as ADM solutions given in 3 .

Convergence of NIM Solutions
In this section, we provide sufficient condition for the convergence of NIM solution series 3.10 .Let f t, x rx 1 − x be defined on a rectangle It is clear that |f| is bounded on R. We have |f t, x | ≤ M, for some real number M on R.

International Journal of Differential Equations
Lemma 4.1.For the initial value problem the m-term approximate NIM solution 3.10 exists on the interval I −χ, χ , where χ min{a, Γ α 1 b/M 1/α }, on rectangle R.Moreover, the points t, u m t ∈ R, for all m 1, 2, . . .whenever t ∈ I.
Proof.We prove the lemma by induction on m.
Clearly, the result is true for m 1.Now, Since the points t, x 0 t ∈ R, for all t ∈ I,
We assume that the result is true for positive integer m − 1 and proves for m.Using NIM solutions, we get

4.5
Induction hypothesis implies that the points t, u m−1 t ∈ R, for all t ∈ I. Hence This implies that the points t, u m t ∈ R, for all m 1, 2, . . .when t ∈ I.
Theorem 4.2.The NIM solution series converges on the interval I −χ, χ to a solution of the IVP 4.2 .
Proof.Note that f t, x rx 1 − x is continuous function defined on a rectangle R and |f| ≤ M on R. Further

4.8
It can be proved by induction that m 1, 2, . ... This shows that the mth term of the NIM solution series is bounded by M/K times the mth term of the Mittag-Leffler function E α K|t| α .

4.10
Thus, the NIM series is convergent under given conditions.

Illustrative Examples
Mathematica 8 has been used for solving the following problems.

International Journal of Differential Equations
Equivalent integral equation of 5.1 is In view of NIM, x 0 0.5, and so on.Few terms of ADM 3 and HPM solution are

5.4
Figure 1 shows the solutions of 5.2 for α 1 by NIM 5-term solution , ADM, and HPM 6-term solution and exact solution x t exp 1/2 t / 1 exp 1/2 t .In Figures 2 and  3 we compare ADM, HPM, and NIM solutions of 5.2 for α 0.5 and α 0.7, respectively.Figure 4 shows NIM solutions for different values of α.
Example 5.2.Consider the following logistic equation of fractional order x 0 1 3 .5.6 Applying I α on both sides of 5.5 and using 5.6 , we get In Figure 5, we compare the solution of 5.7 for α 1 by NIM 5-term solution , ADM 3 , and HPM 6-term solution and exact solution exp 1/4 t / 2 exp 1/4 t .Figure 6 shows NIM solutions of 5.7 for different values of α.

Conclusions
Fractional-order logistic equation is solved using a new iterative method NIM .From Figures 1-6, it is clear that NIM solutions are more stable than ADM and HPM.Further the condition for the convergence of NIM solution series is also provided.It can be concluded that NIM is a useful technique for solving nonlinear problems.