`Discrete Dynamics in Nature and SocietyVolume 2012 (2012), Article ID 365792, 18 pageshttp://dx.doi.org/10.1155/2012/365792`
Research Article

## Homotopy Analysis Method for Second-Order Boundary Value Problems of Integrodifferential Equations

1Department of Mathematics, Faculty of Science, Al Balqa Applied University, Salt 19117, Jordan
2Department of Mathematics, Faculty of Science, University of Jordan, Amman 11942, Jordan

Received 1 April 2012; Accepted 19 June 2012

#### Abstract

In this paper, series solution of second-order integrodifferential equations with boundary conditions of the Fredholm and Volterra types by means of the homotopy analysis method is considered. The new approach provides the solution in the form of a rapidly convergent series with easily computable components using symbolic computation software. The homotopy analysis method provides us with a simple way to adjust and control the convergence region of the infinite series solution by introducing an auxiliary parameter. The proposed technique is applied to a few test examples to illustrate the accuracy, efficiency, and applicability of the method. The results reveal that the method is very effective, straightforward, and simple.

#### 1. Introduction

Integrodifferential equations (IDEs) are often involved in the mathematical formulation of physical and engineering phenomena. IDEs can be encountered in various fields of science such as physics, chemistry, biology, and engineering. These kinds of equations can also be found in numerous applications, such as electromagnetic, plasma physics, elasticity, fluid dynamics, oscillation theory, polymer rheology, chemical kinetics, biomechanics, and control theory [15]. Since it is usually impossible to obtain the closed-form solutions to second-order boundary value problems of IDEs met in practice, these problems must be attacked by various approximate and numerical methods.

The purpose of this paper is to extend the application of the homotopy analysis method (HAM) to provide symbolic approximate solution for the second-order boundary value problems of IDEs of the following Fredholm type: and the following Volterra type: subject to the boundary conditions where are real finite constants, is an unknown function to be determined, and are continuous functions on and , respectively, and are nonlinear continuous functions of . For details about the existence and uniqueness of solutions for such problems, see [68]. However, we assume that (1.1) and (1.2) subject to the boundary conditions (1.3) have a unique analytic solution on the given interval.

The numerical solvability of second-order IDEs with boundary conditions of the Fredholm and Volterra types and other related equations has been pursued by several authors. To mention a few, in [9], the authors have discussed the Legendre polynomials method for solving Fredholm equation . Furthermore, the compact finite difference method is carried out in [10] for the Volterra IDE and the Fredholm IDE . Recently, the monotone iterative sequences method for solving Volterra equation is proposed in [11].

The HAM, which is proposed by Liao [1217], is effectively and easily used to solve some classes of linear and nonlinear problems without linearization, perturbation, or discretization. The HAM is based on the homotopy, a basic concept in topology. The auxiliary parameter is introduced to construct the so-called zero-order deformation equation. Thus, unlike all previous analytic techniques, the HAM provides us with a family of solution expressions in auxiliary parameter . As a result, the convergence region and rate of solution series are dependent upon the auxiliary parameter and thus can be greatly enlarged by means of choosing a proper value of . This provides us with a convenient way to adjust and control convergence region and rate of solution series given by the HAM.

In the recent years, extensive work has been done using HAM, which provides analytical approximations for linear and nonlinear equations. This method has been implemented in several functional equations, such as nonlinear water waves [15], unsteady boundary-layer flows [16], solitary waves with discontinuity [17], Klein-Gordon equation [18], projectile motion with the quadratic resistance law [19], systems of fractional differential equations [20], nonlinear fractional differential equations [21, 22], systems of fractional algebraic-differential equations [23], fractional SIR model [24], singular Volterra integral equation [25], charged particle motion for certain configurations of oscillating magnetic fields [26], Volterra population growth model [27], systems of fractional IDEs [28], and high-order IDEs [29].

The outline of the paper is as follows. In the next section, the basic idea of the HAM is introduced. In Section 3, we utilize the statement of the HAM for solving second-order IDEs with boundary conditions of the Fredholm and Volterra types. In Section 4, numerical examples are given to illustrate the capability of HAM. This paper ends in Section 5 with some concluding remarks.

#### 2. Homotopy Analysis Method

The principles of the HAM and its applicability for various kinds of functional equations are given in [1228]. For convenience of the reader, we will present a review of the HAM [1217] then we will extend the HAM to construct a symbolic approximate solution to second-order IDEs with boundary conditions of the Fredholm and Volterra types. To achieve our goal, we consider the following nonlinear functional equation: where N is a nonlinear operator and is an unknown function of independent variable .

Let denote an initial guess of the exact solution of (2.1), an auxiliary parameter, an auxiliary function, and an auxiliary linear operator with the property when . The auxiliary parameter , the auxiliary function , and the auxiliary linear operator play important roles within the HAM [12]. As we will see later, choosing suitable values of the auxiliary parameter will help us to adjust and control the convergence region of the series solution.

Liao [1217] constructs, using as an embedding parameter, the so-called zero-order deformation equation as where is the solution of (2.1) which depends on , , , , and . When , the zero-order deformation equation (2.2) becomes and when , since and , the zero-order deformation equation (2.2) reduces to

So, is exactly the solution of the nonlinear equation (2.1). Thus, according to (2.3) and (2.4), as increases from 0 to 1, the solution varies continuously from the initial approximation to the exact solution .

Define the so-called th-order deformation derivatives expanding in Taylor series with respect to the embedding parameter , using (2.3) and (2.5), we obtain

Assume that the auxiliary parameter, the auxiliary function , the initial approximation , and the auxiliary linear operator are properly chosen so that the power series (2.6) of converges at . Then, we have under these assumptions the series solution .

Define the vector . Differentiating equation (2.2) -times with respect to embedding parameter , then setting and dividing them by . Using (2.5), we have the so-called th-order deformation equation as where

For any given nonlinear operator , the term can be easily expressed by (2.8). Thus, we can gain by means of solving the linear high-order deformation equation (2.7) one after the other in order. The th-order approximation of is given by .

#### 3. Solving Second-Order Boundary Value Problems of IDEs by HAM

Throughout this section, we will utilize the construction for the IDEs of the Fredholm type in order not to increase the length of the paper without the loss of generality for the remaining type. However, similar construction can be implemented for the Volterra type.

Accordingly, we extend the application of the HAM to solve the second-order IDEs of the Fredholm type subject to the boundary conditions (1.3).

First of all, we assume that (3.1) satisfies the initial condition ,  where the unknown constant can be later determined by substituting the boundary condition into the obtained solution. After that, we chose the initial guess such that the initial conditions and are satisfied. That is, the initial guess will be of the form . Finally, we take the auxiliary linear operator which satisfies the property , where are constants of integrations.

The so-called zero-order deformation equation will be defined as where when , (3.2) reduces to . In this case, the so-called th-order deformation equation can be constructed as Operating the operator , the inverse operator of to both sides of (3.4), then, the th-order deformation equation will have the following form: where and the th-order approximation of can be given as

If we substitute the boundary condition into (3.7), then we obtain an equation of the variable . And so, the approximate solution for (3.1) and (1.3) will be completely constructed. For simplicity, in this paper, we assume that the value of the auxiliary function is kept at .

#### 4. Numerical Results and Discussion

HAM provides an analytical approximate solution in terms of an infinite power series. However, there is a practical need to evaluate this solution and to obtain numerical values from the infinite power series. The consequent series truncation and the practical procedure are conducted to accomplish this task. In this section, we consider four examples to demonstrate the performance and efficiency of the present technique. Through this paper all the symbolic and numerical computations are performed by using Mathematica 7.0 software package.

To show the accuracy of the present method for our problems, we report three types of error. The first one is the residual error, , defined as for the Fredholm type and similarly for the Volterra type, while the exact, , and relative, , errors are defined, respectively, by where , is the th-order approximation of obtained by the HAM, and is the exact solution.

Example 4.1. Consider the following linear Fredholm IDE: subject to the boundary conditions The exact solution is , where , , and .
According to (3.6) and (2.5), we have
Assuming the initial approximation of (4.3) has the form , which satisfies the initial conditions and . Consequently, the first few terms of the HAM series solution for (4.3) according to these initial conditions are as follows: To determine the value of , we must chose a value to the auxiliary parameter . It was proved that if we set , then we have the Adomian decomposition solution which is a special case of the HAM solution [30, 31]. Therefore, is available value for. Now, substitute and the boundary condition at into the 10th order approximation of to get .
Figure 1 shows the 10th-order approximate solution of (4.3) and (4.4) obtained by HAM at together with the exact solution. It is clear from the figure that the approximate solution is in high agreement with the exact solution.
Since the value of has been identified previously, the valid region for the values of can be obtained by the so-called -curve which correspond to the horizontal line as shown in Figure 2. Thus, the valid region of is .

Figure 1: Solutions of (4.3) and (4.4). Dotted line: the 10th-order approximate HAM solution at and solid line: the exact solution.
Figure 2: The -curves of and which are corresponding to the 10th-order approximate HAM solution of (4.3) and (4.4). Dotted line: and solid line: .

Figure 3 shows a comparison between the exact and residual errors obtained from the 10th-order approximation of for (4.3) and (4.4) at , while Figures 4 and 5 show respectively how the exact and residual errors change with varying the values of .

Figure 3: Errors of (4.3) and (4.4) at . Dotted line: the exact error and solid line: the residual error.
Figure 4: Exact error of (4.3) and (4.4) at different values of . Dash-dotted line: , dotted line: , and solid line: .
Figure 5: Residual error of (4.3) and (4.4) at different values of . Dash-dotted line: , dotted line: , and solid line: .

Example 4.2. Consider the following nonlinear Fredholm IDE: subject to the boundary conditions The exact solution is .
According to (3.6) and (2.5), we have Choose the initial guess approximation as , which satisfies the initial conditions and . Then, we have the following approximations terms of the HAM series solution for (4.8): As in the previous example, we determine an introductory value of the constant by substituting and the boundary condition at into the obtained HAM solution. So, the 5th-order approximation of gives . Figure 6 shows the 5th-order approximate solution of (4.8) and (4.9) obtained by HAM at together with the exact solution, while Figure 7 shows the -curve corresponding to the 5th-order approximate HAM solution. It is clear from Figure 7 that the valid region of is .
Our next goal is to show how the auxiliary parameter affects the approximate solutions. However, it is evident from Figure 6 that both curves corresponding to the th-order approximate HAM solution and the exact solution are not coinciding in general at . To improve the solution, let us compute other values of at different values of. Table 1 shows the values of the constant at different values of the auxiliary parameter.
Representation the exact solution and the 5th-order approximate HAM solution of (4.8) and (4.9) at different values of mentioned in Table 1 shows that all these solutions are approximately identical. However, Figure 8 shows the correspondence between the exact solution and the 5th-order approximate HAM solution of (4.8) and (4.9) at . From the last mentioned figure, we see that we can achieve a good approximation with the exact solution by using a few terms in HAM.
Figure 9 shows a comparison between the exact errors of (4.8) and (4.9) at different values of corresponding to the 5th-order approximate HAM solution. It is to be noted that the accuracy of certain node at is more convincing than the other values of .

Table 1: The values of at different values of .
Figure 6: Solutions of (4.8) and (4.9). Dotted line: the 5th-order approximate HAM solution at and solid line: the exact solution.
Figure 7: The -curves of and which are corresponding to the 5th-order approximate HAM solution of (4.8) and (4.9). Dotted line: and solid line: .
Figure 8: Solutions of (4.8) and (4.9). Dotted line: the 5th-order approximate HAM solution at and solid line: the exact solution.
Figure 9: Exact error of (4.8) and (4.9) at different values of . Dash-dotted line: , dotted line: , and solid line: .

Example 4.3. Consider the following linear Volterra IDE: subject to the boundary conditions The exact solution is .
According to (3.6) and (2.5), we have Choose the initial guess approximation as , which satisfies the initial conditions and . Using the iteration formula (3.5), we can directly obtain the following approximations terms of the HAM series solution for (4.13): The rest of components of the iteration formula (3.5) can be obtained in a similar way.
To find the value of , substitute and the boundary condition at into the 5th-order approximation of to get . After that, substituting the value of in the 5th-order approximation of and plotting the -curve, we can obtain the valid region of which is as shown in Figure 10.

Figure 10: The -curves of and which are corresponding to the 5th-order approximate HAM solution of (4.13) and (4.14). Dotted line: and solid line: .

The graph of the 5th-order approximate HAM solution of (4.13) and (4.14) at together with the exact solution are depicted in Figure 11. It is evident from the figure that the 5th-order approximate HAM solution agrees with the exact solution in the interval .

Figure 11: Solutions of (4.13) and (4.14). Dotted line: the 5th-order approximate HAM solution at and solid line: the exact solution.

In Table 2, the exact, relative, and residual errors have been calculated for various in to measure the extent of agreement between the 5th-order approximate HAM solution at and the exact solution. From the table, it can be seen that the HAM provides us with the accurate approximate solution of (4.13) and (4.14).

Table 2: Exact, relative, and residual errors of (4.13) and (4.14) at and.

Example 4.4. Consider the following nonlinear Volterra IDE: subject to the boundary conditions The exact solution is .
According to (3.6) and (2.5), we have Assume that the initial approximation of (4.18) has the form , which satisfies the initial conditions and . Using the iteration formula (3.5), we can directly obtain the following approximations terms of the HAM series solution for (4.18) subject to these initial conditions: Similarly to the previous discussion, one can substitute and the boundary condition at into the 5th-order approximation of to get . Figure 12 shows the agreement between the 5th-order approximate HAM solution of (4.18) and (4.19) at and the exact solution.

Figure 12: Solutions of (4.18) and (4.19). Dotted line: the 5th-order approximate HAM solution at and solid line: the exact solution.

Since the value of has been identified previously, the valid region for the values of can be obtained by the so-called -curve which correspond to the horizontal line as shown in Figure 13. Thus, it is clear from the figure that the valid region of is .

Figure 13: The -curves of and which are corresponding to the 5th-order approximate HAM solution of (4.18) and (4.19). Dotted line: and solid line: .

The detailed data of for (4.18) and (4.19) that includes the exact nodal values, the 5th-order approximate HAM solution nodal values, and the exact error are given in Table 3. It is clear that the accuracy obtained using HAM is in advanced by using only five terms approximations.

Table 3: Numerical results of (4.18) and (4.19) at and .

#### 5. Conclusion

The main concern of this work has been to propose an efficient algorithm for the solution of second-order IDEs with boundary conditions of the Fredholm and Volterra types. The goal has been achieved by extending the HAM to solve this class of boundary value problems. We can conclude that the HAM is a powerful and efficient technique in finding approximate solutions for linear and nonlinear boundary value problems of second-order IDEs of different types. The proposed algorithm produced a rapidly convergent series by choosing suitable values of the auxiliary parameter. After computing several approximants and using the boundary conditions at the boundary points, we can easily determine the approximate solution.

There are two important points to make here. First, HAM provides us with a simple way to adjust and control the convergence region of the series solution by introducing the auxiliary parameter. Second, the results obtained by HAM are very effective and convenient in linear and nonlinear cases with less computational work. This confirms our belief that the efficiency of our technique gives it much wider applicability for general classes of linear and nonlinear problems.

#### References

1. R. P. Kanwal, Linear Integral Differential Equations: Theory and Technique, Birkhauser, Boston, Ga, USA, 2nd edition, 1996.
2. F. Bloom, “Asymptotic bounds for solutions to a system of damped integro-differential equations of electromagnetic theory,” Journal of Mathematical Analysis and Applications, vol. 73, no. 2, pp. 524–542, 1980.
3. K. Holmåker, “Global asymptotic stability for a stationary solution of a system of integro-differential equations describing the formation of liver zones,” SIAM Journal on Mathematical Analysis, vol. 24, no. 1, pp. 116–128, 1993.
4. L. K. Forbes, S. Crozier, and D. M. Doddrell, “Calculating current densities and fields produced by shielded magnetic resonance imaging probes,” SIAM Journal on Applied Mathematics, vol. 57, no. 2, pp. 401–425, 1997.
5. A. J. Jerri, Introduction to Integral Equations with Applications, John Wiley & Sons, New York, NY, USA, 2nd edition, 1999.
6. R. P. Agarwal, “Boundary value problems for higher order integro-differential equations,” Nonlinear Analysis, vol. 7, no. 3, pp. 259–270, 1983.
7. R. P. Agarwal, Boundary Value Problems for High Ordinary Differential Equations, World Scientific, Singapore, 1986.
8. Z. Wang, L. Liu, and Y. Wu, “The unique solution of boundary value problems for nonlinear second-order integral-differential equations of mixed type in Banach spaces,” Computers & Mathematics with Applications, vol. 54, no. 9-10, pp. 1293–1301, 2007.
9. A. Saadatmandi and M. Dehghan, “Numerical solution of the higher-order linear Fredholm integro-differential-difference equation with variable coefficients,” Computers & Mathematics with Applications, vol. 59, no. 8, pp. 2996–3004, 2010.
10. J. Zhao and R. M. Corless, “Compact finite difference method for integro-differential equations,” Applied Mathematics & Computation, vol. 177, no. 1, pp. 271–288, 2006.
11. Q. M. Al-Mdallal, “Monotone iterative sequences for nonlinear integro-differential equations of second order,” Nonlinear Analysis, vol. 12, no. 6, pp. 3665–3673, 2011.
12. S. J. Liao, Beyond Perturbation: Introduction to the Homotopy Analysis Methods, Chapman & Hall, Boca Raton, Fla, USA, 2003.
13. S. Liao, “On the homotopy analysis method for nonlinear problems,” Applied Mathematics & Computation, vol. 147, no. 2, pp. 499–513, 2004.
14. S. J. Liao, “Homotopy analysis method: a new analytic method for nonlinear problems,” Applied Mathematics & Mechanics, vol. 19, no. 10, pp. 957–962, 1998.
15. S. J. Liao and K. F. Cheung, “Homotopy analysis of nonlinear progressive waves in deep water,” Journal of Engineering Mathematics, vol. 45, no. 2, pp. 105–116, 2003.
16. S. Liao, “Series solutions of unsteady boundary-layer flows over a stretching flat plate,” Studies in Applied Mathematics, vol. 117, no. 3, pp. 239–263, 2006.
17. W. Wu and S. J. Liao, “Solving solitary waves with discontinuity by means of the homotopy analysis method,” Chaos, Solitons and Fractals, vol. 26, no. 1, pp. 177–185, 2005.
18. Q. Sun, “Solving the Klein-Gordon equation by means of the homotopy analysis method,” Applied Mathematics & Computation, vol. 169, no. 1, pp. 355–365, 2005.
19. K. Yabushita, M. Yamashita, and K. Tsuboi, “An analytic solution of projectile motion with the quadratic resistance law using the homotopy analysis method,” Journal of Physics A, vol. 40, no. 29, pp. 8403–8416, 2007.
20. M. Zurigat, S. Momani, Z. Odibat, and A. Alawneh, “The homotopy analysis method for handling systems of fractional differential equations,” Applied Mathematical Modelling, vol. 34, no. 1, pp. 24–35, 2010.
21. Z. Odibat, S. Momani, and H. Xu, “A reliable algorithm of homotopy analysis method for solving nonlinear fractional differential equations,” Applied Mathematical Modelling, vol. 34, no. 3, pp. 593–600, 2010.
22. I. Hashim, O. Abdulaziz, and S. Momani, “Homotopy analysis method for fractional IVPs,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 3, pp. 674–684, 2009.
23. M. Zurigat, S. Momani, and A. Alawneh, “Analytical approximate solutions of systems of fractional algebraic-differential equations by homotopy analysis method,” Computers & Mathematics with Applications, vol. 59, no. 3, pp. 1227–1235, 2010.
24. O. Abu Arqub and A. El-Ajou, “Solution of the fractional epidemic model by homotopy analysis method,” Journal of King Saud University. In press.
25. 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.
26. H. S. Nik, S. Effati, and R. Buzhabadi, “Analytic-approximate solution for an integro-differential equation arising in oscillating magnetic fields using homotopy analysis method,” Iranian Journal of Optimization, vol. 2, no. 3, pp. 518–535, 2010.
27. N. A. Khan, A. Ara, and M. Jamil, “Approximations of the nonlinear Volterra's population model by an efficient numerical method,” Mathematical Methods in the Applied Sciences, vol. 34, no. 14, pp. 1733–1738, 2011.
28. M. Zurigat, S. Momani, and A. Alawneh, “Homotopy analysis method for systems of fractional integro-differential equations,” Neural, Parallel & Scientific Computations, vol. 17, no. 2, pp. 169–186, 2009.
29. H. M. Jaradat, F. Awawdeh, and O. Alsayyed, “Series solution to the high-order integro-differential equations,” Analele Universităţii din Oradea—Fascicola Matematică, vol. 16, pp. 247–257, 2009, Tom XVI.
30. J. Cang, Y. Tan, H. Xu, and S. J. Liao, “Series solutions of non-linear Riccati differential equations with fractional order,” Chaos, Solitons and Fractals, vol. 40, no. 1, pp. 1–9, 2009.
31. F. M. Allan, “Derivation of the Adomian decomposition method using the homotopy analysis method,” Applied Mathematics & Computation, vol. 190, no. 1, pp. 6–14, 2007.