Analytical and Numerical Approaches for Complicated Nonlinear Equations
View this Special IssueResearch Article  Open Access
Bezier Curves Method for FourthOrder Integrodifferential Equations
Abstract
The Bezier curves method is applied to solve both linear and nonlinear BVPs for fourthorder integrodifferential equations. Also, the presented method is developed for solving BVPs which arise from the problems in calculus of variation. These BVPs result from the EulerLagrange equations which are the necessary conditions of the extremums of problems in calculus of variation. Some numerical examples demonstrate the validity and applicability of the technique.
1. Introduction
Recently, there has been much attention devoted to the search for reliable and more efficient solution methods for equations modelling physical phenomena in various fields of engineering (see [1, 2]). One of the methods which has received much concern is the Adomian decomposition method (ADM) (see [3, 4]). The ADM has been employed to solve various scientific models. In [5], Wazwaz's main objective was to obtain the exact solutions to two fourthorder integrodifferential equations. Hashim [4] determined the accuracy and efficiency of the ADM in solving integrodifferential equations. In the present work, we suggest a technique similar to the one which was used in [6–8] for solving both linear and nonlinear boundary value problems (BVPs) for fourthorder integrodifferential equations.
Now, we consider the following class of twopoint BVPs for fourthorder integrodifferential equations where is a real nonlinear continuous function, , , , , and are real constants, and , , and are given and can be approximated by Taylor polynomials. The conditions for existence and uniqueness of solutions of (1) are given in [9].
The rest of this paper is organized as follows. In Section 2, we review the Beizer curves method. Several illustrative examples are provided in Section 3 for confirming the effectiveness of the presented method. Section 4 contains some conclusions and notations about the future works.
2. The Beizer Curves Method
Consider the problem (1). Divide the interval into a set of grid points such that where and is a positive integer. Let for . Then, for , the problem (1) can be decomposed to the following problems: where is considered in . Let where is the characteristic function of for . It is trivial that .
Our strategy is using Bezier curves to approximate the solutions by where is given below. Individual Bezier curves that are defined over the subintervals are joined together to form the Bezier spline curves. For , define the Bezier polynomials of degree that approximate the action of over the interval as follows: where is the Bernstein polynomial of degree over the interval and is the control points (see [7]). By substituting (4) in (3), one can define for as Let where is the characteristic function of for . Besides the boundary conditions on , at each node, we need to impose continuity condition on each successive pair of to guarantee the smoothness. Since the differential equation is of first order, the continuity of (or ) and its first derivative gives where is the th derivative with respect to at . Thus, the vector of control points () must satisfy (see [7]) According to the definition of the , we get that . Therefore Ghomanjani et al. [7] proved the convergence of this method where .
Now, the residual function can be defined in as follows: where is the Euclidean norm. Our aim is solving the following problem over : The mathematical programming problem (11) can be solved by many subroutine algorithms. Here, Maple can solve this optimization problem.
Remark 1. Consider the following functional (see [10]):
to find the extreme value of , the boundary points of the admissible curves are known as
To extremize the necessary condition, is that it should satisfy the EulerLagrange equations
with boundary conditions given in (13). The system of BVPs (14) does not always have a solution, and if the solution exists, it cannot be unique. Note that in many variational problems, the existence of a solution is obvious from the geometrical or physical meaning of the problem, and if the solution of the Euler equation satisfies the boundary conditions, it is unique. Also this unique extremal will be the solution of the given variational problem. Thus, another approach for solving the variational problem (12) is finding the solution of the system of ordinary differential equations (ODEs) (14) which satisfies the boundary conditions in (13) which were called systems of BVPs. The simplest form of the variational problem (12) is
with the given boundary conditions
For the extremum of the functional (15), the necessary condition is to satisfy the following secondorder differential equation:
with boundary conditions given in (16).
One can emphasize that our method is to solve BVPs such as (14) and (17).
3. Numerical Results and Discussion
To demonstrate the accuracy of the decomposition method, we consider some examples with known exact solutions.
Example 1. First, we consider the linear fourthorder integrodifferential equation as in (1) with , , , and ; that is, where the exact solution is with and , and one can find the following approximated solution where the maximum absolute errors in [4] and presented method are, respectively, and . The graphs of approximated and exact solution are shown in Figure 1.
Example 2. Now, consider the nonlinear fourthorder BVP (1) with , , , and , and , where the exact solution is with and , and one can find the following approximated solution: where the maximum absolute errors in [4] and presented method are, respectively, and . The graphs of approximated and exact solution are shown in Figure 2.
Example 3. Consider the problem of finding the minimum of the integral with the boundary conditions where the exact solution is According to (17), the associated EulerLagrange equation is as follows: with and , and the computational results are shown in Table 1, where the computed values are compared with the values obtained from the analytical solution and Legendre polynomials in [11].

4. Conclusions
In this sequel, the Bezier curves method was employed to solve linear and nonlinear BVPs for fourthorder integrodifferential equations. The presented algorithm produced results which are of reasonable accuracy. Numerical examples show that the proposed method is efficient and very easy to use.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The authors are very grateful to the referees for their valuable suggestions and comments that improved the paper, and they would like to express their sincere gratitude to Doctor Y. Damchi. The third author gratefully acknowledges the partial support from University Putra Malaysia.
References
 A. Atangana and D. Baleanu, “Nonlinear fractional JaulentMiodek and WhithamBroerKaup equations within Sumudu transform,” Abstract and Applied Analysis, vol. 2013, Article ID 160681, 8 pages, 2013. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 A. Atangana and A. Secer, “The timefractional coupledKortewegdeVries equations,” Abstract and Applied Analysis, vol. 2013, Article ID 947986, 8 pages, 2013. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 G. Adomian, Solving Frontier Problems of Physics: The Decomposition Method, vol. 60 of Fundamental Theories of Physics, Kluwer Academic Publishers Group, Dordrecht, The Netherlands, 1994, With a preface by Yves Cherruault. View at: MathSciNet
 I. Hashim, “Adomian decomposition method for solving BVPs for fourthorder integrodifferential equations,” Journal of Computational and Applied Mathematics, vol. 193, no. 2, pp. 658–664, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 A.M. Wazwaz, “A reliable algorithm for solving boundary value problems for higherorder integrodifferentiable equations,” Applied Mathematics and Computation, vol. 118, no. 23, pp. 327–342, 2001. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 F. Ghomanjani and M. H. Farahi, “The Bezier control points method for solving delay differential equation,” Intelligent Control and Automation, vol. 3, no. 2, pp. 188–196, 2012. View at: Google Scholar
 F. Ghomanjani, M. H. Farahi, and M. Gachpazan, “Bézier control points method to solve constrained quadratic optimal control of time varying linear systems,” Computational & Applied Mathematics, vol. 31, no. 3, pp. 433–456, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 F. Ghomanjani, A. Kılıçman, and S. Effati, “Numerical solution for IVP in Volterra type linear integrodifferential equations system,” Abstract and Applied Analysis, vol. 2013, Article ID 490689, 4 pages, 2013. View at: Publisher Site  Google Scholar
 R. P. Agarwal, Boundary Value Problems for Higher Order Differential Equations, World Scientific, Singapore, 1986. View at: MathSciNet
 E. Tohidi and A. Kılıçman, “Collocation method based on the bernoulli operational matrix for solving nonlinear BVPs which arise from the problems in calculus of variation,” Mathematical Problems in Engineering, vol. 2013, Article ID 757206, 9 pages, 2013. View at: Publisher Site  Google Scholar
 R. Y. Chang and M. L. Wang, “Shifted Legendre direct method for variational problems,” Journal of Optimization Theory and Applications, vol. 39, no. 2, pp. 299–307, 1983. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
Copyright
Copyright © 2013 F. Ghomanjani 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.