We solve some higher-order boundary value problems by the optimal homotopy asymptotic method (OHAM). The proposed method is capable to handle a wide variety of linear and nonlinear problems effectively. The numerical results given by OHAM are compared with the exact solutions and the solutions obtained by Adomian decomposition (ADM), variational iteration (VIM), homotopy perturbation (HPM), and variational iteration decomposition method (VIDM). The results show that the proposed method is more effective and reliable.

1. Introduction

In this paper, we consider a well-posed š‘›th-order problem of the formš‘¢(š‘›)ī€·=šœ“š‘¢,š‘¢ī…ž,ā€¦,š‘¢(š‘›āˆ’1)ī€ø+šœ™(š‘Ÿ),š‘Ž<š‘Ÿ<š‘,(1.1) with boundary conditions: š‘¢(š‘˜)(š‘Ž)=š›¼š‘– and š‘¢(š‘˜)(š‘)=š›½š‘–, where š‘˜(<š‘›) is a nonnegative integer, š›¼š‘– and š›½š‘– are real finite constants, and šœ™(š‘Ÿ) is a continuous function on [š‘Ž,š‘].

Such types of problems have been investigated by many authors [1, 2] due to their mathematical importance and the potential for applications in hydrodynamic and hydromagnetic stability. Fifth-order boundary value problems arise in the mathematical modeling of viscoelastic flows. Sixth- and eighth-order differential equation govern physics of some hydrodynamic stability problems. When an infinite horizontal layer of fluid is heated from below and is subject to the action of rotation, instability sets in. When this instability is as ordinary convection, the ordinary differential equation is sixth order; when the instability sets in as overstability, it is modeled by an eighth-order ordinary differential equation. If an infinite horizontal layer of fluid is heated from below, with the supposition that a uniform magnetic field is also applied across the fluid in the same direction as gravity and the fluid is subject to the action of rotation, instability sets in. When instability sets in as ordinary convection, it is modeled by tenth-order boundary value problem.

So for the solution of these problems is concerned, many methods appeared in literature. The recent analytic methods are Adomian decomposition method (ADM) [3ā€“5], variational iteration method (VIM) [6], homotopy perturbation method (HPM) [7ā€“9], homotopy analysis method (HAM) [10, 11], differential transform method (DTM) [12], and so forth. Classical perturbation methods are based on the assumptions of small or large parameters, and they cannot produce a general form of an approximate solution. The nonperturbation methods like ADM and DTM can deal strongly with nonlinear problems, but the convergence region of their series solution is generally small. The HPM, which is an elegant combination of homotopy and perturbation technique, overcomes the restrictions of small or large parameters in the problems. It deals with a wide variety of nonlinear problems effectively. Recently, Marinca et al. [13ā€“17] introduced OHAM for approximate solution of nonlinear problems of thin film flow of a fourth-grade fluid down a vertical cylinder. In their work, they have used this method to understand the behavior of nonlinear mechanical vibration of electrical machine. They also used the same method for the solution of nonlinear equations arising in the steady-state flow of a fourth-grade fluid past a porous plate and for the solution of nonlinear equation arising in heat transfer. This method is straight forward, reliable, and explicitly defined. It provides a convenient way to control the convergence of the series solution and allows adjustment of convergence region where it is needed.

Fifth- and sixth-order linear and nonlinear problems were solved by Wazwaz [18, 19], while using decomposition method. Noor et al. [20ā€“25] investigated these type of problems using variational iteration method (VIM), homotopy perturbation method (HPM), and variational iteration decomposition method (VIDM). Modified variational iteration method (MVIM) and iterative method (ITM) were used by Mohyud-Din et al. [26, 27] for such type of problems. Kasi Viswanadham and Murali Krishna [28] used Quintic B-Spline Galerkin method for fifth-order boundary value problems. Siraj-ul-Islam et al. [29, 30] used numerical scheme for the solution of fifth- and sixth-order boundary value problems.

Recently, Ali et al. [31, 32] used OHAM for the solution of multipoint boundary value problems and twelfth-order boundary value problems. We use OHAM to find the approximate analytic solution of some higher-order BVPs. The results of OHAM are compared with those of exact solution, and the errors are compared with the existing results. This paper is organized as follows: Section 2 is devoted to the analysis of the proposed method. Some numerical examples are presented in Section 3. In Sectionā€‰ā€‰4, we concluded by discussing results of the numerical simulation using Mathematica.

2. Method Analysis for Two-Point Boundary Value Problems

Consider the differential equation ā„’š‘¢=š’©š‘¢+šœ™,(2.1) along with boundary conditions: ā„¬ī‚€š‘¢,šœ•š‘¢ī‚šœ•š‘Ÿ=0,(2.2) where ā„’ is linear, š’© is a nonlinear, and ā„¬ is a boundary operator. šœ™ is a known function which is continues for š‘Ÿāˆ¶š‘ŸāˆˆĪ©. According to OHAM, we can construct a homotopy defined by[](1āˆ’š‘)ā„’(š‘¢(š‘Ÿ;š‘))āˆ’š’½(š‘)ā„’(š‘¢(š‘Ÿ;š‘))āˆ’š’©(š‘¢(š‘Ÿ;š‘))āˆ’šœ™(š‘Ÿ)=0,(2.3) where š‘āˆˆ[0,1] is an embedding parameter, and š’½(š‘) is a nonzero auxiliary function for š‘ā‰ 0 and ā„Ž(0)=0. Equation (2.3) satisfies ā„’š‘¢=0,forš‘=0,ā„’š‘¢=š’©š‘¢+šœ™,forš‘=1.(2.4) The solution, š‘¢(š‘Ÿ,0)=š‘£0(š‘Ÿ), of ā„’š‘¢=0 traces the solution curve š‘¢(š‘Ÿ) of (2.1), continuously as š‘ approaches to 1, where š‘£0 is the solution of the zeroth-order problem, that will come in the next few lines.

The auxiliary function š’½(š‘) is chosen in the form (it is a commonly used form) š’½(š‘)=š‘šī“š‘–=1š‘š‘–š¶š‘–,(2.5) where š¶š‘–:š‘–=1,2,ā€¦,š‘š are the convergence controlling constants which are to be determined. We will use this function unless otherwise stated. The auxiliary function can be chosen in a variety of ways, as reported by Marinca et al. [13ā€“17]. We will use some other forms of š’½(š‘) as well.

To get an approximate solution, we expand š‘¢(š‘Ÿ,š‘) in Taylorā€™s series about š‘ in the following manner:š‘¢(š‘Ÿ;š‘)=š‘£0(š‘Ÿ)+āˆžī“š‘š=1š‘£š‘šī€·š‘Ÿ,š¶1,š¶2,ā€¦,š¶š‘šī€øš‘š‘š.(2.6) Substituting (2.5) and (2.6) into (2.3) and equating the coefficient of like powers of š‘, we obtain the following linear equations which can be integrated directly.

Zeroth-order problem: ā„’š‘£0ī‚µš‘£=0,ā„¬0,šœ•š‘£0ī‚¶šœ•š‘›=0.(2.7)

First-order problem: ā„’š‘£1=ī€·1+š¶1ī€øā„’š‘£0āˆ’š¶1ī€·š’©0š‘£0ī€øī‚µš‘£āˆ’šœ™,ā„¬1,šœ•š‘£1ī‚¶šœ•š‘›=0.(2.8)

Second-order problem: ā„’š‘£2=ī€·1+š¶1ī€øā„’š‘£1āˆ’š¶1š’©1ī€·š‘£0,š‘£1ī€ø+š¶2ī€·ā„’š‘£0āˆ’š’©0š‘£0ī€øī‚µš‘£āˆ’šœ™,ā„¬2,šœ•š‘£2ī‚¶šœ•š‘›=0.(2.9) Though we can construct higher-order problems easily, solutions upto the second-order problems are enough to produce excellent results.

If the series (2.6) is convergent at š‘=1, then the approximate solution in our case is,Ģƒš‘¢(š‘Ÿ)=š‘£(š‘Ÿ)=š‘£0(š‘Ÿ)+š‘£1ī€·š‘Ÿ,š¶1ī€ø+š‘£2ī€·š‘Ÿ,š¶1,š¶2ī€ø.(2.10) By substituting (2.10) into (2.1), the resulting residual isā„›ī€·š‘Ÿ,š¶1,š¶2ī€ø=ā„’(Ģƒš‘¢(š‘Ÿ))āˆ’š’©(Ģƒš‘¢(š‘Ÿ))āˆ’šœ™(š‘Ÿ).(2.11) If ā„›=0, Ģƒš‘¢ will be the exact solution. Otherwise, we minimize ā„› over domain of the problem. To find the optimal values of š¶š‘–which minimizes ā„›, many methods can be applied [13ā€“17]. We follow two methods: the method of least squares and the Galerekinā€™s method. According to the method of least squares, we first construct the functionalš’„ī€·š¶1,š¶2ī€ø=ī€œš‘š‘Žā„›2š‘‘š‘Ÿ,(2.12) and then minimizing it, we havešœ•š’„šœ•š¶1=šœ•š’„šœ•š¶2=0.(2.13) According to the Galerekinā€™s method, we solve the following system for š¶1 andš¶2:ī€œš‘š‘Žā„›šœ•Ģƒš‘¢šœ•š¶1ī€œš‘‘š‘Ÿ=0,š‘š‘Žā„›šœ•Ģƒš‘¢šœ•š¶2š‘‘š‘Ÿ=0.(2.14) Knowing š¶1 and š¶2, the approximate solution is well determined.

3. Some Numerical Examples

Example 3.1 (fifth-order linear). Consider the following problem: š‘¦(š‘£)(š‘„)=š‘¦āˆ’15š‘’š‘„āˆ’10š‘„š‘’š‘„,0<š‘„<1,(3.1) with boundary conditions š‘¦(0)=0,š‘¦ī…ž(0)=1,š‘¦ī…žī…ž(0)=0,š‘¦(1)=0,š‘¦ī…ž(1)=āˆ’š‘’.(3.2) The exact solution of this problem is š‘¦(š‘„)=š‘„(1āˆ’š‘„)š‘’š‘„.
We choose the auxiliary function as ā„Ž(š‘)=š‘(š¶1+š¶2š‘„). Plugging in this value in (2.3) of Section 2, we obtain the following linear problems which can be integrated directly.
Zeroth-order problem: š‘¦0(5)š‘¦(š‘„)=0,0(0)=0,š‘¦ī…ž0(0)=0,š‘¦0ī…žī…ž(0)=0,š‘¦0(1)=0,š‘¦ī…ž0(1)=āˆ’š‘’.(3.3)
First-order problem: š‘¦1(5)ī€·š‘„,š¶1,š¶2ī€ø=5š‘’š‘„ī€·š¶(3+2š‘„)1+š¶2š‘„ī€øāˆ’ī€·š¶1+š¶2š‘„ī€øš‘¦0š‘¦(š‘„),1(0)=0,š‘¦ī…ž1(0)=0,š‘¦1ī…žī…ž(0)=0,š‘¦1(1)=0,š‘¦ī…ž1(1)=0.(3.4)
Second-order problem: š‘¦2(5)ī€·š‘„,š¶1,š¶2ī€ø=ī€·1+š¶1+š¶2š‘„ī€øš‘¦1(5)ī€·š‘„,š¶1,š¶2ī€øāˆ’ī€·š¶1+š¶2š‘„ī€øš‘¦1ī€·š‘„,š¶1,š¶2ī€ø,š‘¦2(0)=0,š‘¦ī…ž2(0)=0,š‘¦2ī…žī…ž(0)=0,š‘¦2(1)=0,š‘¦ī…ž2(1)=0.(3.5) Adding up the solutions of these problems, the second-order approximate solution, Ģƒš‘¦(š‘„)=š‘¦0(š‘„)+š‘¦1ī€·š‘„,š¶1,š¶2ī€ø+š‘¦2ī€·š‘„,š¶1,š¶2ī€øī€·š‘„+š‘‚15ī€ø,(3.6) is determined by knowing the optimal values of the auxiliary constants, š¶1 and š¶2. Using Galerkinā€™s method, we obtain š¶1=āˆ’1.000245451, š¶2=0.000124615.
By considering these values, (3.6) becomes Ģƒš‘¦(š‘„)=š‘„āˆ’0.5š‘„3āˆ’0.333333š‘„4āˆ’0.125š‘„5āˆ’0.0333333š‘„6āˆ’0.00694444š‘„7āˆ’0.00119049š‘„8āˆ’0.000173601š‘„9āˆ’0.0000220495š‘„10āˆ’2.48013Ɨ10āˆ’6š‘„11āˆ’2.50501Ɨ10āˆ’7š‘„12āˆ’2.27501Ɨ10āˆ’8š‘„13āˆ’2.0326Ɨ10āˆ’9š‘„14ī€·š‘„+š‘‚15ī€ø.(3.7) Numerical results of the solution (3.7) are displayed in Table 1.

Example 3.2 (another fifth-order linear). Consider the following problem: š‘¦(5)(š‘„)+š‘„š‘¦(š‘„)=19š‘„Cos(š‘„)āˆ’2š‘„3Cos(š‘„)āˆ’41Sin(š‘„)+2š‘„2Sin(š‘„),(3.8) with boundary conditions š‘¦š‘¦(āˆ’1)=š‘¦(1)=Cos(1),ī…ž(āˆ’1)=š‘¦ī…žš‘¦(1)=āˆ’4Cos(1)+Sin(1),ī…žī…ž(āˆ’1)=3Cos(1)āˆ’8Sin(1).(3.9) Exact solution of this problem is š‘¦(š‘„)=(2š‘„2āˆ’1)Cos(š‘„).
Considering the second-order solution Ģƒš‘¦(š‘„)=š‘¦0(š‘„)+š‘¦1(š‘„,š¶1)+š‘¦2(š‘„,š¶1,š¶2)+š‘‚(š‘„13), we use the method of least squares to obtain š¶1=0.9940605306, š¶2=āˆ’3.9762851376.
Having these values, our solution in this case is Ģƒš‘¦(š‘„)=āˆ’0.999978+2.49992š‘„2āˆ’1.04155š‘„4+0.0846365š‘„6āˆ’0.00276866š‘„8+0.0000420239š‘„10+3.38286Ɨ10āˆ’7š‘„12ī€·š‘„+š‘‚13ī€ø.(3.10) Numerical results of the solution (3.10) are displayed in Table 2.

Example 3.3 ([33] fifth-order nonlinear). Consider the following problem: š‘¦(š‘£)(š‘„)=š‘¦3(š‘„)š‘’āˆ’š‘„,0<š‘„<1,(3.11) with boundary conditions š‘¦(0)=1,š‘¦ī…ž(0)=1/2,š‘¦ī…žī…ž(0)=1/4,š‘¦(1)=š‘’1/2,š‘¦ī…ž1(1)=2š‘’1/2.(3.12) The exact solution for this problem is š‘¦(š‘„)=š‘’š‘„/2.
We consider the second-order solution, Ģƒš‘¦(š‘„)=š‘¦0(š‘„)+š‘¦1(š‘„,š¶1)+š‘¦2(š‘„,š¶1,š¶2)+š‘‚(š‘„15).
Using Galerkinā€™s procedure in Section 2, we obtain the following values: š¶1=0.010868466,š¶2=āˆ’0.029423113.(3.13) The second-order approximate solution is š‘„Ģƒš‘¦(š‘„)=1+2+š‘„28+0.0205993š‘„3+0.0030533š‘„4+0.0000630671š‘„5+5.2556Ɨ10āˆ’6š‘„6+3.754Ɨ10āˆ’7š‘„7+2.29401Ɨ10āˆ’8š‘„8+1.74892Ɨ10āˆ’9š‘„9āˆ’3.2692Ɨ10āˆ’11š‘„10āˆ’1.48596Ɨ10āˆ’12š‘„11āˆ’6.20243Ɨ10āˆ’14š‘„12āˆ’2.58495Ɨ10āˆ’15š‘„13+2.75287Ɨ10āˆ’17š‘„14ī€·š‘„+š‘‚15ī€ø.(3.14) Numerical results of the solution (3.14) are displayed in Table 3.

Example 3.4 (sixth-order nonlinear). Consider the following problem: š‘¦(š‘£š‘–)(š‘„)=š‘’š‘„š‘¦2(š‘„),0<š‘„<1,(3.15) with boundary conditions š‘¦(0)=1,š‘¦ī…ž(0)=1,š‘¦ī…žī…ž(0)=1,š‘¦(1)=š‘’āˆ’1,š‘¦ī…ž(1)=āˆ’š‘’āˆ’1,š‘¦ī…žī…ž(1)=š‘’āˆ’1.(3.16) The exact solution is š‘¦(š‘„)=š‘’āˆ’š‘„.
For this problem, we take the auxiliary function ā„Ž(š‘)=š‘(š¶1+š¶2š‘’āˆ’š‘„), Ģƒš‘¦(š‘„)=š‘¦0(š‘„)+š‘¦1ī€·š‘„,š¶1,š¶2ī€ø+š‘¦2ī€·š‘„,š¶1,š¶2ī€øī€·š‘„+š‘‚13ī€ø.(3.17) Using Galerkinā€™s method, we obtain š¶1=0.41243798998, š¶2=0.0014069149.
OHAM solution in this case is š‘„Ģƒš‘¦(š‘„)=1āˆ’0.999999994š‘„+22āˆ’0.166666775š‘„3+š‘„424āˆ’0.008332465š‘„5+0.001387441š‘„6āˆ’0.000197784š‘„7+0.000025071š‘„8āˆ’3.002Ɨ10āˆ’6š‘„9+2.918Ɨ10āˆ’7š‘„10āˆ’6.7Ɨ10āˆ’9š‘„11āˆ’2.257Ɨ10āˆ’9š‘„12+1.806Ɨ10āˆ’10š‘„13āˆ’1.974Ɨ10āˆ’11š‘„14ī€·š‘„+š‘‚15ī€ø.(3.18) Numerical results of the solution (3.18) are displayed in Table 4.

Example 3.5 (eighth-order nonlinear). Consider the following problem: š‘¦(š‘£š‘–š‘–š‘–)(š‘„)=š‘’āˆ’š‘„š‘¦2(š‘„),0<š‘„<1,(3.19) with boundary conditions š‘¦ī…ž(0)=1,š‘¦ī…ž(0)=1,š‘¦ī…žī…ž(0)=1,š‘¦ī…žī…žī…ž(0)=1,š‘¦(1)=š‘’,š‘¦ī…ž(1)=š‘’,š‘¦ī…žī…ž(1)=š‘’.(3.20) Considering the second-order solution Ģƒš‘¦(š‘„)=š‘¦0(š‘„)+š‘¦1(š‘„,š¶1)+š‘¦2(š‘„,š¶1,š¶2)+š‘‚(š‘„13), the following values of the convergence controlling constants are obtained by using Galerkinā€™s method: š¶1=āˆ’1.451894673Ɨ10āˆ’11,š¶2=āˆ’0.000647581.(3.21) The approximate solution in this case is Ģƒš‘¦(š‘„)=1+š‘„+0.5š‘„2+0.166667š‘„3+0.0416275š‘„4+0.00884857š‘„5+0.00117027š‘„6+0.00033161š‘„7+1.6061Ɨ10āˆ’8š‘„8+1.78456Ɨ10āˆ’9š‘„9+1.78456Ɨ10āˆ’10š‘„10+1.62233Ɨ10āˆ’11š‘„11+1.3494Ɨ10āˆ’12š‘„12ī€·š‘„+š‘‚13ī€ø.(3.22) If the method of least squares is used to determine š¶ā€™s, we have then š¶1=1.793Ɨ10āˆ’8,š¶2=āˆ’1.001347284.(3.23) The approximate solution in this case is Ģƒš‘¦(š‘„)=1+š‘„+0.5š‘„2+0.166667š‘„3+0.0416667š‘„4+0.00833313š‘„5+0.00138918š‘„6+0.000198233š‘„7+0.000024835š‘„8+2.75944Ɨ10āˆ’6š‘„9+2.75944Ɨ10āˆ’7š‘„10+2.50859Ɨ10āˆ’8š‘„11+2.08656Ɨ10āˆ’9š‘„12ī€·š‘„+š‘‚13ī€ø.(3.24) Let us use the auxiliary function ā„Ž(š‘)=š‘(š¶1+š¶2š‘’āˆ’š‘„) and consider the second-order solution Ģƒš‘¦(š‘„)=š‘¦0(š‘„)+š‘¦1ī€·š‘„,š¶1,š¶2ī€ø+š‘¦2ī€·š‘„,š¶1,š¶2ī€øī€·š‘„+š‘‚13ī€ø.(3.25) Using Galerkinā€™s method, we obtain š¶1=āˆ’0.9993171458, š¶2=āˆ’0.0012314995.
The OHAM solution in this case is š‘„Ģƒš‘¦(š‘„)=1+š‘„+2+š‘„2!33!+0.041666667š‘„4+0.008333333š‘„5+0.001388889š‘„6+1.984Ɨ10āˆ’4š‘„7+2.480Ɨ10āˆ’5š‘„8+2.756Ɨ10āˆ’6š‘„9+2.756Ɨ10āˆ’7š‘„10+2.505Ɨ10āˆ’8š‘„11+2.088Ɨ10āˆ’9š‘„12ī€·š‘„+š‘‚13ī€ø.(3.26) Numerical results of the solutions (3.22), (3.24), and (3.26) are displayed in Table 5.

Example 3.6 (nineth-order linear). Consider the following problem: š‘¦(9)(š‘„)=š‘¦(š‘„)āˆ’9š‘’š‘„,(3.27) with boundary conditions š‘¦(0)=1,š‘¦ī…ž(0)=0,š‘¦ī…žī…ž(0)=āˆ’1,š‘¦ī…žī…žī…ž(0)=āˆ’2,š‘¦ī…žī…žī…žī…ž(0)=āˆ’3,š‘¦(1)=0,š‘¦ī…ž(1)=āˆ’š‘’,š‘¦ī…žī…ž(1)=āˆ’2š‘’,š‘¦ī…žī…žī…ž(1)=āˆ’3š‘’.(3.28) Exact solution is š‘¦(š‘„)=(1āˆ’š‘„)š‘’š‘„.
For this linear problem, we take ā„Ž(š‘)=š‘(š¶1+š¶2š‘„), and according to the rest of the procedure of OHAM, the second-order solution, Ģƒš‘¦(š‘„)=š‘¦0(š‘„)+š‘¦1(š‘„,š¶1,š¶2)+š‘¦2(š‘„,š¶1,š¶2)+š‘‚(š‘„13), is determined by the values of š¶š‘–:š‘–=1,2. Following the Galerkinā€™s method, we obtain š¶1=āˆ’1, š¶2=0, for š‘Ž=0 and š‘=1.
The second-order approximate solution is š‘„š‘¦(š‘„)=1āˆ’22āˆ’š‘„33āˆ’š‘„44āˆ’0.03333333333š‘„5āˆ’0.0069444444š‘„6āˆ’0.0011904762š‘„7āˆ’0.0001736111š‘„8āˆ’0.00002204586š‘„9āˆ’2.48016Ɨ10āˆ’6š‘„10āˆ’2.50521Ɨ10āˆ’7š‘„11āˆ’2.29644Ɨ10āˆ’8š‘„12ī€·š‘„+š‘‚13ī€ø.(3.29) Numerical results of the solution (3.29) are displayed in Table 6.

Example 3.7 (tenth-order nonlinear). Consider the following problem:š‘¦(š‘‹)(š‘„)=š‘’āˆ’š‘„š‘¦2(š‘„),0<š‘„<1,š‘¦(0)=1,š‘¦ī…ž(0)=1,š‘¦ī…žī…ž(0)=1,š‘¦ī…žī…žī…ž(0)=1,š‘¦(š‘–š‘£)(0)=1,š‘¦(1)=š‘’,š‘¦ī…ž(1)=š‘’,š‘¦ī…žī…ž(1)=š‘’,š‘¦ī…žī…žī…ž(1)=š‘’,š‘¦(š‘–š‘£)(1)=š‘’.(3.30) We consider the second-order solution Ģƒš‘¦(š‘„)=š‘¦0(š‘„)+š‘¦1(š‘„,š¶1)+š‘¦2(š‘„,š¶1,š¶2)+š‘‚(š‘„13).
To find the values of š¶š‘–, we apply the Galarkinā€™s method. So solving the system ī€œš‘š‘Žš‘…šœ•Ģƒš‘£šœ•š¶1ī€œš‘‘š‘Ÿ=0,š‘š‘Žš‘…šœ•Ģƒš‘£šœ•š¶2š‘‘š‘Ÿ=0,(3.31) we have š¶1=0, š¶2=āˆ’1.023966086.
In this case, the approximate solution is š‘„Ģƒš‘¦(š‘„)=1+š‘„+22+š‘„36+š‘„424+0.008333323š‘„5+0.00138894š‘„6+0.000198312š‘„7+0.000024898š‘„8+2.712Ɨ10āˆ’6š‘„9+2.8218Ɨ10āˆ’7š‘„10+2.5652Ɨ10āˆ’8š‘„11+2.1377Ɨ10āˆ’9š‘„12ī€·š‘„+š‘‚13ī€ø.(3.32) Numerical results of the solution (3.32) are displayed in Table 7.

4. Conclusions

In this paper, we have used OHAM to find the approximate analytic solution to higher-order two-point boundary value problems in finite domain. It is observed that the method is explicit, effective, and reliable. It works well for higher-order problems and represents the fastest convergence as well as a remarkable low error. The OHAM also provides us with a very simple way to control and adjust the convergence of the series solution using the auxiliary constants š¶š‘–ā€™s which are optimally determined. Furthermore, by using different forms of the auxiliary function, more accuracy can be obtained. It has been also observed that for determining the optimal values of š¶ā€™s, the performance of both the least squares and the Galerkinā€™s method is problem dependent. One can select one of these two which best suits the problem solution.