On the Approximate Solution of Fractional Logistic Differential Equation Using Operational Matrices of Bernstein Polynomials

In this paper, operational matrices of Bernstein polynomials (BPs) are presented for solving the non-linear fractional Logistic differential equation (FLDE). The fractional derivative is described in the Riemann-Liouville sense. The operational matrices for the fractional integration in the Riemann-Liouville sense and the product are used to reduce FLDE to the solution of non-linear system of algebraic equations using Newton iteration method. Numerical results are introduced to satisfy the accuracy and the applicability of the proposed method.


Introduction
It is well-known that the fractional differential equations (FDEs) have been the focus of many studies due to their frequent appearance in various applications, such as in fluid mechanics, viscoelasticity, biology, physics and engineering applications, for more details see for example ( [1] [2]).Consequently, considerable attention has been given to the efficient numerical solutions of FDEs of physical interest, because it is difficult to find exact solutions.Different numerical methods have been proposed in the literature for solving FDEs ([3]- [6]).Recently, several numerical and approximate methods to solve FDEs have been given, such as variational iteration method [7], homotopy perturbation method [7] and collocation method ([8]- [13]).
The fractional Logistic model can obtain by applying the fractional derivative operator on the Logistic equation.The model is initially published by Pierre Verhulst in 1838 ( [14] [15]).The continuous Logistic model is described by first order ordinary differential equation.The discrete Logistic model is simple iterative equation that reveals the chaotic property in certain regions [16].There are many variations of the population modeling [17].The Verhulst model is the classic example to illustrate the periodic doubling and chaotic behavior in dynamical system [16].The model is described the population growth may be limited by certain factors like population density ( [15] [17]).
The solution of Logistic equation explains the constant population growth rate which does not include the limitation on food supply or spread of diseases [15].The solution curve of the model increases exponentially from the multiplication factor up to saturation limit which is maximum carrying capacity [15], d 1 d where N is the population with respect to time, ρ is the rate of maximum population growth and K is the car- rying capacity.The solution of continuous Logistic equation is in the form of constant growth rate as in formula where 0 N is the initial population [18].In this article, we consider FLDE of the form the parameter 0 1 α < ≤ refers to the fractional order derivative We also assume an initial condition The exact solution to this problem at 1 The existence and the uniqueness of the proposed problem (1) are introduced in details in ( [19] [20]).Khader and Hendy [21] introduced a new approximate formula of the fractional derivative using Legendre series expansion and used it to solve numerically the fractional delay equation.In this article, we extended this work to study the numerical solution of the non-linear FLDE.An approximate formula of the fractional derivative is presented.Special attention is given to study the convergence analysis and estimate an upper bound of the error of the introduced formula.

Preliminaries and Notations
In this section, we present some necessary definitions and mathematical preliminaries of the fractional calculus theory and the Bernstein polynomials that will be required in the present paper.

The Fractional Integral and Derivative Operators
We present some necessary definitions and mathematical preliminaries of the fractional calculus theory that will be required in the present paper.
The Riemann-Liouville fractional integral operator a J α of order α is defined on [ ] The Riemann-Liouville fractional derivative operator a D α of order α ( )

Definition 3.
The Caputo fractional derivative operator c a D α of order α ( ) Similar to integer-order differentiation, Caputo fractional derivative operator is a linear operation where λ and µ are constants.For the Caputo's derivative we have [2] 0, is a constant, ( ) ( ) 0 0 0 for and ; 1 , for and . 1 We use the ceiling function α     to denote the smallest integer greater than or equal to α , and Recall that for α ∈  , the Caputo differential operator coincides with the usual differential op- erator of integer order.
For more details on fractional derivatives definitions and its properties see The ( )

Bernstein Polynomials and Their Properties
Bernstein polynomials of degree n are defined on the interval [ ] 0,1 as [22] ( ) ( ) , 1 , 0,1, , , where is a binomial coefficient.The first few Bernstein basis polynomials are: The Bernstein polynomials have the following properties 1) where δ is the Kronecker delta function; 3) ; We can write , where A is an upper triangular matrix, ( ) . For more details about the definition, properties and the convergence analysis of Bernstein polynomials [23].

BPs Operational Matrix of Riemann-Liouville Fractional Integration
Theorem 1. [23] The Bernstein polynomials operational matrix F α from order ( ) ( ) for the Riemann-Liouville fractional integral is defined as follows

Definition 5.
We can define the dual matrix ( ) ( ) Q + × + on the basis of Bernstein polynomials of mth degree as follows ( ) ( ) ( ) where Lemma 2. [24] Let [ ] 2 0,1 L be a Hilbert space with the inner product ( ) ( ) we can find the unique vector is the best approximation of ( )

Let ( )
u x be a continuous function on the interval [ ] 0,1 .Then we can approximate it in the following polynomial in Bernstein form of degree n

0,1 u x C
∈ and any 0 δ > , there exists an integer N such that The Bernstein polynomials operational matrix are used for solving many class of fractional differential equations, they used to solve numerically the fractional heat-and wave-like equations [25] and the multi-term orders fractional differential equations [26] and others [27].

Implementation of Bernstein Polynomials Operational Matrix for Solving FLDE
In this section, we introduce a numerical algorithm using Bernstein polynomials operational matrix method for solving the fractional Logistic differential equation of the form (1).
) By solving this system we can obtain the vector C.Then, we can get The numerical results of the proposed problem (1) are given in Figure 1 and From these figures we can conclude that the obtained numerical solutions are in excellent agreement with the exact solution.

Conclusion and Remarks
In this article, we used operational matrices of the Riemann-Liouville fractional integral and the product by Bernstein polynomials for solving the fractional Logistic differential equation.The properties of these operational matrices are used to reduce FLDE to a non-linear system of algebraic equations which solved by Newton iteration method.From the obtained numerical results, we can conclude that this method gives results with an excellent agreement with the exact solution.All numerical results are obtained using Matlab program 8.
basis, so, we can write any polynomial ( ) u x of degree m in terms of linear combination of

Figure 2 Figure 1 ,
with different values of α in the interval [ ] we presented a comparison between the behavior of the exact solution and the approximate solution using the introduced technique at 1 α = (left), and the behavior of the approximate solution using the proposed method at 0.9 α = (right).But, in Figure2we presented the behavior of the approximate solution with different values of α (

Figure 1 .Figure 2 .
Figure 1.A comparison between the approximate solution and the exact solution at 1α = (left).The behavior of the ap- proximate solution using the proposed method at 0.9 α = (right).