The Scientific World Journal

Volume 2014 (2014), Article ID 465951, 7 pages

http://dx.doi.org/10.1155/2014/465951

Research Article

The Convergence Study of the Homotopy Analysis Method for Solving Nonlinear Volterra-Fredholm Integrodifferential Equations

Behzad Ghanbari

Department of Basic Sciences, Kermanshah University of Technology, P.O. Box 63766-67178, Kermanshah, Iran

Received 13 August 2013; Accepted 19 October 2013; Published 12 January 2014

Academic Editors: Y.-M. Chu, Z. Mukandavire, and K. Zhu

Copyright © 2014 Behzad Ghanbari. 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 aim to study the convergence of the homotopy analysis method (HAM in short) for solving special nonlinear Volterra-Fredholm integrodifferential equations. The sufficient condition for the convergence of the method is briefly addressed. Some illustrative examples are also presented to demonstrate the validity and applicability of the technique. Comparison of the obtained results HAM with exact solution shows that the method is reliable and capable of providing analytic treatment for solving such equations.

1. Introduction

In recent years, there has been a growing interest in the nonlinear Volterra–Fredholm integrodifferential equations, NVFID, in short, which are a combination of differential and Fredholm-Volterra integral equations. These equations occur frequently as a model in mathematical biology and the physical sciences. Volterra–Fredholm integrodifferential equations is sometimes difficult to solve analytically. Therefore, finding either the analytical approximation or numerical solution of such equations is of great interest. The interested reader can refer to [1–11] (and the references therein) for more research works.

We aim in this work to study the convergence of the homotopy analysis method, for solving NVFID equations, of the form:

where is a real number, the kernel , are known continuous functions over , and , are given continuous function defined over .

The homotopy analysis method (HAM) [12] has been proved to be one of the useful techniques to solve numerous linear and nonlinear functional equations [13–17].

As mentioned in [13, 14], unlike all previous analytic techniques, the homotopy analysis method provides great freedom to express solutions of a given nonlinear problem by means of different base functions. Also this method provides a way to adjust and control the convergence region and the rate of convergence of solution series, by introducing the auxiliary parameter .

By properly choosing the base functions, initial approximations, auxiliary linear operators, and auxiliary parameter , HAM gives rapidly convergent successive approximations of the exact solution. A systematic description of this analytic technique, for nonlinear problems, can be found in [13].

This paper is organized as follows. In Section 2, a short description of the basic ideas of the homotopy analysis method will be stated. In Section 3, the homotopy analysis method is applied to construct approximate solution of (1). Section 4 is devoted to the convergence analysis of the method. In Section 5, our numerical findings are reported and demonstrate the accuracy of the proposed scheme, by considering three numerical examples. Finally, conclusions are stated in the last section.

2. Basic Idea of HAM

Let us assume the following nonlinear differential:

where is a nonlinear operator, is an independent variable, and is the solution of the equation.

By utilizing the concept of homotopy in topology, Liao [13] constructs the so-called zero-order deformation equation

where is the embedding parameter, is a nonzero auxiliary parameter, is an auxiliary linear operator, is an initial guess of , and is an unknown function, respectively.

In the view of HAM the solution of original equation is assumed to be expanded in terms of embedding parameter as

where ’s are obtained as follows.

Appling recently proposed “th-order homotopy-derivative operator” [14]:

to both sides of (3), one reads

subject to initial condition

where

In this way ’s, for , are dependent upon not only but also the auxiliary parameter .

As proved in [13], whenever (4) is a convergent series at , its limit must satisfy the original equation . To find a proper convergence-control parameter , to get a convergent series solution or to accelerate its convergence rate, there is a classic way of plotting the so-called ‘‘-curves” or ‘‘curves for convergence-control parameter”. Such a region can be found, although approximately, by plotting the curves of these unknown quantities versus .

However, it is a pity that curves for convergence-control parameter (i.e., -curves) give us only a graphical region and cannot tell us which value of gives the fastest convergent series. Moreover, recently in [18] a misinterpreted usage of -curves has been reported.

3. HAM for NVFID

To investigate the exact solution of (1), by means of the homotopy analysis method, let us consider the nonlinear operator

and linear operator

with the following property:

where is constant.

The initial guess is chosen such that it satisfies the initial condition of problem; that is,

from (8) and by using (9), we have

From (13), it should be noted that the right-hand side of (6) is only dependent upon . Thus, , will be computed recursively, by means of solving the linear high-order deformation equation (6) subject to (7). So, the th-order approximation of is given by

4. Convergence Analysis

In this section, some theorems and conditions in the framework of convergence of the homotopy analysis method are stated.

Theorem 1. Whereas the series

converges, where ’s are resulted from (6), (7), and (13), the limit of the series is an exact solution of (1).

Proof. Since, by hypothesis, the series is convergent, it holds

So, the necessary condition for the convergence of the series is valid; that is,

Using (6) and (17), we have

Since , we must have

On the other hand, we have

So, from (19), we obtain

Also, from initial conditions (7) and (12), the following holds:

Since satisfies (21) and (22), we conclude that it is an exact solution of (1). This completes the proof.

Theorem 2. Suppose that is a Banach space with a suitable norm , say , over which the sequence is defined for a prescribed value of . Assume also that the initial approximation remains inside the domain of the solution . Taking as a constant, the following statements holds.

(i)If there exists some , such that for all we have , then the series solution converges absolutely at , over the domain of .(ii)If there exists some , such that for all we have , then the series solution diverges, at , over the domain of definition of .

Proof. Indeed, this is a special case of Banach' fix point theorem, which can be found in standard texts on real analysis such as in [19].

Theorems 1 and 2 state that the homotopy series solution , of nonlinear problem (1), converges to an exact solution , under the condition that such that , for some . In other words, if, for every , we consider the parameters ’s, as were defined in [20],

for , then the series solution of problem (1) converges to an exact solution , when .

5. Numerical Results

To illustrate the convergence study of HAM, two examples are presented.

Example 1. Let us first consider the following linear Fredholm integral equation:

where , with the initial condition . The exact solution is [6].

To solve (24) by means of HAM and according to the initial condition, we can choose the initial guess of as follows:

The linear operator and the zero-order deformation equations are as defined by (10) and (3), respectively.

Therefore, the th-order deformation equations (6) read as

subject to initial condition

where

Now, we can successively obtain the solution of the th-order deformation equation (26) and (27), for .

Table 1 shows some values of ’s, defined as in (23), which were obtained from the HAM series solution (14) by using different values of . Form Table 1, since for , , and , we can conclude that the HAM approach converges to the exact solution for , , and . Also, one can observe that ’s are not less than one for . So, the HAM approach may be divergent for .

In Table 2 the relative errors of the terms approximation of HAM, which is defined as

for different values of , are presented. It is evident that the auxiliary parameter can also be effectively implemented to adjust and control the rate of convergence of the series solutions by HAM.

Table 1: Numerical values of ’s for different values of in Example 1.

0.675127819
0.558404720
0.371647762
1.53698244

0.250131472
0.119092776
0.360689399
1.540085431

0.466457456
1.052930969
0.097317362
1.543843327

0.971262707
0.366763445
0.669779432
1.548761762

0.378085555
0.182601166
0.068939396
1.555540384

0.305191320
1.321695721
0.700358533
1.565066508

Table 2: Comparison of relative errors for Example 1.

93.053

45.719

29.642

21.377

16.237

12.660

9.975

7.846

6.080

4.555

Example 2. Consider now the following equation [7, 11]:

with initial condition . The exact solution is .

To obtain the approximate solution of this example, if we choose , then the solutions of th-order deformation equations can be resulted in (6) and (7), where

Since HAM provides a family of solution expression in the auxiliary parameter , the convergence region of solution series depends upon the value of .

Table 3 shows the computed values of ’s for different values of , in this example. Also, for different values of , the relative errors of approximate analytic solutions by HAM solution are reported in Table 4.

Clearly, one can observe that the approximate solution for , is more accurate than the approximate solutions obtained for , , and . Therefore, it can be claimed that the auxiliary parameter plays an important role to adjust and control the convergence of the series solution. It can be seen from Tables 3 and 4 that, for smaller values of , our approximate solutions are in very excellent agreement with the exact solutions.

Table 3: Numerical values of ’s for different values of in Example 2.

0.2161299265
0.2360588266
0.335703327
0.4021329943

0.9046884934
0.7386201225
0.182536164
0.0517028491

0.1557524739
0.0828773249
0.015625179
1.2627041992

0.3573468358
0.82893171716
0.342118434
0.3785790309

0.5384820933
0.3132925258
0.4796520139
0.1375116383

0.1126211304
0.0312756855
0.0847887058
0.1544221928

Table 4: Comparison of relative errors for Example 2.

6. Conclusion

In this paper, we presented the application of the homotopy analysis method (HAM) for solving a special form of nonlinear Volterra-Fredholm integro-differential equation.

The sufficient condition for the convergence of the method is illustrated and then verified for two examples.

As we can see in Tables 2 and 4, HAM solutions have a good agreement with the numerical results provided that appropriate values for the convergence control parameter are chosen.

The ability of the HAM is mainly due to the fact that the method provides a way to ensure the convergence of series solution.

Convergence conditions in applying HAM for other equations, and systems of ordinary and partial differential equations, integral equations, and integro-differential equations are also under study in our research team.

Conflict of Interests

The author declare that there is no conflict of interests regarding the publication of this paper.

References

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.

E. Babolian, Z. Masouri, and S. Hatamzadeh-Varmazyar, “Numerical solution of nonlinear Volterra-Fredholm integro-differential equations via direct method using triangular functions,” Computers and Mathematics with Applications, vol. 58, no. 2, pp. 239–247, 2009.

J. Biazar and H. Ghazvini, “Numerical solution for special non-linear Fredholm integral equation by HPM,” Applied Mathematics and Computation, vol. 195, no. 2, pp. 681–687, 2008.

J. Biazar, M. Gholami Porshokouhi, B. Ghanbari, and M. Gholami Porshokouhi, “Numerical solution of functional integral equations by the variational iteration method,” Journal of Computational and Applied Mathematics, vol. 235, no. 8, pp. 2581–2585, 2011.

A. Yildirim, “Solution of BVPs for fourth-order integro-differential equations by using homotopy perturbation method,” Computers and Mathematics with Applications, vol. 56, no. 12, pp. 3175–3180, 2008.

E. Babolian, Z. Masouri, and S. Hatamzadeh-Varmazyar, “Numerical solution of nonlinear Volterra-Fredholm integro-differential equations via direct method using triangular functions,” Computers and Mathematics with Applications, vol. 58, no. 2, pp. 239–247, 2009.

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.

K. Maleknejad and M. T. Kajani, “Solving linear integro-differential equation system by Galerkin methods with hybrid functions,” Applied Mathematics and Computation, vol. 159, no. 3, pp. 603–612, 2004.

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.

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

K. Maleknejad, B. Basirat, and E. Hashemizadeh, “Hybrid Legendre polynomials and Block-Pulse functions approach for nonlinear VolterraFredholm integro-differential equations,” Computers and Mathematics with Applications, vol. 61, no. 9, pp. 2821–2828, 2011.

S. J. Liao, Proposed homotopy analysis techniques for the solution of non-linear problem [Ph.D. dissertation], Shanghai Jiao Tong University, 1992.

S. J. Liao, Beyond Perturbation: Introduction to HAM, Chapman & Hall, CRC Press, Boca Raton, Fla, USA, 2003.

S. J. Liao, “Notes on the homotopy analysis method: some definitions and theorems,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 4, pp. 983–997, 2009.

B. Ghanbari, L. Rada, and K. Chen, “A restarted iterative homotopy analysis method for two nonlinear models from image processing,” International Journal of Computer Mathematics.

J. Biazar and B. Ghanbari, “The homotopy perturbation method for solving neutral functional-differential equations with proportional delays,” Journal of King Saud University, vol. 24, no. 1, pp. 33–37, 2012.

J. Biazar and B. Ghanbari, “HAM solution of some initial value problems arising in heat radiation equations,” Journal of King Saud University, vol. 24, no. 2, pp. 161–165, 2012.

A. Mehmood, S. Munawar, and A. Ali, “Comments to: “Homotopy analysis method for solving the MHD flow over a non-linear stretching sheet (Communications in Nonlinear Science and Numerical Simulation, vol. 14, pp. 2653–2663, 2009)”,” Communications in Nonlinear Science and Numerical Simulation, vol. 15, no. 12, pp. 4233–4240, 2010.

W. Rudin, Principles of Mathematical Analysis, McGraw-Hill Book, New York, NY, USA, 1953.

Z. M. Odibat, “A study on the convergence of homotopy analysis method,” Applied Mathematics and Computation, vol. 217, no. 2, pp. 782–789, 2010.