`Abstract and Applied AnalysisVolume 2012, Article ID 960289, 16 pageshttp://dx.doi.org/10.1155/2012/960289`
Research Article

## On Spectral Homotopy Analysis Method for Solving Linear Volterra and Fredholm Integrodifferential Equations

Department of Mathematics, Universiti Putra Malaysia (UPM), Selangor, 43400 Serdang, Malaysia

Received 12 June 2012; Revised 19 September 2012; Accepted 4 October 2012

#### Abstract

A modification of homotopy analysis method (HAM) known as spectral homotopy analysis method (SHAM) is proposed to solve linear Volterra integrodifferential equations. Some examples are given in order to test the efficiency and the accuracy of the proposed method. The SHAM results show that the proposed approach is quite reasonable when compared to SHAM results and exact solutions.

#### 1. Introduction

Functional equations such as partial differential equations, integral and integrodifferential equations and others are widely employed to model complex real-life problems. Many physical events in different fields of sciences and engineering have been formulated using integrodifferential equations [1, 2]. Integrodifferential equations are usually difficult to solve analytically so there is a need to obtain an efficient approximate solution [3, 4]. Homotopy analysis method (HAM) was first proposed by Liao in [5]. it is based on the homotopy concept in topology for solving nonlinear differential equations. Unlike the traditional perturbation methods like Lyapunov's artificial small parameter method [6], Adomian decomposition method [710] and the -expansion method [11], which are the special cases of HAM, this approach does not need a small perturbation parameter. In the HAM, original nonlinear problem is converted to infinite number of linear problems without using the perturbation techniques [12]. Homotopy analysis method is more powerful than traditional perturbation methods since it is applicable for solving problems with strong nonlinearity [1315], even if they do not have any small or large parameters.

We can adjust the convergence region and the rate of approximation series to give us freedom to use different base function to approximate a nonlinear problem. He in [16] proposed the so called Homotopy perturbation method (HPM). Later it was pointed by Hayat et al. [17] that HPM is only a special case of the HAM . Many nonlinear problems were solved by Homotopy perturbation method [18, 19]. The HPM does not provide us with convergent series solution for strongly nonlinear problems as was indicated by Abbasbandy in [20] using a simple problem in which the physical parameters were large. In 2010, Motsa et al. [21] employed Chebyshev polynomials to solve higher-order deformation equation and called his approach spectral homotopy analysis method (SHAM). The spectral homotopy analysis method has been used only for solving partial and ordinary differential equations [22, 23]. Motsa et al. [2123] found that the Spectral homotopy analysis method is more efficient than the homotopy analysis method as it does not depend on the rule of solution expression and the rule of ergodicity. This method is more flexible than homotopy analysis method since it allows for a wider range of linear and nonlinear operators and one is not restricted to use the method of higher order differential mapping for solving boundary value problems in bounded domains, unlike the homotopy analysis method. The range of admissible values is much wider in spectral homotopy analysis method than in homotopy analysis method. The main restriction of HAM in solving integral equations is to choose the best initial guess as the series solution be convergent. In SHAM the initial approximation is taken to be the solution of the nonhomogeneous linear part of the given equation.

In this paper, we propose spectral homotopy analysis method (SHAM) to solve linear Volterra and Fredholm type of integrodifferential equations. Volterra integrodifferential equation is given by

The paper is organized in the following way. Section 2 included a brief introduction in homotopy analysis method. Spectral homotopy analysis method for solving integral equations is presented in Section 3. We then propose the way to calculate which is needed to obtain in Section 4. Numerical examples are presented in Section 5. In Section 6, concluding remarks are given.

#### 2. Homotopy Analysis Solution

In this section, we give a brief introduction to HAM. We consider the following differential equation in a general form: where is nonlinear operator, denotes independent variables and is an unknown function, respectively. For simplicity we disregard all initial and all boundary conditions which can be dealt in similar way. The so-called zero-order deformation equation was constructed by Liao as follows: where is the embedding parameter and, is a nonzero convergence-parameter, is an auxiliary function and is called an initial guess of and is an unknown function. In addition, is an auxiliary linear operator and is nonlinear operator as follows: with the property where and are constants and is a nonlinear operator. Obviously, when and , it holds and . In this way, as increase from to , alter from initial guess to the solution and is expanded in Taylor series with respect to : where

The series (2.5) converges at if the auxiliary linear operator, the initial guess, the convergence parameter, and the auxiliary function are properly selected:

The admissible and valid values of the convergence-parameter are found from the horizontal portion of the -curves. Liao proved that is one of the solutions of original nonlinear equation. As so (2.2) becomes

Define the vector . Operating on both sides with , we have the so called th-order deformation equation where for is governed by the linear equation (2.10) and can be solved by symbolic computation software such as Maple, Matlab, and so on.

#### 3. Spectral Homotopy Analysis Solution

Consider the second-order Volterra integrodifferential equation

To obtain initial approximation we solve the following two-point boundary value problem: subject to boundary conditions We use the Chebyshev pseudospectral method to solve (3.2). At first we approximate by a truncated series of Chebyshev polynomial of the form where is the th Chebyshev polynomials, , are coefficients and are Gauss-Lobatto points which are the extrema of the th-order Chebyshev polynomial, defined by

Derivatives of the functions at the collocation points are represented as where is the order of differentiation and is the Chebyshev spectral differentiation matrix. By following [24] we express the entries of the differentiation matrix as Substituting (3.4)–(3.6) into (3.2) results in subject to the boundary conditions where

To satisfy the boundary conditions we delete the first and the last rows and columns of and the first and the last rows of and . The values of , are determined from the equation which is the initial approximation for the SHAM solution of the governing equation (3.1). The zeroth-order deformation equation is given by where is the nonzero convergence controlling auxiliary parameter, and is a nonlinear operator given by

Differentiating (3.2) times with respect to the embedding parameter and then setting and finally dividing them by , we have the so-called th-order deformation equation subject to boundary conditions where

Applying the Chebyshev pseudospectral transformation on (3.15)–(3.17) results in subject to the boundary conditions where and are as defined in (3.10) and To implement the boundary conditions we delete the first and the last rows of and the first and the last rows and columns of .

Finally, this recursive formula can be written as following for : so, starting from the initial approximation. We can obtain higher-order approximation for , recursively.

#### 4. Calculation of

Consider the well-known Chebyshev polynomials of first kind of degree , defined by [25]

Also they are derived from the following recursive formula:

In this section we approximate the integrand in (3.15) by letting where , are the basis functions.

Clearly, , can be obtained by solving the following matrix equation: where , are Chebyshev collocation points defined in (2.6). We can obtain by solving the above system so the integral in (3.15) can be rewritten as following: where where as in [26], that is the operational matrix for integration as follows:

#### 5. Numerical Examples

In this section we apply the technique described in Section 3 to some illustrative examples of second-order Volterra integrodifferential equations.

Example 5.1. Consider the second-order Volterra integrodifferential equation subject to the boundary conditions with the exact solution .

We employ HAM and SHAM to solve this example. In Table 1, there is a comparison of the numerical result against the HAM and SHAM approximation solutions at different orders. It is worth noting that the SHAM results become very highly accurate only with a few iterations and fourth-order results become very close to the exact solution. Comparison of the numerical solution with the 4th-order SHAM solution for is made in Figure 1. It can be seen from the -curves (Figures 2 and 3) that if and then the SHAM an HAM approximate solutions are in well agreement with the exact solution. In Table 2, we present that the SHAM results for another value of in (), are also very highly accurate, and in sixth order we obtain the exact numerical solutions. In HAM we choose as initial guess.

Table 1: Numerical result of Example 5.1 against the HAM and the SHAM solutions whit .
Table 2: Numerical result of Example 5.1 against the SHAM solutions whit .
Figure 1: Comparison of the numerical solution with the 4th-order SHAM solution for .
Figure 2: 4th-order SHAM -curve.
Figure 3: 11th-order HAM -curve.

Example 5.2. Consider the second-order Volterra integrodifferential equation subject to the boundary conditions with the exact solution .

Similar to the previous example, we employ HAM and SHAM to solve this example. In Table 3, there is a comparison of the numerical results against the HAM and SHAM approximation solutions (5.2) at different orders. It is worth noting that the rate of convergence in SHAM results is higher than the HAM results. Comparison of the numerical solution with the 16th-order SHAM solution for is made in Figure 4. It can be seen from the -curves (Figures 5 and 6) that if and then the approximate solution is in well agreement with the exact solution. In HAM, we choose as initial guess.

Table 3: Numerical result of Example 5.2 against the HAM and the SHAM solutions whit .
Figure 4: Comparison of the numerical solution with the 16th-order SHAM solution for .
Figure 5: 16th-order SHAM -curve3.
Figure 6: 23th-order HAM -curve.

Example 5.3. Consider the first-order Fredholm integrodifferential equation [27] subject to the boundary conditions .

In this example, by a change of variables, the Fredholm-integro differential equation with nonhomogeneous conditions are transformed to the equation with homogeneous conditions and then apply the spectral homotopy analysis method for this problem. According to governing equation and the boundary condition it is reasonable to set the transformation to make the governing boundary conditions homogeneous. Substituting (5.4) in (5.3) and boundary conditions gives subject to the boundary condition

Comparison between absolute errors in solutions by SHAM and Legendre collocation matrix method (LCMM) is tabulated in Table 4. It is worth noting that the SHAM results become very highly accurate only with four iterations.

Table 4: Absolute errors for Example 5.3.

Example 5.4. Consider the following second-order Fredholm integrodifferential equation [28]: where , .

In this example, by a change of variables, the Fredholm-integro differential equation with nonhomogeneous conditions are transformed to the equation with homogeneous conditions and then we apply the spectral homotopy analysis method for this problem. According to governing equation and the boundary condition it is reasonable to set the transformation to make the governing boundary conditions homogeneous. Substituting (5.8) in (5.7) and boundary conditions gives subject to boundary condition

Comparison between absolute errors in solutions by SHAM and HAM is tabulated in Table 5. It is worth noting that the SHAM results become very highly accurate only with four iterations.

Table 5: Absolute errors for Example 5.4.

#### 6. Conclusion

In this paper for the first time in the literature we presented the application of spectral homotopy analysis method (SHAM) in solving linear integro differential equations. A comparison was made between exact analytic solutions and numerical results from the spectral homotopy analysis method and homotopy analysis method (HAM) solutions. The numerical results indicate that in SHAM the rate of convergence is faster than HAM. For example, in Example 5.1 we found that the fourth-order SHAM approximation sufficiently gives a match with the numerical results up to twelve decimal places.In contrast, HAM solutions have a good agreement with the numerical results in 11th-order. In this paper we choose the admissible value of from -curve for and . As it is shown in Figures 2 and 3 the range of admissible values is much wider in SHAM than in HAM. In using HAM there is restriction on choosing an initial approximation that the higher-order deformation can be integrated whereas in SHAM the initial approximation is the solution of the nonhomogeneous linear part of (3.1). The spectral homotopy analysis method is more efficient as it does not depend on the rule of coefficient ergodicity unlike the HAM. With some changes we can apply this method on Fredholm integrodifferential equations. In Example 5.3 [27], we compare the SHAM solution with Legendre collocation matrix method (LCCM) to solve linear Fredholm integrodifferential equation. As can be seen in Table 4 SHAM is more accurate and efficient than LCCM. In this work we obtained the numerical results up to twelve decimal places having well agreement with exact solution but reported them only up to six decimal places.

In this paper, we described the spectral homotopy analysis method to solve linear Volterra and Fredholm integrodifferential equations; however, it remains to be generalized and verified for more complicated integral equations that we consider it as future works.

#### Acknowledgment

The authors gratefully acknowledge that this research was partially supported by the Universiti Putra Malaysia under the ERGS Grant Scheme having project number 5527068.

#### References

1. 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.
2. K. Parand, S. Abbasbandy, S. Kazem, and J. A. Rad, “A novel application of radial basis functions for solving a model of first-order integro-ordinary differential equation,” Communications in Nonlinear Science and Numerical Simulation, vol. 16, no. 11, pp. 4250–4258, 2011.
3. P. Darania and A. Ebadian, “A method for the numerical solution of the integro-differential equations,” Applied Mathematics and Computation, vol. 188, no. 1, pp. 657–668, 2007.
4. A. Karamete and M. Sezer, “A Taylor collocation method for the solution of linear integro-differential equations,” International Journal of Computer Mathematics, vol. 79, no. 9, pp. 987–1000, 2002.
5. S. J. Liao, The proposed homotopy analysis technique for the solution of nonlinear problems [Ph.D. thesis], Shanghai Jiao Tong University, 1992.
6. A. M. Lyapunov, The General Problem of the Stability of Motion, Taylor & Francis, London, UK, 1992.
7. G. Adomian, “A review of the decomposition method and some recent results for nonlinear equations,” Mathematical and Computer Modelling, vol. 13, no. 7, pp. 17–43, 1990.
8. G. Adomian and R. Rach, “Noise terms in decomposition solution series,” Computers & Mathematics with Applications, vol. 24, no. 11, pp. 61–64, 1992.
9. G. Adomian and R. Rach, “Analytic solution of nonlinear boundary value problems in several dimensions by decomposition,” Journal of Mathematical Analysis and Applications, vol. 174, no. 1, pp. 118–137, 1993.
10. G. Adomian, Solving Frontier Problems of Physics: The Decomposition Method, vol. 60, Kluwer Academic, Boston, Mass, USA, 1994.
11. P. K. Bera and J. Datta, “Linear delta expansion technique for the solution of anharmonic oscillations,” Pramana Journal of Physics, vol. 68, no. 1, pp. 117–122, 2007.
12. S. Liao, “On the homotopy analysis method for nonlinear problems,” Applied Mathematics and Computation, vol. 147, no. 2, pp. 499–513, 2004.
13. J. H. He, “The homotopy perturbation method for nonlinear oscillators with discontinuities,” Applied Mathematics and Computation, vol. 151, no. 1, pp. 287–292, 2004.
14. G. Domairry, A. Mohsenzadeh, and M. Famouri, “The application of homotopy analysis method to solve nonlinear differential equation governing Jeffery-Hamel flow,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 1, pp. 85–95, 2009.
15. Z. Ziabakhsh and G. Domairry, “Solution of the laminar viscous flow in a semi-porous channel in the presence of a uniform magnetic field by using the homotopy analysis method,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 4, pp. 1284–1294, 2009.
16. J.-H. He, “Homotopy perturbation technique,” Computer Methods in Applied Mechanics and Engineering, vol. 178, no. 3-4, pp. 257–262, 1999.
17. T. Hayat, S. Noreen, and M. Sajid, “Heat transfer analysis of the steady flow of a fourth grade fluid,” International Journal of Thermal Sciences, vol. 47, no. 5, pp. 591–599, 2008.
18. M. Jalaal, D. D. Ganji, and G. Ahmadi, “Analytical investigation on acceleration motion of a vertically falling spherical particle in incompressible Newtonian media,” Advanced Powder Technology, vol. 21, no. 3, pp. 298–304, 2010.
19. D. D. Ganji, G. A. Afrouzi, and R. A. Talarposhti, “Application of variational iteration method and homotopy-perturbation method for nonlinear heat diffusion and heat transfer equations,” Physics Letters, Section A, vol. 368, no. 6, pp. 450–457, 2007.
20. S. Abbasbandy, “Modified homotopy perturbation method for nonlinear equations and comparison with Adomian decomposition method,” Applied Mathematics and Computation, vol. 172, no. 1, pp. 431–438, 2006.
21. S. S. Motsa, P. Sibanda, and S. Shateyi, “A new spectral-homotopy analysis method for solving a nonlinear second order BVP,” Communications in Nonlinear Science and Numerical Simulation, vol. 15, no. 9, pp. 2293–2302, 2010.
22. S. S. Motsa, P. Sibanda, F. G. Awad, and S. Shateyi, “A new spectral-homotopy analysis method for the MHD Jeffery-Hamel problem,” Computers & Fluids, vol. 39, no. 7, pp. 1219–1225, 2010.
23. S. S. Motsa and P. Sibanda, “A new algorithm for solving singular IVPs of Lane-Emden type,” in Proceedings of the 4th International Conference on Applied Mathematics, Simulation, Modelling (ASM '10), pp. 176–180, Corfu Island, Greece, July 2010.
24. W. S. Don and A. Solomonoff, “Accuracy and speed in computing the Chebyshev collocation derivative,” SIAM Journal on Scientific Computing, vol. 16, no. 6, pp. 1253–1268, 1995.
25. L. M. Delves and J. L. Mohamed, Computational Methods for Integral Equations, Cambridge University Press, Cambridge, UK, 1985.
26. K. Maleknejad, S. Sohrabi, and Y. Rostami, “Numerical solution of nonlinear Volterra integral equations of the second kind by using Chebyshev polynomials,” Applied Mathematics and Computation, vol. 188, no. 1, pp. 123–128, 2007.
27. S. Yalçinbaş, M. Sezer, and H. H. Sorkun, “Legendre polynomial solutions of high-order linear Fredholm integro-differential equations,” Applied Mathematics and Computation, vol. 210, no. 2, pp. 334–349, 2009.
28. H. M. Jaradat, F. Awawdeh, and O. Alsayyed, “Series solution to the high-order integro-differential equations,” Fascicola Matematica, vol. 16, pp. 247–257, 2009.