Journal of Applied Mathematics

Volume 2013 (2013), Article ID 147921, 15 pages

http://dx.doi.org/10.1155/2013/147921

## Second-Order Systems of ODEs Admitting Three-Dimensional Lie Algebras and Integrability

^{1}Department of Mathematics, Quaid-i-Azam University, Islamabad 45320, Pakistan^{2}Centre for Differential Equations, Continuum Mechanics and Applications, School of Computational and Applied Mathematics, University of the Witwatersrand, Wits 2050, South Africa

Received 27 November 2012; Accepted 20 January 2013

Academic Editor: Mehmet Pakdemirli

Copyright © 2013 Muhammad Ayub 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

We present a systematic procedure for the determination of a complete set of *k*th-order
() differential invariants corresponding to vector fields in three variables for three-dimensional Lie algebras. In addition, we give a procedure for the construction of a system of two *k*th-order ODEs admitting three-dimensional Lie algebras from the associated complete set of invariants and show that there are 29 classes for the case of *k* = 2 and 31 classes for
the case of . We discuss the singular invariant representations of canonical forms for systems of two second-order ODEs admitting three-dimensional Lie algebras. Furthermore, we give an integration procedure for canonical forms for systems of two second-order ODEs admitting three-dimensional Lie algebras which comprises of two approaches, namely, division into four types I, II, III, and IV and that of integrability of the invariant representations. We prove that if a system of two second-order ODEs has a three-dimensional solvable Lie algebra, then, its general solution can be obtained from a partially linear, partially coupled or reduced invariantly represented system of equations. A natural extension of this result is provided for a system of two *k*th-order () ODEs. We present illustrative examples of familiar integrable physical systems which admit three-dimensional Lie algebras such as the classical Kepler problem and the generalized Ermakov systems that give rise to closed trajectories.

#### 1. Introduction

Realizations of Lie algebras in terms of vector fields are applied, in particular, for the integration and classification of ordinary differential equations (see, e.g. [1–5]). In spite of the importance for applications, the problem of complete description of realizations has not been solved even for cases when either the dimension of the algebras or the dimension of the realization space is a fixed small integer [6]. The classification of finite-dimensional Lie algebras was initiated by Lie [7] who also presented realizations in terms of vector fields in (1+1)- and (1+2)-dimensional spaces. In the case of vector fields in the two-dimensional space, he gave a complete classification in the complex domain. For realizations of vectors fields in the three-dimensional space, Lie presented a partial classification which was recently completed. All the possible complex Lie algebras of dimension less than four were already obtained by Lie himself [8, 9]. Bianchi [10] investigated three-dimensional real Lie algebras. Complete and correct classification of four-dimensional real Lie algebras was obtained by Mubarakzjanov [11]. Wafo Soh and Mahomed [12] used the results of Bianchi and Mubarakzjanov to classify realizations of three- and four-dimensional real Lie algebras in the space of three variables, and they also described systems of two second-order ODEs admitting real four-dimensional Lie algebras.

The realizations of real Lie algebras of dimension no greater than four of vector fields in the space of an arbitrary (finite) number of variables were given in [6]. Popovych et al. [6] used the results of Bianchi and Mubarakzjanov with a different approach. They found that the work [12] does not provide the complete list of realizations of real Lie algebras of dimensions three and four.

Systems of two second-order ODEs arise frequently in several applications in classical, fluid, and quantum mechanics as well as in general relativity. Hence, these have been studied vigorously over the years. A sample of the literature as far as symmetry Lie algebras of ODEs are concerned can be found in [13–15]. In particular, the classification of systems of two second-order ODEs admitting real four-dimensional Lie algebras was attempted in [12]. Gaponova and Nesterenko [16] gave a classification of system of two second-order ODEs admitting real three- and four-dimensional Lie algebras by using the results of [6].

In the this paper, by utilizing realizations given in [6], a systematic procedure for finding a complete set of th-order () differential invariants including basis of invariants corresponding to vector fields in three variables for three-dimensional Lie algebras is presented. In addition, we give a procedure for the construction of two th-order ODEs possessing three-dimensional Lie algebras from the associated complete set of invariants. We show that there are classes for and classes for . Moreover, we mention those cases in which regular systems of two second-order ODEs are not obtained for three-dimensional Lie algebras. We present a brief discussion on the singular representation of canonical forms for systems of two second-order ODEs admitting three-dimensional Lie algebras. Furthermore, we provide an integration procedure for canonical forms for systems of two second-order ODEs admitting three-dimensional Lie algebras. This procedure is composed of two approaches, namely, one which depends on general observations of canonical forms which admit three-dimensional Lie algebras; the second depends on a complete set of differential invariants and is applicable only for those cases which admit three-dimensional solvable Lie algebras. We show by familiar physical examples how these integration approaches work. In the second approach, for solvable Lie algebras, we analyze integrability in terms of integration of the invariant representations. We prove that if a system of two second-order ODEs admits a solvable three-dimensional Lie algebra, then, its general solution can be deduced from a partially linear, partially coupled or invariant reduced system of equations. A natural extension of this result is stated for a system of two th-order () ODEs.

In Table 1, for ease of reference, we list the realizations of three-dimensional Lie algebras as given in [6]. In the second section, for the Lie algebras of vector fields considered, we discuss and deduce a complete set of differential invariants including basis of invariants with the construction of the corresponding systems of two th-order ODEs. Moreover, for systems of two second-order ODEs admitting three-dimensional Lie algebras, we discuss the singular representation and those cases in which we do not obtain systems. In the end, in Table 2, we give a complete list of second-order differential invariants corresponding to vector fields in three variables of three-dimensional real Lie algebras and the associated complete list of canonical forms for systems of two second-order ODEs. In Section 3, we provide an integration procedure for the canonical forms of systems of two th-order () ODEs admitting three-dimensional Lie algebras. This is stated as theorems for systems admitting solvable algebras. Reduction of canonical forms for systems of two second-order ODEs into invariant systems of first-order and algebraic equations are given in Table 2. The paper ends with a brief conclusion.

#### 2. Invariants and Systems of ODEs

When one calculates the symmetries of a given differential equation, one finds the generators in the form of vector fields and then computes the Lie brackets to get the structure constants of the particular Lie algebra one has found. We can start from a given Lie algebra with a set of structure constants and ask which vector fields in at most three variables satisfy the given Lie bracket relations with none of the vector fields vanishing. We thus ask for possible realizations of the given Lie algebra. We make use of the realizations of three-dimensional Lie algebras in the (1+2)-dimensional space given in [6]. We utilize these realizations in order to obtain differential invariants and invariant representations for system of two second-order and higher-order ODEs. Extensions to higher order invariants are given to understand the behavior of such invariants. This is the inverse problem in symmetry analysis and was initiated by Lie.

For completeness and ease of reference, we list the realizations of three-dimensional real algebras as given in [6] in Table 1. However, we mostly keep the format of [12].

Differential invariants play an important role in the structure and reduction of differential equations. Differential invariants can be obtained by various approaches and we refer, for example, to [16] for realizations in (1+2)-dimensional space. Differential invariants obtained by various approaches corresponding to three-dimensional Lie algebras have very interesting properties. In the investigation of the invariant approach for a system of two second-order ODEs admitting three-dimensional Lie algebras, it is found that, for a regular system of ODEs, there are four independent invariants, three of which play the role of basis and the fourth one is derived from the remaining three.

For the case of a regular system of two second-order ODEs admitting three-dimensional Lie algebras, there are four functionally independent invariants. For example, consider the system of two second-order ODEs admitting the three-dimensional Lie algebra . We have the four independent invariants Here, the invariants and are of order two, that is, equal to the order of the system of ODEs that is constructible. Similarly for the regular system of two third-order ODEs admitting the same Lie algebra ; we have the six invariants from which we determine the third-order system. Here, the invariants and are of order three, that is, equal to the order of the considered system of ODEs. By proceeding in this way for regular systems of two th-order ODEs admitting three-dimensional Lie algebras, the number of order invariants is two.

To study differential invariants further, we consider the following example. Consider the system admitting the three-dimensional Lie algebra . We arrive at the following set of invariants This set involves three invariants , and which play the role of basis of invariants for the considered system of second-order ODEs admitting the three-dimensional Lie algebra . This set is a complete set of invariants for the system. In the complete set, one of the invariants is obtained by differentiation. We state a result on invariant differentiation of invariants corresponding to a Lie algebra admitted by any system of ODEs which is already well-known for scalar ODEs in the literature.

Proposition 1. * If and are invariants of a Lie algebra admitted by any system of ODEs, then, is also its invariant. *

The proof is identical (except here we have many dependent variables) to the scalar case and uses the operator identity (see, e.g., [17]).

One can write where the invariant differentiation operator is There arises a natural question. Can we construct a th-order system from a basis of invariants corresponding to the symmetry vectors of a second-order system? We answer this question.

Suppose that we have a complete set of invariants representing a second-order system of ODEs admitting a three-dimensional Lie algebra, say for the above case . Now if we differentiate both of the highest-order invariants from this set with respect to another basic invariant, we get two new invariants and the total number of invariants become six for this case, including the previous four. These two new invariants play the role of third-order invariants representing the third-order system of ODEs admitting the three-dimensional Lie algebra . We therefore have for this case where and and are third-order invariants. We can thus construct the following third-order system We have constructed a system of two third-order ODEs admitting the three-dimensional Lie algebra from their complete set of invariants. This system of ODEs is represented by a complete set of invariants which contains six functionally independent invariants, namely, , and , where three of them and are a basis of invariants and two of them are third-order invariants.

Now if we repeat this process on the third-order invariants, we obtain eight invariants including the previous ones which form a complete set of invariants for the fourth prolonged group. From these, we can construct a system of two fourth-order ODEs admitting the three-dimensional Lie algebra . We can generalize this to a system of two th-order ODEs. However, there are some nonstandard cases to look at first. These behave differently. There are four cases , and in which two of them and yield singular invariants and the other two and give a single equation and do not form a system.

##### 2.1. Singular Case

In both of the cases and , we have an extra condition, due to which the rank of the coefficient matrix on the solution manifold is less than the rank of the coefficient matrix on the generic manifold. So singular invariants are obtained in both of these cases.

Consider . We find the following set of singular invariants where and the corresponding system are In both of the cases, the singularity of the invariants is removed in the third- and higher-order systems. In the case of , we get the following complete set of invariants: and the corresponding third-order system is In this case, and form a basis of invariants, and here the third-order invariants are and , where . For the other case , it is easy to see that where again . For each of the remaining cases and , we obtain the same set of second-order invariants for which the highest-order invariant is only one and . So a system of second-order ODEs cannot be constructed. Only a single scalar equation is obtained, namely However, in both of these cases, this behavior changes for higher orders. In fact for , we have the third-order invariants and for the third-order invariants are They give rise to a regular system of third-order ODE with a complete set of invariants having the same format as the singular cases similar to (13).

We find that, for a system of two second-order ODEs admitting three-dimensional Lie algebras, there are 29 classes which are listed in Table 2 together with their second-order differential invariants.

In the case of a system of two third- and higher-order ODEs admitting three-dimensional Lie algebras, there are 31 canonical forms, and these can be obtained by invariant differentiation as described above. Note here that and are admitted by a system of two third- and higher-order ODEs whereas they do not form a system of two second-order ODEs.

We state these results in the form of a theorem.

Theorem 2. *A regular system of two th-order ODEs admitting three-dimensional Lie algebras can be represented by a complete set of invariants which contains complete functionally independent invariants in one of the following forms:
**
In both and , , and form a basis of invariants, where the highest order invariants in are and in they are .**The total number of systems of two th-order ODEs having the above structure is thirty-one, in which form corresponds to the cases , and and to the remainder.*

#### 3. Integrability

There are basically two approaches to integrability using Lie point symmetries. One is that of successive reduction of order using differential invariants for scalar ODEs. The other is that of canonical forms which may apply to scalar as well as system of ODEs, although the latter is not straightforward. However, the canonical forms approach only applies when equations are classified according to the symmetry Lie algebras. In Table 2, we have given the complete set of second-order differential invariants of three-dimensional Lie algebras admitted by systems of two second-order ODEs and also classified systems of two second-order ODEs for their point symmetries together with their canonical forms. We provide two approaches for the integration procedure on the basis of these canonical forms. The first is a general approach depending on the subdivision of the canonical forms based on general observations. The second is the differential invariants approach based on the representation of the canonical forms in terms of first-order differential and algebraic invariants.

##### 3.1. General Integration Approach

In this approach, we judiciously subdivide all the canonical forms into four types. They are as follows.

*Type I*:

*Type II*:

*Type III*:

*Type IV*: Remaining canonical forms except those given above.

We notice that the Types I and II systems comprise linear equations. The Type I system is completely integrable. Also that the Type II system is uncoupled and depends on two quadratures of the equation for its solution. Further, we observe that the Type III has equations that comprise scalar second-order ODEs in the dependent variable . In view of these observations, we are led to the following definitions.

*Definition 3. *A system of two second-order ODEs is said to be partially uncoupled if one of its equations constitutes a scalar second-order equation in one of the dependent variables and the other has both.

It is easily seen that the Types I to III equations are partially uncoupled.

*Definition 4. *A partially uncoupled system of two second-order ODEs is partially linear if one of its equations forms a scalar linear second-order equation in its own right.

We see that the Types I and II systems are partially linear.

*Definition 5. *A system of two second-order ODEs is said to be partially linearizable by point transformation if it can be transformed via an invertible transformation to a partially linear form.

The following system is partially linearizable by a point transformation as the equation in is the well-known Painlevé-Ince equation which is linearizable via a point transformation.

Let be a basis of a three-dimensional algebra and let the rank of the coefficient matrix associated with these operators be .

We formulate the following theorems the proofs of which are evident.

Theorem 6. *A system of two second-order ODEs is partially linearizable by point transformation if and only if it admits the Lie algebra or with rank in both cases. *

Theorem 7. *A system of two second-order ODEs can be reduced to a partially uncoupled system by point transformation if and only if it admits one of the Lie algebras , , and with rank in all cases. *

The proofs follow from the above listed Types I to III algebras and representative equations.

Generally, we can integrate Type I completely. Types II and III can be integrated when the associated second-order equation in the system is integrable. Moreover, only particular cases of Type IV can be integrated depending on the given Lie algebra and corresponding system of ODEs. Some of these difficulties were also pointed out for systems of two second-order ODEs that admit four symmetries [12]. Here, we further discuss all the cases in detail before embarking on a further systematic study that captures the main results.

Type I is trivially integrable by quadratures.

The Type II system is integrable if the associated second-order ODE in admits a two-dimensional Lie algebra and Lie's four canonical forms for scalar second-order ODEs become applicable. So, we have Lie reducibility here.

In the case of Type III, one requires one more symmetry for the associated equation in as it already has translations in symmetry. As an example we consider the system The system (25) has symmetry generators and its Lie algebra is three-dimensional. The associated equation in is . Its Lie algebra is two-dimensional spanned by . By using the transformation the associated equation in transformed form becomes . We have , which is Type III in Lie's classification of scalar second-order ODEs. Solving it in transformed form by Lie's method and reverting back to the original variables, we obtain , where and are constants. By using the value of in system (25) and then solving for , we deduce In general, the Type III systems of second-order ODEs are integrable if the associated second-order ODEs in have another symmetry generator besides translations in , so that Lie's classification of scalar second-order ODEs admitting two-dimensional algebras is applicable here as well.

Theorem 8. *If a singular system of two second-order ODEs admits the solvable symmetry algebra , then, it is partially linearizable via point transformation and its general solution can be trivially obtained by quadratures. *

The proof of this follows from the Type I canonical form.

Theorem 9. *If a singular system of two second-order ODEs admits the three-dimensional solvable symmetry algebra , then, it can be reduced to a partially uncoupled system and its general solution can be obtained by quadratures from the general solution of the associated transformed scalar second-order ODE. *

The proof here follows from the Type II canonical form.

For systems that fall in the Type IV category one requires more insight and the examples given in [12] for systems that admit four symmetries become relevant. Their integrability needs more than just the reductions via invariants. This category also has the nonsolvable algebras. We consider this in detail next.

##### 3.2. Differential Invariant Approach

In this approach, we use the basis of invariants representing all the canonical forms which we have given in Table 2. For each regular canonical form, we find the reduced system in terms of a first-order equation and an algebraic equation based on the basis of invariants. These can easily be seen from Table 2.

We have subdivided all the reduced systems in two main categories, namely, *Category A* for those admitting three-dimensional solvable Lie algebras and *Category B* for those admitting three-dimensional nonsolvable Lie algebras. The Types I to III of the previous subsection all fall in *Category A*. The systems admitting nonsolvable algebras comprise *Category B* with the rest firmly falling in *Category A*. It is interesting that, in both categories, uncoupled cases arise when the choice of taking one of the second-order invariant is unique.

We present the integration strategy for *Category A* for regular systems. The systems of *Category B* are not amenable to such integration strategies. This is quite straightforward to work out if one commences with the simplest system in this category. We further subdivide *Category A* into coupled and uncoupled systems. All partially linear cases fall in the uncoupled systems. Here, we give details of integration for each case of *Category A*, which depends on the assumption that we can solve in both the coupled and the uncoupled systems. At the end of this section, we present some physical examples with detailed calculations.

Note here that, in both the uncoupled and coupled cases, for ease and simplicity of the calculations, we only mention the relevant algebra of the considered system and variables and as well as in precise form. The details are given in Table 2. In the following, we assume that is the explicit solution of each of the first order ODEs listed in Table 2. The argument presented below is similar to that used in the discussion of reduction to quadrature for the scalar ODEs (see [1–3]). Moreover, , and are taken as arbitrary functions and and as well as are arbitrary integration constants.

###### 3.2.1. Uncoupled Cases

. Let . Then by using the values of and in Table 2, we have Thus, we can obtain a further quadrature in . Therefore, the equation is solvable by two quadratures upon insertion of the solution.

. Let . Then by invoking the values of and as listed in Table 2, we obtain After some manipulations with , we find By using (30) in (29), we determine which implies that Utilizing (32) in (30), after some further calculations, we deduce Hence, (32) and (33) are the solutions of the system of ODEs admitting .

. Let . By using the values of and as listed in Table 2, we arrive at After calculations with , we find Invoking (35) in (34) and solving for , we obtain which implies that After some computations and using (37) in (35), we determine It is hence the case that (37) and (38) are the solutions of the system of ODEs corresponding to .

. Let . Using the values of and , we obtain which implies that Now utilizing (39) and (40) in , we deduce that which results in Now (40) and (42) are the solutions of the system of ODEs possessing .

. Let . Utilizing the values of and as listed, we obtain After some manipulations with , we find Inverting, we have By using (45) in (43), we have that Thus, (45) and (46) are the solutions of the system of ODEs corresponding to .

For , and , the integration procedure works the same as in the case .

###### 3.2.2. Coupled Cases

. Let . Using the values of and as listed in Table 2, we obtain After calculations with , we arrive at By utilizing (48) in (47), in the same manner as in the previous cases, we find that Hence (48) and (49) are the solutions of the system of ODEs corresponding to .

. The procedure of integration is the same as in the case .

. Let . Using the values of and as listed in Table 2, we determine After some manipulations with , we find From the second equation of the system in Table 2 and some manipulations with (50), we find that Some further computations give which implies that Using (53) and (54) in (51), we find Therefore, (54) and (55) are the solutions of the system of ODEs corresponding to .

. The procedure of integration is the same as the case .

. The procedure of integration is similar to the case .

. Let . Using the values of and as listed in Table 2, we obtain Using the chain rule and some manipulations with , we arrive at Further simplifications give Utilizing (57) and (58) in (56), we find that which implies that Invoking (60) in (58), we have Thus, (60) and (61) are the solutions of the system of ODEs possessing .

. The procedure of integration is similar to the case .

. The procedure of integration is the same as the case .

. Let . Utilizing the values of and as listed in Table 2, we obtain By the chain rule and some simplifications with , we find Substituting (63) in (62) and further manipulations give Invoking this in (63), we deduce that Now after substitution of (65), (64) takes the form Here, (65) and (66) are the solutions of the system of ODEs corresponding to .

. The procedure of integration is the same as the case .

. The procedure of integration is again similar to the case .

. The procedure of integration corresponds to the case and some details are given in the case .

. Let . Using the values of and as listed in Table 2, we obtain After manipulations with , we find Invoking (68) in (67), we determine We have In (68)–(70), Also, we find that Using (69) and (72) in (70), we deduce which shows that By invoking (74) in (68), we have that Thus, (74) and (75) are the solutions of the system of ODEs corresponding to .

In the view of the preceding reductions, we have the following theorem.

Theorem 10. *If a regular system of two second-order ODEs admits a three-dimensional solvable symmetry algebra, then, the general solution of the regular system can be obtained by quadratures from the general solution of the invariant representation of the regular system given by the first-order ordinary differential equation and algebraic equation. *

Theorem 11. *If a system of two th-order ODEs possess a three-dimensional solvable symmetry algebra, then, the general solution of the system can be deduced by quadratures from the general solution of the invariant representation given by one of the following systems:
*

The proof follows from Theorem 2 and the discussion preceding Theorem 10.

###### 3.2.3. Illustrative Examples

Here, we present some familiar physical examples to illustrate the above general integration procedure which we have given in detail for each of the case of system of ODEs admitting solvable three-dimensional Lie algebras. We consider two different cases, namely, one from the uncoupled cases and one from the coupled cases. We also present a case of nonsolvable Lie algebras. Here, we observe interesting results rather than integration.

(i) Consider the system (25) once again from a new angle. We have The system (25) has symmetry generators and its Lie algebra is three-dimensional which is in our classification. From Table 2, the first-order differential invariants for this case are and the reduced system of first-order equation and algebraic equation is The solution of in system (79), after substituting the value of from (25), gives where is constant. Reverting back to the original variables and after solving for , we obtain where and are constants. Substituting the values in and solving for , we find Both (81) and (82) provide the solution of system (25).

(ii) The classical Kepler problem in reduced Cartesian coordinates is where is a constant. In polar coordinates, they are It has the three-dimensional Lie algebra of symmetries [18] spanned by By using the transformation the system (84) transforms into with associated transformed Lie algebra of symmetries which is in our classification. From Table 2, the first-order differential invariants for this case are and the reduced system of first-order equation and algebraic equation is By utilizing the values of and from (87), we have that Now the solution of in the system (91) is In and variables, after some simplification, (92) and from (91) give where and are constants. Reverting back to the original variables, we deduce The term on the left-hand side of (95) is the angular momentum and here it is constant from the calculations and denoted by . Performing a simple integration procedure, we arrive at the equation of the orbit, namely, where and . This corresponds to that presented in [19].

(iii) Consider a case of nonsolvable Lie algebras, namely, the generalized Ermakov system in polar coordinates of the form which was considered in [20]. This system has the three symmetries and forms the Lie algebra sl(2,) which is in our classification. The first-order differential invariants for this case are , and the reduced system of first-order equation and an algebraic equation is . By using these and considering the angular momentum as constant, , we find that . In the original variables, it is where . Here is the Ermakov invariant. The classical integrability was discussed in [20].

#### 4. Conclusion

The integrability of systems of ordinary differential equations (ODEs) which admit symmetry algebras has been of interest in the literature. The traditional symmetry reductions do not work for systems of ODEs as in the case of scalar ODEs (see, e.g., [5, 12, 17]). Here, we have focused on systems of two th-order () ODEs that admit three-dimensional Lie algebras.

By utilizing the realizations given recently, we have provided a systematic approach for finding a complete set of differential invariants including basis of invariants corresponding to vector fields in three variables for systems of two th-order () ODEs which admit three-dimensional Lie algebras. In particular, we presented a complete list of second-order differential invariants for systems of two second-order ODEs that admit three-dimensional Lie algebras. We have shown how invariant differentiations provide higher-order differential invariants. We have given a complete list of invariant characterizations of canonical forms for systems of two second-order ODEs that admit real Lie algebras of dimension three. These are new in the literature. In addition, we have given a procedure for the construction of two th-order ODEs from their associated complete set of invariants. We have shown that there are classes for the case of and classes for the case of of such ODEs. Furthermore, we have given a discussion of those cases in which regular systems of two second-order ODEs are not obtained for three-dimensional Lie algebras. We also have given a brief discussion on the singularity for canonical forms for a system of two second-order ODEs admitting-three dimensional Lie algebras and found two cases of singular invariants, namely, those involving the algebras and .

We have provided an integration procedure for the canonical forms for a system of two second-order ODEs which possess three-dimensional Lie algebras. This procedure is composed of two approaches, namely, one which depends on general observations of canonical forms; the other depends on basis of differential invariants and the solution procedure is applicable only for those cases which admit three-dimensional solvable Lie algebras. The latter was extended in a natural way to th-order () systems that admit three-dimensional Lie algebras. Specifically in the second approach, we obtain the reducibility for canonical forms for systems of two second-order ODEs which have three-dimensional Lie algebras in terms of basis of invariants; these are given in Table 2. Then, we gave details of integration for each case of solvable Lie algebras in a general manner.

We also have obtained two cases of partial linearization, namely, Types I and II. It was noted that the Type I system is trivially integrable. We have further shown how the associated scalar second-order ODEs in the Types II to III can be solved in terms of quadratures resorting to the classical Lie table for integrability of scalar ODEs. This was illustrated by means of an example. There are certain Type IV systems which have complications in their integrability in that they admit nonsolvable algebras.

We have presented familiar physical systems taken from the literature, namely, the classical Kepler problem and the generalized Ermakov systems that give rise to closed trajectories.

It would be of interest to further pursue integrability properties of Type IV systems having nonsolvable Lie algebras. Also worthy of investigation are systems that admit higher number of symmetries such as those studied in [12, 16]. We return to these in the future.

#### Remarks

(i)In the tables, is used as a place holder for the relevant real Lie algebras the th algebra of dimension with superscripts , if any, indicates parameters on which the Lie algebra depends; also the column in the tables gives the algebra realizations; the realization is referred to by a superscript . Further we use the notation . Finally, the elements of a basis of a given Lie algebra are named , where is less than or equal to the dimension of the underlying real Lie algebra. Here, , and are arbitrary functions.(ii)For the analysis of the invariant representations, it is observed that for each of the cases and , with the same set of invariants, a single equation is determined. We find the differential invariants in both cases. However, the canonical form is not obtained. Consider, for example, . For this, we find the following invariants and the corresponding single equation, namely, . For each of the cases and , due to the extra condition, the rank of the coefficient matrix on the solution manifold is less than the rank of generic coefficient matrix. So only a set of singular invariants are obtained for these cases. Consider, for example, . For this, we find the following set of singular invariants (see Table 2). The two cases and are not given in the invariant representations listed in Table 2.(iii)The algebras which are referred to by * give rise to singular invariants when associated with systems of two second-order ODEs as given in Table 2. We do however arrive at regular invariants associated with systems of two third-order ODEs for these cases.(iv)The algebras do not correspond to systems of two second-order ODEs and so are not presented in Table 2.#### Acknowledgments

M. Ayub and M. Khan are grateful to the Higher Education Commission (HEC) of Pakistan for financial support (via Indigenous Scholarship and IRSIP program). Also F. M. Mahomed thanks the HEC for granting him a visiting professorship and the Department of Mathematics at the Quaid-i-Azam University, Islamabad, for its hospitality during the time this work was initiated. M. Ayub also thanks the Centre for Differential Equations, Continuum Mechanics and Applications, School of Computational and Applied Mathematics, University of the Witwatersrand, South Africa, for its hospitality during the time this project was completed. Finally F. M. Mahomed is grateful to the NRF of South Africa for research funding. Parts of this work was presented at the SDEA2012 Jan Conference.

#### References

- S. Lie, “Klassiffikation und Integration von gewönlichen Differentialgleichungen zwischen
*x, y*, die eine Gruppe von Transformationen gestaten,”*Archiv der Mathematik*, vol. 8, no. 9, p. 187, 1883. View at Google Scholar - S. Lie,
*Vorlesungen Über Differentialgleichungen Mit Bekannten Infinitesimalen Transformationen*, B. G. Teubner, Leipzig, Germany, 1891. - F. M. Mahomed and P. G. L. Leach, “Lie algebras associated with scalar second-order ordinary differential equations,”
*Journal of Mathematical Physics*, vol. 30, no. 12, pp. 2770–2777, 1989. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - F. M. Mahomed and P. G. L. Leach, “Symmetry Lie algebras of $n$th order ordinary differential equations,”
*Journal of Mathematical Analysis and Applications*, vol. 151, no. 1, pp. 80–107, 1990. View at Publisher · View at Google Scholar · View at MathSciNet - C. Wafo Soh and F. M. Mahomed, “Reduction of order for systems of ordinary differential equations,”
*Journal of Nonlinear Mathematical Physics*, vol. 11, no. 1, pp. 13–20, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - R. O. Popovych, V. M. Boyko, M. O. Nesterenko, and M. W. Lutfullin, “Realizations of real low-dimensional Lie algebras,”
*Journal of Physics*, vol. 36, no. 26, pp. 7337–7360, 2003. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. Lie,
*Theories der Transformation Gruppen*, Chelsea Publishing, New York, NY, USA, 2nd edition, 1970, Originally published in Leipzig, 1888, 1890, 1893. - S. Lie, “Theorie der transformationsgruppen,”
*Mathematische Annalen*, vol. 16, no. 4, pp. 441–528, 1880. View at Publisher · View at Google Scholar · View at MathSciNet - S. Lie, “Theorie der transformationsgruppen,” in
*Gesammetle Abhandlungen*, vol. 6, pp. 1–94, B. G. Teubner, Leipzig, Germany, 1927. View at Google Scholar - L. Bianchi,
*Lezioni Sulla Teoria dei Gruppi Continui Finiti di Transformazioni*, Enirco Spoerri, Pisa, Italy, 1918. - G. M. Mubarakzjanov, “On solvable Lie algebras,”
*Izvestija Vysših Učebnyh Zavedeniĭ Matematika*, vol. 1, no. 32, pp. 114–123, 1963 (Russian). View at Google Scholar · View at MathSciNet - C. Wafo Soh and F. M. Mahomed, “Canonical forms for systems of two second-order ordinary differential equations,”
*Journal of Physics*, vol. 34, no. 13, pp. 2883–2911, 2001. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. V. Aminova and N. A. M. Aminov, “Projective geometry of systems of differential equations: general conceptions,”
*Tensor*, vol. 62, no. 1, pp. 65–86, 2000. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - P. G. L. Leach, “A further note on the Hénon-Heiles problem,”
*Journal of Mathematical Physics*, vol. 22, no. 4, pp. 679–682, 1981. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - P. G. L. Leach and V. M. Gorringe, “The relationship between the symmetries of and the existence of conserved vectors for the equation $\ddot{x}+f(r)L+g(r)\widehat{r}$ = 0,”
*Journal of Physics*, vol. 23, no. 13, pp. 2765–2774, 1990. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - O. Gaponova and M. Nesterenko, “Systems of second-order ODEs invariant with respect to low-dimensional Lie algebras,”
*Physics AUC*, vol. 16, part 2, pp. 238–256, 2006. View at Google Scholar - P. J. Olver,
*Applications of Lie Groups to Differential Equations*, vol. 107 of*Graduate Texts in Mathematics*, Springer, New York, NY, USA, 2nd edition, 1993. View at Publisher · View at Google Scholar · View at MathSciNet - V. Gorringe, V. M, and P. G. L. Leach, “Kepler's third law and the oscillator isochronism,”
*American Journal of Physics*, vol. 61, pp. 991–995, 1993. View at Google Scholar - N. H. Ibragimov, F. M. Mahomed, D. P. Mason, and D. Sherwell,
*Di Erential Equations and Chaos*, New Age International, New Delhi, India, 1996. - K. S. Govinder and P. G. L. Leach, “Integrability of generalized Ermakov systems,”
*Journal of Physics*, vol. 27, no. 12, pp. 4153–4156, 1994. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet