Abstract
The main objective of this paper is to examine the stability and convergence of the Laplace-Adomian algorithm to approximate solutions of the pantograph-type differential equations with multiple delays. This is done by comparatively investigating it with other methods.
1. Introduction
Delay differential equations (DDEs) are a large and important class of differential equations in which the derivative of the unknown function at a certain time is given in terms of the values of the function at previous times. They often arise in wide and diverse range of applications in population studies [1, 2], economics [3], medical biology [4, 5], controls of mechanical systems [6], and so forth. The pantograph equation is one of the most important kinds of DDEs. The name pantograph was used by Iserles and Liu [7] to study how the electric current is collected by the pantograph of an electric locomotive, from where it gets its name.
In recent years, some promising approximate analytical solutions have been proposed, such as the Taylor collocation method [8], Bernstein polynomials [9], spline method [10, 11], and other methods are reviewed in [2, 12, 13].
In this work we consider the following problems
Problem 1.
Problem 2.
where is an analytical function, , and ; , .
The aim of this paper is to employ and examine the stability of the Laplace-Adomian algorithm (LAA) for solving Problems 1 and 2. This method was first proposed by Khuri [14], who applied the scheme to a class of nonlinear differential equations. In this method the solution is given as an infinite series usually converging very rapidly to the exact solution of the problem.
A major advantage of this method is that it is free from round-off errors and without any discretization or restrictive assumptions. Therefore, results obtained by LAA are more accurate and efficient. LAA has been shown to accurately and easily approximate solutions of large class of linear and nonlinear ODEs and PDEs [14–16]; for example, Ongun [17] employed LAA to give an approximate solution of nonlinear ordinary differential equation systems, such as a model for HIV infection of T cells, Wazwaz [18] also used this method for handling the nonlinear Volterra integro-differential equations, Khan and Faraz [19] modified LAA to obtain series solutions of the boundary layer equation, and Yusufoglu [20] adapted LAA to solve the Duffing equation.
The numerical technique of LAA basically illustrates how the Laplace transforms are used to approximate the solution of the nonlinear differential equations by manipulating the decomposition method that was first introduced by Adomian [21, 22].
2. Laplace-Adomian Algorithm
To illustrate the basic idea of the Laplace-Adomian algorithm, we consider the following nonlinear operator: with the initial condition where is a linear operator, is a nonlinear operator, and is a given analytical function.
The technique consists first of applying the Laplace transform (denoted throughout this paper by ) to both sides of (3), to get Applying the formulas of the Laplace transform, we obtain where Suppose the answer to (3) is as follows: where the terms are to be recursively computed, the nonlinear operator is decomposed as follows: where is an infinite series of the Adomian polynomials , calculated by the formula [21]
Substituting (8) and (9) into (6) leads to
Using the linearity of the Laplace transform gives (11) as
Matching both sides of (11) yields
Generally
Applying the inverse Laplace transform to (13) gives the initial approximation
Substituting this value of into the inverse Laplace transform of (14) gives . The other terms can be obtained recursively in similar fashion from
Using (18), we can rewrite them as follows to obtain :
Subsiutiting the values of into (8) gives
The sum of terms in (8), that is, , where , tends to . This means that is an appropriate approximation of . The terms in the above series soon tend to zero where has been the coefficients of calculations derived from the operation , is the number of terms, and is the order of the operation derivation. Therefore, it has a rapid convergence.
3. Test Problems
In this section, we will apply Laplace-Adomian algorithm to solve the pantograph delay equation. The absolute errors in Tables 1–4 are the values of at selected points. All iterates are calculated by using Matlab 7.
3.1. Stability of the Laplace-Adomian Algorithm
A method is said to be stable when the obtained solution undergoes small variations as there are slight variations in inputs and parameters and when probable perturbations in parameters that are effective in equations and conditions prevailing them do not introduce, in comparison to the physical reality of the problem, any perturbations in what is returned. We propose here to compare the Laplace-Adomian algorithm with other numerical methods by offering examples and examining the stability of the Laplace-Adomian algorithm.
Example 1 (Evans and Raslan [13]). Consider the following pantograph equation:
The exact solution is . Applying the result of (17) gives us
The iteration formula (18) for this example is
In Table 1 we make a comparison between the absolute errors obtained by the Bernstien series [9] in column 2, together with the spline method [10, 11] in columns 3 and 4, and finally, the Laplace-Adomian algorithm with , and in the last three columns. Given the following absolute errors that were obtained from these methods, we can conclude that the rate of convergence of LAA is higher than the other methods. However, taking a closer look at the errors of LAA shows that the error increases significantly more than the other methods in the bottom of the columns. This means that the stability of LAA is decreasing more than the other methods.
Example 2 (Muroya et al. [23]). Consider the following pantograph equation:
The exact solution is . Applying the result of (17), gives us
The other terms of the sequences can be obtained directly from (18) as in general
Table 2 compares the results of the Laplace-Adomian algorithm and the Taylor collocation method [8]. It can be seen that the Laplace-Adomian algorithm is weaker in the stability than Taylor collocation method and yet is stronger in the convergence.
Example 3 (Liu and Li [12]). Consider the multipantograph delay equation with variable coefficients Here , and . The initial approximation for this example is And from (18), we obtain Comparisons of approximate solutions for few terms with exact solution are illustrated in Table 3.
Example 4. Consider the pantograph equation of third order
Let us start with an initial approximation
The other terms of the sequences can be obtained directly from (18)
where is the Adomian polynomials that represent the nonlinear term . The Adomian polynomials for are given by
Table 4 shows the absolute error of the numerical approximation by using a few iterations.
The results for Examples 3–4 in Figures 1 and 2, respectively, show that the approximations converge rapidly and only a few iterations are sufficient to obtain accurate solutions.
Example 5. Consider the multipantograph equation
Following the procedures in the previous examples gives us
According to (18), we obtain
Thus
which is the exact solution.
4. Conclusion
In this paper, the Laplace-Adomian algorithm has been successfully applied to the generalized pantograph equations with multiple delays to obtain high approximate solutions or exact solutions with little iterations used. Moreover, the Laplace-Adomian algorithm is simple and easy to use and stronger in convergence compared with other methods. Despite these advantages, it is seen that its stability is lower than other numerical methods.