Mathematical Problems in Engineering
Volume 2017 (2017), Article ID 7498136, 10 pages
https://doi.org/10.1155/2017/7498136

Research Article

Application of Reproducing Kernel Hilbert Space Method for Solving a Class of Nonlinear Integral Equations

Sedigheh Farzaneh Javan,¹ Saeid Abbasbandy,¹ and M. Ali Fariborzi Araghi²

¹Department of Mathematics, Science and Research Branch, Islamic Azad University, Tehran 14515, Iran
²Department of Mathematics, Central Tehran Branch, Islamic Azad University, Tehran, Iran

Correspondence should be addressed to M. Ali Fariborzi Araghi; fariborzi.araghi@gmail.com

Received 11 November 2016; Revised 7 February 2017; Accepted 8 February 2017; Published 16 March 2017

Academic Editor: Haipeng Peng

Copyright © 2017 Sedigheh Farzaneh Javan 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.

A new approach based on the Reproducing Kernel Hilbert Space Method is proposed to approximate the solution of the second-kind nonlinear integral equations. In this case, the Gram-Schmidt process is substituted by another process so that a satisfactory result is obtained. In this method, the solution is expressed in the form of a series. Furthermore, the convergence of the proposed technique is proved. In order to illustrate the effectiveness and efficiency of the method, four sample integral equations arising in electromagnetics are solved via the given algorithm.

1. Introduction

Electromagnetics is the phenomenon associated with electric and magnetic fields and their interactions which is generally one of the most important sciences. Exterior calculus is given in [1, 2] inside some textbooks. A way to teach electromagnetics can be approached via the use of differential forms which is given in [3]. According to electromagnetic field problems from many years ago, some solutions via linear and nonlinear integral equations (NIE) have been given which can be useful in the field. In those methods like block-pulse functions (BPFs), Galerkin, and collocation, the most important ways are basic functions and appropriate projection. Based on the Reproducing Kernel Hilbert Space method, an approach has been found to solve some electromagnetic issues.

Nonlinear integral equations are encountered in different fields of science and numerous applications as elasticity, plasticity, heat and mass transfer, oscillation theory, fluid dynamics, filtration theory, electrostatics, electrodynamics, biomechanics, game theory, control, queuing theory, electrical engineering, economics, and medicine, among others. There are different types of NIE usually which cannot be worked out explicitly, so it should be approached approximately.

Therefore, many researchers studied and focused on different numerical techniques which can work out these integral equations. For instance, in [4, 5], the authors presented the homotopy analysis method to solve the second kind of nonlinear Fredholm and Volterra integral equations. The linear multistep techniques were applied in [6], to obtain the numerical solution of a singular nonlinear Volterra integral equation. In [7], an asymptotic technique to approach numerically the nonlinear Abel-Volterra integral equation was applied.

Reproducing Kernel Hilbert Space (RKHS) was introduced by Minggen et al. [8, 9], and it was developed in different areas, including approximation theory, statistics, machine learning theory, group representation theory, and various areas of complex analysis. Reproducing Kernel Hilbert Space Method (RKHSM) is a kernel based approximation method which was applied for solving nonlinear boundary value problems [7–12], generalized singular nonlinear Lane-Emden type equations [13], integrodifferential equations [14–16], integrodifferential fractional equations [17], Bratus Problem [18], and so forth.

Consider the following nonlinear integral equation:

where , are real constants, is an unknown function which can be determined, is a continuous function on , is a continuous function on , is a continuous term in as , , and is Reproducing Kernel Space. Equation (1) has a continuous solution on [19]. The existence and uniqueness conditions of the solution for (1) were discussed in [19–23]. We assume that the solution of (1) is unique.

Over several decades, numerical methods in electromagnetic problems have been one of the most important subjects of extensive researches [1–4]. On the other hand, many problems in electromagnetics can be modeled by integral equations mentioned in [24–26], for example, electric field integral equation (EFIE) and magnetic field integral equation (MFIE). In recent years, several numerical methods for solving linear and nonlinear integral equations have been presented. Applicable equations of electromagnetics have been implied in the presented paper.

In previous works like [13–15], the Gram-Schmidt orthogonalization process has been considered to implement RKHSM. Since this process is unstable numerically and it may take a lot of time to run the algorithm, here, we put away this process and act with another way. Our approach combines the methods mentioned in [13–17]. More specifically, on the contrary to [13–15], without use of the orthogonalization process, the RKSHM is applied successfully to solve the nonlinear problem (1).

The structure of this paper would be described as follows. In Section 2, the basic definitions, assumptions, and preliminaries of RKHS are described. The main idea and convergence of the proposed scheme are discussed in Section 3. Section 4 contains the numerical experiments. Finally, Section 5 is dedicated to a brief conclusion.

2. Preliminaries

In this section, some basic definitions and important properties of Reproducing Kernel Hilbert Spaces (RKHS) are mentioned [8, 9, 27–29].

Definition 1. A Hilbert Space is an inner product space that is complete and separable with respect to the norm defined by the inner product. Completeness of the space holds provided that every Cauchy sequence of points in that has a limit that is also in and separable of admits a countable orthonormal basis of it.

Definition 2. For an abstract set , let be a Hilbert Space of real or complex-valued functions on set . We say is a Reproducing Kernel Hilbert Space if there exist a linear and bounded evaluation functional over , or, equivalently,

Riesz Representation Theorem implies that for all in there exists a unique function of with the reproducing property,

for each where represent the inner product of the Hilbert Space .

Definition 3. The space is interpreted as

The inner product and the norm in are of forms

Lemma 4 (see [9, 29]). Functional space is inner space.

Theorem 5 (see [9, 29]). Functional space is a Hilbert Space.

Theorem 6 (see [9, 29]). Functional space is Reproducing Kernel Hilbert Space.

Now, it is taken away that expression form of the Reproducing Kernel function .

Based on essay, it is easy to prove that is the answer of the following generalized differential equation [9, 29]:

where is Dirac’s delta function. While , is the answer of the following constant linear homogeneous differential equation with order:

with the boundary conditions

Equation (7) is characteristic . Then the general solution of Equation (7) is

where coefficients and , , could be calculated by solving the following linear equations:

Subsequently, the representation of the Reproducing Kernel of is provided by

3. Main Idea and Theoretical Discussion

The uniqueness conditions for nonlinear problems exist in [21–23]. The unique solution of (1) is assumed in this paper. The solution of (1) is given in space. We consider (1) as

where . It is obvious that is the bounded linear operator of to . Put and , where is the adjoint operator of . In fact, for , we have and .

The orthonormal system of from the space can be derived from Gram-Schmidt orthogonal process of :

Definition 7. In a topological space , a subset of is called dense in if .

Theorem 8. If is dense on then is the complete function system of the space and , where the subscript t in the operator indicates that the operator applies to the function of .

Proof. We have

Clearly . For each fixed , let , , which means

Assume that is dense on and so . It follows that from the existence of . Now, the theorem is proved.

Theorem 9. If is dense on and the solution of (12) is unique, then the solution of (12) is

Proof. Using (13), we have

On the other hand, and , , are the Fourier series expansion about normal orthogonal system and is the Hilbert Space. Thus the series is convergent in the sense of and the proof would be complete.

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

Theorem 10. If then exists such that .

Proof. We have , for any . We know . Thus .

Theorem 11. The approximate solution is uniformly convergent.

Proof. Assuming , by Theorems 9 and 10, it can be proved that

In the sequel, a new iterative method to achieve the solution of (12) is presented. If

then (16) can be written as

Now suppose, for some , is known. There is no problem if we assume . We put and define the -term approximation to by

where

In the following, it would be proven that the approximate solution in the iterative (22) is convergent to the exact solution of (12) uniformly.

Theorem 12. Suppose that is bounded in (22). If is dense on then -term approximate solution in the iterative (22) converges to the exact solution of (12) and , whereas is given by (23).

Proof. First of all, the convergence of from (22) would be proven. We infer

Subsequence is orthogonal, and it yields that

It is obvious that the sequence is monotonically increasing. Because is bounded and is convergent, then is bounded and this implies that .
If then

So . Consequently as . To prove the completeness of it requires , where that as . Now we can prove is the solution of (12).
If we take limit from (22), we will have , so . Let , and then

From (23) and (29), it is concluded that .
is dense on . For each , subsequence exists that as . Hence, when , we have which indicates that is the solution of (12).

The mentioned scheme above is an efficient method of solving nonlinear equations [31–33]. However, in implementing this algorithm on a computer, is not quite orthogonal, due to rounding errors. In other words, Gram-Schmidt process is numerically unstable and the computational cost of the algorithm is high. Therefore, the following process is suggested similar for linear problems in [20, 34]. This is the subject of the next theorem, where the following notations are used:

Theorem 13 (). The approximate solution obtained from (22) is found as follows:

where

Proof. Suppose that . Since ,

and () are linear independence, and therefore

Equation (12) implies . For we have

Both sides of (33) provide

From (32) and (35) the following equation can be reached:

Equation (32) implies . So

which proves the theorem.

Algorithm 14. The following steps exist for approximating the solution without applying Gram-Schmidt orthogonal process:

 Step  1. Fix and . If , set . Else set . Step  2. For set . Set . Step  3. Set . Step  4. For set . Step  5. . Step  6. Set . Step  7. Set . Step  8. If then set and go to step 6. Else stop.

Algorithm 15. The following steps exist for approximating the solution by applying Gram-Schmidt orthogonal process:

 Step  1. Fix and . If , set . Else set . Step  2. For set . Set . Step  3. For and , if then set and . Else . Else . Step  4. For set . Step  5. Set . Step  6. Set . Step  7. Set . Step  8. Set . Step  9. If then set and go to step 7. Else stop.

4. Numerical Experiments

In this part, four numerical examples are solved for potency and utility of the present method. All computations are performed by MAPLE package. Results which are taken by this method show a proper agreement with the exact solution. A comprehensive applicability of this method is given the stability and consistence of the presented method. The reliability of the method and increasing of the accuracy cause this method to be more applicable.

Example 1. For first applicable instance, we offer nonlinear Fredholm integral equations [26, 30]:

The exact solution of this equation is . According to (38) we can assume an initial approximation . The numerical results are given in Table 1 by taking , and . In Table 1 a comparison between the absolute errors of our method and the Haar wavelet method [30] is given. Figure 1 shows the approximate solution and its errors.

Table 1: The absolute errors comparison between the proposed approach and method [30].

Exact solution
Approximate solution
Method in [30]
Absolute error
Absolute error in method [30]

0.1
1.105170918
1.105154993
1.114627560
0.000015925
0.009456642

0.2
1.221402757
1.221389903
1.226758840
0.000012855
0.005356083

0.3
1.349858806
1.349845490
1.345458179
0.000013317
0.004400627

0.4
1.491824696
1.491803913
1.480202470
0.000020784
0.011622226

0.5
1.648721268
1.648680858
1.671769819
0.000040412
0.023048551

0.6
1.822118797
1.822039541
1.838854903
0.000079259
0.016736106

0.7
2.013752703
2.013608464
2.022118086
0.000144243
0.008365383

0.8
2.225540923
2.225304288
2.223139077
0.000236640
0.002401846

0.9
2.459603104
2.459253211
2.443684898
0.000349899
0.015918206

Figure 1: The absolute errors comparison between the proposed approach and method [30].

Example 2. For second applicable example, an electromagnetic problem is solved via the presented method. It is simulated to nonlinear Volterra integral equations model [26, 30]:

The exact solution of this equation is . According to (39) we can assume an initial approximation . Numerical results are given in Table 2 by taking , and . In Table 2 a comparison between the absolute errors of the proposed method and the BPFs method [26] is given. Figure 2 shows the approximate solution and its errors. In Table 3, a comparison execution time between Algorithms 14 and 15 is given.

Table 2: The absolute errors comparison between the proposed approach and method [26].

Exact solution
Approximate solution
Method in [26]
Absolute error
Absolute error in method [26]

0.1
0.099833
0.098524
0.101388
0.001308
0.00155

0.2
0.198669
0.199766
0.194073
0.001097
0.004596

0.3
0.295520
0.294857
0.299995
0.000662
0.004475

0.4
0.389418
0.389658
0.387978
0.000240
0.001440

0.5
0.479426
0.479585
0.486266
0.000160
0.006840

0.6
0.564642
0.563993
0.565930
0.000648
0.001288

0.7
0.644218
0.645477
0.640624
0.001260
0.003594

0.8
0.717356
0.715373
0.720611
0.001982
0.003255

0.9
0.783327
0.785653
0.782351
0.002326
0.000976

Table 3: The execution time (seconds) comparison between Algorithms 14 and 15.

Algorithm 14 (Sec)
Algorithm 15 (Sec)

4
3.47
2.87

8
9.52
4.82

16
59.75
31.77

Figure 2: The absolute errors comparison between the proposed approach and method [26].

Example 3. An electromagnetic problem is solved via our method for another applicable example. It is simulated to nonlinear Fredholm integral equations model [26, 30]:

The exact solution of this equation is . According to (40), the initial approximation is chosen. Numerical results are given in Table 4 by taking , and . In Table 4 a comparison between the absolute errors of the proposed method and the BPFs method [26] is given. Figure 3 shows the approximate solution and its error. In Table 5 a comparison execution time between Algorithms 14 and 15 is given.

Table 4: The absolute errors comparison between the proposed approach and method [26].

Exact solution
Approximate solution
Method in [26]
Absolute error
Absolute error in method [26]

0.1
0.010000
0.010762
0.010308
0.000762
0.000308

0.2
0.040000
0.040655
0.038140
0.000655
0.001860

0.3
0.090000
0.090444
0.092828
0.000444
0.002828

0.4
0.160000
0.16023
0.158746
0.000230
0.001254

0.5
0.250000
0.250001
0.257867
0.000001
0.007867

0.6
0.360000
0.359606
0.361871
0.000394
0.001871

0.7
0.490000
0.489448
0.483453
0.000552
0.006547

0.8
0.640000
0.639391
0.647515
0.000609
0.007515

0.9
0.810000
0.811638
0.807183
0.001638
0.002817

Table 5: The execution time (seconds) comparison between Algorithms 14 and 15.

Algorithm 14 (Sec)
Algorithm 15 (Sec)

4
3.1
2.59

8
8.46
5.09

16
34.2
28.89

Figure 3: The absolute errors comparison between the proposed approach and method [26].

Example 4. A nonlinear Fredholm integral problem [26, 30] is solved via our method for this applicable example:

The exact solution of this equation is . According to (41) we consider an initial approximation as . Numerical results are given in Table 6 by taking , and . In Table 6 a comparison between the absolute errors of the proposed method and the Haar wavelet method [30] is given. Figure 4 shows the approximate solution and its error.

Table 6: The absolute errors comparison between the proposed approach and method [30].

Exact solution
Approximate solution
Method in [30]
Absolute error
Absolute error in method [30]

0.1
1.105170918
1.105495796
1.111099473
0.000324878
0.005928555

0.2
1.221402757
1.220931139
1.220296615
0.000471618
0.001106142

0.3
1.349858806
1.349502887
1.339817084
0.000355919
0.010041722

0.4
1.491824696
1.491740269
1.471965091
0.000084427
0.019859605

0.5
1.648721268
1.648649846
1.457127691
0.000071422
0.191593577

0.6
1.822118797
1.822243106
1.831888927
0.000124309
0.009770130

0.7
2.013752703
2.013650232
2.011936556
0.000102471
0.001816147

0.8
2.225540923
2.224684503
2.209678381
0.000856420
0.015862542

0.9
2.459603104
2.458626573
2.426854676
0.000976531
0.032748428

Figure 4: The absolute errors comparison between the proposed approach and method [30].

5. Conclusion

According to this essay, supplementary of iterative Reproducing Kernel Hilbert Space Method was introduced and applied to acquire the approximate solution of some nonlinear integral equation. In this method, unlike other similar methods, orthogonal process is not used. However the time is increasing, and accuracy is also increasing. The main point which is mentioned in this paper is that Algorithm 14 has higher execution time in comparison with Algorithm 15, but the approximate solution in Algorithm 14 is more accurate than Algorithm 15. Current uniform convergence method is stated and proved. The obtained numerical results confirm that it is a good candidate for solving the nonlinear integral equation.

Conflicts of Interest

The authors declare no conflicts of interest.

Acknowledgments

The authors would like to thank Professor Haipeng Peng for instructive comments and recommendations to improve the quality of this work and also the Islamic Azad University, Science and Research Branch, Tehran, for supporting this project.

References

K. Meetz and W. Engi, Electromagnetische Felder, Springer, Berlin, Germany, 1979.

W. Thirring, Lehrbuch der Mathematischen Physik, vol. 2, Springer, Vienna, Austria, 1978.

K. F. Warnick, R. H. Selfridge, and D. V. Arnold, “Teaching electromagnetic field theory using differential forms,” IEEE Transactions on Education, vol. 40, no. 1, pp. 53–68, 1997.

S. Abbasbandy, “Homotopy analysis method for heat radiation equations,” International Communications in Heat and Mass Transfer, vol. 34, no. 3, pp. 380–387, 2007.

H. Vosughi, E. Shivanian, and S. Abbasbandy, “A new analytical technique to solve Volterra's integral equations,” Mathematical Methods in the Applied Sciences, vol. 34, no. 10, pp. 1243–1253, 2011.

J.-P. Kauthen, “A survey of singularly perturbed Volterra equations,” Applied Numerical Mathematics. An IMACS Journal, vol. 24, no. 2-3, pp. 95–114, 1997.

A. A. Kilbas and M. Saigo, “On solution of nonlinear Abel-Volterra integral equation,” Journal of Mathematical Analysis and Applications, vol. 229, no. 1, pp. 41–60, 1999.

C. Minggen and D. Zhongxing, “On the best operator of interpolation,” Mathematica Numerica Sinica, vol. 8, no. 2, pp. 209–216, 1986.

C. Minggen and L. Yingzhen, Nonlinear Numerical Analysis in the Reproducing Kernel Space, Nova Science, New York, NY, USA, 2009.

S. Abbasbandy, B. Azarnavid, and M. S. Alhuthali, “A shooting reproducing kernel Hilbert space method for multiple solutions of nonlinear boundary value problems,” Journal of Computational and Applied Mathematics, vol. 279, pp. 293–305, 2015.

N. Shawagfeh, O. Abu Arqub, and S. Momani, “Analytical solution of nonlinear second-order periodic boundary value problem using reproducing kernel method,” Journal of Computational Analysis and Applications, vol. 16, no. 4, pp. 750–762, 2014.

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.

B. Azarnavid, F. Parvaneh, and S. Abbasbandy, “Picard-reproducing kernel Hilbert space method for solving generalized singular nonlinear Lane-Emden type equations,” Mathematical Modelling and Analysis, vol. 20, no. 6, pp. 754–767, 2015.

O. A. Arqub, M. Al-Smadi, and N. Shawagfeh, “Solving Fredholm integro-differential equations using reproducing kernel Hilbert space method,” Applied Mathematics and Computation, vol. 219, no. 17, pp. 8938–8948, 2013.

L. Yang and M. Cui, “New algorithm for a class of nonlinear integro-differential equations in the reproducing kernel space,” Applied Mathematics and Computation, vol. 174, no. 2, pp. 942–960, 2006.

O. Abu Arqub, M. Al-Smadi, and S. Momani, “Application of reproducing kernel method for solving nonlinear Fredholm-Volterra integrodifferential equations,” Abstract and Applied Analysis, vol. 2012, Article ID 839836, 16 pages, 2012.

S. Bushnaq, B. Maayah, S. Momani, and A. Alsaedi, “A reproducing kernel Hilbert space method for solving systems of fractional integrodifferential equations,” Abstract and Applied Analysis, vol. 2014, Article ID 103016, 6 pages, 2014.

M. Inc, A. Akgül, and F. Geng, “Reproducing kernel Hilbert space method for solving Bratu's problem,” Bulletin of the Malaysian Mathematical Sciences Society, vol. 38, no. 1, pp. 271–287, 2015.

A.-M. Wazwaz, A First Course in Integral Equations, World Scientific, 1997.

S. Javadi, E. Babolian, and E. Moradi, “New implementation of reproducing kernel Hilbert space method for solving a class of functional integral equations,” Communications in Numerical Analysis, vol. 2014, Article ID cna-00205, 7 pages, 2014.

W. Jiang and Z. Chen, “Solving a system of linear Volterra integral equations using the new reproducing kernel method,” Applied Mathematics and Computation, vol. 219, no. 20, pp. 10225–10230, 2013.

A.-M. Wazwaz, Linear and Nonlinear Integral Equations, Methods and Applications, Springer, Berlin, Germany, 2011.

I. L. El-Kalla, “Convergence of the Adomian method applied to a class of nonlinear integral equations,” Applied Mathematics Letters. An International Journal of Rapid Publication, vol. 21, no. 4, pp. 372–376, 2008.

G. Gumah, K. Moaddy, M. Al-Smadi, and I. Hashim, “Solutions to uncertain Volterra integral equations by fitted reproducing kernel Hilbert space method,” Journal of Function Spaces, vol. 2016, Article ID 2920463, 11 pages, 2016.

V. Sizikov and D. Sidorov, “Generalized quadrature for solving singular integral equations of Abel type in application to infrared tomography,” Applied Numerical Mathematics. An IMACS Journal, vol. 106, pp. 69–78, 2016.

E. Babolian, Z. Masouri, and S. Hatamzadeh-Varmazyar, “New direct method to solve nonlinear Volterra-Fredholm integral and integro-differential equations using operational matrix with block-pulse functions,” in Progress in Electromagnetics Research B, vol. 8, pp. 59–76, EMW Publishing, Cambridge, Mass, USA, 2008.

D. Alpay, Ed., Reproducing kernel spaces and applications, vol. 143 of Operator Theory: Advances and Applications, Birkhäuser Basel, 2003.

S. Saitoh, D. Alpay, J. A. Ball, and T. Ohsawa, Reproducing Kernels and Their Applications, vol. 11, Springer Science & Business Media, Berlin, Germany, 2013.

S. Saitoh, Integral Transforms, Reproducing Kernels and Their Applications, vol. 369 of Pitman Research Notes in Mathematics Series, Longman, Harlow, UK, 1997.

E. Babolian and A. Shahsavaran, “Numerical solution of nonlinear Fredholm integral equations of the second kind using Haar wavelets,” Journal of Computational and Applied Mathematics, vol. 225, no. 1, pp. 87–95, 2009.

R. Ketabchi, R. Mokhtari, and E. Babolian, “Some error estimates for solving Volterra integral equations by using the reproducing kernel method,” Journal of Computational and Applied Mathematics, vol. 273, pp. 245–250, 2015.

A. Alvandi, M. Paripour, and Z. Roshani, “Reproducing kernel method for solving a class of Fredholm integro-differential equations,” in Proceedings of the 46th Annual Iranian Mathematics Conference (AIMC 46), p. 505, Yazd University, 2015.

I. Komashynska and M. Al-Smadi, “Iterative reproducing kernel method for solving second-order integrodifferential equations of fredholm type,” Journal of Applied Mathematics, vol. 2014, Article ID 459509, 11 pages, 2014.

E. Babolian, S. Javadi, and E. Moradi, “Error analysis of reproducing kernel Hilbert space method for solving functional integral equations,” Journal of Computational and Applied Mathematics, vol. 300, pp. 300–311, 2016.