Journal of Applied Mathematics

Volume 2014 (2014), Article ID 678989, 10 pages

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

## Variational Iteration Method for Volterra Functional Integrodifferential Equations with Vanishing Linear Delays

Department of Mathematics, Faculty of Arts and Sciences, Celal Bayar University, Muradiye Campus, Yunus Emre, 45140 Manisa, Turkey

Received 30 April 2014; Revised 29 June 2014; Accepted 29 June 2014; Published 15 July 2014

Academic Editor: Luigi Muglia

Copyright © 2014 Ali Konuralp and H. Hilmi Sorkun. 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

Application process of variational iteration method is presented in order to solve the Volterra functional integrodifferential equations which have multi terms and vanishing delays where the delay function vanishes inside the integral limits such that for , . Either the approximate solutions that are converging to the exact solutions or the exact solutions of three test problems are obtained by using this presented process. The numerical solutions and the absolute errors are shown in figures and tables.

#### 1. Introduction

Nowadays, understanding also from the work of Brunner [1], we are faced with some important problems including the numerical analysis of Volterra functional equations with vanishing delays. He exposed open problems about numerical analysis of th-order Volterra functional integrodifferential equations (VFIDE): where for , as well as the multidelay pantograph-type VFIDE in [1]. During 1990s, Brunner et al. [2] and Hu [3] introduced geometric mesh concept in collocation methods in order to obtain the collocation solutions of the problems. On the other hand, the pantograph differential equations are employed for their numerical solutions by using various methods such as Taylor matrix method [4], variational iteration method [5, 6], differential transform method [7], and methods in other papers [8–12]. In [5], the process of VIM is given for the first order multipantograph equations. In [6], the variational iteration method is applied to some examples in order to obtain the numerical or exact solutions of the multipantograph equation where the coefficient of is a constant and additionally the Lagrange multiplier is given for (1) without terms (2) and (3), that is, only for the extended multipantograph equation.

The variational iteration method which obtains the analytical or numerical solutions of a wide spectrum of differential equations, as well as integral equations, was proposed in the late 90s by He [13–15] and has been used in hundreds of papers by many authors in order to solve the well-known famous equations and to show the effectiveness, straightness, and convergence of that powerful method [16–20]. Furthermore, since it is a useful mathematical tool that ensures its reliability, the method has been developed according to the needs [21–24] and has also been extended for fractional differential and fractional integrodifferential equations [25–28].

The basic idea behind variational iteration method is to construct an iteration formula for the considered equation. After finding optimized Lagrange multiplier for constructed correction functional (i.e., iteration formula), the method takes into account that corrected iteration formula and starts to iterate with an initial function. In most cases the method provides exact solutions or the series form of exact solutions.

In addition to [28], in this paper, the procedure of the variational iteration method is presented for the Volterra functional integrodifferential equations with vanishing delays (1), where the Volterra integral terms are as in (2) and the delayed Volterra integral terms are as in (3); then this extended scheme is applied to three test problems for showing the applicability of the procedure and the convergent numerical solutions to the exact solutions. The numerical data for different parameter values are also given by tables and figures.

#### 2. The Application Process of Variational Iteration Method to VFIDE

Now the process how to apply the method for (1) is given. Firstly, for the basic idea of the method, let us consider the following nonlinear equation: where is a linear operator, is a nonlinear operator, and is a known analytical function. According to variational iteration method, we can construct the following correction functional: where is a general Lagrange multiplier which can be identified by variational theory, is an initial approximation with possible unknowns, and is considered as restricted variation, that is, [14].

Now the above idea can be extended as follows. Considering th-order Volterra functional integrodifferential equation (1) with (2) and (3), the correction functional according to relation (5) can be written as where except the name of the variable involved for and are as in (2) and (3), respectively. In order to specify the iteration, the Lagrange multiplier has to be found. The form of the Lagrange multiplier that will be determined by transferring the derivation from to can be either a linear or a nonlinear function of and . Supposing all the functions on the right hand side of the th-order derivation in (6) as a function such that the correction functional (6) can be written as for , . Making the above correction functional stationary and noticing that which represents all variables to be restricted, we obtain In order to find the conditions on , highest order derivative appearing in the integrand in (10) is transferred from to and so (10) becomes After repeating this process times, we will eventually obtain an expression which () can be decomposed from and their coefficients yield the stationary conditions. Thus, as it is indicated by using the analogue way in [22], the Lagrange multiplier is calculated as the polynomial type Therefore, on account of (12) the iteration functional is where , , , and .

Nevertheless, for the correction functional becomes and Lagrange multiplier is evaluated as from (12). But it can also be chosen as a nonlinear function in order to accelerate the convergent rate of correction functional (14) so that the function can be modified as by excepting the function with its coefficients. Now we have Making (15) stationary and noticing that , it is obviously obtained that and then we have the following stationary conditions: which are the same conditions obtained in [5], so the Lagrange multiplier is found: Accordingly, the correction functional for is By using formulae (13) and (19), with the given initial term , the sequence is identified and the th term of the sequence should be the approximate solution of problem (1).

#### 3. Numerical Examples

In this section we show how the method can be applied to such problems so we give some examples that are the modifications of (1) for the numerical verification of the presented method in Section 2.

*Example 1. *Firstly, we have
where the functions are considered as , , , , , , in (1) with (2) and (3).

In order to solve the problem (20) by means of VIM, we use the proposed procedure in Section 2. Considering the nonlinear Lagrange multiplier (18), since , the Lagrange multiplier is directly calculated so that the correction functional is Starting with and since the initial function is , from iteration formula (21) The other terms of the sequences can be found by using the iteration formula with the previous terms and for the large values of , because the exact solution of (20) is , for have more terms than previous ones and are not necessary to write here, but we can give these limitations It is obviously seen that using this iteration formula (21), the approximate solution of the problem (20) that is convergent to the exact solution even for the large values of is found in the beginning terms of the sequence . Tables 1 and 2 show the values of solutions for comparison purposes and Figures 1 and 2 support the efficiency and the accuracy of the method.

*Example 2. *Now, we have
where the functions are considered as , , , , , , in (1) with (2) and (3).

To solve the problem (24) by means of VIM, we consider the linear Lagrange multiplier evaluated from (12) for this problem. Since , the Lagrange multiplier is directly calculated so that the correction functional is as in (13). To be more accurate in finding the solution, it is obvious to start the initial function as a polynomial type of order two because of the structure of the equation in (24). Starting with and since the initial function is , from iteration formula (25) is found. From the condition , so (25) is and because this must correspond to the initial polynomial function, , . Thus the first iteration solution of the problem (24) is (27) with substituting , which is the exact solution.

*Example 3. *Finally, we have
where the functions are considered as , , , for , , , in (1) with (2) and (3).

This example is different from other two examples and now it is not important which formula of the Lagrange multiplier will be used, because the coefficient is zero and both of the linear form (12) and nonlinear form (18) give the desired multiplier as so that the correction functional is Starting with and since the initial function is , from iteration formula (29) are obtained. The other terms of the sequences can be found by using the iteration formula with substituting the previous terms. The fourth iteration solution coincides with the exact solution of (28), and the approximate solution is changed infinitesimally by the parameter () in the . For example, for the fourth iteration solution is So it is clearly seen from Figure 3 that the fourth iteration solution is the approximate solution with the errors indicated in Tables 3 and 4.

#### 4. Conclusion

In this study, the process of variational iteration method for the Volterra functional integrodifferential equations with vanishing delays (1), where the Volterra integral terms are as in (2), the delayed Volterra integral terms are as in (3), and for , is the linear delay function, is constructed and it is applied to the problems that are the different types of problem (1). In Section 2, two types of Lagrange multiplier are given, that is, linear one and nonlinear one. From also the previous papers [19, 26], it is understood that sometimes the nonlinear multiplier yields the more accurate approach than the linear one. The method is applicable also in the pantograph-type differential equations and Volterra integrodifferential equations with linear delay functions.

#### Conflict of Interests

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

#### References

- H. Brunner, “Current work and open problems in the numerical analysis of Volterra functional equations with vanishing delays,”
*Frontiers of Mathematics in China*, vol. 4, no. 1, pp. 3–22, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - H. Brunner, Q. Hu, and Q. Lin, “Geometric meshes in collocation methods for Volterra integral equations with proportional time delays,”
*Seminar für Angewandte mathematik*1999-25, Eidgenössische Technische Hochschule, Zürich, Switzerland, 1999. View at Google Scholar - Q. Hu, “Geometric meshes and their application to Volterra integro-differential equations with singularities,”
*IMA Journal of Numerical Analysis*, vol. 18, no. 1, pp. 151–164, 1998. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - M. Sezer and S. Yalçinbaş, “Approximate solution of multi-pantograph equation with variable coefficients,”
*Journal of Computational and Applied Mathematics*, vol. 214, no. 2, pp. 406–416, 2008. View at Publisher · View at Google Scholar · View at MathSciNet - Z.-H. Yu, “Variational iteration method for solving the multi-pantograph delay equation,”
*Physics Letters A*, vol. 372, no. 43, pp. 6475–6479, 2008. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - A. Saadatmandi and M. Dehghan, “Variational iteration method for solving a generalized pantograph equation,”
*Computers & Mathematics with Applications*, vol. 58, no. 11-12, pp. 2190–2196, 2009. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - Y. Keskin, A. Kurnaz, M. E. Kiris, and G. Oturanc, “Approximate solutions of generalized pantograph equations by the differential transform method,”
*International Journal of Nonlinear Sciences and Numerical Simulation*, vol. 8, no. 2, pp. 159–164, 2007. View at Google Scholar · View at Scopus - Y. Liu, “Numerical investigation of the pantograph equation,”
*Applied Numerical Mathematics*, vol. 24, no. 2-3, pp. 309–317, 1997. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - H. Brunner, “The numerical analysis of pantograph-type Volterra functional integral equations,” in
*Proceedings of the HKBU-NSYSU Workshop in Computational and Applied Mathematics*, p. 6, Department of Mathematics, National Sun Yat-sen University, Hong Kong, April 2007. - M. Z. Liu and D. Li, “Properties of analytic solution and numerical solution of multi-pantograph equation,”
*Applied Mathematics and Computation*, vol. 155, no. 3, pp. 853–871, 2004. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - G. Derfel and A. Iserles, “The pantograph equation in the complex plane,”
*Journal of Mathematical Analysis and Applications*, vol. 213, no. 1, pp. 117–132, 1997. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - S. Sedaghat, Y. Ordokhani, and M. Dehghan, “Numerical solution of the delay differential equations of pantograph type via Chebyshev polynomials,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 17, no. 12, pp. 4815–4830, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - J. He, “A new approach to nonlinear partial differential equations,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 2, no. 4, pp. 230–235, 1997. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - J. He, “Approximate solution of nonlinear differential equations with convolution product nonlinearities,”
*Computer Methods in Applied Mechanics and Engineering*, vol. 167, no. 1-2, pp. 69–73, 1998. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - J. He, “Variational iteration method for delay differential equations,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 2, no. 4, pp. 235–236, 1997. View at Publisher · View at Google Scholar · View at Scopus - M. A. Abdou and A. A. Soliman, “Variational iteration method for solving BURger's and coupled BURger's equations,”
*Journal of Computational and Applied Mathematics*, vol. 181, no. 2, pp. 245–251, 2005. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - J. H. He, “Some asymptotic methods for strongly nonlinear equations,”
*International Journal of Modern Physics B*, vol. 20, no. 10, pp. 1141–1199, 2006. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - A.-M. Wazwaz, “The variational iteration method for exact solutions of Laplace equation,”
*Physics Letters A*, vol. 363, no. 4, pp. 260–262, 2007. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - A. Konuralp, “The steady temperature distributions with different types of nonlinearities,”
*Computers & Mathematics with Applications*, vol. 58, no. 11-12, pp. 2152–2159, 2009. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - A. Barari, H. D. Kaliji, M. Ghadim, and G. Domairry, “Non-linear vibration of Euler-Bernoulli beams,”
*Latin American Journal of Solids and Structures*, vol. 8, no. 2, pp. 139–148, 2011. View at Google Scholar · View at Scopus - M. Tatari and M. Dehghan, “On the convergence of He's variational iteration method,”
*Journal of Computational and Applied Mathematics*, vol. 207, no. 1, pp. 121–128, 2007. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - J. H. He and X. H. Wu, “Variational iteration method: new development and applications,”
*Computers & Mathematics with Applications*, vol. 54, no. 7-8, pp. 881–894, 2007. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - J.-H. He, “Variational iteration method—some recent results and new interpretations,”
*Journal of Computational and Applied Mathematics*, vol. 207, no. 1, pp. 3–17, 2007. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - L. Xu, “Variational iteration method for solving integral equations,”
*Computers & Mathematics with Applications*, vol. 54, no. 7-8, pp. 1071–1078, 2007. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - Z. M. Odibat and S. Momani, “Application of variational iteration method to nonlinear differential equations of fractional order,”
*International Journal of Nonlinear Sciences and Numerical Simulation*, vol. 7, no. 1, pp. 27–34, 2006. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - A. Konuralp, Ç. Konuralp, A. Yildirim, and A. Yıldırım, “Numerical solution to the Van der Pol equation with fractional damping,”
*Physica Scripta*, vol. 2009, Article ID 014034, 2009. View at Google Scholar - X.-J. Yang and D. Baleanu, “Fractal heat conduction problem solved by local fractional variation iteration method,”
*Thermal Science*, vol. 17, no. 2, pp. 625–628, 2013. View at Publisher · View at Google Scholar · View at Scopus - A. Konuralp and H. H. Sorkun, “The numerical solutions of volterra functional integro-differential equations having vanishing linear delays,” in
*Proceedings of the 4th International Conference on Mathematical and Computational Applications*, Manisa, Turkey, June 2013.