Laplace Discrete Adomian Decomposition Method for Solving Nonlinear Integro Differential Equations ()
1. Introduction
One of the issues in finding the solutions of integro-differential equations is the choices of approaches used that are either analytical or numerical. Previous literature suggested that analytical solutions of integro-differential equations are not usually obtainable especially when the equations are entangled with nonlinear terms. Therefore, it is required to find an efficient approximate solution numerically. Much attention has been devoted for the search of better and more efficient methods in recent years with the introduction of several numerical methods such as the Galerkin methods [1] , Runge-Kutta methods [2] , Chebyshev collocation method [3] , Taylor collection method [4] , rationalized Haar functions method [5] , Galerkin methods with hybrid functions [6] , and Adomian Decomposition Method (ADM) [7] - [12] . In addition to these numerical methods, Khuri [13] [14] used Laplace transform numerical scheme, based on the ADM to solve nonlinear differential equations and Bratu’s problem, respectively. This method is popularly known as Laplace Adomian Decomposition Method (LADM) where the ADM and Laplace transforms are combined.
The LADM is known for its rapid convergence in solution and also uses only little iteration as successfully applied in Kiymaz [15] . Furthermore, several modifications of the ADM and LADM methods can be seen in [16] - [22] with wide applications ranging from differential equations, partial differential equations, integral equations and integro-differential equations among others. It seems that the LADM method is always open for further modifications especially on discretizing the Adomian decomposition. In this paper, we aim at extending the Laplace Adomian decomposition method for finding the solution of nonlinear integro-differential equations by firstly discretizing the Adomian decomposition method, followed by coupling some numerical integration schemes or quadrature rules. These quadrature rules are used to approximate the definite integrals which are analytically impossible [23] [24] [25] [26] as the solution is given at the nodes. For convenience, the proposed extension is known as Laplace Discrete Adomian Decomposition Method (LDADM). The rest of this paper is arranged as follows. Section 2 describes the proposed LDADM. Section 3 presents numerical results with four examples and finally Section 4 concludes.
2. Laplace Discrete Adomian Decomposition Method
In this section, we present the Laplace discrete Adomian decomposition method over the integro-differential equation of the form
(2.1)
Subject to the following initial conditions
To solve the nonlinear integro-differential Equation (2.1) using the Laplace transform method, we recall the Laplace transform of the second derivative of
, that is
. Thus on applying the Laplace transform to both sides of Equation (2.1) we obtain
(2.2)
We thus obtain the following equation after using the above prescribed initial conditions
(2.3)
The decomposition method represents the solution
as a series of the form:
(2.4)
and the nonlinear term
is decomposed into an infinite series of the form:
(2.5)
where
’s are the Adomian polynomials of
given by the formula
(2.6)
On substituting Equation’s (2.4) and (2.5) in Equation (2.3) and making comparison between the right and left hand sides, we thus obtain:
(2.7)
Furthermore, if the evaluation of the integrals in Equation (2.7) is analytically possible, the ADM can be applied directly. However, in the case where the evaluation of integrals is analytically impossible we consider the numerical integration scheme given by the following formula:
(2.8)
where
is continuous function on
,
are the nodes of the numerical integration,
is the fixed step length and
are the weights functions. Now, applying the formula given in Equation (2.8) on Equation (2.7) to obtain
(2.9)
Finally, on applying the inverse Laplace transform to the first part of Equation (2.9) gives
, and, consequently will define
. Also, using
enables us to evaluate
. The determination of
and
leads to the determination of
that will allow us to determine
, and so on. This in turn will lead to the complete determination of the components of
upon using the second part of Equation (2.9).
(2.10)
From this recursive relation in Equation (2.10),
can be calculated. The solution of Equation (2.1) is now determined in Equation (2.10). However, in practice the series
must be truncated to the series:
with
(2.11)
2.1. Trapezoidal Method (TR)
We couple the trapezoidal method (TR) to Equation (2.10) to obtain:
(2.12)
2.2. Simpson’s Method (SM)
Also, on coupling the Simpson’s method (SM) to Equation (2.10), we get:
(2.13)
3. Computational Results and Analysis
In this section, we consider several nonlinear integro-differential equations as examples in order to show the efficiency and the simplicity of the proposed method. We start with the nonlinear Volterra and Fredholm integro-differential equations down to their systems, respectively.
Example 3.1
Consider the nonlinear Volterra integro-differential equation
(3.1)
with the exact solution given by
.
In order to use the quadrature rule for Equation (3.1), let
, we get
Taking the Laplace transform of both sides of the above equation gives
So that,
or equivalently
1) Trapezoidal Method (TM)
We divide the interval (0, 1) into subinterval of equal lengths
and denote
.
The recursive relation is given by
or
Taking the inverse Laplace transform of both sides of the first part of recursive relation, and using this recursive relation gives
The series solution is obtained by summing the following iterates
.
The results produced by the present method with only few components (m = 5) are in a very good agreement with the exact solution results as shown in Table 1 and illustrated graphically in Figure 1.
2) Simpson’s Method (SM)
We divide the interval (0, 1) into subinterval of equal lengths
and denote
. The recursive relation is given by
Table 1. Comparison between the exact solution
and approximate solution using LDADM based on TM in EX. 3.1.
Figure 1. Curves of the exact solution
and approximate solution using LDADM based on TM in EX. 3.1.
Or
Taking the inverse Laplace transform of both sides of the first part of recursive relation, and using this recursive relation gives
The series solution is obtained by summing the following iterates
.
The results produced by the present method with only few components (m = 5) are in a very good agreement with the exact solution results as shown in Table 2 and illustrated graphically in Figure 2.
Example 3.2
Consider the nonlinear Fredholm integro-differential equation
Table 2. Comparison between the exact solution
and approximate solution using LDADM based on SM in EX. 3.1.
Figure 2. Curves of the exact solution
and approximate solution using LDADM based on SM in EX. 3.1.
with the exact solution is
.
Taking the Laplace transform of both sides of the given equation gives
So that
or equivalently
1) Trapezoidal Method (TM)
Let
, the recursive relation is given by
Taking the inverse Laplace transform of both sides of the first part of recursive relation, and using this recursive relation gives
The series solution is obtained by summing the following iterates
.
The results produced by the present method with only few components (m = 5) are in a very good agreement with the exact solution results as shown in Table 3 and illustrated graphically in Figure 3.
2) Simpson’s Method (SM)
Let
, the recursive relation is given by
Table 3. Comparison between the exact solution
and Approximate solution using LDADM based on RM in EX. 3.2.
Figure 3. Curves of the exact solution
and approximate solution using LDADM based on TM in EX. 3.2.
Taking the inverse Laplace transform of both sides of the first part of recursive relation, and using this recursive relation gives
The series solution is obtained by summing the following iterates
.
The results produced by the present method with only few components (m = 5) are in a very good agreement with the exact solution results as shown in Table 4 and illustrated graphically in Figure 4.
Example 3.3
Consider the system of nonlinear Volterra integro differential equation
(3.2)
with the exact solution
.
In order to use the quadrature rule for Equation (3.2), let
, we get
Table 4. Comparison between the exact solution
and approximate solution using LDADM based on SM in EX. 3.2.
Figure 4. Curves of the exact solution
and approximate solution using LDADM based on SM in EX. 3.2.
Taking the Laplace transform of both sides of the above equation gives
So that
1) Trapezoidal Method (TM)
We divide the interval (0, 1) into subinterval of equal lengths
and denote
. The recursive relations as
Taking the inverse Laplace transform of both sides of the first part of recursive relation, and using this recursive relation gives
The series solutions are
The results produced by the method with only few components (m = 5) are in a very good agreement with the exact solution results as shown in Table 5 and Table 6 and illustrated graphically given in Figure 5 and Figure 6.
2) Simpson’s Method (SM)
We divide the interval (0, 1) into subinterval with equal length
and denote
. The recursive relation is given by
Taking the inverse Laplace transform of both sides of the first part of recursive relation, and using this recursive relation gives
Table 5. Comparison between the exact solution
and approximate solution using LDADM based on TM in EX. 3.3.
Table 6. Comparison between the exact solution
and approximate solution using LDADM based on TM.
Figure 5. Curves of the exact solution
and approximate solution using LDADM based on TM in EX. 3.3.
Figure 6. Curves of the exact solution
and approximate solution using LDADM based on TM in EX. 3.3.
The series solutions are
The results produced by the method with only few components (m = 5) are in a very good agreement with the exact solution results as shown in Table 7 and Table 8 and illustrated graphically given in Figure 7 and Figure 8.
Table 7. Comparison between the exact solution
and approximate solution using LDADM based on SM in EX. 3.3.
Table 8. Comparison between the exact solution
and approximate solution using LDADM based on SM in EX. 3.3.
Figure 7. Curves of the exact solution
and approximate solution using LDADM based on SM in EX. 3.3.
Figure 8. Curves of the exact solution
and approximate solution using LDADM based on SM in EX. 3.3.
Example 3.4
Consider the system of nonlinear Fredholm integro differential equation
with the exact solution
.
As usual, on taking the Laplace transform of both sides of the above equation gives
So that
1) Trapezoidal Method (TM)
We divide the interval (0, 1) into subinterval of equal lengths
and denote
. The recursive relations are
Taking the inverse Laplace transform of both sides of the first part of recursive relation, and using this recursive relation gives
The series solutions are
The results produced by the method with only few components (m = 5) are in a very good agreement with the exact solution results as shown in Table 9 and Table 10 and illustrated graphically given in Figure 9 and Figure 10.
2) Simpson’s Method (SM)
We divide the interval (0, 1) into subinterval of equal lengths
and denote
. The recursive relation is given by
Taking the inverse Laplace transform of both sides of the first part of recursive relation, and using this recursive relation gives
Table 9. Comparison between the exact solution
and approximate solution using LDADM based on TM in EX. 3.4.
Table 10. Comparison between the exact solution
and approximate solution using LDADM based on TM in EX. 3.4.
Figure 9. Curves of the exact solution
and approximate solution using LDADM based on TM in EX. 3.4.
Figure 10. Curves of the exact solution
and approximate solution using LDADM based on TM in EX. 3.4.
The series solutions are
The results produced by the method with only few components (m = 5) are in a very good agreement with the exact solution results as shown in Table 11 and Table 12 and illustrated graphically given in Figure 11 and Figure 12.
Table 11. Comparison between the exact solution
and approximate solution using LDADM based on SM in EX. 3.4.
Table 12. Comparison between the exact solution
and approximate solution using LDADM based on SM in EX. 3.4.
Figure 11. Curves of the exact solution
and approximate solution using LDADM based on SM in EX. 3.4.
Figure 12. Curves of the exact solution
and approximate solution using LDADM based on SM in EX. 3.4.
4. Conclusion
In this paper, a modification of the Laplace Adomian decomposition method inspired by property of discretization is proposed. We developed a new Laplace Discrete Adomian decomposition method (LDADM) in which has been successfully applied to finding efficient numerical solutions of integro-differential equations featuring both nonlinear Volterra and Fredholm integrals. The method was based on the well-known Adomian decomposition method coupled with some numerical integration schemes (quadrature rules) alongside utilizing the famous and most used Laplace transform. The method gives approximate solutions iteratively with less number of computational steps. The results reveal that the proposed method is simple to execute and effective. Thus, many highly nonlinear integro-differential equations can be solved using the proposed method.