Abstract

We present geometric based methods for solving systems of discrete or difference equations and introduce a technique for finding conservation laws for such systems.

1. Introduction

Depending on a model being studied, some physical laws may be described by differential equations (DEs). Lie group theory provides us with powerful tools for obtaining analytical solutions of such equations [1]. Over the last 30 years, a considerable amount of work has been invested into applying Lie’s theory to solve and classify difference equations (see [26] and references therein).

The use of symmetry methods for ordinary difference equations has been introduced by Maeda [2]. He showed that the resulting linearized symmetry condition (LSC) amounts to a set of functional equations which is hard to solve in general. Hydon introduced a technique to solve LSC and obtain symmetries in closed form by repetitive differentiations [5]. For example, a full classification of second-order according to their point symmetries exists in the literature [6]. However, for systems of difference equations , most of results are based on induction methods (see [79] and references therein). In this paper, we present a method for solving using their underlying symmetry.

It has been proved in [5] that every second-order linear homogeneous has an eight-dimension Lie algebra isomorphic to . This is not valid for second-order system of difference equations . We shall prove this in Section 3.

2. Groundwork

Let us consider an -th order system of We assume that for each there exists at least one such that .

Consider a point transformationwhere are continuous variables. will be called one-parameter Lie group of transformations if it satisfies the following properties:(i) is the identity map, i.e., for (ii) for every and close to 0.(iii)Each can be expanded as a Taylor series in a neighbourhood of .

Therefore, we havewhere are continuous functions which we shall refer to as characteristics, , and is the “shift” operator. It is defined as follows:We define the discrete differentiation operator as follows:where is the identity operator.

The symmetry condition for the (1) iswhenever (1) holds.

Lie symmetries are obtained by linearizing the symmetry condition (6) about the identity. We have the following system of linearized symmetry condition (SLSC):where the symmetry generator is given by

Definition 1. A function is invariant function under the Lie group of transformations if

where can be found by solving the characteristic equation

Theorem 2. The discrete differential operator in (5) and the generator of symmetry in (8) commute.

Proof. We prove the theorem for any generalisation is straightforward.

Corollary 3. For each invariant , is also an invariant.

Proof. We have .
Equivalently,
or .

We shall use this corollary for reductions in Section 3.

A first integral for the system (1) is a quantity such thatwhenever (1) holds.

In Section 5, we shall use the condition (12) to develop a constructive technique for obtaining first integrals.

Remark 4. In this paper, we shall consider Lie point symmetry; i.e, the characteristics are given by .

We refer the reader to [1] for more information on symmetry methods for differential equations.

3. Symmetries and Reductions

3.1. Finding Characteristics

Consider a second-order system of 2 EsWe assume that or and or , so the system is of second order.

The SLCS (7) reduces towhere , and .

The functional equations (14) and (15) contain functions and with different pairs of arguments making them difficult to solve. For concreteness, if, for instance, the discrete variable stands for “state” in physics, and belong to two different states.

To overcome this, we proceed as follows.

Step 1 (elimination of and ). We differentiate (total differentiation) (14) and (15) with respect to and , respectively, keeping and fixed. Here, we take as function of and as function of .

The total derivative operators are given byIn this case, this is simplified toSo we apply the operator (17) to (14) and (18) to (15) keeping and fixed. This leads to the determining system

Step 2 (elimination of and ). We now differentiate (19) and (20) with respect to and , respectively, keeping and fixed. This means that we apply the operator on (19) and on (20). For a second-order we need at most to differentiate four times. After separating with respect to and the resulting equations, we obtain a system of DEs in and which is solvable by hand or by using a computer algebra package.

Step 3 (explicit form of constant of integration). When integrating in Step 2 to obtain the characteristics and , we have constant of integration which appears to be functions of . To obtain their explicit form, we need to substitute the results obtained in Step 2 in (19) and (20). If we do not succeed in obtaining all the constant of integration, we need further substitution in the SLSC (14) and (15).

3.2. Reductions

Consider a second-order and its symmetry generatorThe method of characteristics for partial differential equations (PDEs)leads to three independent constants of integration . Each invariant under is function of those constant. .

For second-order systems, two invariants suffice to do reduction of the systems.

Letbe the invariants functions under . We choose them in a way that the Jacobian is nonzero.That is, (24) can be inverted as follows:By Corollary 3, and are also invariant functions. Therefore, the solution of (21) satisfies(27) is a first-order SEs which can be solved by further reductions or by using computer algebra software (maple, Mathematica,…) for linear systems. Note that there exist some first-order systems which cannot be solved analytically.

The general solution isfor some constant .

So the second-order system (21) is equivalent to the first-order system obtained by substituting (28) in (26):(29) also admits the symmetries generated by . The best way to integrate any first-order analytic E is to use its canonical coordinates [10],which satisfyThe obvious choice of canonical coordinates is (see [10])

3.3. Applications
3.3.1. Example 1

Consider the most general homogeneous second-order linear system of difference equationswhere are arbitrary functions.

One can readily verify that the determining system (19) and (20) amounts toThereforewhere are constants.

The characteristics in (35) must satisfy the SLSC (14) and (15). Hence, we haveandSo, (35) is simplified toThe first generator of symmetry for a second-order homogeneous linear system (36) is the scaling symmetry given byThe system (36), which governs the remaining generators of the Lie point symmetry for the system (33), is of second order in and . Its general solution iswhere are constants.

So, the most large Lie algebra of symmetry generators which can be obtained from a homogeneous second-order system of 2 difference equations has dimension five.

For clarification let us consider

and . The system (33) becomes The system which governs the remaining generators of the Lie point symmetry in this case is given byThe general solutions for this system will beTherefore we have 5 generators of the Lie point symmetry spanned by

3.3.2. Example 2

Consider the system(45) is a special case of systems investigated in [11], where the author looked at the stability of the systems.

We choose the ansatz .

The determining system (19) and (20) amounts toDifferentiating twice (46) with respect to and twice (47) with respect to keeping and fixed we obtain, after separating with respect to and , the following system of Des:whose most general solutions areTo obtain the nature of functions we substitute (49) in (46) and (47). After separating with respect to and we get the following SEs:whose solutions areThe remaining unknown functions and are determined by substituting (51) and (49) into the SLSC (14) and (15). This leads to the SEsThe general solutions to (52) are given bywhere are arbitrary constants. It follows that the characteristics are given byTherefore, we have six generators of Lie point symmetryEach generator in (55) can be used to reduce the order of (45).

Let us consider . By the characteristic method for Partial Differential Equations, the invariants are given by following equation:We getwhere are constants.

If we choose , we haveand if we choose , we haveFrom (58) and (59)

we deduceLet us now consider the generator . The resulting invariants arewhereNote also the relationship between themFrom (61), we deduce the following relation:One can readily check that the general solution to (64) is given bywhere and are defined in (62).

From (61) we obtainwhich is a first-order system after substitution of by the results given in (65). Its solutions can be obtained by using the following canonical coordinates:This leads to the the following linear system with variable coefficients:where and .

The latter is a linear first-order system with variable coefficients. Its general solution isThe general solution of (45) is obtained by substituting (69) into (67).

4. Conservation Laws

In Section 2, we have defined a first integral associated with a second-oreder SEs. It is given by (12)LetBy differentiating (70) with respect to and we obtainandThe substitution of (72) in (73) leads to the following second-order system of functional equations:As for SLSC, we differentiate repeatedly to obtain a system of DEs for and . Given the solutions of (74) we easily construct . For consistency of our solutions, we must check the integrability conditionsandThe first integral is then given byThe constant of integration which is a function depending on is determined by substituting (77) in (70).

4.1. Applications

Let us consider the second-order SEsBy choosing the ansatz and one can readily check that the determining system (74) is simplified towhere and denote the right-hand side of (78).

Differentiating (79) with respect to and leads toThus, we have from (72)Substituting (80) in (79) and separating with respect to and we obtain the systemThe solutions to (82) will provide us with the explicit form of

The first integral is then given byfor some constants .

For clarification, let us consider , that is,The solutions to (82) will bewhere are constants. We have twelve solutions for and . That is,(1)If then(2)If then(3)If then(4)Ifthen(5)Ifthen(6)If then(7)If then(8)If then(9)Ifthen(10)Ifthen(11)Ifthen(12)IfthenTherefore, we obtain twelve conservation laws for the system (84). They are given by

5. Conclusion and Discussions

We have presented a method for obtaining nontrivial symmetries and how to use them for solving a second-order SEs. Each symmetry can be used to reduce the order. However, different symmetries lead to different reductions (see (60) and (64)), but the same solution. We also proposed a technique to construct first integral associated to second-order systems of difference equations.

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.