Abstract and Applied Analysis

Volume 2014, Article ID 978636, 15 pages

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

## Comparison of Different Approaches to Construct First Integrals for Ordinary Differential Equations

^{1}Centre for Mathematics and Statistical Sciences, Lahore School of Economics, Lahore 53200, Pakistan^{2}Centro de Matemática, Computação e Cognição, Universidade Federal do ABC, Rua Santa Adélia 166, Bairro Bangu, 09.210-170 Santo André, SP, Brazil^{3}Department of Mathematics, School of Science and Engineering, Lahore University of Management Sciences, Lahore Cantt 54792, Pakistan

Received 15 December 2013; Accepted 2 March 2014; Published 7 May 2014

Academic Editor: Mariano Torrisi

Copyright © 2014 Rehana Naz 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

Different approaches to construct first integrals for ordinary differential equations and systems of ordinary differential equations are studied here. These approaches can be grouped into three categories: direct methods, Lagrangian or partial Lagrangian formulations, and characteristic (multipliers) approaches. The direct method and symmetry conditions on the first integrals correspond to first category. The Lagrangian and partial Lagrangian include three approaches: Noether’s theorem, the partial Noether approach, and the Noether approach for the equation and its adjoint as a system. The characteristic method, the multiplier approaches, and the direct construction formula approach require the integrating factors or characteristics or multipliers. The Hamiltonian version of Noether’s theorem is presented to derive first integrals. We apply these different approaches to derive the first integrals of the harmonic oscillator equation. We also study first integrals for some physical models. The first integrals for nonlinear jerk equation and the free oscillations of a two-degree-of-freedom gyroscopic system with quadratic nonlinearities are derived. Moreover, solutions via first integrals are also constructed.

#### 1. Introduction

The study of conserved quantities plays a great role in mathematical physics and in applied mathematics. For instance, a considerable number of phenomena have some kind of “conservation.” Examples can be easily found from the hydrodynamics, electrodynamics, shallow water phenomena, and so forth. One can also mention the celebrated Kepler's third law or the conservation of energy in the classical mechanics, particularly the one-dimensional harmonic oscillator. In regard to these last two phenomena, the conserved quantity is called first integral, which is the analogous of conservation laws for ordinary differential equations models.

In a recent paper, Naz et al. [1] studied different approaches to construct conservation laws for partial differential equations. The purpose of this paper is to analyze all different approaches for construction of first integrals for ordinary differential equations. In fact, different approaches to derive first integrals can be grouped into three categories: direct methods, Lagrangian or partial Lagrangian formulations, and characteristic (multipliers) approaches. In 1798, Laplace [2] developed a method known as the direct method for the construction of first integrals. Although such a method does not originally require any symmetry of the considered equation, Kara and Mahomed [3] developed a relationship between symmetries and conservation laws. The joint conditions of symmetry and direct method are used to construct the first integrals.

Noether's theorem [4] is a powerful technique to derive first integrals for the differential equations having Lagrangian formulations using its symmetries, although it requires a suitable Lagrangian. Kara et al. [5] developed the partial Noether approach. The partial Noether approach is applicable to differential equations with or without a Lagrangian. The interested readers are referred to [6–11] for discussions on first integrals by the Noether approach and partial Noether approach. Ibragimov [12] introduced the concept of formal Lagrangian formulation for differential equations and its adjoint as a system. Atherton and Homsy [13] introduced the adjoint variational principle for such systems. Then Ibragimov [12] incorporated symmetry considerations and provided formulas to construct the first integrals similar to those provided by Noether's approach. The concepts of self-adjointness [14, 15], weak self-adjointness [16], nonlinear self-adjointness [17–20], and quasi-self-adjointness [21, 22] are used to construct first integrals by this approach.

The characteristic, multiplier, or integrating factor methods are also very powerful and elegant methods for construction of the first integrals. There are four different approaches based on the knowledge of the characteristics. The first method developed by Steudel [23] in 1962 involves writing a first integral in the characteristic form. The characteristics or integrating factors are the multipliers of the differential equations that makes them exact. To derive the first integrals by this method first of all the characteristics need to be determined. The second method is based on the first method and it involves the variational derivative (see Proposition 5.49 in Olver [24]). The reader is referred to [25–27] for a good account of understanding how to compute multipliers and first integrals. In the third approach, we compute the variational derivatives on the solution space of given differential equations and these characteristics sometimes may correspond to an adjoint symmetry not to first integral. The fourth approach according to Anco and Bluman [28] provides formulae for finding first integrals. In the last few decades, the researchers focused on the development of symbolic computational packages based on different approaches of first integrals and these packages are well documented in [29] and references therein.

The well-known Noether identity can be expressed in terms of Hamiltonian function and symmetry operators (see, e.g., [30, 31]). This is a simple and elegant way to construct first integrals of Hamiltonian equations consisting of a system of first-order differential equations. No integration is required here to construct solutions.

Lie approach as described, for example, by Ibragimov and Nucci [32] and Mahomed [33], is successfully applied to differential equations to derive the exact solutions. On the other hand, the knowledge of first integrals enables one to reduce or completely solve an ordinary differential equation. In fact, if one considers an th-order scalar ordinary differential equation having independent first integrals, one can obtain from those first integrals the general solution of the considered equation possessing constants. Kara et al. [34] explored the solutions of differential equations using the Noether symmetries of a Lagrangian associated with the first integrals. Using the relationship between Noether symmetries and first integrals [35] the reductions and exact solutions of differential equations were derived. The generalization of this idea is the association of Lie-Bäcklund symmetries [3] and nonlocal symmetries [36, 37] with a first integral and it led to the development of the double reduction theory to find reductions and solutions [38–43].

The paper is organized in the following manner. The fundamental relations are defined in Section 2. We present the main ideas behind the mentioned methods in Section 3. Then, in Section 4, we apply these different approaches to the harmonic oscillator equation. In Section 5 some solutions are obtained via first integrals. Relations between Hamiltonian functions and first integrals are discussed in Section 6. In Section 7, the first integrals for nonlinear jerk equation are computed by the Noether approach for the equation and its adjoint as a system and by the multiplier approach. The exact solutions of jerk equation for different cases are also established via first integrals. The first integrals for the free oscillations of a two-degree-of-freedom gyroscopic system with quadratic nonlinearities are also derived. Finally, concluding remarks are presented in Section 8.

#### 2. Fundamental Relations

The following definitions are taken from the literature (see, e.g., [44, 45]).

Consider a th-order ordinary differential equation system where is the independent variable and , are the dependent variables. We will adopt the summation convention and there is summation over repeated upper and lower indices.

The total derivative with respect to is The following are the basic operators defined in , the vector space of differential functions.

The Lie-Bcklund operator is defined as where in which .

The Euler operator is given by

The characteristic form of Lie-Bcklund operator (3) is
in which is the* Lie characteristic function* defined by

The* Noether operator* associated with a Lie-Bäcklund operator is
where

*First Integral.* A first integral of system (1) is a differential function , such that
for every solution of (1).

#### 3. Approaches to Construct First Integrals

Now we present various approaches to construct first integrals taken from the literature.

##### 3.1. Direct Method

The direct method was first used by Laplace [2] in 1798 to construct all local first integrals. The determining equations for the first integrals for the direct method are

##### 3.2. Symmetry and First Integral Relation

Kara and Mahomed [3] added a symmetry condition to the direct method. The Lie-Bäcklund symmetry generator and the first integral are associated with the following equation:

The first integrals are computed by the joint conditions (11) and (12).

##### 3.3. Noether’s Approach

In 1918, Noether developed a new approach to construct first integrals [4] and it is currently known as Noether's approach.

*Euler-Lagrange Differential Equations.* If there exists a function satisfying
then is called a Lagrangian of (1) and relationship (13) yields the Euler-Lagrange differential equations. Equations (1) and (13) are equivalent.

*Noether Symmetry Generator.* A Lie-Bäcklund operator is a Noether symmetry generator associated with a Lagrangian of the Euler-Lagrange differential equations (13) if there exists a function such that

*Noether First Integral.* For each Noether symmetry generator associated with a given Lagrangian corresponding to the Euler-Lagrange differential equations, there corresponds a function known as a first integral and is defined by
or
where are the characteristics of the first integral.

In the Noether approach we need to construct Lagrangian . The Noether symmetries are then computed from (14) and finally (16) provides the first integrals corresponding to each Noether symmetry. The reader is guided to [46] for further discussions about this technique and its relations with the so-called Noether symmetries.

##### 3.4. Partial Noether Approach

The partial Noether approach for construction of first integrals was introduced by Kara et al. [5] and it can be useful for constructing first integrals when the differential equation does not have a known Lagrangian.

*Partial Lagrangian.* Suppose that the th-order differential system (1) can be expressed as
A function , is known as a partial Lagrangian of system (17) if
provided for some . Here is an invertible matrix.

*Partial Noether Operator.* The operator satisfying
is a partial Noether operator corresponding to the partial Lagrangian .

The first integrals of the system (1) associated with a partial Noether operator corresponding to the partial Lagrangian are determined from (16).

We can also use the partial Noether approach for equations arising from the variational principal and have the Lagrangian.

##### 3.5. Noether Approach for a System and Its Adjoint

*Adjoint Equations.* Let be the new dependent variables. The system of adjoint equations to the system of th-order differential equations (1) is defined by (Atherton and Homsy [13], Ibragimov [12])
where

*Symmetries of Adjoint Equations.* Suppose system (1) has the generator
Ibragimov [12] showed that the following operator is a Lie point symmetry for the system (1) and (20):
The operator (23) is an extension of (22) to the variable and
yields .

*Conservation Theorem.* Every Lie point, Lie Bcklund, and nonlocal symmetry of the system of th-order differential equations (1) yields a first integral for the system consisting of (1) and the adjoint equations (21). Let be the Lagrangian defined by

Then the first integrals are given from the formula where , are the coefficient functions of the generator (22). The first integrals constructed from (26) contain the arbitrary solutions of adjoint equation (21) and, thus, for each solution one has first integrals.

The dependence on the nonlocal variable provides a nonlocal first integral. One can eliminate such variable if the original system of ODEs is nonlinearly self-adjoint [14, 17] and to the equations admitting this remarkable property one can establish a first integral for the original system.

##### 3.6. Characteristic Method

According to Steudel [23] and Olver [24], the first integral can be expressed in the characteristic form as where are the characteristics or multipliers.

##### 3.7. Variational Approach

The variational approach was developed by Olver [24]. The variational derivative of (27) yields all the multipliers for which the equation can be expressed as a local first integral. The multipliers determining equation is and it holds for arbitrary functions .

##### 3.8. Variational Approach on Solution Space of the Differential Equation

In this approach, the multiplier determining equation is obtained by taking the variational derivative of (27) on the solution space of the differential equation; that is, Equations (29) are less overdetermined than (28). Sometimes the characteristics constructed from (29) may correspond to adjoint symmetries (see [28]) and not to a first integral.

##### 3.9. Integrating Factor Method for First Integrals

Consider the system (1) and let and such that

On the solutions of system (1), it is concluded that , which means that is a first integral of (1) and the functions are integrating factors; see [28] for further details and deeper discussion.

The linearized system to (1) is given by where

The adjoint of the linearized system (31) is given by

Moreover, the operators and satisfy the identity where If satisfy the condition where then the first integral is In (38) and are any fixed functions such that the function is finite, while

We finish with the following definition.

*Definition 1. *The system (1) is said to be self-adjoint if and only if .

#### 4. First Integrals of Simple Harmonic Oscillator

We compute the first integrals of simple harmonic oscillator by utilizing different approaches. Consider

##### 4.1. Direct Method

Equation (11) with becomes where

Equation (43) after expansion results in or If we further restrict to be then (46) becomes Splitting (47) according to derivatives of , we obtain The system of (49) is solved for , and to obtain

##### 4.2. Symmetry and First Integral Relation

The first-order prolongation of the Lie point symmetry generators of (42) is where

The first integrals are computed by the joint conditions (11) and (12).

The second important aspect of this approach is that we can associate a symmetry with a first integral. The relationship (12) holds for symmetry and first integral and, thus, symmetry is associated with . Similarly is associated with . This association of symmetries with a first integral helps in finding a solution via double reduction theory [38–43].

##### 4.3. Noether's Approach

Equation (42) admits the standard Lagrangian which satisfies the Euler Lagrange equation . Now we show how to compute the Noether symmetries corresponding to a Lagrangian (53).

The Noether symmetry determining (14) results in where , , lie symmetry operators and is the gauge terms. Expansion of (54) gives Noether symmetry determining equation

The separation of (55) with respect to powers of derivatives of gives rise to

The solution of system (56) is

Formula (16) with , , and from (57) yields the first integrals (50).

##### 4.4. Partial Noether Approach

Equation (42) admits the partial Lagrangian and the corresponding partial Euler-Lagrange equation is where

The partial Noether operators corresponding to satisfy (19); that is, The usual separation by derivatives of gives

System (61) yields

Formula (16) with , , and from (62) yields the first integrals (50). Hence the first integrals in each case are with respective characteristic . Here the partial Noether's approach yields all nontrivial first integrals as obtained by Noether's approach. The difference here lies in the forms of and which are distinct from the ones used in the Noether approach.

##### 4.5. Noether Approach for a System and Its Adjoint

The adjoint equation for (42) is and this yields

Let , then , and . Substituting these expressions of and into (64) and equating the coefficients of and to zero, one obtains . Thus, (42) is nonlinearly self-adjoint.

The Lagrangian for the system consisting of (42) and (64) is The Lagrangian satisfies The formula for first integrals from (26) is Equation (67) with from (65) results in According to the conservation theorem, every Lie point, Lie-Bäcklund, and nonlocal symmetry of the system of second-order differential equation (42) yields a first integral for the system consisting of (42) and the adjoint equation (64). For the Lie symmetry , the first integrals are where is the solution of adjoint equation (64). Note that yields and , can be obtained from the substitution and , respectively. One can use the other Lie symmetries given in (51) to derive the first integrals but one requires the solution of adjoint equation to construct first integrals. In order to illustrate this fact, let us consider the generator It is easy to check that (70) is not a Noether symmetry generator. Substituting and into (68), one arrives at The substitution in (71) provides the trivial first integral . However, setting and , respectively, into (71), then the first integrals and are obtained again.

##### 4.6. Characteristic Method

For (42) assume the characteristics and first integrals of form ; then formula (27) yields Equation (72), after separating, with respect to gives Equations (73) finally result in (46) and, hence, following the same procedure as we did for the direct method, five first integrals (50) are obtained.

##### 4.7. Variational Approach

For the variational approach with multiplier of form , we have Equation (74), after expansion, takes the following form: Separation of (75) with respect to yields

The solution of (76) gives rise to the following characteristics: A multiplier has the property Equation (78) with multipliers from (77) and yields the first integrals as given in (50).

##### 4.8. Variational Approach on Space of Solutions of the Differential Equation

For (42), condition (29) results in Expanding (79), we have and this yields where are constants. The multipliers with respect to constants are the same as obtained in Section 4.7 and yield the first integrals obtained in (50). The multiplier associated with is which does not correspond to any first integral. It might correspond to an adjoint symmetry.

##### 4.9. Integrating Factor Method for First Integrals

Applying the results of Section 3.9 to our considered equation means to follow closely the first example presented in [28].

Substituting into (33) and (36), it is, respectively, obtained

Noticing that is a solution of both (82) and (83), one can take into (41). Therefore, from (39) and (38), it is obtained and

#### 5. Exact Solutions via First Integrals

The Noether symmetries associated with the first integrals can be utilized to derive the exact solutions of ordinary differential equations [34].

If is a Noether symmetry and is a first integral of (1) corresponding to a first-order Lagrangian , then the following properties are satisfied [34]: where

Proposition 2. *Suppose is a symmetry of , where ; then it satisfies
*

Proposition 3. *In (87) if and , then is a point symmetry of reduced equation , in which is an arbitrary constant.*

Now we will compute the exact solutions of (42) using its first integrals which are reduced forms of the equation under consideration. The first integrals reduce an th-order ODE to th-order ODE. For scalar first-order ODE, the first integrals transform to quadrature whereas for scalar second-order ODE the first integrals result in the first-order ODEs. Some of these reduced forms (first integrals) can be solved directly. The other reduced form can be transformed to quadrature by using the Noether symmetries with its associated first integrals which yield the exact solutions. The first three integrals , , and of (42) yield a solution directly. Since which implies , the first integral in (50) can be written as which can also be expressed as Equation (153) is a variable separable and yields and this comprises the exact solution of ODE (42).

A similar procedure is adapted to get the following exact solution of (42) using or :

Now we show how one can find the exact solution of (42) using Noether symmetries associated with the first integral. The Noether symmetry is associated with the first integral in (50). The induced equation can be expressed as Using (92) one can easily find the invariant and it reduces (93) to Equation (95) is expressible as a variable separable and it finally yields or The solution in (94) is an exact solution of (42) with which can be determined from (96) or (97).

Similarly, the Noether symmetry associated with in (50) provides the exact solution where satisfies or

#### 6. Hamiltonian Functions and First Integrals

Suppose is the independent variable and are the phase space coordinates. The derivatives of , with respect to are given by where is known as the total derivative operator with respect to . Here we present the basic operators needed in the sequel after introducing the necessary notations.

The Euler operator, for each , is and the associated variational operator is Applying operators (103) and (104) on equated to zero results in the following canonical Hamilton equations: Equations (106) are obtained using and . Equation (105) is the well-known Legendre transformation which relates the Hamiltonian and Lagrangian, where and .

Let be the operator in the space . The operator in (107) is a generator of a point symmetry of the canonical Hamiltonian system (106) if it satisfies [24] on the system (106).

In [24], the authors have studied the Hamiltonian symmetries in evolutionary or canonical form. The symmetry properties of the Hamiltonian action have been considered in the space in [30, 31]. They presented the Hamiltonian version of Noether's theorem considering the general form of the symmetries (107).

The following important results which are analogs of Noether symmetries and the Noether theorem (see [24, 30, 44, 47] for a discussion) were established.

Theorem 4 (Hamilton action symmetries). *A Hamiltonian action
**
is said to be invariant up to gauge associated with a group generated by (107) if
*

Theorem 5 (Hamiltonian version of Noether's theorem). *The canonical Hamilton system (106) which is invariant has the first integral
**
for some gauge function if and only if the Hamiltonian action is invariant up to divergence with respect to the operator given in (107) on the solutions to (106).*

##### 6.1. First Integrals of Harmonic Oscillator in Hamiltonian Framework

Let us transfer the preceding example into the Hamiltonian framework and define The Hamiltonian function for this problem is The canonical Hamiltonian equations (106) for Hamiltonian function (113) result in The Hamiltonian operator determining equation (110), after expansion, yields in which we assume that , , and . One can also assume these functions to be dependent on . We have chosen dependence to simplify the calculations here and this leads to at least one Hamiltonian Noether operator. Equation (115) with the help of (114) can be written as

One can separate (116) with respect to powers of derivatives of and finally arrive at the following Hamiltonian Noether operators and the gauge terms:

The first integrals from formula (111) are It is worthy to notice here that no integration is required to derive solutions of (114).

#### 7. Applications to Some Models from Real World

In this section we apply the considered techniques to some equations arising from concrete problems, namely, the jerk equation and free oscillations with two-degree-of-freedom gyroscopic system with quadratic nonlinearities.

##### 7.1. Jerk Equation

According to Gottlieb [48, 49], the most general nonlinear jerk equation is where the prime denotes differentiation with respect to and are constants. In (119), at least one of should be different from zero and if , then so that the jerk equation is not a derivative of a second-order ODE [50].

###### 7.1.1. First Integrals for Nonlinear Jerk Equation by Noether Approach for a System and Its Adjoint

Let us look for first integrals for (119) using Noether approach for a system and its adjoint. The adjoint equation for (119) is

The Lagrangian for system and is . The system and possesses a first integral where is any Lie point symmetry of (119). Now we use strictly self-adjointness for eliminating the nonlocal variable in the first integral (121).

Using the self-adjoint condition we conclude that and . Then we conclude that the strictly self-adjoint subclass of (119) is given by the family

The only admitted Lie point symmetry generator of (124) is If we take , (119) admits not only (125), but also

The Lie point symmetry generator (125) provides the trivial first integral . However, from (126) one can construct a nontrivial one. In fact, using (126), we obtain . Substituting this expression for , , and into (121) and setting , after reckoning, we have

###### 7.1.2. First Integrals for Nonlinear Jerk Equation by Multipliers Approach

Now we will derive first integrals for nonlinear jerk equation by multipliers approach. Assume multipliers of form . The multipliers determining equation (29) becomes After expansion of (128), the multipliers and first integrals are computed for specific values of parameters.

*Case 1 (). *
Equation (119) reduces to
The multipliers and first integrals are

*Case 2 (). *
Equation (119) reduces to
The multipliers and first integrals are

*Case 3 (). *
Equation (119) reduces to
The multipliers and first integrals are

*Case 4 (). *
Equation (119) reduces to
Equation (128) yields two multipliers
The first integrals are given by

*Case 5 (). *
Equation (119) reduces to
The multipliers and first integrals are

*Case 6 (). *
The multipliers and first integrals for this case are
and is the same as the first integral (127) derived by Noether approach for system and its adjoint.

###### 7.1.3. Reduction of Order and Implicit Solution to Jerk Equation

We now utilize first integrals to compute the exact solution of the jerk equation (119).

We firstly consider the first integrals and obtained in Case 2. Setting and , straightforward calculations yield The routine calculations show that (141) finally results in which satisfies , , and

Now we obtain an implicit solution to using the first integrals given in Case 1. Again, assuming that and , where and are arbitrary constants and and are the first integrals established in Case 1, we finally arrive at

From (1