# Second-order systems of ODEs admitting three-dimensional lie algebras and integrability.

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 kth-order (k [greater than or equal to] 2) 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 kth-order ODEs possessing three-dimensional Lie algebras from the associated complete set of invariants. We show that there are 29 classes for k = 2 and 31 classes for k [greater than or equal to] 3. 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 kth-order (k [greater than or equal to] 3) 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 kth-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 kth-order (k [greater than or equal to] 2) 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

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (1)

admitting the three-dimensional Lie algebra [A.sup.3.sub.3.2]. We have the four independent invariants

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (2)

Here, the invariants [w.sub.3] and [w.sub.4] 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

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (3)

admitting the same Lie algebra [A.sup.3.sub.3,2]; we have the six invariants

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (4)

from which we determine the third-order system. Here, the invariants [w.sub.5] and [w.sub.6] 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 kth-order ODEs admitting three-dimensional Lie algebras, the number of order k invariants is two.

To study differential invariants further, we consider the following example. Consider the system

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (5)

admitting the three-dimensional Lie algebra [A.sup.a,2.sub.3,6], [absolute value of a] [less than or equal to] 1, a [not equal to] 0,1. We arrive at the following set of invariants

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (6)

This set involves three invariants u, v, and w which play the role of basis of invariants for the considered system of second-order ODEs admitting the three-dimensional Lie algebra [A.sup.a,2.sub.3,6], [absolute value of a] [less than or equal to] 1, a [not equal to] 0,1. 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 u and v are invariants of a Lie algebra admitted by any system of ODEs, then, dv/du is also its invariant.

The proof is identical (except here we have many dependent variables) to the scalar case and uses the operator identity [XD.sub.t] - [D.sub.t]X = -([D.sub.t][xi])[D.sub.t] (see, e.g., [17]).

One can write

[d.sub.v]/du = [D.sub.t]v/[D.sub.t]u = Dv, (7)

where the invariant differentiation operator is

D = [([D.sub.t]u).sup.-1] [D.sub.t]. (8)

There arises a natural question. Can we construct a kth-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 [A.sup.a,2.sub.3,6], [absolute value of a] [less than or equal to] 1, a [not equal to] 0,1. 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 [A.sup.a,2.sub.3,6], [absolute value of a] [less than or equal to] 1, a [not equal to] 0,1. We therefore have for this case

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (9)

where [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] and Dw and [D.sup.2]v are third-order invariants. We can thus construct the following third-order system

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (10)

We have constructed a system of two third-order ODEs admitting the three-dimensional Lie algebra [A.sup.a,2.sub.3,6], [absolute value of a] [less than or equal to] 1, a [not equal to] 0,1 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, u, v, w, Dv, Dw, and [D.sup.2]v, where three of them u, v, and w are a basis of invariants and two of them Dw, [D.sup.2]v 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 [A.sup.a,2.sub.3,6], [absolute value of a] [less than or equal to] 1, a [not equal to] 0,1. We can generalize this to a system of two kth-order ODEs. However, there are some nonstandard cases to look at first. These behave differently. There are four cases [A.sup.3.sub.3,1], [A.sup.3.sub.3,5], [A.sup.4.sub.3,1], and [A.sup.1.sub.3,8] in which two of them [A.sup.3.sub.3,1] and [A.sup.3.sub.3,5] yield singular invariants and the other two [A.sup.4.sub.3,1] and [A.sup.1.sub.3,8] give a single equation and do not form a system.

2.1. Singular Case. In both of the cases [A.sup.3.sub.3,1] and [A.sup.3.sub.3,5], 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 [A.sup.3.sub.3,1] [equivalent] [[partial derivative].sub.x], t[[partial derivative].sub.x], y[[partial derivative].sub.x]. We find the following set of singular invariants

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (11)

where D = [D.sub.t] and the corresponding system are

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (12)

In both of the cases, the singularity of the invariants is removed in the third- and higher-order systems. In the case of [A.sup.3.sub.3,1], we get the following complete set of invariants:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (13)

and the corresponding third-order system is

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (14)

In this case, u, v, and w form a basis of invariants, and here the third-order invariants are w and [D.sup.3]v, where D = [D.sub.t]. For the other case [A.sup.3.sub.3,5], it is easy to see that

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (15)

where again D = [D.sub.t]. For each of the remaining cases [A.sup.4.sub.3,1] and [A.sup.1.sub.3,8], we obtain the same set of second-order invariants

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (16)

for which the highest-order invariant is only one and D = [[??].sup.-1][D.sub.t]. So a system of second-order ODEs cannot be constructed. Only a single scalar equation is obtained, namely

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (17)

However, in both of these cases, this behavior changes for higher orders. In fact for [A.sup.4.sub.3,1], we have the third-order invariants

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (18)

and for [A.sup.1.sub.3,8] the third-order invariants are

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (19)

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 [A.sup.4.sub.3,1] and [A.sup.1.sub.3,8] 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 kth-order ODEs (k [greater than or equal to] 3) admitting three-dimensional Lie algebras can be represented by a complete set of invariants which contains 2k complete functionally independent invariants in one of the following forms:

(a) u, v, w, [dv/du], [dw/du], [[d.sup.2]v/d[u.sup.2]], ..., [[d.sup.k-2]w/d[u.sup.k-2]], [[d.sup.k-1]w/d[u.sup.k-1]],

(b) u, v, [dv/du], [[d.sup.2]v/d[u.sup.3]], w, [[d.sup.3]v/d[u.sup.3]], [dw/du], [[d.sup.4]v/d[u.sup.4]], ..., [d.sup.k-3]w/d[u.sup.k-3]], [[d.sup.k]v/d[u.sup.k]]. (20)

In both (a) and (b), u, v, and w form a basis of invariants, where the highest order invariants in (a) are [d.sup.k-2]w/d[u.sup.k-2], [d.sup.k-1]/d[u.sup.k-1] and in (b) they are [d.sup.k-3]w/d[u.sup.k-3], [d.sup.k]v/d[u.sup.k].

The total number of systems of two kth-order ODEs having the above structure is thirty-one, in which form (b) corresponds to the cases [A.sup.3.sub.3,1], [A.sup.3.sub.3,5], [A.sup.4.sub.3,1], and [A.sup.1.sub.3,8] and (a) 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:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (21)

Type II:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (22)

Type III:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (23)

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 y equation for its solution. Further, we observe that the Type III has equations that comprise scalar second-order ODEs in the dependent variable y. 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

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (24)

is partially linearizable by a point transformation as the equation in y is the well-known Painleve-Ince equation which is linearizable via a point transformation.

Let {[X.sub.i]} be a basis of a three-dimensional algebra and let the rank of the coefficient matrix associated with these operators be r.

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 [A.sub.3,1] or [A.sub.3,5] with rank r = 1 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 [A.sup.4.sub.3,2], [A.sup.3.sub.3,4], [A.sup.4.sub.3,5], [A.sup.a,3.sub.3,6], and [A.sup.a,3.sub.3,7] with rank r = 2 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 y 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 y as it already has translations in t symmetry. As an example we consider the system

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (25)

The system (25) has symmetry generators [X.sub.1] = [[partial derivative].sub.x], [X.sub.2] = t[[partial derivative].sub.x], [X.sub.3] = -[[partial derivative].sub.t] + x[[partial derivative].sub.x] and its Lie algebra is three-dimensional. The associated equation in y is [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. Its Lie algebra is two-dimensional spanned by [Y.sub.1] = [[partial derivative].sub.t], [Y.sub.2] = t[[partial derivative].sub.t] - 2y[[partial derivative].sub.y]. By using the transformation

[bar.t] = y, [bar.x] = x, [bar.y] = t, (26)

the associated equation in transformed form becomes [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. We have [[bar.Y].sub.1] = [[partial derivative].sub.[bar.y]], [[bar.Y].sub.2] = [bar.y][[partial derivative].sub.[bar.y]] - 2[bar.t][[partial derivative].sub.[bar.t]], 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 y = 2[c.sub.1.sup.2] [sech.sup.2] [c.sub.1](t - [c.sub.2]), where [c.sub.1] and [c.sub.2] are constants. By using the value of y in system (25) and then solving for x, we deduce

x = 2048 [c.sub.1.sup.18][integral]{[integral] [e.sup.-t] [sech.sup.14] ([c.sub.1] (t - [c.sub.2])) x [tanh.sup.4] ([c.sub.1] (t - [c.sub.2])) dt} dt + [c.sub.3]t + [c.sub.4]. (27)

In general, the Type III systems of second-order ODEs are integrable if the associated second-order ODEs in y have another symmetry generator besides translations in t, 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 [A.sup.3.sub.3,1], 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 systemo f two second-order ODEs admits the three-dimensional solvable symmetry algebra [A.sup.3.sub.3,5], 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 TypeII 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 w 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 dv/du = h 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 u, v, and w as well as dv/du in precise form. The details are given in Table 2. In the following, we assume that v = [lambda](u, [c.sub.1]) is the explicit solution of each of the first order ODEs dv/du = h(u, v) 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, f, g, h, k, [phi], [psi], and [lambda] are taken as arbitrary functions and [c.sub.1], [c.sub.2], and [c.sub.3] as well as [c.sub.4] are arbitrary integration constants.

3.2.1. Uncoupled Cases

[A.sup.4.sub.3,2]. Let v = [lambda](u, [c.sub.1]). Then by using the values of u and v in Table 2, we have

t[??] = [lambda](y, [c.sub.1]). (28)

Thus, we can obtain a further quadrature in y. Therefore, the equation [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] is solvable by two quadratures upon insertion of the y solution.

[A.sup.1.sub.3,3]. Let v = [lambda](u, [c.sub.1]). Then by invoking the values of u and v as listed in Table 2, we obtain

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (29)

After some manipulations with w = f(u, v), we find

[??] = -1/[integral] [lambda](y, [c.sub.1])f(y, [lambda](y, [c.sub.1]))dy + [c.sub.2]. (30)

By using (30) in 29), we determine

[??][lambda] (y, [c.sub.1]){[integral] [lambda](y, [c.sub.1]) f (y, [lambda] (y, [c.sub.1])) dy + [c.sub.2]} = -1, (31)

which implies that

y = [phi](t, [c.sub.1], [c.sub.2], [c.sub.3]). (32)

Utilizing (32) in (30), after some further calculations, we deduce

x = [psi](t, [c.sub.1], [c.sub.2], [c.sub.3], [c.sub.4]). (33)

Hence, (32) and 33) are the solutions of the system of ODEs admitting [A.sup.1.sub.3,3].

[A.sup.2.sub.3,3]. Let v = [lambda](u, [c.sub.1]). By using the values of u and v as listed in Table 2, we arrive at

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (34)

After calculations with w = f(u, v), we find

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (35)

Invoking (35) in (34) and solving for [??], we obtain

[??] = 2/[([integral] f(y, [lambda](y, [c.sub.1]))dy + [c.sub.2]).sup.2] - 2[lambda] (y, [c.sub.1]), (36)

which implies that

y = [phi](t, [c.sub.1], [c.sub.2], [c.sub.3]). (37)

After some computations and using (37) in (35), we determine

x = [psi](t, [c.sub.1], [c.sub.2], [c.sub.3], [c.sub.4]). (38)

It is hence the case that (37) and 38) are the solutions of the system of ODEs corresponding to [A.sup.2.sub.3,3].

[A.sup.3.sub.3,4]. Let v = [lambda](u, [c.sub.1]). Using the values of u and v, we obtain

[??] = [lambda](y, [c.sub.1]), (39)

which implies that

y = [phi](t, [c.sub.1], [c.sub.2]). (40)

Now utilizing (39) and (40) in w = f(u, v), we deduce that

x = [integral]{[integral][e.sup.-t]f ([phi],[[phi].sub.t]) dt}dt + [c.sub.3]t + [c.sub.4], (41)

which results in

x = [psi](t, [c.sub.1], [c.sub.2], [c.sub.3], [c.sub.4]). (42)

Now (40) and (42) are the solutions of the system of ODEs possessing [A.sup.3.sub.3,4].

[A.sup.1.sub.3,5]. Let v = [lambda](u, [c.sub.1]). Utilizing the values of u and v as listed, we obtain

[??] = [lambda](y, [c.sub.1]). (43)

After some manipulations with w = f(u, v), we find

[integral][dy/[c.sub.2] exp{[integral] f(y, [lambda] (u, [c.sub.1]))dy} = t + [c.sub.3]. (44)

Inverting, we have

y = [phi](t, [c.sub.1], [c.sub.2], [c.sub.3]). (45)

By using (45) in (43), we have that

x = [psi](t, [c.sub.1], [c.sub.2], [c.sub.3], [c.sub.4]). (46)

Thus, (45) and (46) are the solutions of the system of ODEs corresponding to [A.sup.1.sub.3,5].

For [A.sup.4.sub.3,5], [A.sup.a,3.sub.3,6], and [A.sup.a,3.sub.3,7], the integration procedure works the same as in the case [A.sup.3.sub.3,4].

3.2.2. Coupled Cases

[A.sup.1.sub.3,1]. Let v = [lambda](u, [c.sub.1]). Using the values of u and v as listed in Table 2, we obtain

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (47)

After calculations with w = f(u, v), we arrive at

x = [phi](t, [c.sub.1], [c.sub.2], [c.sub.3]). (48)

By utilizing (48) in (47), in the same manner as in the previous cases, we find that

y = [psi](t, [c.sub.1], [c.sub.1], [c.sub.3], [c.sub.4]). (49)

Hence (48) and 49) are the solutions of the system of ODEs corresponding to [A.sup.1.sub.3,1].

[A.sup.2.sub.3,1]. The procedure of integration is the same as in the case [A.sup.1.sub.3,1].

[A.sup.1.sub.3,2]. Let v = [lambda](u, [c.sub.1]). Using the values of u and v as listed in Table 2, we determine

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (50)

After some manipulations with w = f(u, v), we find

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (51)

From the second equation of the system in Table 2 and some manipulations with (50), we find that

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (52)

Some further computations give

[??][e.sup.y] = [phi](y + [c.sub.2]), (53)

which implies that

y = h(t, [c.sub.1], [c.sub.2], [c.sub.3]). (54)

Using (53) and (54) in (51), we find

x = k(t, [c.sub.1], [c.sub.2], [c.sub.3], [c.sub.4]). (55)

Therefore, (54) and (55) are the solutions of the system of ODEs corresponding to [A.sup.1.sub.3,2].

[A.sup.2.sub.3,2]. The procedure of integration is the same as the case [A.sup.1.sub.3,5].

[A.sup.3.sub.3,2]. The procedure of integration is similar to the case [A.sup.1.sub.3,3].

[A.sup.3.sub.3,3]. Let v = [lambda](u, [c.sub.1]). Using the values of u and v as listed in Table 2, we obtain

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (56)

Using the chain rule and some manipulations with w = f(u, v), we arrive at

y = [phi](y - [1/[??]], [c.sub.1], [c.sub.2]). (57)

Further simplifications give

[??] = 1/y - [[phi].sup.-1](y, [c.sub.1], [c.sub.2]). (58)

Utilizing (57) and (58) in (56), we find that

[??](y - [[phi].sup.-1] (y, [c.sub.1], [c.sub.2])) [lambda]([[phi].sup.-1](y, [c.sub.1], [c.sub.2]), [c.sub.1]) = 1, (59)

which implies that

y = h(t, [c.sub.1], [c.sub.2], [c.sub.3]). (60)

Invoking (60) in 58), we have

x = k(t, [c.sub.1], [c.sub.2], [c.sub.3], [c.sub.1]). (61)

Thus, (60) and (61) are the solutions of the system of ODEs possessing [A.sup.3.sub.3,3].

[A.sup.1.sub.3,4]. The procedure of integration is similar to the case [A.sup.1.sub.3,3].

[A.sup.2.sub.3,4]. The procedure of integration is the same as the case [A.sup.3.sub.3,3].

[A.sup.2.sub.3,5]. Let v = [lambda](u, [c.sub.1]). Utilizing the values of u and v as listed in Table 2, we obtain

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (62)

By the chain rule and some simplifications with w = f(u, v), we find

x = [psi](y, [c.sub.1], [c.sub.2]). (63)

Substituting (63) in 62) and further manipulations give

y = [psi](x, [c.sub.1], [c.sub.2], [c.sub.3]). (64)

Invoking this in (63), we deduce that

x = h(t, [c.sub.1], [c.sub.2], [c.sub.3], [c.sub.4]). (65)

Now after substitution of (65), (64) takes the form

y = [psi](t, [c.sub.1], [c.sub.2], [c.sub.3], [c.sub.4]). (66)

Here, (65) and (66) are the solutions of the system of ODEs corresponding to [A.sup.2.sub.3,5].

[A.sup.a,1.sub.3,6]. The procedure of integration is the same as the case [A.sup.1.sub.3,3].

[A.sup.a,1.sub.3,6]. The procedure of integration is again similar to the case [A.sup.a,1.sub.3,6].

[A.sup.a,1.sub.3,7]. The procedure of integration corresponds to the case [A.sup.1.sub.3,3] and some details are given in the case [A.sup.a,2.sub.3,7].

[A.sup.a,2.sub.3,7]. Let v = [lambda](u, [c.sub.1]). Using the values of u and v as listed in Table 2, we obtain

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (67)

After manipulations with w = f(u, v), we find

x = tan {[integral] [f(u, [lambda](u, [c.sub.1]))du/[lambda](u, [c.sub.1]) + f(u, [lambda](u, [c.sub.1]))] + [c.sub.2]}. (68)

Invoking (68) in (67), we determine

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (69)

We have

du/dt = [??]([lambda](u, [c.sub.1]) + f(u, [lambda](u, [c.sub.1])))/[lambda](u, [c.sub.1]). (70)

In (68)-(70),

u = y + arctan [??]. (71)

Also, we find that

u = [[phi].sup.-1] (t + [c.sub.3]),

y + arctan [??] = [[phi].sup.-1] (t + [c.sub.3]). (72)

Using (69) and (72) in (70), we deduce

u = [psi](t, [c.sub.1], [c.sub.2], [c.sub.3]), (73)

which shows that

y = h(t, [c.sub.1], [c.sub.2], [c.sub.3]). (74)

By invoking (74) in 68), we have that

x = k(t, [c.sub.1], [c.sub.2], [c.sub.3], [c.sub.4]). (75)

Thus, (74) and (75) are the solutions of the system of ODEs corresponding to [A.sup.a,2.sub.3,7].

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 kth-order (k [greater than or equal to] 3) 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:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (76)

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 anew angle. We have

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (77)

The system (25) has symmetry generators [X.sub.1] = [[partial derivative].sub.x], [X.sub.2] = t[[partial derivative].sub.x], [X.sub.3] = -[[partial derivative].sub.t] + x[[partial derivative].sub.x] and its Lie algebra is three-dimensional which is [A.sup.3.sub.3,4] in our classification. From Table 2, the first-order differential invariants for this case are

u = y, v = [??], (78)

and the reduced system of first-order equation and algebraic equation is

dv/du = q(u, v), w = f(u, v). (79)

The solution of dv/du in system (79), after substituting the value of g from (25), gives

v = [+ or -]u[square root of 2([c.sup.*.sub.1] - u)], (80)

where [c.sup.*.sub.1] is constant. Reverting back to the original variables and after solving for y, we obtain

y = 2[c.sub.1.sup.2] [sech.sup.2][c.sub.1] (t - [c.sub.2]), (81)

where [c.sub.1] = [square root of [c.sup.*.sub.1]/2] and [c.sub.2] are constants. Substituting the values in w and solving for x, we find

x = 2048 [c.sub.1.sup.18] [integral]{[integral] [e.sup.-t] [sech.sup.14] ([c.sub.1] (t - [c.sub.2])) x [tanh.sup.4] ([c.sub.1] (t - [c.sub.2])) dt} dt + [c.sub.3]t + [c.sub.4]. (82)

Both (81) and (82) provide the solution of system (25).

(ii) The classical Kepler problem in reduced Cartesian coordinates is

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (83)

where [mu] is a constant. In polar coordinates, they are

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (84)

It has the three-dimensional Lie algebra of symmetries [18] spanned by

[X.sub.1] = [[partial derivative].sub.t], [X.sub.2] = [[partial derivative].sub.[theta]], [X.sub.3] = t[[partial derivative].sub.t] - [2/3] r[[partial derivative].sub.r]. (85)

By using the transformation

[bar.t] = t, [bar.x] = [theta], [bar.y] = [3/2] ln r (86)

the system (84) transforms into

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (87)

with associated transformed Lie algebra of symmetries

[X.sub.1] = [[partial derivative].sub.[bar.t]], [X.sub.2] = [[partial derivative].sub.[bar.x]], [X.sub.3] = [bar.t] [[partial derivative].sub.[bar.t]] + [[partial derivative].sub.[bar.y]], (88)

which is [A.sup.1.sub.3.2] in our classification. From Table 2, the first-order differential invariants for this case are

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (89)

and the reduced system of first-order equation and algebraic equation is

dv/du = g(u, v), w = f(u, v). (90)

By utilizing the values of f and g from (87), we have that

dv/du = -4v + 9([mu]v/[u.sup.2]) - 9[v.sup.3]/2u - 9([mu]/u) + 9[v.sup.2]u, w = -4/3v. (91)

Now the solution of dv/du in the system (91) is

[u.sup.4][v.sup.2] - [9/2][mu][u.sup.2][v.sup.2] + [9/4] [u.sup.4][v.sup.4] = h. (92)

In [bar.x] and [bar.y] variables, after some simplification, (92) and w from (91) give

d([e.sup.(2/3)[bar.y]])/dt = [+ or -][square root of [4h/9[k.sup.2]] + 2[mu][e.sup.(2/3)[bar.y]] - [e.sup.(4/3)[bar.y]]], [([e.sup.(2/3)[bar.y]).sup.2] [??] = k, (93)

where h and k are constants. Reverting back to the original variables, we deduce

dr/dt = [+ or -][square root of (4h/9[k.sup.2]] + 2[mu]/r - 1/[r.sup.2])],

[r.sup.2][theta] = k. (95)

The term on the left-hand side of (95) is the angular momentum and here it is constant from the calculations and denoted by L/m. Performing a simple integration procedure, we arrive at the equation of the orbit, namely,

d[theta] = ds/[- or +][square root of ([h.sup.*] + 2[m.sup.2][mu]s/[L.sup.2] - [s.sup.2])], (96)

where s([theta]) = 1/r([theta]) and [h.sup.*] = 4h/9[k.sup.4]. 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

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (97)

which was considered in [20]. This system has the three symmetries

[X.sub.1] = [[partial derivative].sub.t], [X.sub.2] = 2t[[partial derivative].sub.t] + r[[partial derivative].sub.t], [X.sub.3] = [t.sup.2][[partial derivative].sub.t] + tr[[partial derivative].sub.r], (98)

and forms the Lie algebra sl(2,R) which is [A.sup.4.sub.3,8] in our classification. The first-order differential invariants for this case are [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], and the reduced system of first-order equation and an algebraic equation is dv/du = g(u, v), w = f(u, v). By using these and considering the angular momentum as constant, L = m[r.sup.2][??], we find that dv/du = -vh(u, L/m)' where g(u, v) = -vh(u, L/m). In the original variables, it is

M = [1/2][([r.sup.2][??]).sup.2] - [integral] h([theta])d[theta], (99)

where h([theta]) = h(u, L/m). Here M 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 kth-order (k [greater than or equal to] 2) 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 kth-order (k [greater than or equal to] 2) 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 kth-order ODEs from their associated complete set of invariants. We have shown that there are 29 classes for the case of k = 2 and 31 classes for the case of k [greater than or equal to] 3 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 [A.sup.3.sub.3,1] and [A.sup.3.sub.3,5].

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 kth-order (k [greater than or equal to] 3)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, A is used as a place holder for the relevant real Lie algebras [A.sup.a,b,n.sub.i,j] ... (the jth algebra of dimension i with superscripts a, b, if any, indicates parameters on which the Lie algebra depends; also the column N in the tables gives the algebra realizations; the realization is referred to by a superscript n, e.g., [A.sup.3.sub.3,1]). Further we use the notation [[partial derivative].sub.t] = [partial derivative]/[partial derivative]t, [[partial derivative].sub.x] = [partial derivative]/[partial derivative]x, [[partial derivative].sub.y] = [partial derivative]/[partial derivative]y. Finally, the elements of a basis of a given Lie algebra are named [X.sub.i], where i is less than or equal to the dimension of the underlying real Lie algebra. Here, f, g, h, k, and [phi] are arbitrary functions.

(ii) For the analysis of the invariant representations, it is observed that for each of the cases [A.sup.4.sub.3,1] and [A.sup.1.sub.3,8], 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, [A.sup.4.sub.3,1] [equivalent] [[partial derivative].sub.t], x[[partial derivative].sub.t], [phi]"(x) [not equal to] 0. For this, we find the following invariants [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] and the corresponding single equation, namely, [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. For each of the cases [A.sup.3.sub.3,1] and [A.sup.3.sub.3,5], 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, [A.sup.3.sub.3,1] [equivalent] [[partial derivative].sub.x], t[[partial derivative].sub.x], y[[partial derivative].sub.x]. For this, we find the following set of singular invariants u = t, v = y, [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], and dv/du = [??] = 0 (see Table 2). The two cases [A.sup.4.sub.3,1] and [A.sup.1.sub.3,8] are not given in the invariant representations listed in Table 2.

(iii) The algebras [A.sup.3.sub.3,1] and [A.sup.3.sub.3,5] 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 [A.sup.4.sub.3,1] and [A.sup.1.sub.3,8] do not correspond to systems of two second-order ODEs and so are not presented in Table 2.

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

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 Quaidi-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

[1] S. Lie, "Klassifikation und Integration von gewonlichen Differentialgleichungen zwischen x, y, die eine Gruppe von Transformationen gestaten," Archiv der Mathematik, vol. 8, no. 9, p. 187, 1883.

[2] S. Lie, Vorlesungen Uber Differentialgleichungen Mit Bekannten Infinitesimalen Transformationen, B. G. Teubner, Leipzig, Germany, 1891.

[3] 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.

[4] F.M. Mahomed and P.G.L. Leach, "Symmetry Lie algebras of nth order ordinary differential equations," Journal of Mathematical Analysis and Applications, vol.151, no.1, pp.80-107,1990.

[5] 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.

[6] 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.

[7] S. Lie, Theories der Transformation Gruppen, Chelsea Publishing, New York, NY, USA, 2nd edition, 1970, Originally published in Leipzig, 1888, 1890, 1893.

[8] S. Lie, "Theorie der transformationsgruppen," Mathematische Annalen, vol.16, no.4, pp.441-528,1880.

[9] S. Lie, "Theorie der transformationsgruppen," in Gesammetle Abhandlungen, vol.6, pp.1-94, B.G. Teubner, Leipzig, Germany, 1927.

[10] L. Bianchi, Lezioni Sulla Teoria dei Gruppi Continui Finiti di Transformazioni, Enirco Spoerri, Pisa, Italy, 1918.

[11] G. M. Mubarakzjanov, "On solvable Lie algebras," Izvestija Vyssih Ucebnyh Zavedenii Matematika, vol.1, no. 32, pp. 114-123, 1963 (Russian).

[12] 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.

[13] 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.

[14] P. G. L. Leach, "A further note on the Heenon-Heiles problem," Journal of Mathematical Physics, vol.22, no.4, pp.679-682, 1981.

[15] P.G.L. Leach and V.M. Gorringe, "The relationship between the symmetries of and the existence of conserved vectors for the equation [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]," Journal of Physics, vol. 23, no.13, pp.2765-2774,1990.

[16] 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.

[17] 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.

[18] 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.

[19] N. H. Ibragimov, F. M. Mahomed, D. P. Mason, and D. Sherwell, Di Erential Equations and Chaos, New Age International, New Delhi, India, 1996.

[20] K. S. Govinder and P. G. L. Leach, "Integrability of generalized Ermakov systems," Journal of Physics, vol.27, no.12, pp.4153-4156, 1994.

Muhammad Ayub, (1,2) Masood Khan, (1) and F.M. Mahomed (2)

(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

Correspondence should be addressed to F. M. Mahomed; fazal.mahomed@wits.ac.za

Received 27 November 2012; Accepted 20 January 2013

```
TABLE 1: Realizations of three-dimensional Lie algebras
in three variables.

Algebra                N   Realization                         Rank

[A.sub.3,1]            1   [[partial derivative].sub.t],        3
[[partial derivative].sub.x],
[[partial derivative].sub.y]

2   [[partial derivative].sub.x],        2
[[partial derivative].sub.y],
t[[partial derivative].sub.x] +
f(t)[[partial
derivative].sub.y],

3   [[partial derivative].sub.x],        1
t[[partial derivative].sub.x],
y[[partial derivative].sub.x]

4   [[partial derivative].sub.t],        1
x[[partial derivative].sub.t],
[phi](x)[[partial
derivative].sub.t], [phi]" (x)
[not equal to] 0

[A.sub.3,2]            1   [[partial derivative].sub.t],        3
[[X.sub.1],                t[[partial derivative].sub.t] +
[X.sub.2]] = X1            [[partial derivative].sub.y],
[[partial derivative].sub.x]

2   [[partial derivative].sub.t],        2
t[[partial derivative].sub.t] +
y[[partial derivative].sub.x],
[[partial derivative].sub.x]

3   [[partial derivative].sub.t],        2
t[[partial derivative].sub.t],
[[partial derivative].sub.x]

4   [[partial derivative].sub.x],        2
t[[partial derivative].sub.t] +
x[[partial derivative].sub.x],
t[[partial derivative].sub.x]

[A.sub.3,3]            1   [[partial derivative].sub.t],        2
[[X.sub.2],                [[partial derivative].sub.x],
[X.sub.3]] =               x[[partial derivative].sub.t]
[X.sub.1]

2   [[partial derivative].sub.t],        2
[[partial derivative].sub.x],
x[[partial derivative].sub.t] +
y[[partial derivative].sub.x]

3   [[partial derivative].sub.t],        3
[[partial derivative].sub.x],
x[[partial derivative].sub.t] +
[[partial derivative].sub.y]

[A.sub.3,4]            1   [[partial derivative].sub.t],        2
[[X.sub.1],                [[partial derivative].sub.x], (t
[X.sub.3]] =               + x)[[partial derivative].sub.t]
[X.sub.1]                  + x[[partial derivative].sub.x]
[[X.sub.2],
[X.sub.3]] =           2   [[partial derivative].sub.t],        3
[X.sub.1] +                [[partial derivative].sub.x], (t
[X.sub.2]                  + x)[[partial derivative].sub.t]
+ x[[partial derivative].sub.x]
+ [[partial derivative].sub.y]

3   [[partial derivative].sub.x],        2
t[[partial derivative].sub.x],
-[[partial derivative].sub.t] +
x[[partial derivative].sub.x]

[A.sub.3,5]            1   [[partial derivative].sub.t],        2
[[X.sub.1],                [[partial derivative].sub.x],
[X.sub.3]] =               t[[partial derivative].sub.t] +
[X.sub.1]                  x[[partial derivative].sub.x]
[[X.sub.2],
[X.sub.3]] =           2   [[partial derivative].sub.t],        3
[X.sub.2]                  [[partial derivative].sub.x],
t[[partial derivative].sub.t] +
x[[partial derivative].sub.x] +
[[partial derivative].sub.y]

3   [[partial derivative].sub.x],        1
t[[partial derivative].sub.x],
x[[partial derivative].sub.x]

4   [[partial derivative].sub.x],        2
t[[partial derivative].sub.x],
x[[partial derivative].sub.x] +
[[partial derivative].sub.y]

[A.sup.a.sub.3,6],     1   [[partial derivative].sub.t],        2
[absolute value of         [[partial derivative].sub.x],
a] [less than or           t[[partial derivative].sub.t] +
equal to] 1, a [not        ax[[partial derivative].sub.x]
equal to] 0, 1
[[X.sub.1],
[X.sub.3]] =           2   [[partial derivative].sub.t],        3
[X.sub.1]                  [[partial derivative].sub.x],
[[X.sub.2],                t[[partial derivative].sub.t] +
[X.sub.3]] =               ax[[partial derivative].sub.x] +
a[X.sub.2]                 [[partial derivative].sub.y]

3   [[partial derivative].sub.x],        2
t[[partial derivative].sub.x],
(1 - a)t[[partial
derivative].sub.t] + x[[partial
derivative].sub.x]

[A.sup.a.sub.3,7], a   1   [[partial derivative].sub.t],        2
[greater than or           [[partial derivative].sub.x],
equal to] 0                (at + x)[[partial
[[X.sub.1],                derivative].sub.t] + (ax -
[X.sub.3]] =               t)[[partial derivative].sub.x]
a[X.sub.1] -
[X.sub.2]              2   [[partial derivative].sub.t],        3
[[X.sub.2],                [[partial derivative].sub.x],
[X.sub.3]] =               (at + x)[[partial
[X.sub.1] +                derivative].sub.t] + (ax -
a[X.sub.2]                 t)[[partial derivative].sub.x] +
[[partial derivative].sub.y]

3   [[partial derivative].sub.x],        2
t[[partial derivative].sub.x],
-(1 + t2)[[partial
derivative].sub.t] + (a - t)
x[[partial derivative].sub.x]

[A.sub.3,8], {s1(2,    1   [[partial derivative].sub.t],        1
R)}                        t[[partial derivative].sub.t],
[[X.sub.1],                [t.sup.2][[partial
[X.sub.2]] =               derivative].sub.t]
[X.sub.1]
[[X.sub.2],            2   [[partial derivative].sub.t],        3
[X.sub.3]] =               t[[partial derivative].sub.t] +
[X.sub.3]                  x[[partial derivative].sub.x],
[[X.sub.3],                [t.sup.2][[partial
[X.sub.1]] = -             derivative].sub.t] +
2[X.sub.2]                 2tx[[partial derivative].sub.x]
+ x[[partial derivative].sub.y]

3   [[partial derivative].sub.t] +       2
[[partial derivative].sub.x],
t[[partial derivative].sub.t] +
x[[partial derivative].sub.x],
[t.sup.2][[partial
derivative].sub.t] +
[x.sup.2][[partial
derivative].sub.x]

4   -t[[partial derivative].sub.x],      2
[1/2] (-t[[partial
derivative].sub.t] + x[[partial
derivative].sub.x]), x[[partial
derivative].sub.t]

5   [[partial derivative].sub.t],        2
t[[partial derivative].sub.t] +
x[[partial derivative].sub.x],
([t.sup.2] - [x.sup.2])
[[partial derivative].sub.t] +
2tx[[partial derivative].sub.x]

[A.sub.3,9], {so(3)}   1   (1 + [t.sup.2])[[partial             2
[[X.sub.1],                derivative].sub.t] + xt[[partial
[X.sub.2]] =               derivative].sub.x], x[[partial
[X.sub.3]                  derivative].sub.t] -t[[partial
[[X.sub.3],                derivative].sub.x], -
[X.sub.1]] =               xt[[partial derivative].sub.t] -
[X.sub.2]                  (1 + [x.sup.2])[[partial
[[X.sub.2],                derivative].sub.x]
[X.sub.3]] =           2   -sin t tan x[[partial                3
[X.sub.1]                  derivative].sub.t] -cos
t[[partial derivative].sub.x] +
sin t sec x[[partial
derivative].sub.y], -cos t tan
x[[partial derivative].sub.t] +
sin t[[partial
derivative].sub.x] + cos t sec
x[[partial derivative].sub.y],
[[partial derivative].sub.t]
```