- About this Journal ·
- Abstracting and Indexing ·
- Advance Access ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents

Abstract and Applied Analysis

Volume 2014 (2014), Article ID 520467, 7 pages

http://dx.doi.org/10.1155/2014/520467

## Solution of the Differential-Difference Equations by Optimal Homotopy Asymptotic Method

Department of Mathematics, Abdul Wali Khan University, Mardan 23200, Pakistan

Received 20 January 2014; Revised 1 May 2014; Accepted 1 May 2014; Published 29 May 2014

Academic Editor: Fuding Xie

Copyright © 2014 H. Ullah et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

We applied a new analytic approximate technique, optimal homotopy asymptotic method (OHAM), for treatment of differential-difference equations (DDEs). To see the efficiency and reliability of the method, we consider Volterra equation in different form. It provides us with a convenient way to control the convergence of approximate solutions when it is compared with other methods of solution found in the literature. The obtained solutions show that OHAM is effective, simpler, easier, and explicit.

#### 1. Introduction

In physical science the DDEs play a vital role in modeling of the complex physical phenomena. The DDEs models are used in vibration of particles in lattices, the flow of current in a network, and the pulses in biological chains. The models containing DDEs have been investigated by numerical techniques such as discretizations in solid state physics and quantum mechanics. In the last decades most of the research work has been done on DDEs. Levi and Yamilov [1, 2] attribute their work to classification of DDEs and connection of integrable partial differential equation (PDEs) and DDEs [3, 4]. In [5–8] the exact solutions of the DDEs have been studied. For the solution of the DDEs, Zou et al. extended the homotopy analysis method (HAM) [9, 10]. Wang et al. extended Adomian decomposition method (ADM) for solving nonlinear difference-differential equations (NDDEs) and got a good accuracy with the analytic solution [11]. Recently, Marinca et al. introduced optimal homotopy asymptotic method (OHAM) [12–16] for the solution of nonlinear problems. The validity of OHAM is independent of whether or not the nonlinear problems contain small parameters.

The motivation of this paper is to extend OHAM for the solution of NDDEs. In [17–20] OHAM has been proved to be useful for obtaining an approximate solution of nonlinear boundary value problems of differential equations. In this work, we have proved that OHAM is also useful and reliable for the solution of NDDEs, hence showing its validity and great potential for the solution of NDDEs phenomenon in science and engineering.

In the succeeding section, the basic idea of OHAM [12–16] is formulated for the solution of boundary value problems. In Section 3, the effectiveness of the modified formulation of OHAM for NDDEs has been studied.

#### 2. Basic Mathematical Theory of OHAM for NDDE

Let us take OHAM to the following differential-difference equation: with boundary conditions: where is a differential operator, is an unknown function, and denote spatial and temporal independent variables, respectively, is the boundary of , and is a known analytic function. can be divided into two parts and such that where is simpler part of the partial differential equation which is easier to solve and contains the remaining part of . So (1) can be written as According to OHAM, one can construct an optimal homotopy satisfying where is an embedding parameter, is an unknown function, and is a nonzero auxiliary function. The auxiliary function is nonzero for and . Equation (5) is the structure of OHAM homotopy. Clearly we have Obviously, when and , we obtain respectively. Thus, as varies from to , the solution varies from to , where is the zeroth order solution and can be obtained from (10).

Expanding , , by Taylors series and choosing as given where , are constants and to be determined and .

Now substituting (9) into (5) and equating the coefficient of like powers of , we obtain We obtained the zeroth order solution and first, second, third, and the general order solutions by solving (10)–(14). In the same way the remaining solutions can be determined. It has been observed that the convergence of the series (9) depends on the auxiliary constants . If it is convergent at , one has Substituting (15) into (1), it results in the following expression for residual: If , then is the exact solution of the problem. Generally it does not happen, especially in nonlinear problems.

For the determinations of auxiliary constants, , , there are different methods like Galerkin’s method, Ritz method, least squares method, and collocation method. One can apply the method of least squares as follows: where and are two values, depending on the nature of the given problem.

The auxiliary constants can be optimally calculated as The order approximate solution can be obtained by these obtained constants.

The convergence of OHAM is directly proportional to the number of optimal constants .

#### 3. Application of OHAM to Difference-Differential Equations

To show the validity and effectiveness of OHAM formulation to difference-differential equation, we apply it to Voltera equations in different form.

*Model 1. *We consider Voltera equation of the form
with initial conditions
The exact solution of (3.1) is given by
and we take
According to (1), we have
Using (19), (20), and (22) into (5) and using the method discussed in Section 2 lead to the following. *Zeroth Order Problem*. Consider
with
from which we obtain
*First Order Problem*. One has
with
Its solution is
*Second Order Problem*. Consider
with
*Third Order Problem*. One has
with
Adding (26), (29), (32), and (35), we obtain
For the computation of the constants , , and applying the method of least square mentioned in (16)–(18), we get
Putting these values in (37), we obtained the approximate solution of the form
The exact solution of the problem (3.1) is [11]
The Adomian solution with Padé approximation of (3.1) is given by [11]

*Model 2. *Now consider another form of Volterra equation
with initial conditions
The exact solution of (42) is given by
We take
According to (1), we have
Using (42), (43), and (45) into (5) and using the method discussed in Section 2 lead to the following. *Zeroth Order Problem*. Consider
with
from which we obtain
*First Order Problem*. One has the following:
with
Its solution is

Second Order Problem. Consider
with
*Third Order Problem*. One has
with
Adding (49), (52), (55), and (58), we obtain
For the computation of the constants , , and applying the method of least square mentioned in (16)–(18), we get
The exact solution of problem (42) is [9]
The HAM solution of (42) is given by [9]

#### 4. Results and Discussions

The formulation presented in Section 2 provides highly accurate solutions for the problems demonstrated in Section 3. We have used Mathematica 7 for most of our computational work. Table 1 and Figure 1 give the comparisons of OHAM with ADM and exact solutions for Model 1. Also the absolute errors at different values of at are given (Figure 4). Figure 3 shows the comparisons of OHAM results with HAM, exact at different values of at for Model 2. The convergence of approximate orders to exact solutions is given in Table 2 for Model 1 and Figure 5 for Model 2. The residuals have been plotted for Model 1 in Figure 2 and for Model 2 in Figure 5. We have concluded that the results obtained by OHAM are identical to the results obtained by ADM, HAM, and exact. OHAM converges rapidly with increasing of the order of approximation.

#### 5. Conclusion

In this work, we have seen the effectiveness of OHAM [12–16] to differential-difference equations. By applying the basic idea of OHAM to differential-difference equations, we found it simpler in applicability, more convenient to control convergence, and involved less computational overhead. Therefore, OHAM shows its validity and great potential for the differential-difference equations arising in science and engineering.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgment

The authors are thankful to the anonymous reviewers for their constructive and valuable suggestions.

#### References

- D. Levi and R. Yamilov, “Conditions for the existence of higher symmetries of evolutionary equations on the lattice,”
*Journal of Mathematical Physics*, vol. 38, no. 12, pp. 6648–6674, 1997. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - R. I. Yamilov, “Construction scheme for discrete Miura transformations,”
*Journal of Physics A: Mathematical and General*, vol. 27, no. 20, pp. 6839–6851, 1994. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - I. Yu. Cherdantsev and R. I. Yamilov, “Master symmetries for differential-difference equations of the Volterra type,”
*Physica D: Nonlinear Phenomena*, vol. 87, no. 1–4, pp. 140–144, 1995. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. B. Shabat and R. I. Yamilov, “Symmetries of nonlinear lattices,”
*Leningrad Mathematical Journal*, vol. 2, p. 377, 1991. View at Google Scholar - D.-J. Zhang, “Singular solutions in Casoratian form for two differential-difference equations,”
*Chaos, Solitons and Fractals*, vol. 23, no. 4, pp. 1333–1350, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - K. Narita, “Soliton solution for a highly nonlinear difference-differential equation,”
*Chaos, Solitons and Fractals*, vol. 3, no. 3, pp. 279–283, 1993. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Z. Wang and H.-Q. Zhang, “New exact solutions to some difference differential equations,”
*Chinese Physics*, vol. 15, no. 10, pp. 2210–2215, 2006. View at Publisher · View at Google Scholar · View at Scopus - W. Zhen and H. Zhang, “A symbolic computational method for constructing exact solutions to difference-differential equations,”
*Applied Mathematics and Computation*, vol. 178, no. 2, pp. 431–440, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - L. Zou, Z. Zong, and G. H. Dong, “Generalizing homotopy analysis method to solve Lotka-Volterra equation,”
*Computers & Mathematics with Applications*, vol. 56, no. 9, pp. 2289–2293, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Z. Wang, L. Zou, and H. Zhang, “Applying homotopy analysis method for solving differential-difference equation,”
*Physics Letters A: General, Atomic and Solid State Physics*, vol. 369, no. 1-2, pp. 77–84, 2007. View at Publisher · View at Google Scholar · View at Scopus - Z. Wang, L. Zou, and Z. Zong, “Adomian decomposition and Padé approximate for solving differential-difference equation,”
*Applied Mathematics and Computation*, vol. 218, no. 4, pp. 1371–1378, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - N. Herişanu and V. Marinca, “Explicit analytical approximation to large-amplitude non-linear oscillations of a uniform cantilever beam carrying an intermediate lumped mass and rotary inertia,”
*Meccanica*, vol. 45, no. 6, pp. 847–855, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - V. Marinca and N. Herişanu, “Application of optimal homotopy asymptotic method for solving nonlinear equations arising in heat transfer,”
*International Communications in Heat and Mass Transfer*, vol. 35, no. 6, pp. 710–715, 2008. View at Publisher · View at Google Scholar · View at Scopus - V. Marinca, N. Herişanu, C. Bota, and B. Marinca, “An optimal homotopy asymptotic method applied to the steady flow of a fourth-grade fluid past a porous plate,”
*Applied Mathematics Letters*, vol. 22, no. 2, pp. 245–251, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - V. Marinca, N. Herisanu, and I. Nemes, “A new analytic approach to nonlinear vibration of an electrical machine,”
*Proceedings of the Romanian Academy*, vol. 9, pp. 229–236, 2008. View at Google Scholar - V. Marinca and N. Herişanu, “Determination of periodic solutions for the motion of a particle on a rotating parabola by means of the optimal homotopy asymptotic method,”
*Journal of Sound and Vibration*, vol. 329, no. 9, pp. 1450–1459, 2010. View at Publisher · View at Google Scholar · View at Scopus - M. Idrees, S. Islam, and A. M. Tirmizi, “Application of optimal homotopy asymptotic method of the Korteqag-de-Varies equation,”
*Computers & Mathematics with Applications*, vol. 63, pp. 695–707, 2012. View at Google Scholar - H. Ullah, S. Islam, M. Idrees, and M. Arif, “Solution of boundary layer problems with heat transfer by optimal homotopy asymptotic method,”
*Abstract and Applied Analysis*, vol. 2013, Article ID 324869, 10 pages, 2013. View at Publisher · View at Google Scholar · View at MathSciNet - R. Nawaz, H. Ullah, S. Islam, and M. Idrees, “Application of optimal homotopy asymptotic method to Burger equations,”
*Journal of Applied Mathematics*, vol. 2013, Article ID 387478, 8 pages, 2013. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - H. Ullah, S. Islam, M. Idrees, and R. Nawaz, “Application of optimal homotopy asymptotic method to doubly wave solutions of the coupled Drinfel’d-Sokolov-Wilson Equations,”
*Mathematical Problems in Engineering*, vol. 2013, Article ID 362816, 8 pages, 2013. View at Publisher · View at Google Scholar