Research Article  Open Access
Reproducing Kernel Method for Solving Nonlinear DifferentialDifference Equations
Abstract
On the basis of reproducing kernel Hilbert spaces theory, an iterative algorithm for solving some nonlinear differentialdifference equations (NDDEs) is presented. The analytical solution is shown in a series form in a reproducing kernel space, and the approximate solution is constructed by truncating the series to terms. The convergence of to the analytical solution is also proved. Results obtained by the proposed method imply that it can be considered as a simple and accurate method for solving such differentialdifference problems.
1. Introduction
Differentialdifference equations play a crucial role in modelling of much physical phenomena such as particle vibrations in lattices, currents in electrical networks, pulses in biological chains, discretization in solid state, quantum physics, textile engineering, stratified hydrostatic flows, and so on. See, for example [1–14].
Recently, there have been lots of efforts in giving exact or approximate solutions of NDDEs. For instance, Zhu [4] extended the Expfunction method for solving NDDEs. In this method, exact solutions are sought in the form of an exponential type rational function in which both the numerator and denominator are polynomials of exponential functions. Mokhtari [5] applied the variational iteration method in which a correction functional is established by a general Lagrange multiplier, which can be identified optimally via the variational theory. The method gives rapidly convergent successive approximations of the exact solution. Qian et al. [6] extended the multilinear variable separation approach to a special differentialdifference equation. Baldwin et al. [7] presented an algorithm to find exact travelling wave solutions of differentialdifference equations in terms of the tanh function and found kinktype solutions in many spatially discrete nonlinear models. Wu et al. [8] and Wang et al. [9] extended the Adomian decomposition method to solve NDDEs. The technique is based on a decomposition of the solution of NDDE in a series of functions. Each term of the series is obtained from a polynomial generated from an expansion of an analytic function into a power series. Many authors [10–14] applied the homotopy perturbation method to solve NDDEs. The method, which is based on homotopy, constructs a continuous mapping of an initial guess through an auxiliary linear operator, and an auxiliary parameter is used to ensure the convergence of solution series.
The theory of reproducing kernels [15] was used for the first time at the beginning of the th century by S. Zaremba in his work on boundary value problems for harmonic and biharmonic functions. This theory has been successfully applied to fractal interpolation [16], solving ordinary differential equations [17–20] and partial differential equations [21, 22]. The book [23] provides excellent overviews of the existing reproducing kernel methods for solving various model problems such as integral and integrodifferential equations.
In this study, a general technique is proposed for solving some NDDEs in the reproducing kernel space. The main idea is to construct the reproducing kernel space satisfying the conditions for determining solution of the NDDEs. The analytical solution is represented in the form of series through the function value at the right side of the equation. For illustration, we apply this method to the Volterra lattice equation, the discretized mKdV lattice equation and the discrete sineGordon equation. The advantages of the approach lie in the following facts. The approximate solution converges uniformly to the analytical solution . The method is meshfree, easily implemented and capable in treating various boundary conditions. Since the method needs no time discretization, there is no matter, in which time the approximate solution is computed, from the both elapsed time and stability problem, point of views. Also we can evaluate the approximate solution for fixed once and use it over and over.
In the next section we describe how to solve a NDDE through the reproducing kernel method and verify convergence of the approximate solution to the exact solution. Several numerical results are presented in Section 3. The last section is a brief conclusion.
2. Construction of the Method
Consider the following NDDE where is a linear differential operator of oneorder derivative in , is a nonlinear differential operator and . Since nonhomogeneous initial condition can be reduced readily to the homogeneous one, we only consider the following homogeneous initial condition We assume that (2.1)(2.2) has a unique solution.
Definition 2.1. are absolutely continuous in , . The inner product and norm in are defined, respectively, by
Remark 2.2. The space is a reproducing kernel space, and its reproducing kernel is given by
Definition 2.3. is absolutely continuous in . The inner product and norm in are defined, respectively, by
Remark 2.4. The space is a reproducing kernel space and its reproducing kernel is given by
Remark 2.5. The definition of the spaces and is suitable for our numerical examples and in general is defined on the basis of the operator , and must be defined based on the operator .
In order to represent the analytical solution of the model problem, we can assume that is an invertible bounded linear operator, choose a countable dense subset in and define where is the adjoint operator of . The orthonormal system of can be derived from the GramSchmidt orthogonalization process of as
Theorem 2.6. Suppose that is dense in , then is a complete system in and .
Proof. See [22].
Theorem 2.7. If is dense in , then the analytical solution of (2.1)(2.2) will be
Proof. Since is a complete system in , we have
Now the approximate solution can be obtained by the term intercept of the analytical solution , that is, Obviously, In the sequel, for evaluating (2.11), we choose nonnegative integer and put the initial function . Then the approximate solution is defined by where
Lemma 2.8. If , is bounded in (2.13), and is continuous, then .
Proof. Note that On the other hand, we have From and former statements, it follows that Therefore, because of the boundedness of . The continuity of implies that
Theorem 2.9. Suppose that is bounded in (2.13), and (2.1)(2.2) has a unique solution. If is dense in , then the term approximate solution derived from the above method converges to the analytical solution of (2.1)(2.2) and where is given by (2.14).
Proof. At first, we prove the convergence of . From (2.13), we infer that The orthonormality of yields that In terms of (2.22), it holds that . Due to the condition that is bounded, is convergent and there exists a constant such that . If , and then On account of , consequently, we have The completeness of shows that as . Now, we prove that is the solution of (2.1)(2.2). Taking limits in (2.13), we get Note that Therefore, consider In view of (2.14), we have Since is dense in , for each , there exist a subsequence such that We have known that Let , by Lemma 2.8 and the continuity of , we have which indicates that satisfies (2.1)(2.2).
3. Numerical Results
To test the accuracy of the proposed method, three examples are treated in this section. The results are compared with the exact solutions. All experiments are done by taking .
Example 3.1. Consider the Volterra equation of the form With the initial condition , it is easy to check that the exact solution will be Taking , relative error (RE) of numerical results for , , and are reported in Table 1. Accuracy of approximate solutions is getting better as increases.

Example 3.2. Consider the discrete mKdV equation
With the initial condition , where and are arbitrary constants, the exact solution will be
Taking , and and choosing and , we have listed relative error of approximate solutions at some nodal points in Table 2. Numerical results are in good agreement with the exact solutions.

Example 3.3. Consider the discrete sineGordon equation of the form
With the initial condition where and are arbitrary constants, the exact solution will be
Taking , and and choosing and , we have given relative error of approximate solutions at some nodal points in Table 3. Numerical results are in good agreement with the exact solutions.

4. Conclusion
In this work, we proposed an algorithm for solving a class of nonlinear differentialdifference equations on the basis of reproducing kernel spaces. Results of numerical examples show that the present method is an accurate and reliable analyticalnumerical technique for solving such differentialdifference equations. The method is shown to be of good convergence, simple in principle, and easy to program.
References
 M. Toda, Theory of Nonlinear Lattices, vol. 20 of Springer Series in SolidState Sciences, Springer, Berlin, Germany, 1981. View at: Publisher Site  Zentralblatt MATH
 G.C. Wu, L. Zhao, and J.H. He, “Differentialdifference model for textile engineering,” Chaos, Solitons and Fractals, vol. 42, no. 1, pp. 352–354, 2009. View at: Publisher Site  Google Scholar
 J.H. He and E. W. M. Lee, “Variational principle for the differentialdifference system arising in stratified hydrostatic flows,” Physics Letters Section A, vol. 373, no. 1819, pp. 1644–1645, 2009. View at: Publisher Site  Google Scholar
 S.D. Zhu, “Discrete $(2+1)$dimensional Toda lattice equation via Expfunction method,” Physics Letters A, vol. 372, no. 5, pp. 654–657, 2008. View at: Publisher Site  Google Scholar
 R. Mokhtari, “Variational iteration method for solving nonlinear differentialdifference equations,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 9, no. 1, pp. 19–23, 2008. View at: Google Scholar
 X.M. Qian, S. Y. Lou, and X.B. Hu, “Variable separation approach for a differentialdifference system: special Toda equation,” Journal of Physics A, vol. 37, no. 6, pp. 2401–2411, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 D. Baldwin, Ü. Göktaş, and W. Hereman, “Symbolic computation of hyperbolic tangent solutions for nonlinear differentialdifference equations,” Computer Physics Communications, vol. 162, no. 3, pp. 203–217, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 L. Wu, L.D. Xie, and J.F. Zhang, “Adomian decomposition method for nonlinear differentialdifference equations,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 1, pp. 12–18, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 Z. Wang, L. Zou, and Z. Zong, “Adomian decomposition and Padé approximate for solving differentialdifference equation,” Applied Mathematics and Computation, vol. 218, no. 4, pp. 1371–1378, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 A. Yildirim, “Exact solutions of nonlinear differentialdifference equations by he's homotopy perturbation method,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 9, no. 2, pp. 111–114, 2008. View at: Google Scholar
 A. Yildirim, “He's homotopy perturbation method for nonlinear differentialdifference equations,” International Journal of Computer Mathematics, vol. 87, no. 5, pp. 992–996, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 M. A. Abdou, “New applications of He's homotopy perturbation method for nonlinear differential difference equations,” Physica Scripta, vol. 81, no. 1, Article ID 015003, 2010. View at: Publisher Site  Google Scholar
 M. M. Mousa and A. Kaltayev, “Homotopy perturbation method for solving nonlinear differentialdifference equations,” Zeitschrift fur Naturforschung Section A, vol. 65, no. 67, pp. 511–517, 2010. View at: Google Scholar
 M. M. Mousa, A. Kaltayev, and H. Bulut, “Extension of the homotopy perturbation method for solving nonlinear differentialdifference equations,” Zeitschrift fur Naturforschung Section A, vol. 65, no. 12, pp. 1060–1064, 2010. View at: Google Scholar
 N. Aronszajn, “Theory of reproducing kernels,” Transactions of the American Mathematical Society, vol. 68, pp. 337–404, 1950. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 P. Bouboulis and M. Mavroforakis, “Reproducing kernel Hilbert spaces and fractal interpolation,” Journal of Computational and Applied Mathematics, vol. 235, no. 12, pp. 3425–3434, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 F. Geng, “Solving singular second order threepoint boundary value problems using reproducing kernel Hilbert space method,” Applied Mathematics and Computation, vol. 215, no. 6, pp. 2095–2102, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 F. Geng, “A novel method for solving a class of singularly perturbed boundary value problems based on reproducing kernel method,” Applied Mathematics and Computation, vol. 218, no. 8, pp. 4211–4215, 2011. View at: Publisher Site  Google Scholar
 F. Geng and M. Cui, “A reproducing kernel method for solving nonlocal fractional boundary value problems,” Applied Mathematics Letters, vol. 25, no. 5, pp. 818–823, 2012. View at: Publisher Site  Google Scholar
 F. Z. Geng and X. M. Li, “A new method for Riccati differential equations based on reproducing kernel and quasilinearization methods,” Abstract and Applied Analysis, vol. 2012, Article ID 603748, 8 pages, 2012. View at: Publisher Site  Google Scholar
 M. Cui and F. Geng, “A computational method for solving onedimensional variablecoefficient Burgers equation,” Applied Mathematics and Computation, vol. 188, no. 2, pp. 1389–1401, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 M. Mohammadi and R. Mokhtari, “Solving the generalized regularized long wave equation on the basis of a reproducing kernel space,” Journal of Computational and Applied Mathematics, vol. 235, no. 14, pp. 4003–4014, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 M. Cui and Y. Lin, Nonlinear Numerical Analysis in the Reproducing Kernel Space, Nova Science, New York, NY, USA, 2009.
Copyright
Copyright © 2012 Reza Mokhtari 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.