`Abstract and Applied AnalysisVolume 2012 (2012), Article ID 258067, 10 pageshttp://dx.doi.org/10.1155/2012/258067`
Research Article

The Combined RKM and ADM for Solving Nonlinear Weakly Singular Volterra Integrodifferential Equations

School of Mathematics and Sciences, Harbin Normal University, Heilongjiang, Harbin, 150025, China

Received 2 July 2012; Accepted 12 November 2012

Copyright © 2012 Xueqin Lv and Sixing Shi. 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

The reproducing kernel method (RKM) and the Adomian decomposition method (ADM) are applied to solve th-order nonlinear weakly singular Volterra integrodifferential equations. The numerical solutions of this class of equations have been a difficult topic to analyze. The aim of this paper is to use Taylor’s approximation and then transform the given th-order nonlinear Volterra integrodifferential equation into an ordinary nonlinear differential equation. Using the RKM and ADM to solve ordinary nonlinear differential equation is an accurate and efficient method. Some examples indicate that this method is an efficient method to solve th-order nonlinear Volterra integro-differential equations.

1. Introduction

In this paper, we consider the following th-order nonlinear weakly singular Volterra integro-differential equation of the following form [14]: where , is a given function, is the unknown function, and is the kernel of the integro equation. We usually assume that the function and are continuous or square integrable on .

Some problems of mathematical physics are described in terms of (1.1) which has been studied by different methods including the spline collocation method [5], piecewise polynomials [6], Haar wavelets [7], the homotopy perturbation method (HPM) [8, 9], the wavelet-Galerkin method [10], Taylor polynomials [11], the Tau method [12], the sinc-collocation method [13], the combined Laplace transform-Adomian decomposition method [14], and the Adomian’s asymptotic decomposition method [15] to determine exact and approximate solutions. But to our knowledge there is still no viable analytic approach for solving weakly singular Volterra integro-differential equations. The present work is motivated by the desire to obtain approximate solution to th-order nonlinear weakly singular Volterra integro-differential equation, where the integrand is weakly singular in the sense that its integral is continuous at the singular point, that is, its kernel is singular as .

Reproducing kernel theory has important application in numerical analysis, differential equation, probability and statistics and so on [16, 17]. And the RKM has been applied successfully to solving linear and nonlinear problems [1820].

The rest of the paper is organized as follows. In the Section 2, transforming (1.1) into an differential equation by Taylor’s approximation. In Section 3, the RKM is introduced. Applying RKM and ADM to solving (1.1) is discussed in Section 4. The numerical examples are presented in Section 5. Finally, a brief conclusion is stated in last section.

2. Taylor’s Approximation

Consider the following th-order nonlinear weakly singular Volterra integro-different equation: We have by setting Rewriting (2.3) as where the solution under the integral has been replaced by . Thus, So that or equivalently We use the following Taylor’s approximation of degree of about : Thus, Substituting the approximate relation (2.9) into the right hand side of (2.7) yields Therefore, (2.1) can be approximated by the th-order nonlinear differential equation (2.10).

3. Analysis of Reproducing Kernel Hilbert Space

Definition 3.1 (reproducing kernel Space , see [17]). and endowed it with the inner product and norm, respectively,

Theorem 3.2. The space is a reproducing kernel space. That is, there exists a function , for each fixed , and for any , satisfying the reproducing kernel can be denoted by where are known coefficients.

Theorem 3.3. Let be a reproducing kernel space and . If converges to in the sense of , then converges to uniformly.

Property 1. If is a reproducing kernel space, the reproducing kernel function in is unique.

Definition 3.4 (reproducing kernel Space , see [17]). and endowed it with the inner product and norm, respectively,

There exists a unique reproducing kernel function , and can be denoted by

The method of obtaining , the coefficients of the reproducing kernel , and the proof of Theorems 3.2 and 3.3 are given in [17].

4.1. Representation of the Inverse Operator

Here, we propose a new differential operator, as follows: letting where , then we convert (2.10) as follows: We now introduce how to determine the inverse operator of . Obviously, is a bounded linear operator.

We choose as any dense set in , and let , where is the conjugate operator of and is given by (3.7). Furthermore, for simplicity, let denote , namely,

Now, several Lemmas are given.

Lemma 4.1. is the complete system of .

Proof. For , let , that is,
Note that is the dense set in , therefore . It follows that from the existence of .

Lemma 4.2. The following formula holds where the subscript of operator indicates that the operator applies to functions of .

Proof. Consider This completes the proof.

The orthonormal system of can be derived from Gram-Schmidt orthogonalization process of , where are orthogonal coefficients.

Theorem 4.3. If the inverse operator , exists and is dense on , then the inverse operator can be determined as where .

Proof. From (4.9), it holds that The proof of the theorem is complete.

From Theorem 4.3, obviously, is determined.

4.2. Decomposition Method

By applying to both sides of (4.3), we have The ADM introduces the solution and the nonlinear function by infinite series where are Adomian polynomials for the nonlinear term and can be found from the following formula: Substituting (4.12) and (4.13) into (4.11) yields According to the ADM, the components can be determined as which gives where , .

From (4.17), we can determine the components , and hence the series solution of in (4.12) can be immediately obtained.

We obtain approximate solution of the following equation:

The following use some examples to demonstrate the effectiveness of the algorithm.

5. Numerical Examples

To illustrate the applicability and effectiveness of our method, we consider the following examples. Symbolic and numerical computations performed by using Mathematica 5.0.

Example 5.1. Consider the following first-order nonlinear weakly singular Volterra integro-differential equation: with , , and .

Let , , .

On select 100 points and get the approximate solution , the results are shown in Figure 1.

Figure 1: The approximate values in comparison with its exact values .

Example 5.2. Consider the following second-order nonlinear weakly singular Volterra integro-differential equation: with , , and .

Let , , , .

On select 100 points and get the approximate solution , the results are shown in Table 1.

Table 1: Comparison of the errors of and .

Example 5.3. Consider the following third-order nonlinear weakly singular Volterra integro-differential equation: with , , , and .

Let , , , , .

On select 100 points and get the approximate solution , the results are shown in Figure 2.

Figure 2: The approximate values in comparison with its exact values .

6. Conclusion

In this paper, we have reduced the solution of nonlinear weakly singular Volterra integro-differential equations to the solution of ordinary nonlinear differential equations by removing the singularity using an appropriate Taylor’s approximation. Then we have demonstrated the solution of these ordinary nonlinear differential equations by RKM and ADM. The ADM is an accurate and efficient method to solve nonlinear weakly singular Volterra integro-differential equations.

Acknowledgments

Research supported by the National Natural Science Foundation of China (61071181), the Educational Department Scientific Technology Program of Heilongjiang Province(12512133 and 12521148) and the Scientific Innovation Project for Graduate of Heilongjiang Province (YJSCX2012-184HLJ).

References

1. J. Morchało, “On two-point boundary value problem for an integro-differential equation of second order,” Polytechnica Posnaniensis, no. 9, pp. 51–66, 1975.
2. R. P. Agarwal, “Boundary value problems for higher order integro-differential equations,” Nonlinear Analysis, vol. 7, no. 3, pp. 259–270, 1983.
3. A.-M. Wazwaz, “A reliable algorithm for solving boundary value problems for higher-order integro-differentiable equations,” Applied Mathematics and Computation, vol. 118, no. 2-3, pp. 327–342, 2001.
4. E. Babolian and A. S. Shamloo, “Numerical solution of Volterra integral and integro-differential equations of convolution type by using operational matrices of piecewise constant orthogonal functions,” Journal of Computational and Applied Mathematics, vol. 214, no. 2, pp. 495–508, 2008.
5. H. Brunner, “On the numerical solution of nonlinear Volterra-Fredholm integral equations by collocation methods,” SIAM Journal on Numerical Analysis, vol. 27, no. 4, pp. 987–1000, 1990.
6. D. Conte and I. D. Prete, “Fast collocation methods for Volterra integral equations of convolution type,” Journal of Computational and Applied Mathematics, vol. 196, no. 2, pp. 652–663, 2006.
7. O. Lepik, “Haar wavelet method for nonlinear integro-differential equations,” Applied Mathematics and Computation, vol. 176, no. 1, pp. 324–333, 2006.
8. M. Ghasemi, M. Tavassoli Kajani, and E. Babolian, “Application of He's homotopy perturbation method to nonlinear integro-differential equations,” Applied Mathematics and Computation, vol. 188, no. 1, pp. 538–548, 2007.
9. J. Saberi-Nadjafi and A. Ghorbani, “He's homotopy perturbation method: an effective tool for solving nonlinear integral and integro-differential equations,” Computers & Mathematics with Applications, vol. 58, no. 11-12, pp. 2379–2390, 2009.
10. Y. Mahmoudi, “Wavelet Galerkin method for numerical solution of nonlinear integral equation,” Applied Mathematics and Computation, vol. 167, no. 2, pp. 1119–1129, 2005.
11. K. Maleknejad and Y. Mahmoudi, “Taylor polynomial solution of high-order nonlinear Volterra-Fredholm integro-differential equations,” Applied Mathematics and Computation, vol. 145, no. 2-3, pp. 641–653, 2003.
12. G. Ebadi, M. Y. Rahimi-Ardabili, and S. Shahmorad, “Numerical solution of the nonlinear Volterra integro-differential equations by the tau method,” Applied Mathematics and Computation, vol. 188, no. 2, pp. 1580–1586, 2007.
13. M. Zarebnia and Z. Nikpour, “Solution of linear Volterra integro-differential equations via Sinc functions,” International Journal of Applied Mathematics and Computer, vol. 2, pp. 1–10, 2010.
14. A.-M. Wazwaz, “The combined Laplace transform-Adomian decomposition method for handling nonlinear Volterra integro-differential equations,” Applied Mathematics and Computation, vol. 216, no. 4, pp. 1304–1309, 2010.
15. L. Bougoffa, R. C. Rach, and A. Mennouni, “An approximate method for solving a class of weakly-singular Volterra integro-differential equations,” Applied Mathematics and Computation, vol. 217, no. 22, pp. 8907–8913, 2011.
16. H. Yao, The reserch of algorithms for some singular differential equations of higher even-order [Ph.D. thesis], Department of Mathematics, Harbin Institute of Technology, 2008.
17. M. Gui and Y. Lin, Nonlinear Numercial Analysis in the Reproducing Kernel Space, Nova Science Publisher, New York, NY, USA, 2008.
18. 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.
19. 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.
20. Z. Chen and Y. Lin, “The exact solution of a linear integral equation with weakly singular kernel,” Journal of Mathematical Analysis and Applications, vol. 344, no. 2, pp. 726–734, 2008.