Laplace transform homotopy perturbation method for the approximation of variational problems.

This article proposes the application of Laplace Transform-Homotopy Perturbation Method and some of its modifications in order to find analytical approximate solutions for the linear and nonlinear differential equations which arise from some variational problems. As case study we will solve four ordinary differential equations, and we will show that the proposed solutions have good accuracy, even we will obtain an exact solution. In the sequel, we will see that the square residual error for the approximate solutions, belongs to the interval [0.001918936920, 0.06334882582], which confirms the accuracy of the proposed methods, taking into account the complexity and difficulty of variational problems.

Background

The calculus of variations is a powerful branch of the analysis with many applications in both pure and practical mathematics (Lanczos 1986). It has also been found that even, the laws of physics can be expressed in a compact and elegant way through variational principles, as occurs with the Lagrange equations of mechanics, which can be deduced from the variational principle of Hamilton. Unlike elementary calculus problems, which seeks to find the points at which a function reaches its maximum and minimum values, the variational calculus considers the problem of some magnitude, whose values depend on a entire curve, throughout an integral (for instance: surface area or descent time). In fact the aim is to find the curve that extremizes the aforementioned quantity in question.

As will be seen in the next section, the procedure that we will follow consist in finding the differential equation for a function, that leads an integral to take an extreme value.

Although some variational problems were solved by using special methods (like Bernoulli’s solution for the brachistochrone problem (see “Basic idea of variational calculus” section) or so-called isoperimetric problem, proposed and solved by the ancient Greeks), it was Euler who presented the variational calculus as a coherent branch of the analysis by discovering the basic differential equation for an extremization curve.

As discussed in the next section, the Euler equation is a nonlinear equation, quite complicated and impossible to solve it in general, and although some applications lead to specific cases which are soluble, it is clear the importance of searching for methods to find approximate solutions to these equations. In fact, nonlinear problems arise in many phenomena, practical or theoretical. For the same reason, several methods focused to find approximate solutions to nonlinear differential equations, algebraic and transcendental equations and integro-differential equations have been reported, such as those based on based on tanh method (Evans and Raslan 2005), exp-function (Xu 2007; Mahmoudi et al. 2008), Adomian’s decomposition method (ADM) (Adomian and Rach 1985, Adomian 1988; Babolian and Biazar 2002; Rach 2012; Fatoorehchi and Abolghasemi 2012, 2013; Fatoorehchi et al. 2014, 2015a, b, c, d; Fatoorehchi and Abolghasemi 2015e; Hashim 2006), parameter expansion (Zhang and Xu 2007), homotopy perturbation method (HPM) (He 1998, 1999, 2003, 2008; Filobello-Nino et al. 2012, 2013c, 2014a, b, 2015b;Ghaderi 2012; Aminikhah and Hemmatnezhad 2012; Aminikhah 2011, 2012; Vazquez-Leal et al. 2012a, b; Fatoorehchi and Abolghasemi 2011), homotopy analysis method (HAM) (Rashidi et al. 2012a, 2012b; Lakestani and Kazemian 2013), series method (Filobello-Nino et al. 2015a), Lie algebra method (Shang 2012, 2013, 2015; Olver 1993; Casas 1996) and perturbation method (PM) (Filobello-Nino et al. 2013a, b) among many others.

A particularly relevant job, which deals with a combined version of the Laplace transform method (LTM) and Adomian’s decomposition method (ADM) (named the H2LTM), it is described in reference (Fatoorehchi and Abolghasemi 2015e). This work introduced a technique for extending the LTM to solve nonlinear differential equations obtaining analytical and accurate solutions.

The rest of this work is organized as follows. “Basic idea of variational calculus” section presents a basic idea of the main results of variational calculus required for this work. In “Standard HPM” section, we introduce the basic idea of standard HPM method. “LT-HPM and some of its modified versions” section presents an introduction of LT-HPM and some of its modified versions useful for this work. Additionally, “Cases study” section presents four cases study of variational problems. In addition a wide discussion on the results is presented in “Discussion” section. Finally, a brief conclusion is given in “Conclusions” section.

Basic idea of variational calculus

This paper deals with the problem of finding a function with continues second derivatives (Elsgolts 1983; Levi 1980; Simmons 1983; Zelikin 2005; Boyce 1998), in order to find an approximate solution to the variational problem of extreme an integral of the formso that its graph join two points, namely and

From variational calculus, we know that if is a function that satisfies the above requirements, then satisfies the Euler equation.

Some interesting cases of the integral (1) are (Simmons 1983; Filobello-Nino et al. 2013d):

The problem of finding the shortest curve joining two given points, requires to minimize:

In the problem of finding the curve that generates the surface of revolution of smallest area when revolves about the axis , we have to minimize:

The problem of brachistochrone, search for the vertical curve without friction that joins two fixed points, through which a particle slides in the shortest possible time (Connor and Robertson 1997, 1998; Filobello-Nino et al. 2013d). This variational problem for the curve of quickest descent, requires minimize:

When the variable y does not appear in the function f, Euler’s equation (2) is reduced immediately to the formintegrating (6), we see that is conserved in according with differential equationwhere c is a constant.

On the other side if f is not an explicit function of x, then Euler’s equation can be expressed in the formwhere

It is customary to refer to the solutions of the Euler equation as stationary functions, and the corresponding value of the integral (1) as stationary value. Besides, the solutions to the above mentioned Euler equation, which are unrestricted by the boundary conditions, are called extremals (Simmons 1983).

Standard HPM

The standard homotopy perturbation method (HPM) was proposed by Ji Huan He, it was introduced to approach various kinds of nonlinear problems. The HPM is considered as a combination of the classical perturbation technique and the homotopy (whose origin is in the topology), but not restricted to small parameters as occur with traditional perturbation methods (He 1998, 1999).

To figure out how HPM works, consider a general nonlinear differential equation in the formwith the following boundary conditionswhere A is a general differential operator, B is a boundary operator, f(x) a known analytical function and Γ is the domain boundary for Ω. A can be divided into two operators L and N, where L is linear and N nonlinear; so that (10) can be rewritten as

Generally, a homotopy can be constructed as (He 1998, 1999).orwhere p is a homotopy parameter, whose values are within range of 0 and 1, is the first approximation for the solution of (12) that satisfies the boundary conditions.

Assuming that solution for (13) or (14) can be written as a power series of as

Substituting (15) into (14) and equating identical powers of terms, there can be found values for the sequence …

When , it yields the approximate solution for (1) in the form

LT-HPM and some of its modified versions

Basic idea of the LT-HPM method

LT-HPM follows the same steps of standard HPM until (14), next we apply Laplace transform on both sides of homotopy equation (14), to obtain (Aminikhah and Hemmatnezhad 2012; Aminikhah 2011, 2012; Filobello-Nino et al. 2013c, 2014a).(employing the freedom of homotopy formulation, it is possible that LT-HPM substitutes for an arbitrary function provided with some unknown parameters to be adequately determined [see cases study 3 and 4)].

Using the differential property (170) of LT, we have (Spiegel 1998):orapplying inverse Laplace transform to both sides of (19), we obtainAssuming that the solutions of (12) can be expressed as a power series of

Then substituting (21) into (20), we getcomparing coefficients of , with the same power leads to

Assuming that the exact solution may be obtained as follows

The above notation has emphasized the fact that the nonlinear operator can be expanded in formal series with respect to the embedding parameter . The same criterion is used throughout “LT-HPM and some of its modified versions” section (Marinca and Herisanu 2011).

Basic idea of non-linearities distribution Laplace transform-homotopy perturbation method (NDLT-HPM)

Vazquez-Leal et al. (2012a) introduced a modified version of homotopy perturbation method, the nonlinearities distribution homotopy perturbation method (NDHPM), which sometimes eases the solutions searching process for (12) and reduces the complexity of solving differential equations in terms of power series.

As first step, we introduced the following homotopy (Vazquez-Leal et al. 2012a).or

It can be noticed that the homotopy function (25) is essentially the same as (13), except for the nonlinear operator and the non homogeneous function , which contain embedded the homotopy parameter . Thus, a way to introducenon-linearities distribution Laplace transform-homotopy perturbation method (NDLT-HPM) is apply directly, Laplace transform on both sides of homotopy equation (36), to obtain (Filobello-Nino et al. 2014b).(more generally, one could substitute in (25) by another function such that, when see cases study above).

Using the differential property of LT, we have (Spiegel 1998).orapplying inverse Laplace transform to both sides of (29), we obtain

Assuming that the solutions of (12) and can be expressed as a power series of andrespectively, where are usually obtained by applying Taylor series method to .

Then substituting (32) and (31) into (30), we getcomparing coefficients of , with the same power leads to

Assuming that the initial approximation has the form: therefore the exact solution may be obtained as follows

Basic idea of Laplace transform homotopy perturbation method with variable coefficients (CVLT-HPM)

To obtain (18), we assumed that the coefficient of is one. In this section we will assume the mentioned coefficient is a positive whole power of , thus we rewrite explicitly (17) for this case as (Filobello-Nino et al. 2015b).where is a positive integer, and employing the versatility and freedom of HPM we have substituted for an arbitrary function . We will choose as a trial function , a polynomial provided with some unknown parameters to be determined [see (44)].

Using the properties (170) and (171), we have (Spiegel 1998).after integrating m times, we obtainor

Assuming also, that the solutions of the ODE to solve, can be expressed as a power series of

Then substituting (40) into (39), we getcomparing coefficients of with the same power leads to

Assuming that then an approximate solution may be obtained as follows

For boundary value problems it is expected that some of the initial conditions above mentioned, be initially unknown; therefore (43) can be expressed as

In order to calculate adequately the values for we will deduce an algebraic system for them, in the following way.

Requiring that (44) satisfies the boundary conditions at the endpoint of the interval.

In order to find the values of the total number of parameters, we require add more algebraic equations to those found in (1), up to get as many equations as parameters to determine. If are required, say j additional equations, then a convenient possibility is incorporating the following j-equations (Marinca and Herisanu 2011). where the residual is defined, by substituting of (44) into (12), to obtain(see 176).

. The above points belong to the interest interval, and it is assumed that is the approximate solution of (12) given by (33).

For the case of second order ODES, where ; it is possible to show the term containing could give rise to an inappropriate term To avoid this potential problem, we should perform, as a work rule, an adequate transformation in order to express our differential equation in its normal form (Simmons 1983). Such as it is known, the normal form does not contain the problematic term , and therefore allows applying VCLT-HPM algorithm in accordance with the above procedure (Filobello-Nino et al. 2015b).

Cases study

This section, applies LT-HPM and some of their modifications aforementioned, in order to find analytical approximate expressions for the stationary functions, solutions of the Euler equations, which turn out from the following cases study.

Case 1

This case study uses basic version of LT-HPM.

We will get an approximate solution for the problem of finding a function that extremizes the value of the integral.

After employing Euler equation (2), for we obtain

As a case study, we consider the following boundary conditions , so that we can express formally our problem in the following way

Unfortunately there is not a systematic way to obtain an exact solution for this boundary value problem.

To obtain an approximate solution for (46), we identify

We construct the following homotopy in accordance with (13)or

Applying Laplace transform to (50) we get

As it is explained in (Spiegel 1998), it is possible to rewrite the above equation, aswhere we have defined .

After applying the initial condition , the last expression can be simplified as followswhere, we defined .

Solving for and applying Laplace inverse transform

Next, we assume the solution for (46) has the formand we chooseas the first approximation for the solution of (46) that satisfies the condition

After substituting (54) and (55) into (53), we get

Equating terms with identical powers of , we obtain

… .

Solving the above equations for we obtain

and so on.

By substituting (60)–(62) into (54) and calculating the limit when , results in a second order approximation

In order to calculate the value of , we require that (63) satisfies the boundary condition , so that we obtain

Thus, substituting (64) into (63), we obtain a handy LT-HPM solution.

Case 2

This case study employs CVLT-HPM method.

We will obtain the exact solution for the problem of finding a function that extremizes the value of the integral.

After employing Euler equation (2), for we obtain

Considering the boundary conditions ,, our purpose will be to get an approximate solution for the following non homogeneous variable coefficients linear equation.

Since it is not known in advance the value of and (67) has a term containing [see “Basic idea of Laplace transform homotopy perturbation method with variable coefficients (CVLT-HPM)” section (2)], then, it is better to write (67) in a simpler way by using the following substitution (Simmons 1983)this transforms an equation in the standard form as (67), in its normal form given bywhere

Besides (69) implies the following boundary conditions for u:since and

Therefore, we can express formally our problem in terms of as follows

In the same way, we note thatimplies the condition

With the purpose of applying CVLT-HPM, we will rewrite (70) as follows

In accordance with “Basic idea of Laplace transform homotopy perturbation method with variable coefficients (CVLT-HPM)” section, it is convenient to construct the following homotopy.where we have chosen as trial function, the constant function .

In order to apply CVLT-HPM, we simplify (76) as follows

After applying Laplace transform to (77) we get

As it is explained in (Spiegel 1998) it is possible to rewrite the above equation, aswhere we have defined and employed the initial conditions , and [see (73) and (74)].

After integrating (78), we obtain

Applying to this last integral expression, yields in

Next, we assume a series solution for , in the form

Substituting (81) into (80), we get

On comparing the coefficients with identical powers of and perform the indicated operations, we have

… .

Solving the above operations for we obtainand so on.

By substituting solutions (88)–(93) into (81) and calculating the limit value , results in.

Clearly, the second term of the right hand side of (94) is a geometric series of ratio . Since , it is well known that the mentioned series converges and the sum of the first n terms is given by (Simmons 1983); therefore (94) becomes in

Assuming that , we get

From the boundary condition we get , therefore

Finally, from (69) we obtain the solution for the boundary value problem (68)

Thus, for this case study CVLT-HPM provides an exact solution for (68).

Case 3

Next, we will analyze the case of

After employing Euler equation (8), we getwhere is a constant.

Later, we will require to express (100) in such a way that it does not involves Therefore we differentiate (100) respect to .

We rewrite (100) as followsorso thatwhere we employed the notation

Next, we propose the following substitution

Thus, (104) adopts the form

We will consider as a case study . It corresponds to the equation (see 100) (Filobello-Nino et al. 2013d).

We note, it is possible to rewrite (108) as

Although (109) can be solved performing the algebraic substitution , the dependent variable remains implicit.

In fact, we propose the following boundary value problem.

In order to get an explicit analytical approximate solution for (110), we will consider (106) and (107) for .

Next, we use Newton’s binomial to transform (112) into the following approximate equation

Such as it will be explained in discussion section, although LT-HPM is an appropriate formulation for this case study, it is required a reformulation because the term , gives rise to an approximate solution inconsistent with the proposed boundary condition (from condition (111) we see that ). In order to delay the effect of the above mention term in the iterative process, we begin proposing the following homotopy, which is an original contribution of this work (see also, “Discussion” section).where we have introduced a trial function which may depend, for our purposes ofand is an integer number .

We rewrite the above equation as follows.

Next, using NDLT-HPM we proposeso that (116) becomes

It is worthwhile to mention, after evaluating the limit , we obtain the equation to solve (116) from (118) as it should be if (118) has to be a valid homotopy.

After considering values as large as we wish (), it is clear that for practical purposes, (118) is equivalent to

Thus, (118) let us to delay sufficiently the role of the term in the iterative process and replacing it for ; it eliminates the mentioned problem because the iterative process resulting will be consistent with . We will show the calculation of parameter will allow recover part of lost information for the omission of term , and the substitution of (113) instead of (112).

Following the usual Laplace transform algorithm, applying to (119) we obtain

As it is explained in (Spiegel 1998) we rewrite (120) as followswhere we have employed the initial condition .

Applying Laplace inverse transform to this last expression, we obtain

Next, we assume a series solution for , in the form

Substituting (123) into (122), we get

After comparing the coefficients of like powers of and perform the indicated operations, we have

Solving the above operations for and we obtainand so on.

By substituting (129)–(132) into (123), calculating the limit when , we obtain the following third order approximation.

From (111) and (133) we obtain

In order to calculate the values of and we will require that (134) satisfies the proposed boundary condition . Also, proceeding as it is indicated in “Basic idea of Laplace transform homotopy perturbation method with variable coefficients (CVLT-HPM)” section, we substitute (134) into (101) for and evaluating the expression so obtained for some value . Choosing , which lies into the interval of interest, (see “Discussion” section below), we obtain an system of algebraic equations for and , whose solution is given by:

After substituting (135) into (134) we obtain the following approximation

Case 4

This case is studied, by using a combination of LT-HPM and NDLT-HPM.

Finally we will consider the problem of extremizing an integral of the formwhere, is an arbitrary function, except having Taylor series.

Defining , it is clear that (f1) is formally equivalent to the case

After employing Euler equation (7), for we obtain the following differential equationwhere is a constant.

Separating variables, we rewrite (139) in terms of

Unfortunately the above integral hardly can be solved in terms of elementary functions for many functions.

In principle it is possible to provide an scheme in order to obtain analytical approximate solutions for problems like (139), by using LT-HPM.

After multiplying by , it is possible to rewrite (139) as

We employ Newton’s binomial to transform (141) into the following approximate equation

Following “Basic idea of non-linearities distribution Laplace transform-homotopy perturbation method (NDLT-HPM)” and “Basic idea of laplace transform homotopy perturbation method with variable coefficients (CVLT-HPM)” sections, it is convenient to construct a homotopy as followswhere in one hand we have chosen as a polynomial trial function (A is a parameter to be determined) and on the other hand we emphasized the possibility of embedded the homotopy parameter into the function (see “Discussion” section).

Next we rewrite the above equation as

Following Laplace transform algorithm, we apply

Such as it was explained in (Spiegel 1998) we rewrite (145) as follows

Applying Laplace inverse transform to the above expression we get

Next, in order to show how to work the methodology indicated by (147), we propose as cases: with two different boundary conditions and , in order to show the flexibility of the proposed methodology. So that for instance, we can express formally (139) for our first case as followsin the same manner after, we will consider the case ,.

Later, we will require expressing (148) in such a way that it does not involves For that purpose we obtain the following equivalent expression for (148)differentiating respect to , in order to eliminate we get

Thus, the proposed analytical approximate solution (147), adopts the form

Next, we assume a series solution for , in the form

In order to keep control of the different order approximations from (150), we introduce the following Taylor series expansions.

Substituting (151) and (152) into (150), we get

On comparing the coefficients of like powers of we have

After solving the above operations for and we obtainand so on.

By substituting (158)–(161) into (151) and calculating the limit when , results in a handy third order approximation

In order to calculate the values of and we will require that (162) satisfies the boundary condition Also, following “Basic idea of Laplace transform homotopy perturbation method with variable coefficients (CVLT-HPM)” section, we substitute (162) into (149) and evaluate the expression so obtained for some value, let us say , which lies into the interval of interest,(see “Discussion” below). After performing the above substitution, we obtain a system of algebraic equations for and , whose solution provides

Thus, substituting (163) into (162), we obtain a handy approximate solution for (158)

Proceeding the same manner for the conditions, , we obtain the values , . Finally, substituting them into (162), we get

Discussion

In order to use a pure numerical solution as reference, we utilized the scheme based on trapezoid combined with Richardson extrapolation (Ascher et al. 1995; Ascher and Petzold 1998) from the built-in numerical routines provided by Maple 15. Moreover, the command was setup with an absolute error (AE) tolerance of .

This work proposed LT-HPM and some its modifications, in order to find analytical approximate expressions for stationary functions solutions of the Euler equations, which turned out from four cases studied. Such as it was already mention, research in this area of knowledge is relevant because calculus of variations is a branch of the analysis with applications for both, pure and applied mathematics and expressing the basic principles of physics with simplicity and elegance among others.

We noted Euler equation (2) which describes these problems is of nonlinear nature and therefore the task of finding analytical solutions for this equation is in general a very difficult task. As a matter of fact, although some cases admit exact solutions, many variational problems admit just numerical solutions (Razzaghi and Elnagar 1994; Razzaghi and Yousefi 2001). Nevertheless this article showed one case study for which was not possible to find a numerical solution. Unlike (Lakestani and Kazemian 2013) where just one particular case of variational problem was solved, the goal of this paper was expose how LT-HPM and some of its modifications are employed in order to analyze a wide variety of cases study. Indeed, cases study 3 and 4 were presented in the most general way possible and to achieve it, we employed some modified versions of LT-HPM.

Thus, case study 1 proposed the basic version of LT-HPM (“Basic idea of the LT-HPM method” section) in order to find an analytical approximate solution for (46); essentially we used the procedure followed in (Filobello-Nino et al. 2013c). The proposed method expresses the problem of solving a nonlinear differential equation in terms of solving an algebraic equation for some unknown initial condition [see (64)]. This formulation has shown dependable performance for the case of differential equations without the presence of singular points and this case satisfies these requirements (Filobello-Nino et al. 2013c). Figure 1 shows the comparison between numerical solution of (46) and LT-HPM second order approximation (65). From this figure is clear the proposed solution provides a good approximation. In more precise terms, to verify the precision of (63) and (65) we calculate the square residual error (SRE) of (65), which is explained in “Appendix 2”. SRE is in general terms a positive number representative of the total error by using an approximate solution for a given problem (Marinca and Herisanu 2011). The resulting value for our case study was 0.001918936920, which confirms the accuracy of the proposed solution. If more accuracy has to be required, one must consider higher order approximations.

Fig. 1

Comparison between numerical solution of (46) (cross) and LT-HPM approximation (65) (line)

Our second variational problem (66) gave rise to the second order non homogeneous linear differential equation with variable coefficients (68). Since the aforementioned equation is something long, and following the indications of “Basic idea of Laplace transform homotopy perturbation method with variable coefficients (CVLT-HPM)” section, we ease the application of the proposed methodology, transforming (67) in its normal form (70), by using the substitution (69) (Simmons 1983). The equivalent equation (73) is also of variable coefficients; therefore we employed the modified version of LT-HPM explained in “Basic idea of Laplace transform homotopy perturbation method with variable coefficients (CVLT-HPM)” section (Filobello-Nino et al. 2015b). CVLT-HPM expresses the problem to be solved in terms of a differential equation for Laplace transform although it easily solvable [see (78)–(80)].

Once was expressed in terms of , we assumed a series solution in the form (81). CVLT-HPM method, calculated the different approximate orders, expressing the n-th iterative process in terms of integrals of the lower order approximations, previously calculated [see (83)–(87)].

Unlike basic LT-HPM algorithm (“Basic idea of the LT-HPM method” section), CVLT-HPM employed its homotopy freedom in order to substitute as it was originally defined by HPM method (“Standard HPM” section), for the polynomial trial function (constant), whose value was finally determined from the boundary condition . It is remarkable to note we obtained an exact solution for (98), through the convergent geometric series (94), which shows that in those cases where a problem has an exact solution, CVLT-HPM method is able to provide it.

Case study 3 was presented in the most general possible manner (99). This case is relevant because, such as it was explained in “Background” section, some values of corresponds to known cases. Thus, defines the problem of brachistochrone (Connor and Robertson 1997, 1998; Filobello-Nino et al. 2013d) and the problem of finding the curve that generates the surface of revolution of smallest area (Simmons 1983).

The corresponding Euler equation (102), presents an additional challenge; the presence of a highly nonlinear term that corresponds to the multiplication of by which difficult the application of LT-HPM algorithm. To overcome these difficulties, we proposed the substitution (106) in order to obtain (107). Despite of this equation involved the derivative of a power of dependent variable the methodology explained following Eqs. (114)–(119), turned out adequate to obtain an accurate analytical approximate solution for the case study [see (110), (111), and (112)]. A relevant point of the proposed method is, although (108) is easily expressed in terms of separated variables, and the resulting integral (109) can be solved performing the algebraic substitution , the dependent variable remains implicit, and therefore Newton–Raphson method would be required to obtain values of , given some values of . In the sequel, this case study not even admits a numerical solution.

On the other hand, such as we mentioned, the formulation of LT-HPM would turn out to be inappropriate for this case, because the term , would give rise to an approximate solution inconsistent with the proposed boundary condition As a matter of fact, if trial function was chosen as then , and the rest of the iterations would vanish. If trial function is (constant), then the subsequent iterative process would give rise to fractional powers of (in fact ). The problem of solutions containing these powers is that, its derivatives becomes infinite at x = 0, unlike the initial conditions of the problem. It is not difficult to verify that a similar problem would occur if the trial function were chosen as one of the following: and so on. Only the case provides a whole powers polynomial function consistent with With the purpose of delay the effect of the problematic term , and introduce in the iterative process the presence of term we proposed employ the homotopy equation (119). In order to motivate it, we proposed the homotopy (114), which allows introduce a trial function, even depending of the homotopy parameter , in order to reach our goals. Nevertheless, we note (114) is a valid homotopy because (113) is recovered in the limit , as expected. In particular by choosing in accordance with (117), let us employed for practical purposes, to Eq. (119) instead of (113) in the limit The above is due to the homotopy freedom of (114), which allowed delay sufficiently the role of the term and essentially substituted it for ; thus the iterative process resulted to be consistent with . The proposed homotopy (114), which is employed in order to delay and substitute a problematic term of a differential equation, including the introduction of a trial function depending of parameter [see (115) and (117)] is an original contribution of this work. A relevant fact of the proposed method is the determination of parameter allowed recovers part of lost information for the omission of term , and the substitution of (113) instead of (112). Then LT-HPM provided a third order approximation for (133), and through Eq. (111) an approximate solution (134) for . In order to calculate the values of and we required one hand, (134) satisfied the proposed boundary condition , and on the other hand we deduced another equation by substituting (134) into (101) for and then, evaluating the expression so obtained for some value which lied into the interval of interest . It is important to notice that unlike the others cases study, this example proposed a nonlinear differential equation defined in an interval relatively large in order to show the flexibility of the proposed method. Further research is required to propose a systematic procedure to obtain suitable cancelling points. For instance, the cancelling point for this case study was arbitrarily selected as after some trial and error attempts. This strategy was effective to obtain a low residual error. The resulting system of algebraic equations was solved to obtain the values of and and finally an approximate solution for [see (136)].

It is worth to note the value of the SRE of (136) is just 0.06334882582; from this we depict the accuracy of our approximate solution (see Fig. 2). If more accuracy has to be required, one must considerer higher order approximations. An outstanding point is that, numerical routines from Maple 15 failed in this case in order to provide a numerical solution as reference. Since the numerical solutions are often the only manner to obtain information from a nonlinear problem, LT-HPM method indeed, is a useful strategy capable of supporting approximate methods, especially for cases like this.

Fig. 2

Combination of NDLT-HPM and LT-HPM to obtain approximation (136)

Finally, our case study 4, proposed the problem of extremizing (137), which was also introduced in the most general way. Such as it was said represents an arbitrary function, although having Taylor series. In the sequel, the mathematical steps (137)–(139) and (141)–(147) were perform without assuming a particular

We note the proposed homotopy (143) chose the polynomial trial function and embedded the homotopy parameter into (see “Basic idea of non-linearities distribution Laplace transform-homotopy perturbation method (NDLT-HPM)” section). (Filobello-Nino et al. 2014b) employed this technique, provided for the freedom of homotopy formulations in order to redistribute the nonlinearities (Vazquez-Leal et al. 2012a) and non-polynomial nonhomogeneous terms, between the successive iterations of LT-HPM method, and thus to ease the search for an approximate solution. The main reason is that NDLT-HPM distributes the contributions of the above mentioned functions along the iterations for the successive parameter powers that allows to NDLT-HPM method being simpler at the different stages of iterations. As a matter of fact, we proposed a combination between LT-HPM and NDLT-HPM method, in order to obtain better results. As a case study, we considered (although for this choosing of (140) is integrable, it is expected the method works even for cases where (140) cannot be expressed in terms of elementary functions). Its Taylor expansions (152) in terms of the homotopy parameter and further substitution in the proposed approximate solution (150) show the manner which a function is included in the iterative process; in the same way we note the trial function was chosen in order to get a whole power series solution, with good precision. With the end of calculating the values of and we required that (162) satisfies the boundary condition Also, following “Basic idea of Laplace transform homotopy perturbation method with variable coefficients (CVLT-HPM)” section, we substituted (162) into (149) and evaluate the expression so obtained for some value, we chose , which lies into . After performing the above substitution, we obtain a system of algebraic equations for and , whose solution provided the approximate solution (164). As a matter of fact, proceeding at the same manner, for the proposed second case study given by the conditions , we obtained the approximate solution (165). For both cases, we obtained handy accurate approximate solutions; the SRE of (164) and (165) were 0.002934863992 and 0.05009724353 respectively (see Figs. 3, 4), from which we infer that the proposed method is potentially useful in the resolution of the kind of problems as (139). It is possible to obtain more accurate solutions, considering higher order approximations, although could be lost the handy character of our approximation. Worth mentioning even in these case studies, the value of , was chosen by trial and error attempts in order to find a good approximate solution.

Fig. 3

Combination of LT-HPM and NDLT-HPM to obtain approximation (164) (line) and numerical solution (cross)

Fig. 4

Combination of LT-HPM and NDLT-HPM to obtain approximation (165) (line) and numerical solution (cross)

It should be mentioned that, variational problems have been successfully studied for other authors, although the most of them have proposed just numerical approximations (Razzaghi and Elnagar 1994; Razzaghi and Yousefi 2001). On the other hand, HAM was employed to obtain an approximate solution for brachistocrone problem (Lakestani and Kazemian 2013). This work obtained precise solutions but long. As it is well known, HAM is accurate and powerful, but sometimes its approximate solutions turn out to be long and cumbersome, and for the same reason, they are not adequate for practical applications; on the other hand we obtained handy accurate approximate solutions with a reduced number of terms. Such as we have already mentioned, unlike (Lakestani and Kazemian 2013) that just studied one variational problem, the goal of this paper was expose in a general manner, how to find analytical solutions for several interesting cases of (2). On the other hand, (Ghaderi 2012) employed homotopy perturbation method (HPM) and Green’s function method in order to find approximate solutions to variational calculus. Nevertheless unlike our work, the mentioned article, dealt with moving boundary and isoperimetric problems. The main advantage of (Ghaderi 2012) is that the mathematical properties of the Green’s function are well established, however its study does not belong to the elementary mathematics; unlike the above, LT-HPM is an iterative method which is based on elementary Laplace transform properties (Spiegel 1998) and simple algebraic steps, making it an ideal technique for practical applications. Finally, future research will be required to improve the present proposal to be applied to infinite intervals.

Conclusions

This work introduced LT-HPM and some of its modifications in order to find analytical approximate solutions for some Euler’s ordinary differential equations, whose solutions extremizing the value of integrals of the form (1). As it is well known, the calculus of variations is a powerful branch of the analysis that can be applied to several problems of pure mathematics, also to express the basic principles of the mathematical physics in an elegant and simple way.

The relevance of this work consisted mainly on two points. On one hand it showed that the proposed methodology is potentially efficient in order to find solutions to the difficult linear and nonlinear problems arising from variational problems (even CVLT-HPM was able to provide exact solutions, see case study 2). On the other hand the article suggested for some cases, appropriate mathematical manipulations of the aforementioned equations, in order to transform them into a more accessible form, to increase the possibility of success (see cases study 3 and 4). The above is even more important, considering that many times even the numerical methods fail to solve these complicated equations, such as happened with our case study 3, and the numerical solutions and qualitative methods are usually the only manner to obtain information from a nonlinear problem. Finally, an outstanding issue for a future work is improving the methodology used to calculate the parameters introduced for our proposed solutions. A forthcoming paper should investigate a systematic manner to find the optimum values for the aforementioned parameters, without using the method of trial and error attempts in its determination.