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.
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.
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.