Abstract
We study higher-order boundary value problems (HOBVP) for higher-order nonlinear differential equation. We make comparison among differential transformation method (DTM), Adomian decomposition method (ADM), and exact solutions. We provide several examples in order to compare our results. We extend and prove a theorem for nonlinear differential equations by using the DTM. The numerical examples show that the DTM is a good method compared to the ADM since it is effective, uses less time in computation, easy to implement and achieve high accuracy. In addition, DTM has many advantages compared to ADM since the calculation of Adomian polynomial is tedious. From the numerical results, DTM is suitable to apply for nonlinear problems.
1. Introduction
Recently, many researchers use ADM to approximate numerical solutions. In [1], Wazwaz proposed a modification of ADM method in series solution to accelerate its rapid convergence, and, in [2], Wazwaz also presented several numerical examples of higher-order boundary value problems for first-order linear equation and second-order nonlinear equation by applying modified decomposition method.
In addition, Wazwaz [3, 4] provided first-order linear and second-order nonlinear problems to solve fifth-order and sixth-order boundary value problems by the modified decomposition method. Later, Meštrović [5] solved eight-order boundary value problems for first-order linear and second-order nonlinear boundary value problems. Similarly, in [6], Hosseini and Jafari used Adomian decomposition method to solve high-order and system of nonlinear differential equations.
The Adomian decomposition method is widely used in applied science to compute the series solution accurately because it provides rapid convergent series to solve the problem. The big application of ADM in research area is stochastic and deterministic problems involving differential, integrodifferential, integral, differential delay, and systems of such equations; for example, see [3, 7]. For example, we use the same method for solving several different problems, such as, in calculus of variations, see [8], for eikonal partial differential equation, see [9], for the Fitzhugh-Nagumo equation which models the transmission of nerve impulses, see [10], for linear and nonlinear systems of Volterra functional equations using Adomian-Pade technique, see [11], for coupled Burgers equations by using Adomian-Pade technique see [12], for solution of a nonlinear time-delay model in biology by using semi-analytical approaches, see [13], for solving the pantograph equation of order , see [14], and for nonclassic problem for one-dimensional hyperbolic equation by using the decomposition procedure, see [15].
Further, Ray and Bera in [16] used ADM to solve analytical solution of a fractional diffusion equation. They performed the explicit solution of the equation in the closed form by using initial value problem.
In [17], Ayaz investigated initial value problem of partial differential equation (PDE) to solve two-dimensional differential transformation method, and we compare the results with Adomian decomposition method. The results show that the solutions of the present method are exactly the same as the decomposition method, but the calculation of DTM is simple and reduces the difficulty of calculations.
Ayaz in [18] performed two- and three-dimensional differential transformation methods to find exact solutions of linear and nonlinear partial differential equations. Results are compared to decomposition method, and DTM has less computational effort. After that, Ertürk and Momani in [19] presented numerical solution by comparing the differential transformation method (DTM) and Adomian decomposition method (ADM) for solving linear and nonlinear fourth-order boundary value problems and proved that DTM is very accurate and efficient in numerical solution.
Recently, Arikoglu and Ozkol solved fractional differential equations by using differential transform method. They applied fractional differential equations to various types of problems such as the Bagley-Torvik, Ricatti, and composite fractional oscillation equations; see [20].
In this study, we make comparison among differential transformation method, Adomian decomposition method and exact solutions. We prove that DTM is more powerful technique than ADM and can be applied to nonlinear problems easily.
2. Differential Transformation Method
Suppose that the function is continuously differentiable in the interval for , then we have the following definition.
Definition 2.1. The differential transform of the function for the th derivative is defined as follows: where is the original function and is the transformed function. The inverse differential transform of is defined as Note that, the substitution of (2.1) into (2.2) yields the following equation: which is the Taylor's series for at .
The following theorems are operations of differential transforms.
Theorem 2.2. If , then .
Theorem 2.3. If , then .
Theorem 2.4. If , then .
Theorem 2.5. If , then .
Theorem 2.6. If , then .
Theorem 2.7. If , then .
Theorem 2.8. If , then , where
Theorem 2.9. If , then .
Theorem 2.10. If , then .
Theorem 2.11. If , then .
Theorem 2.12. If , then .
See the details in [21].
3. General Differential Transformation
Now we prove the following generalized theorem.
Theorem 3.1. The general differential transformation for nonlinear th-order BVPs, , is given by
Proof. We prove this theorem by induction method. First of all, we prove the generalization of the differential equation . Let , then
For , we have
By definition, we have
From (3.1), for , we have
This implies that the theorem holds for . Now assume that for , the DT is as follows:
Note that we have
For , we have
Thus, for , we have
For , we have
By definition, we have
From (3.1), for , we have
This implies that the theorem holds for . Now, we prove the generalization of the BVPs,
For that purpose, we fixed .
For and , we have
By definition, we have
From (3.1), for , we have
Now that is the required result; thus, the theorem holds for .
Now assume that for , the DT is as follows:
For , we have
Thus, for , we have
For , we have
From (3.1), for , we have
Note that the theorem holds for .
Example 3.2. By using DTM, we solve the following nonlinear equation of order three for fifth-order BVP:
subject to the boundary conditions
By using Theorem 3.1, (3.22) is transformed to the following equation:
By using (2.1) at , we obtain the following transformed boundary conditions:
where and .
By using the transformed equation (3.24) together with the transformed boundary conditions (3.25), we can easily solve , for . The constants and can be evaluated by using boundary conditions (3.23) at :
These equations give and .
Finally, the following series solution can be formed by applying the inverse transformation in (2.2) up to :
Example 3.3. Next, we solve the following fourth-order nonlinear for fifth-order BVP by using DTM:
subject to the boundary conditions
By using Theorem 3.1, (3.28) is transformed to the following equation:
By using (2.1) at , we obtain the following transformed boundary conditions:
where and .
By using the transformed equation (3.30) together with the transformed boundary conditions (3.31), we can easily solve , for . The constants and can be evaluated by using boundary conditions (3.29) at :
These equations give and .
Finally, the following series solution can be formed by applying the inverse transformation in (2.2) up to :
Example 3.4. We perform the following third order of nonlinear equation for sixth-order BVP by using DTM:
subject to the boundary conditions
By using Theorem 3.1, (3.34) is transformed to the following equation:
By using (2.1) at , we obtain the following transformed boundary conditions:
where , , and .
By using the transformed equation (3.36) together with the transformed boundary conditions (3.37), we can easily solve , for . The constants and can be evaluated by using boundary conditions (3.35) at :
These equations give , , and .
Consequently, the following series solution can be formed by applying the inverse transformation equation in (2.2) up to :
4. Adomian Decomposition Method
Consider the differential equation where is invertible and is the highest-order derivative, is a linear differential operator, where the order of must be greater than , is a nonlinear operator, and is the source term. By using the given conditions and applying the inverse operator to both sides of (4.1), we get the following equation: where the function is arising from integrating the source term and from applying the given conditions which are prescribed.
The series solution of by the standard Adomian method is given as follows: On the other hand, the nonlinear function by an infinite series of polynomial is given as follows: The components are determined by using the following relation: for .
From the above equations, we observe that the component is identified by the function . We can apply modification by assuming that the function is divisible into two parts, and . Let the function be as follows: Under this assumption, we have a slight variation for components and , where assigned to and is combined with the other terms in (4.5) to assign . The modified recursive algorithm is as follows: for . There are several rules that are needed to follow for Adomian polynomials of nonlinear operator : and so on; see [3].
For comparison purpose, we solve the boundary value problems in Examples 3.2, 3.3, and 3.4 by using the Adomian decomposition method.
Example 4.1. By using ADM, we solve the following nonlinear equation of order three for fifth-order BVP: subject to the boundary conditions Equation (4.10) can be rewritten in operator form as follows: Operating with fivefold integral operator on (4.12) and using the boundary conditions at , we obtain the following equation: Then, determine the constants , . Substituting the decomposition series (4.3) for and the series of polynomials (4.4) into (4.13) yields Then, we split the terms into two parts which are assigned to and that are not included under in (4.14). We can obtain the following recursive relation: To determine the constants and , we use the boundary conditions in (4.11) at on the four-term approximant , where Then, solving the above equation yields Then we get the series solution as follows:
Example 4.2. Now, we solve the following fourth-order nonlinear for fifth-order boundary value problems by using ADM: subject to the boundary conditions Equation (4.19) can be rewritten in operator form as follows: Operating with fivefold integral operator on (4.21) and using the boundary conditions at , we obtain the following equation: Then, determine the constants , . Substituting the decomposition series (4.3) for and the series of polynomials (4.4) into (4.22) yields Then, we split the terms into two parts which are assigned to and that are not included under in (4.23). We can obtain the following recursive relation: To determine the constants and , we use the boundary conditions in (4.20) at on the four-term approximant , where Then, solving the above equation yields Finally, we get the series solution as follows:
Example 4.3. Finally, we perform third order of nonlinear function for sixth-order BVP for ADM as follows: subject to the boundary conditions Equation (4.28) can be rewritten in operator form as follows: Operating with sixfold integral operator on (4.30) and using the boundary conditions at , we obtain the following equation: Then, determine the constants , , and . Substituting the decomposition series (4.3) for and the series of polynomials (4.4) into (4.31) yields Then similarly, we split the terms into two parts which are assigned to and that are not included under in (4.32). We can obtain the following recursive relation: To determine the constants , , and , we use the boundary conditions in (4.29) at on the four-term approximant , where Then, solving the above equation yields Finally, we get the series solution as follows:
5. Results
We provide the results of the given examples in Tables 1, 2, and 3.
6. Conclusion
From the results, the proposed method, as well as the Differential Transformation Method, is more accurate than the Adomian Decomposition Method. The errors between the solutions of Differential Transformation Method and the exact solutions are smaller compared to the errors between the solutions of Adomian Decomposition Method and the exact solutions. In addition, Differential Transformation Method also shows less computational effort because it needs less time in calculation. Besides that, it is hard to calculate Adomian polynomials. From the results we obtained, it can reinforce conclusion made by many researchers that Differential Transformation Method is more efficient and accurate than Adomian Decomposition Method. Therefore, we can conclude that Differential Transformation Method is applicable for such problems in the bounded domains. The computations in all examples were performed by using Maple 13.
Acknowledgments
The authors gratefully acknowledge that this research was partially supported by the University Putra Malaysia under the Research University Grant scheme 05-01-09-0720RU and Fundamental Research Grant scheme 01-11-09-723FR. The authors also thank the referee(s) for very constructive comments and suggestions that improved the paper.