Journal of Function Spaces and Applications

Volume 2012 (2012), Article ID 984682, 23 pages

doi:10.1155/2012/984682

Research Article

Explicit Solution of Telegraph Equation Based on Reproducing Kernel Method

Mustafa Inc,¹ Ali Akgül,² and Adem Kiliçman³

¹Department of Mathematics, Science Faculty, Firat University, 23119 Elazig, Turkey

²Department of Mathematics, Education Faculty, Dicle University, 21280 Diyarbakir, Turkey

³Department of Mathematics and Institute for Mathematical Research, University Putra Malaysia, 43400 Serdang, Selangor, Malaysia

Received 5 August 2012; Revised 21 November 2012; Accepted 21 November 2012

Academic Editor: Gabriel N. Gatica

Copyright © 2012 Mustafa Inc 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 propose a reproducing kernel method for solving the telegraph equation with initial conditions based on the reproducing kernel theory. The exact solution is represented in the form of series, and some numerical examples have been studied in order to demonstrate the validity and applicability of the technique. The method shows that the implement seems easy and produces accurate results.

1. Introduction

In this paper, we consider the telegraph equation of the following form:

over a region and , are known constant coefficients with initial conditions as

where can be voltage or current through the wire at position and time . In (1.1), we have

where is conductance of resistor, is resistance of resistor, is inductance of capacitor, is capacitance of capacitor, and can be considered as a function depending on distance and time , and constants are depending on a given problem and , and are known continuous functions.

The hyperbolic partial differential equations model the vibrations of structures (e.g., buildings, beams, and machines) and are the basis for fundamental equations of atomic physics. Equation (1.1), referred to as the second-order telegraph equation with constant coefficients, models a mixture between diffusion and wave propagation by introducing a term that accounts for effects of finite velocity to the standard heat or mass transport equation [1]. However, (1.1) is commonly used in signal analysis for transmission and propagation of electrical signals [2, 3].

In recent years, much attention has been given in the literature to the development, analysis, and implementation of stable methods for the numerical solution of second-order hyperbolic equations, see, for example, [4–11]. These methods are conditionally stable. In [12], Mohanty carried over a new technique to solve (1.1), which is unconditionally stable and is of second-order accuracy in both the time and space components. Mohebbi and Dehghan [13] presented a high-order accurate method for solving one-space-dimensional linear hyperbolic equations and proved the high-order accuracy due to the fourth-order discretization of spatial derivative and unconditional stability. A compact finite difference approximation was presented in [14] by using the fourth order discretizing spatial derivatives of the linear hyperbolic equation and collocation method for the time component. Another solution is approximated by suing a polynomial at each grid point such that its coefficients were determined by solving a linear system of equations [15]. By using collocation points and approximating the solution by using a thin plate splines radial basis function was presented in [16].

In [17], the author used the Chebyshev cardinal functions. Lakestani and Saray[18] used interpolating scaling function. Ding et al. [19] constructed a class of new difference scheme based on a new nonpolynomial spline method to solve (1.1) and (1.2). Lakoud and Belakroum [20] studied the existence and uniqueness of the solution with integral condition by using the Ro the time discretization method. Dehghan et al. [21] used to compute the solution for the linear, variable coefficient, fractional derivative, and multispace telegraph equations by using the variational iteration method. Further, Biazar et al. [22] obtained an approximate solution by using the variational iteration method. Recently, Yao and Lin [23] investigated a nonlinear hyperbolic telegraph equation with an integral condition by reproducing kernel space at in (1.1). In [24], Yousefi presented a numerical method by using Legendre multiwavelet Galerkin method.

In this paper, the RKHSM [25–47] will be used to investigate the telegraph equation (1.1). Several researches have been devoted to the application of RKHSM to a wide class of stochastic and deterministic problems involving fractional differential equation, nonlinear oscillator with discontinuity, singular nonlinear two-point periodic boundary value problems, integral equations, and nonlinear partial differential equations [27–41]. The method is well suited to physical problems.

The efficiency of the method was used by many authors to investigate several scientific applications. Geng and Cui [27] applied the RKHSM to handle the second-order boundary value problems. Yao and Cui [28] and Wang et al. [29] investigated a class of singular boundary value problems by this method and the obtained results were good. Zhou et al. [30] used the RKHSM effectively to solve second-order boundary value problems. In [31], the method was used to solve nonlinear infinite-delay-differential equations. Wang and Chao [32], Li and Cui [33], and Zhou and Cui [34] independently employed the RKHSM to variable-coefficient partial differential equations. Geng and Cui [35] and Du and Cui [36] investigated to the approximate solution of the forced Duffing equation with integral boundary conditions by combining the homotopy perturbation method and the RKHSM. Lv and Cui [37] presented a new algorithm to solve linear fifth-order boundary value problems. In [38, 39], authors developed a new existence proof of solutions for nonlinear boundary value problems. Cui and Du [40] obtained the representation of the exact solution for the nonlinear Volterra-Fredholm integral equations by using the reproducing kernel space. Wu and Li [41] applied iterative reproducing kernel method to obtain the analytical approximate solution of a nonlinear oscillator with discontinuities. Recently, the method was applied the fractional partial differential equations and the multipoint boundary value problems [42–45]. For more details about RKHSM and the modified forms and their effectiveness, see [25–47].

In the present work, we use the following equation:

by transformation for homogeneous initial conditions of (1.1) and (1.2), we get as follows (1.5):

where

The paper is organized as follows. Section 2 is devoted to several reproducing kernel spaces and a linear operator is introduced. The solution representation in has been presented in Section 3. We prove that the approximate solution converges to the exact solution uniformly. Some numerical examples are illustrated in Section 4. We provide some conclusions in the last sections.

2. Preliminaries

Hilbert spaces can be completely classified: there is a unique Hilbert space up to isomorphism for every cardinality of the base. Since finite-dimensional Hilbert spaces are fully understood in linear algebra, and since morphisms of Hilbert spaces can always be divided into morphisms of spaces with Aleph-null () dimensionality, functional analysis of Hilbert spaces mostly deals with the unique Hilbert space of dimensionality Aleph-null and its morphisms. One of the open problems in functional analysis is to prove that every bounded linear operator on a Hilbert space has a proper invariant subspace. Many special cases of this invariant subspace problem have already been proven [48].

2.1. Reproducing Kernel Spaces

In this section, we define some useful reproducing kernel spaces.

Definition 2.1 (reproducing kernel). Let be a nonempty set. A function is called a reproducing kernel of the Hilbert space if and only if

(a) for all ,(b) for all and all .

The last condition is called “the reproducing property” as the value of the function at the point is reproduced by the inner product of with .

Then we need some notation that we use in the development of the paper. In the next we define several spaces with inner product over those spaces. Thus the space defined as

is a Hilbert space. The inner product and the norm in are defined by

respectively. Thus the space is a reproducing kernel space, that is, for each fixed and any , there exists a function such that

and similarly we define the space

The inner product and the norm in are defined by

respectively. The space is a reproducing kernel Hilbert space and its reproducing kernel function is given by [26] as

and the space

is a Hilbert space where the inner product and the norm in are defined by

respectively. The space is a reproducing kernel space and its reproducing kernel function is given by [26] as

Similarly, the space defined by

is a Hilbert space and then inner product and the norm in are defined by

respectively. The space is a reproducing kernel space and its reproducing kernel function is given by [26] as

Now we have the following theorem.

Theorem 2.2. The space is a complete reproducing kernel space whose reproducing kernel is given by

where

Proof. Since

through iterative integrations by parts for (2.15), we have

Note the property of the reproducing kernel as

If

Then by (2.16) we obtain

when ,

therefore

Since

we have

From (2.18) and (2.23), the unknown coefficients and can be obtained. Thus is given by

Now we note that the space given in [26] as

is a binary reproducing kernel Hilbert space. The inner product and the norm in are defined by

respectively.

Theorem 2.3. The is a reproducing kernel space and its reproducing kernel function is

such that for any ,

Similarly, the space

is a binary reproducing kernel Hilbert space. The inner product and the norm in are defined by [26] as

respectively. is a reproducing kernel space and its reproducing kernel function is

3. Solution Representation in W(Ω)

In this section, the solution of (1.1) is given in the reproducing kernel space . We define the linear operator as

Model problem (1.1) changes to the following problem:

Lemma 3.1. The operator is a bounded linear operator.

Proof. Since

by continuity of , we have

Now similarly for , we obtain

and then

Therefore we conclude

Now, if we choose a countable dense subset in and define

where is the adjoint operator of , then the orthonormal system of can be derived from the process of Gram-Schmidt orthogonalization of as

Then we have the following theorem.

Theorem 3.2. Suppose that is dense in , then is complete system in and

Proof. We have

That is clearly . For each fixed , if

then

Note that is dense in , hence, . It follows that from the existence of . So the proof is complete.

Theorem 3.3. If is dense in , then the solution of (1.5) is given as

Proof. Since is complete system in , we have

Now the approximate solution can be obtained from the term intercept of the exact solution and

Of course it is also easy to show that

3.1. Convergence Analysis

We assume that is dense in . We discuss the convergence of the approximate solutions constructed in Section 3. Let be the exact solution of (1.1) and the term approximation solution of (1.1). Then we have the following theorem.

Theorem 3.4. If , then

Moreover a sequence is monotonically decreasing in .

Proof. From (3.14) and (3.16), it follows that

Thus

In addition,

Then clearly, is monotonically decreasing in .

4. Experimental Results for the Telegraph Equation

In this section, three numerical examples are provided to show the accuracy of the present method. All the computations were performed by Maple 13. Since, the RKHSM does not require discretization of the variables, that is, time and space, it is also not effected by computation round-off errors and no need to face with necessity of large computer memory and time. The accuracy of the RKHSM for the problem (1.1) is controllable and absolute errors are small with present choice of and (see Tables 1, 2, 3, 4, 5, and 6). Thus the numerical results which we obtain justify the advantage of this methodology.

Table 1: The exact solution, approximate solution, absolute error, and relative error of Example 4.1 for initial conditions at .

Exact solution
Approximate solution
Absolute error
Relative error

0.2
0.5877852524
0.5877852576

0.4
0.9510565165
0.951056525

0.6
0.9510565163
0.95105657

0.8
0.5877852522
0.587785294

Table 2: The absolute error ofExample 7for difference schemes and our scheme at .

Difference scheme [13]
Difference scheme [21]
Present scheme

0.2

0.4

0.6

0.8

Table 3: The exact solution, approximate solution, absolute error, and relative error of Example 4.2 for initial conditions at .

Exact solution
Approximate solution
Absolute error
Relative error

0.2
0.04761212907
0.047612168

0.4
0.07703804312
0.07703791

0.000001727977433

0.6
0.07703804310
0.077037983

0.8
0.04761212906
0.04761256075

0.000009066807314

Table 4: The absolute error of Example 4.2 for difference schemes and our scheme at .

Difference scheme [13]
Difference scheme [21]
Our scheme

0.2

0.4

0.6

0.8

Table 5: The exact solution, approximate solution, absolute error, and relative error of Example 4.3 for initial conditions at .

Exact solution
Approximate solution
Absolute error
Relative error

0.2
0.004842491805
0.00484215485

0.00006958297785

0.4
0.02451511476
0.024515542

0.00001742761575

0.6
0.02451511476
0.02451574136

0.00002555974166

0.8
0.004842491805
0.004843445

0.0001968397756

Table 6: The absolute error of Example 4.3 for difference schemes and our scheme at .

Difference scheme [13]
Difference scheme [21]
Our scheme

0.2

0.4

0.6

0.8

Note that the solutions are very rapidly convergent by utilizing the RKHSM. Further, the series solution methodology can be applied to various types of linear or nonlinear system of partial differential equations and single partial differential equations, see, for example, [25–30].

Using our method we choose points in . In Tables 2, 4 and 6 we compute the absolute errors at the following points:

In Table 7, we compute the following relative errors:

at the points . It is possible to refine the result by increasing the intensive points.

Table 7: The relative errors and computational times for Examples 4.1–4.3 at different values of and .

Example 5.1
Time
Example 5.2
Time
Example 5.3
Time

0.1
0.1

1.685

1.809

1.716

0.2
0.2

1.857

1.732

1.872

0.3
0.3

1.841

1.934

1.669

0.4
0.4

1.857

1.872

1.716

0.5
0.5

1.622

1.841

1.716

0.6
0.6

1.826

1.654

1.701

0.7
0.7

1.731

1.560

1.638

0.8
0.8

1.747

1.591

1.763

0.9
0.9

1.701

1.653

1.669

1.0
1.0

1.623

1.513

1.732

We constructed the figures for different values of and . It can be concluded by figures that the speed of convergence is decreasing by increasing the values of and .

Example 4.1. Consider the following telegraph equation with initial conditions:

Then the exact solution is given as

If we apply to (4.3), then we obtain the following equation:

Then we have the estimation in Table 1.

Now if we compare [13, 21] and the present scheme we have Table 2.

We have Figures 1, 2, and 3 for this example, where ES = exact solution and AS = approximate solution.

Figure 1: The compared results of the analytical and numerical solutions of Example 4.1 for .

Figure 2: The compared results of the analytical and numerical solutions of Example 4.1 for .

Figure 3: The compared results of the analytical and numerical solutions of Example 4.1 for .

Example 4.2. Consider the following telegraph equation with initial conditions:

The exact solution . If we apply to (4.6), then we obtain

then similarly, in Table 3 we have the estimation among the exact and approximate solutions and the error terms.

Then the comparison yields in Table 4.

We have Figures 4, 5, and 6 for this example.

Figure 4: The compared results of the analytical and numerical solutions of Example 4.2 for .

Figure 5: The compared results of the analytical and numerical solutions of Example 4.2 for .

Figure 6: The compared results of the analytical and numerical solutions of Example 4.2 for .

Example 4.3. Consider the following telegraph equation with initial conditions:

The exact solution . If we apply

to (4.8) then (4.10) is obtained as

We have Figures 7, 8, 9, 10, and 11 for this example.

Figure 7: The compared results of the analytical and numerical solutions of Example 4.3 for .

Figure 8: The compared results of the analytical and numerical solutions of Example 4.3 for .

Figure 9: The compared results of the analytical and numerical solutions of Example 4.3 for .

Figure 10: The compared results of the analytical and numerical solutions of Examples 4.1–4.3.

Figure 11: The compared results of the relative errors for Examples 4.1–4.3.

Remark 4.4. Ding et al. [19] has solved Examples 4.1–4.3 by using new polynomial spline methods. In our method the numerical solutions are in good agreement with analytical solutions (see Figures 1–11). In addition, the figure of the relative error is drawn for the value of and also for three examples (see Figure 11).

One may view Tables 1, 2, 3, 4, 5 and 6 for the confidence of the method and the comparison with the other methods. In Table 7, computing time with relative error is also given for each example.

5. Conclusion

In this paper, the RKHSM was used for the telegraph equation with initial conditions. The approximate solutions to the equations have been calculated by using the RKHSM without any need to transformation techniques and linearization or perturbation of the equations. In closing, the RKHSM avoids the difficulties and massive computational work by determining the analytic solutions. We compare our solutions with the exact solutions and the results of [19].

A clear conclusion can be drawn from the numerical results as the RKHSM algorithm provides highly accurate numerical solutions without spatial discretizations for the nonlinear partial differential equations. It is also worth noting that the advantage of this methodology displays a fast convergence of the solutions. The illustrations show the dependence of the rapid convergence depends on the character and behavior of the solutions just as in a closed form solutions.

Acknowledgments

The authors express their sincere thanks to the referees for the careful and detailed reading of the paper and the very helpful suggestions that improved the paper substantially. The third author gratefully acknowledges that this paper was partially supported by the University Putra Malaysia under the ERGS Grant Scheme having project no. 5527068.

References

M. S. El-Azab and M. El-Gamel, “A numerical algorithm for the solution of telegraph equations,” Applied Mathematics and Computation, vol. 190, no. 1, pp. 757–764, 2007.

A. C. Metaxas and R. J. Meredith, Industrial Microwave, Heating, Peter Peregrinus, London, UK, 1993.

G. Roussy and J. A. Pearcy, Foundations and Industrial Applications of Microwaves and Radio Frequency Fields, John Wiley, New York, NY, USA, 1995.

M. Dehghan, “On the solution of an initial-boundary value problem that combines Neumann and integral condition for the wave equation,” Numerical Methods for Partial Differential Equations, vol. 21, no. 1, pp. 24–40, 2005.

R. K. Mohanty, M. K. Jain, and K. George, “On the use of high order difference methods for the system of one space second order nonlinear hyperbolic equations with variable coefficients,” Journal of Computational and Applied Mathematics, vol. 72, no. 2, pp. 421–431, 1996.

E. H. Twizell, “An explicit difference method for the wave equation with extended stability range,” Nordisk Tidskrift for Informationsbehandling (BIT), vol. 19, no. 3, pp. 378–383, 1979.

F. Gao and C. Chi, “Unconditionally stable difference schemes for a one-space-dimensional linear hyperbolic equation,” Applied Mathematics and Computation, vol. 187, no. 2, pp. 1272–1276, 2007.

A. Saadatmandi and M. Dehghan, “Numerical solution of the one-dimensional wave equation with an integral condition,” Numerical Methods for Partial Differential Equations, vol. 23, no. 2, pp. 282–292, 2007.

A. Borhanifar and R. Abazari, “An unconditionally stable parallel difference scheme for telegraph equation,” Mathematical Problems in Engineering, vol. 2009, Article ID 969610, 17 pages, 2009.

A. Saadatmandi and M. Dehghan, “Numerical solution of hyperbolic telegraph equation using the Chebyshev tau method,” Numerical Methods for Partial Differential Equations, vol. 26, no. 1, pp. 239–252, 2010.

M. A. Abdou, “Adomian decomposition method for solving the telegraph equation in charged particle transport,” Journal of Quantitative Spectroscopy and Radiative Transfer, vol. 95, no. 3, pp. 407–414, 2005.

R. K. Mohanty, “An unconditionally stable difference scheme for the one-space-dimensional linear hyperbolic equation,” Applied Mathematics Letters, vol. 17, no. 1, pp. 101–105, 2004.

A. Mohebbi and M. Dehghan, “High order compact solution of the one-space-dimensional linear hyperbolic equation,” Numerical Methods for Partial Differential Equations, vol. 24, no. 5, pp. 1222–1235, 2008.

M. Dehghan, “Finite difference procedures for solving a problem arising in modeling and design of certain optoelectronic devices,” Mathematics and Computers in Simulation, vol. 71, no. 1, pp. 16–30, 2006.

M. Dehghan, “Implicit collocation technique for heat equation with non-classic initial condition,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 7, no. 4, pp. 461–466, 2006.

M. Dehghan and A. Shokri, “A numerical method for solving the hyperbolic telegraph equation,” Numerical Methods for Partial Differential Equations, vol. 24, no. 4, pp. 1080–1093, 2008.

M. Dehghan and M. Lakestani, “The use of Chebyshev cardinal functions for solution of the second-order one-dimensional telegraph equation,” Numerical Methods for Partial Differential Equations, vol. 25, no. 4, pp. 931–938, 2009.

M. Lakestani and B. N. Saray, “Numerical solution of telegraph equation using interpolating scaling functions,” Computers & Mathematics with Applications, vol. 60, no. 7, pp. 1964–1972, 2010.

H.-F. Ding, Y.-X. Zhang, J.-X. Cao, and J.-H. Tian, “A class of difference scheme for solving telegraph equation by new non-polynomial spline methods,” Applied Mathematics and Computation, vol. 218, no. 9, pp. 4671–4683, 2012.

A. G. Lakoud and D. Belakroum, “Rothe's method for a telegraph equation with integral conditions,” Nonlinear Analysis. Theory, Methods & Applications, vol. 70, no. 11, pp. 3842–3853, 2009.

M. Dehghan, S. A. Yousefi, and A. Lotfi, “The use of He's variational iteration method for solving the telegraph and fractional telegraph equations,” International Journal for Numerical Methods in Biomedical Engineering, vol. 27, no. 2, pp. 219–231, 2011.

J. Biazar, H. Ebrahimi, and Z. Ayati, “An approximation to the solution of telegraph equation by variational iteration method,” Numerical Methods for Partial Differential Equations, vol. 25, no. 4, pp. 797–801, 2009.

H. Yao and Y. Lin, “New algorithm for solving a nonlinear hyperbolic telegraph equation with an integral condition,” International Journal for Numerical Methods in Biomedical Engineering, vol. 27, no. 10, pp. 1558–1568, 2011.

S. A. Yousefi, “Legendre multiwavelet Galerkin method for solving the hyperbolic telegraph equation,” Numerical Methods for Partial Differential Equations, vol. 26, no. 3, pp. 535–543, 2010.

M. Cui and Z. Deng, “Solutions to the define solution problem of differential equations in space ,” Advances in Mathematics, vol. 17, pp. 327–328, 1986.

M. Cui and Y. Lin, Nonlinear Numerical Analysis in the Reproducing Kernel Space, Nova Science Publishers, New York, NY, USA, 2009.

F. Geng and M. Cui, “Solving a nonlinear system of second order boundary value problems,” Journal of Mathematical Analysis and Applications, vol. 327, no. 2, pp. 1167–1181, 2007.

H. Yao and M. Cui, “A new algorithm for a class of singular boundary value problems,” Applied Mathematics and Computation, vol. 186, no. 2, pp. 1183–1191, 2007.

W. Wang, M. Cui, and B. Han, “A new method for solving a class of singular two-point boundary value problems,” Applied Mathematics and Computation, vol. 206, no. 2, pp. 721–727, 2008.

Y. Zhou, Y. Lin, and M. Cui, “An efficient computational method for second order boundary value problems of nonlinear differential equations,” Applied Mathematics and Computation, vol. 194, no. 2, pp. 354–365, 2007.

X. Lü and M. Cui, “Analytic solutions to a class of nonlinear infinite-delay-differential equations,” Journal of Mathematical Analysis and Applications, vol. 343, no. 2, pp. 724–732, 2008.

Y.-L. Wang and L. Chao, “Using reproducing kernel for solving a class of partial differential equation with variable-coefficients,” Applied Mathematics and Mechanics, vol. 29, no. 1, pp. 129–137, 2008.

F. Li and M. Cui, “A best approximation for the solution of one-dimensional variable-coefficient Burgers' equation,” Numerical Methods for Partial Differential Equations, vol. 25, no. 6, pp. 1353–1365, 2009.

S. Zhou and M. Cui, “Approximate solution for a variable-coefficient semilinear heat equation with nonlocal boundary conditions,” International Journal of Computer Mathematics, vol. 86, no. 12, pp. 2248–2258, 2009.

F. Geng and M. Cui, “New method based on the HPM and RKHSM for solving forced Duffing equations with integral boundary conditions,” Journal of Computational and Applied Mathematics, vol. 233, no. 2, pp. 165–172, 2009.

J. Du and M. Cui, “Solving the forced Duffing equation with integral boundary conditions in the reproducing kernel space,” International Journal of Computer Mathematics, vol. 87, no. 9, pp. 2088–2100, 2010.

X. Lv and M. Cui, “An efficient computational method for linear fifth-order two-point boundary value problems,” Journal of Computational and Applied Mathematics, vol. 234, no. 5, pp. 1551–1558, 2010.

W. Jiang and M. Cui, “Constructive proof for existence of nonlinear two-point boundary value problems,” Applied Mathematics and Computation, vol. 215, no. 5, pp. 1937–1948, 2009.

J. Du and M. Cui, “Constructive proof of existence for a class of fourth-order nonlinear BVPs,” Computers & Mathematics with Applications, vol. 59, no. 2, pp. 903–911, 2010.

M. Cui and H. Du, “Representation of exact solution for the nonlinear Volterra-Fredholm integral equations,” Applied Mathematics and Computation, vol. 182, no. 2, pp. 1795–1802, 2006.

B. Wu and X. Li, “Iterative reproducing kernel method for nonlinear oscillator with discontinuity,” Applied Mathematics Letters, vol. 23, no. 10, pp. 1301–1304, 2010.

W. Jiang and Y. Lin, “Representation of exact solution for the time-fractional telegraph equation in the reproducing kernel space,” Communications in Nonlinear Science and Numerical Simulation, vol. 16, no. 9, pp. 3639–3645, 2011.

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.

Y. Lin and M. Cui, “A numerical solution to nonlinear multi-point boundary value problems in the reproducing kernel space,” Mathematical Methods in the Applied Sciences, vol. 34, no. 1, pp. 44–47, 2011.

F. Z. Geng, “A numerical algorithm for nonlinear multi-point boundary value problems,” Journal of Computational and Applied Mathematics, vol. 236, no. 7, pp. 1789–1794, 2012.

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.

B. Y. Wu and X. Y. Li, “A new algorithm for a class of linear nonlocal boundary value problems based on the reproducing kernel method,” Applied Mathematics Letters, vol. 24, no. 2, pp. 156–159, 2011.

L. M. Surhone, M. T. Tennoe, and S. F. Henssonow, Reproducing Kernel Hilbert Space, Betascript Publishing, Berlin, Germany, 2010.