Mathematical Problems in Engineering

Volume 2013 (2013), Article ID 261645, 7 pages

http://dx.doi.org/10.1155/2013/261645

Research Article

The Approximate Solutions of Fredholm Integrodifferential-Difference Equations with Variable Coefficients via Homotopy Analysis Method

Seydi Battal Gazi Karakoç, Aytekin Eryılmaz, and Musa Başbük

Department of Mathematics, Nevşehir University, 50300 Nevşehir, Turkey

Received 8 January 2013; Accepted 25 February 2013

Academic Editor: Safa Bozkurt Coskun

Copyright © 2013 Seydi Battal Gazi Karakoç et al. 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.

Numerical solutions of linear and nonlinear integrodifferential-difference equations are presented using homotopy analysis method. The aim of the paper is to present an efficient numerical procedure for solving linear and nonlinear integrodifferential-difference equations. The reliability and efficiency of the proposed scheme are demonstrated by some numerical experiments and performed on the computer algebraic system.

1. Introduction

Currently, considerable interest in mixed integrodifferential-difference equations (IDDEs) has been stimulated due to their numerous applications in the areas of engineering, mechanics, physics, chemistry, astronomy, biology, economics, potential theory, and electrostatics [1–7]. In integrodifferential-difference equations, the unknown function appears to be under integration sign, and it may also include the derivatives and functional arguments of the unknown function. This type of equations can be grouped into Fredholm integrodifferential-difference equations and Volterra integrodifferential-difference equations. The upper bound of the integral part of Volterra type is variable, while it is a fixed number for that of Fredholm type [1]. Since IDDEs are usually difficult to solve in an analytical manner, or to obtain closed form solution, a numerical method is needed. Many numerical methods were applied such as Taylor matrix method [2], Chebyshev finite difference method [8], Legendre tau method [4], Bessel matrix method [5], and variational iteration method [9].

In this paper, we propose an efficient method, namely homotopy analysis method, to obtain the numerical solution of th-order linear Fredholm integrodifferential-difference equation with variable coefficients. Consider

with mixed conditions

where , and are functions defined on ; the real coefficients , and are appropriate constants.

Homotopy analysis method (HAM) was first introduced by Liao [10] to obtain series solutions of various linear and nonlinear problems. HAM is a promising method that gives us acceptable analytical results with convenient convergence [10]. In contrast to perturbation techniques, this approach is independent of any small parameters. In addition, different from all other analytic techniques, HAM provides us with a simple procedure to obtain the convergence of series of solutions, so that one can obtain accurate enough approximations by auxiliary convergence-controller parameter . In his papers, Liao used this method to solve many nonlinear problems. In [10], especially, he pointed out the basic ideas of the HAM and gave the details in theory and he has successfully applied it to many nonlinear problems [10–15]. Recently, this technique has successfully been applied to several nonlinear problems, such as the viscous flows of non-Newtonian fluids [16, 17], nonlinear heat transfer [18], nonlinear Fredholm integral equations [19], the KdV-type equations [20], differential-difference equations [21], time-dependent Emden-Fowler-type equations [22], Laplace equation with Dirichlet and Neumann conditions [23], and pantograph equations [24].

2. Homotopy Analysis Method

The purpose of the study is to apply the homotopy analysis method to the integrodifferential-difference equation. Now suppose that denotes an initial approximation guess of exact solution of . The auxiliary parameter , which is a convergence-controller parameter, the auxiliary function , and the auxiliary linear operator play important roles within the homotopy analysis method to adjust and control the convergence region of solution series. Liao constructs, using as an embedding parameter, the so-called zero-order deformation equation

where is unknown function to be determined and is given by

When , the zero-order deformation equation (3) becomes , and when , then the zero-order deformation equation (3) becomes . As increases from 0 to 1, the solution varies from the initial guess to the solution . So is exactly the solution of nonlinear equation (1). Expanding in Taylor series with respect to q, we have

where

If the auxiliary linear operator, the initial guess, the auxiliary parameter which will be determined, later and the auxiliary function are properly chosen, the power series of converges at , then the following series solution can be obtained:

where the terms can be determined by the so-called high-order deformation equations. Now the vector is defined as , and differentiating (3) equation times with respect to embedding parameter q, then setting , and finally dividing by , we obtain the so-called mth-order deformation equation in the following form:

where

In order to obey both the first rule of solution expression and the rule of the coefficient ergodicity, the corresponding auxiliary function is determined uniquely by [10].

For any given nonlinear operator N and the term can be easily expressed by (10). So, we can obtain by means of solving the linear high-order deformation equation (8). The mth-order approximation of is given by

The foregoing approximate solution consists of , which is a cornerstone of the HAM in determining convergence of series solution rapidly. We may adjust and control the convergence region and rate of the solution series (11) by means of the auxiliary parameter . To obtain valid region of , we first plot the so-called curves of , , and , where and so on. According to these curves, it is easy to discover the valid region of , which corresponds to the line segments nearly parallel to the horizontal axis.

Theorem 1 (convergence theorem). As long as the series (8) converges to , where is governed by the high-order deformation equation (9) under the definitions (10) and (11), it must be the exact solution of (3).

3. Numerical Examples and Error Analysis

In the present paper, the HAM technique is applied to the integrodifferential-difference equations, and how one can control the convergence of approximate solution and make the convergence fast is shown. The method is applied to the following problems.

Example 2. Firstly, let us consider the second-order integrodifferential-difference equation with variable coefficients [4]

with the exact solution and initial conditions

To solve (12) by means of homotopy analysis method, we choose as initial approximation which satisfies (13) and let us define as follows:

Hence, the high-order deformation equation (8)–(10) is as follows:

Therefore, by starting with , we successively obtain , by (15),

In this way we can get

by choosing proper convergence parameter with . To find the valid region of on the convergence of as goes to infinity, we plot the curves of and as shown in Figure 1 which clearly indicates that the valid region of is about −0,6 < < −0,15.

Figure 1: The curves of and .

Some numerical results of , , and for = −0,368 are presented in Table 1, and the absolute errors of , , and for = −0,368 are given in Table 2. We also give the graphs of the exact solution and 20th order approximation for = −0,368 in Figure 2 and the absolute error of the 20th order approximation for = −0,368 in Figure 3.

Table 1: Comparison of the exact solution and approximations of (12).

−1.0
2.7182818
2.7182784
2.7182076
2.7163599

−0.8
2.2255409
2.2255388
1.8221033
1.8216473

−0.6
1.8221188
1.8221182
1.8221033
1.8216473

−0.4
1.4918246
1.4918250
1.4918242
1.4916964

−0.2
1.2214027
1.2214031
1.2214053
1.2213995

1.0000000
1.0000000
1.0000000
1.0000000

Table 2: The absolute errors approximations of (12).

−1.0

0.000074218443
0.0019219267

−0.8

0.000040710174
0.0010634492

−0.6

0.000015492072
0.00047149553

−0.4

0.00012819928

−0.2

0
0
0

Figure 2: Exact solution and 20th order approximation for = −0,368.

Figure 3: Absolute error of the 20th order approximation for = −0,368.

Example 3. Now, consider the second-order integrodifferential-difference equation with variable coefficients [1]

with exact solution and subjects to initial conditions

In this case, we determine the high-order deformation equation (8)–(10) for (18) and (19)

Choosing as initial approximation which satisfies (19), we successively obtain , by (20), and we can get

The curves of and have been shown in Figure 4, and the valid region of is about −0,7 < < −0,2.

Figure 4: The curves of and .

Finally, we have solutions of (15) and (18) for = −0,37 as shown in Table 3, and the absolute errors of , , and for = −0,37 are given in Table 4. We also give the graphs of the exact solution and 20th order approximation for = −0,37 in Figure 5 and the absolute error of the 20th order approximation for = −0,37 in Figure 6.

Table 3: Comparison of the exact solution and approximations of (18).

0.3678794
0.3678819
0.3678847
0.3679715

0.4493289
0.4493295
0.4493312
0.4493805

0.5488116
0.5488117
0.5488126
0.5488364

0.6703200
0.6703200
0.6703203
0.6703286

−0.2
0.8187307
0.8187307
0.8187308
0.8187320

1.0000000
1.0000000
1.0000000
1.0000000

1.2214027
1.2214027
1.2214026
1.2214008

1.4918246
1.4918246
1.4918236
1.4918054

1.8221188
1.8221185
1.8221138
1.8220361

2.2255409
2.2255399
2.2255246
2.2252881

2.7182818
2.7182787
2.7182364
2.7176359

Table 4: The absolute errors approximations of (18).

−1.0

−0.2

0
0
0

Figure 5: Exact solution and 20th order approximation for = −0,37.

Figure 6: Absolute error of the 20th order approximation for = −0,37.

Example 4. Finally, consider the third-order nonlinear integrodifferential-difference equation [9],

subject to initial conditions

with exact solution . Now, the high-order deformation equation (8)–(10) for (22) and (23) is

We choose as initial guess which satisfies (23) and we get

In Figure 7, the curves of and have been plotted, and the valid region of is determined as −0,7 < < 0.

Figure 7: The graphic of and .

Then, we also have solutions of (22) for = −0,64 as shown in Table 5, and the absolute errors of , and , for = −0,64 are given in Table 6. We also give the graphs of the exact solution and 20th order approximation for = −0,64 in Figure 8 and the absolute error of the 20th order approximation for = −0,64 in Figure 9.

Table 5: Comparison of the exact solution and approximations of (22).

−1.0

1.4918246
1.4908049
1.4896912
1.4882483

1.3498588
1.3494086
1.3489259
1.3482953

1.2214027
1.2212638
1.2211163
1.2209227

−0.2
1.1051709
1.1051529
1.1051338
1.1051088

1.0000000
1.0000000
1.0000000
1.0000000

Table 6: The absolute errors approximations of (22).

−1.0

−0.2

0
0
0

Figure 8: Exact solution and 20th order approximation for = −0,64.

Figure 9: Absolute error of the 20th order approximation for = −0,64.

4. Conclusions

In this work, a reliable algorithm based on the HAM to solve various integrodifferential-difference equations with constant and variable coefficients is presented. Some examples are given to illustrate the validity and accuracy of this procedure. The series solutions of (12), (18), and (22) obtained by HAM contain the auxiliary parameter . In general, by means of the so-called curve, it is straightforward to choose a proper value of which ensures that the series solution is convergent. Figures 1, 4, and 7 show the curves obtained from the mth-order HAM approximation solutions. From these figures, the valid regions of correspond to the line segments nearly parallel to the horizontal axis. Finally, it can be concluded that the homotopy analysis method is a promising tool for both linear and nonlinear IDDEs.

References

M. Gülsu and Y. Öztürk, “A new collocation method for solution of mixed linear integro-differential-difference equations,” Applied Mathematics and Computation, vol. 216, no. 7, pp. 2183–2198, 2010.

M. Gülsu and M. Sezer, “Approximations to the solution of linear Fredholm integrodifferential-difference equation of high order,” Journal of the Franklin Institute, vol. 343, no. 7, pp. 720–737, 2006.

K. Maleknejad and Y. Mahmoudi, “Numerical solution of linear Fredholm integral equation by using hybrid Taylor and block-pulse functions,” Applied Mathematics and Computation, vol. 149, no. 3, pp. 799–806, 2004.

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.

N. şahin, S. Yüzbaşı, and M. Sezer, “A Bessel polynomial approach for solving general linear Fredholm integro-differential-difference equations,” International Journal of Computer Mathematics, vol. 88, no. 14, pp. 3093–3111, 2011.

W. Wang, “An algorithm for solving the high-order nonlinear Volterra-Fredholm integro-differential equation with mechanization,” Applied Mathematics and Computation, vol. 172, no. 1, pp. 1–23, 2006.

W. Wang and C. Lin, “A new algorithm for integral of trigonometric functions with mechanization,” Applied Mathematics and Computation, vol. 164, no. 1, pp. 71–82, 2005.

M. Dehghan and A. Saadatmandi, “Chebyshev finite difference method for Fredholm integro-differential equation,” International Journal of Computer Mathematics, vol. 85, no. 1, pp. 123–130, 2008.

J. Biazar and M. Gholami Porshokouhi, “Application of variational iteration method for linear and nonlinear integrodifferential-difference equations,” International Mathematical Forum, vol. 5, no. 65–68, pp. 3335–3341, 2010.

S. Liao, Beyond Perturbation: Introduction to the Homotopy Analysis Method, vol. 2 of Modern Mechanics and Mathematics, CRC press, Boca Raton, Fla, USA, 2004.

S. J. Liao, Proposed homotopy analysis techniques for the solution of nonlinear problems [Ph.D. thesis], Dissertion, Shanghai Jiao Tong University, 1992.

S. J. Liao, “Homotopy analysis method: a new analytic method for nonlinear problems,” Applied Mathematics and Mechanics, vol. 19, no. 10, pp. 885–890, 1998.

S. J. Liao, “A uniformly valid analytic solution of two-dimensional viscous flow over a semi-infinite flat plate,” Journal of Fluid Mechanics, vol. 385, pp. 101–128, 1999.

S. Liao, “An explicit analytic solution to the Thomas-Fermi equation,” Applied Mathematics and Computation, vol. 144, no. 2-3, pp. 495–506, 2003.

S. Liao, “On the homotopy analysis method for nonlinear problems,” Applied Mathematics and Computation, vol. 147, no. 2, pp. 499–513, 2004.

T. Hayat, T. Javed, and M. Sajid, “Analytic solution for rotating flow and heat transfer analysis of a third-grade fluid,” Acta Mechanica, vol. 191, no. 3-4, pp. 219–229, 2007.

T. Hayat, M. Khan, M. Sajid, and S. Asghar, “Rotating flow of a third grade fluid in a porous space with hall current,” Nonlinear Dynamics, vol. 49, no. 1-2, pp. 83–91, 2007.

S. Abbasbandy, “Homotopy analysis method for heat radiation equations,” International Communications in Heat and Mass Transfer, vol. 34, no. 3, pp. 380–387, 2007.

S. Abbasbandy and E. Shivanian, “A new analytical technique to solve Fredholm's integral equations,” Numerical Algorithms, vol. 56, no. 1, pp. 27–43, 2011.

S. Abbasbandy and F. S. Zakaria, “Soliton solutions for the 5th-order KdV equation with the homotopy analysis method,” Nonlinear Dynamics, vol. 51, no. 1-2, pp. 83–87, 2008.

Z. Wang, L. Zou, and H. Zhang, “Applying homotopy analysis method for solving differential difference equation,” Physics Letters A, vol. 369, no. 1-2, pp. 77–84, 2007.

A. S. Bataineh, M. S. M. Noorani, and I. Hashim, “Solutions of time-dependent Emden-Fowler type equations by homotopy analysis method,” Physics Letters A, vol. 371, no. 1-2, pp. 72–82, 2007.

M. Inc, “On exact solution of laplace equation with dirichlet and neumann boundary conditions by the homotopy analysis method,” Physics Letters A, vol. 365, no. 5-6, pp. 412–415, 2007.

F. Awawdeh, A. Adawi, and S. Al-Shara', “Analytic solution of multipantograph equation,” Journal of Applied Mathematics and Decision Sciences, vol. 2008, Article ID 605064, 11 pages, 2008.