Analytical Solution of Nonlinear System of Fractional Differential Equations ()
1. Introduction
This paper is concerned with the analytical solution of a nonlinear system of fractional differential equations. Systems of fractional differential equations (FDEs) have many applications in engineering and science, including electrical networks, fluid flow, control theory, fractals theory, electromagnetic theory, viscoelasticity, potential theory, chemistry, biology, optical and neural network systems ( [1] - [16] ). We use Adomain decomposition method ( [17] - [24] ) for solving this type of equations. The existence and uniqueness of the solution are proved, the convergence of ADM series solution is discussed and the error analysis is given. This method has many advantages; it is efficiently working with different types of linear and nonlinear equations in deterministic or stochastic fields and gives an analytic solution for all these types of equations without linearization or discretization.
2. Formulation of the Problem
Consider the following nonlinear system of FDEs:
(1)
subject to the initial conditions,
(2)
where,
where
is bounded
,
,
,
are finite constants and
satisfy Lipschitz condition with Lipschitz constants
such as,
(3)
and the fractional derivative in this system is of sequential Caputo sense which defined as
In the applications, the Caputo sense is preferred to be used because the initial conditions of
and its derivatives will be of integer orders and have a physical meaning.
Now performing subsequently the fractional integration of order
, this reduces the system (1)-(2) to the system of FIEs,
(4)
and has Adomian polynomials representation,
(5)
where,
(6)
Substitute from Equation (5) into Equation (3), we get
(7)
Let
in (7) we get,
(8)
(9)
Finally, the solution is,
(10)
3. Analysis of Convergence
3.1. The Uniqueness of Solution
In the previous section, we find the series solution (10) of the system (1)-(2) and here we want to prove the existence and uniqueness of this series solution.
Theorem 1. If
where
, then the series (10) is the solution of the system (1)-(2) and this solution is unique, where
,
.
Proof. For existence,
then the Adomian’s series solution satisfy Equation (4) which is the reduced system of FIEs to the system (1)-(2).
For uniqueness of the solution: Assume that
and
are two different solutions to the system (1)-(2) and hence,
Let
where,
then,
but,
and since,
then,
this implies that,
and this completes the proof.
3.2. Proof of Convergence
Theorem 2. The series solution (10) of the system (1)-(2) using ADM converges if
and
,
, where
,
.
Proof. Define the Banach space
, the space of all continuous functions on J with the norm
and a sequence
such that,
. We have,
Let,
and
be two arbitrary partial sums with
. Now, we are going to prove that
is a Cauchy sequence in this Banach space.
Let
then,
Using the triangle inequality,
Since,
, and
then,
. Consequently,
but,
and as
then,
and hence,
is a Cauchy sequence in this Banach space so, the series
converges and the proof is complete.
3.3. Error Analysis
For ADM, we can estimate the maximum absolute truncated error of the Adomian’s series solution in the following theorem.
Theorem 3. The maximum absolute truncation error of the series solution (10) to the system (1)-(2) is estimated to be,
Proof. From Theorem 2 we have,
But,
as then,
so,
So, the maximum absolute truncation error in the interval J is,
and this completes the proof.
4. Numerical Examples
Example 1. Consider the following nonlinear system of FDEs,
(11)
subject to the initial conditions,
where
.
This system was discussed before in [25], it is solved by using the iterative method. Now, we will solve it by using ADM. Applying ADM to system (11) leads to the following recursive relations,
(12)
(13)
(14)
where
and
represent the Adomian polynomials of the nonlinear terms
and
respectively.
Using the relations (12)-(14), the first three terms of the series solution when
are,
(15)
while for
are,
(16)
and for
are,
(17)
Figures 1(a)-(c) show ADM solution of
and
at different values of
(
).
Example 2. Consider the following nonlinear system of FDEs,
(18)
Figure 1. (a) ADM Sol. [n = 5]; (b) ADM Sol. [n = 5]; (c) ADM Sol [n = 5].
subject to the initial conditions,
which has the exact solution
,
and
.
Using ADM to system (18) leads to the following scheme,
(19)
(20)
(21)
where
and
represent the Adomian polynomials of the nonlinear terms
and
respectively.
Using relations (19)-(21), the first few terms of the series solution are,
(22)
A comparison between ADM solution and exact solution of
and
is given in Figures 2(a)-(c) (
).
Figure 2. (a) ADM and Exact Sol; (b) ADM and Exact Sol; (c) ADM and Exact Sol.
Example 3. Consider the following nonlinear system of FDEs [26],
(23)
subject to the initial conditions,
where
.
Using ADM to system (23) leads to the following scheme,
(24)
(25)
where
and
represent the Adomian polynomials of the nonlinear terms
and
respectively.
Using the relations (24)-(25), the first four terms of the series solution are,
(26)
(27)
Figure 3(a) and Figure 3(b) show ADM solution of
and
at different
Figure 3. (a) ADM Sol. [n = 5]; (b) ADM Sol. [n = 5].
values of
(
).
Example 4. Consider the nonlinear system of FDEs,
(28)
subject to the initial conditions,
Using ADM to the system (28), we get
(29)
From the relations (29), the first two terms of the series solution are,
(30)
Figure 4(a) and Figure 4(b) show ADM solution of
and
(
).
Now, we will use Theorem 3 to evaluate the maximum absolute truncated error of the series solution (30). So, we evaluate the following values,
·
:
.
·
:
.
·
:
.
·
:
.
·
:
.
·
,
.
The maximum error of
:
,
· For
:
,
· For
:
,
· For
:
,
· For
:
.
The maximum error of
:
,
· For
:
,
· For
:
,
· For
:
,
· For
:
.
5. Application: On Fractional-Order Rabies Model
The fractional-order rabies model,
(31)
subject to the initial conditions,
was discussed before in [27], it was solved by using Adams-type predictor-corrector method. Now, we will solve it by using ADM.
Applying ADM to the system (31) leads to the following scheme,
(32)
(33)
where
represents the Adomian polynomials of the nonlinear term
.
Using the relations (32)-(33) and taking
, the first five-terms of the series solution are,
Figure 5. (a) ADM Sol. [α = 0.9 ]; (b) ADM Sol. [α = 0.9]; (c) ADM Sol. [α = 0.8 ]; (d) ADM Sol. [α = 0.8 ].
(34)
while for
, the series solution will be,
(35)
Figures 5(a)-(d) show ADM solution of
and
(
) at different values of
.
6. Conclusion
In this paper, we use a simple method to solve nonlinear system of FDEs, this method gives a good approximate analytical solution of this type of equation as we compare ADM solution with the exact solution and also by evaluating the maximum absolute error which results from using partial sum of the series ADM solution.