Iterative Solution and Convergence of Nonlinear Volterra Integral Equations of The Second Kind using HPM

Download Full-Text PDF Cite this Publication

Text Only Version


Iterative Solution and Convergence of Nonlinear Volterra Integral Equations of The Second Kind using HPM

Charles Ekene Chika Department of Mathematical Sciences, University of Nigeria, Nsukka, Nigeria.


Chinedu Junior Okere

China University of Petroleum, Beijing, China.

Abstract:- An iterative method is developed to solve a class of nonlinear Volterra integral equations. This method uses the concept of homotopy perturbation to approximate the exact solution of the integral equation. The convergence is discussed and illustrated with examples. Examples are used to show the validity of the algorithm and theorems presented which shows that the method is effective and easy to implement.

Keywords: Homotopy perturbation method; Nonlinear Volterra integral; Convergence, Numerical Method

Subject Class: Numerical Method


    Solution to integral equations is of high interest in science since these equations appear in many applied problems; which in many cases takes the form of the so called Fredholm or Volterra Integral Equations. This equations result from direct model or from solving another problem like the one encountered by Chika and Hooshyar [5]. Di_erent methods have been developed to solve these integral equations, some which are presented in [1]. These methods have some assumptions like requiring separable kernel, linearity etc. hence, there is need to develop more methods that will solve di_erent classes of these equations.

    Homotopy Perturbation Method (HPM) has been good in solving nonlinear equations ever since it was introduced by He [7], and this has been applied to solve linear integral equations like the one by Abbasbandy [8]. The convergence of this method has been studied as well like the ones presented in [4] and [6]. J. Biazar and H. Ghazvini [2] used this method to solve nonlinear Fredholm integral which worked well. The aim of this work is to apply this method in nonlinear Volterra integral and study it’s convergence.

    Numerical (or iterative) algorithm developed using homotopy perturbation as well as examples showing the validity of the algorithm is presented in section 2. Section 3 discusses the convergence and section 4 has the conclusion and future works.

  2. Numerical Algorithm and Examples

    Let’s consider the nonlinear Volterra Integral Equation:

    x [a; b] and m is a positive integer. This type of equations occurs in applied problems like those arising from solitons, fractals, etc.




    To get the solution to Eq(1), consider the expansion () = 1 () where p [0; 1] such that the solution u(x) to Eq(1)

    is u(x) = lim 1 (x)



    Using the homotopy equation




    Substituting () = 1 () in Eq(2) we have

    Expanding the series in Eq(3) and equating like powers of p gives

    It follows that

    Therefore, the numerical algorithm for solving Eq(1) can be written as, () = ()

    for n > 0


    In practice, we take u(x) ()

    is the so called partial sum.

    To check the feasibility of the algorithm, we look at some examples:




    Example 1: Find the function u(x) such that u(x) = + 1

    (1 3) + 1

    3(), [0, ] is satisfied.


    Using the algorithm for n > 1


    100 0 2

    We use MATLAB to compute the above quantity for di_erent values of x and N = 4. The exact solution to this problem is u(x) = ex. As can be seen in Figure 1 and Table 1, there is good agreement between the exact solution and the numerical solution.

    Figure 1: The graph of exact and approximate solution by HPM in Example 1

    Example 2: Find u(x) that satisfied u(x) = 1 + 1 2 1

    4 1

    6 + 1 ( ) 2() u(x). The exact solution to this

    160 480


    80 0

    equation is u(x) = 1 + x2. Applying the algorithm, we have; (x) = 1 + 1

    2 1

    4 1


    for n > 1






    () = 80 ( ) (()

    0 =0

    Looking at Figure 2 and Table 2, the exact solution and the numerical solution are in good agreement.

    Example 3: Consider the equation () = cos() + 1

    cos(2) 1 2 1

    + 1 ( )2 (), x [0; 1] whose solution

    is u(x) = cos (x).


    80 160

    20 0

    Table 1: Numerical output for example 1

    Table 2: Numerical output for example 2

    Applying the algorithm

    Figure 3 and Table 3 shows the comparison of the exact and the numerical solutions. It can be seen that the numerical approximation performed well.

  3. Convergence and Error Estimation

    In this section, we discuss the convergence of the algorithm presented above. We are going to use Abel’s Theorem [4], i.e if =




    0 has radius of

    Figure 2: The graph of exact and approximate solution by HPM in Example 2

    Figure 3: The graph of exact and approximate solution by HPM in Example 3 Table 3: Numerical output for example 3




    Convergence not less than one and the series 0 is absolutely convergent, then,


    u(x) = lim =



    For this first work, we take m=2 to investigate the convergence.

    Theorem 1: for m = 2, let K (x, t) and f(x) be continuous in the regions = [a; b] × [a; b] and = [a; b] respectively. K (x, t) and

    f(x) is bounded such that | (, )| and |()|

    for all x, t [a, b] and C = NkN0. If || < 1 , then the algorithm


    above is uniformly convergent in [a, b] for each p [0, 1].

    Proof for m=2

    Using Eq(6), for p [0, 1]

    Since m=2

    This is geometric series, for convergence;


    or more strongly

    Theorem 2: The Nth order approximation error of the solution to Eq(1) () =

    (||())+1 for all x [a, b] or more strongly


    (||())+1 1 ||()

    1|| ()


    Note: This estimates the error between partial sum and infinite sum using HPM (not exact solution). Proof

    Using Eq(6) and m = 2

    From Theorem 1, if other quantities are constant; smaller , gives better convergence and from Theorem 2 reduces the error; this is illustrated in the next example.

    Example 4: Using equation in example 3 and changing from 1 to 1

    , the equation becomes

    20 100

    This implies that || < 1


    , therefore, the algorithm converges.

    Comparing values in Table 3 and Table 4 as well as Figure 3 and Figure 4, HPM gives a better approximation in example 4 than example 3, which confirms the validity of the Theorems.

    Table 4: Numerical output for example 4

  4. Conclusion

The algorithm developed using homotopy perturbation method performs well as shown in the examples, the percentage error in the shown examples are all less than 0.65% except in example 3 which has up to 2.65%. This is explained by the convergence theorems. Since 4th partial sum is used in all the computations, better approximation can be gotten by adding more terms. This algorithm is easy to implement and gives good approximation. Given the impressive performance, future communication will involve more analysis and extending the method to more general class of nonlinear Volterra integral.

Figure 4: The graph of exact and approximate solution by HPM in Example 4


  1. Abdul-Majid Wazwaz , ” Linear and Nonlinear Integral Equations,” Higher Education Press; Springer, Beijing (2011).
  2. J. Biazar, H. Ghazvini , “Numerical solution for special non-linear Fredholm integral equation by HPM,” Appl. Math. Comput. 195 (2008) 681-687.
  3. S. Abbasbandy , “Numerical solutions of the integral equations: Homotopy perturbation method and Adomian’s decomposition method ,” Appl. Math. Comput. 173 (2006) 493-500.
  4. Edyta Hetmaniok, Damian Slota, Roman Witula , “Convergence and error estimation of homotopy perturbation method for Fredholm and Volterra integral equations,” Appl. Math. Comput. 218 (2012) 10717-10725.
  5. Chika CE, Hooshyar MA. , “An electromagnetic inverse scattering problem for dielectrics that depend on two spatial variables via Eikonal approximation,” Microw Opt Technol Lett., 2019;1{7.
  6. Behzad Ghanbari, “On the Convergence of the Homotopy Analysis Method for Solving Fredholm Integral Equations,” Walailak J Sci & Tech, 2013; 10(4):395{403.
  7. J.H. He, “Homotopy perturbation Technique,” Computer Methods in Applied Mechanics and Engineering, 178 (1999) 257-262.
  8. S. Abbasbandy, “Numerical solutions of the integral equations: Homotopy perturbation method and Adomain’s decomposition method,” Applied Mathematics and Computation, 173 (2006) 493-500.


Leave a Reply

Your email address will not be published. Required fields are marked *