Darboux Integrability of Trapezoidal $H^{4}$ and $H^{6}$ Families of Lattice Equations II: General Solutions

In this paper we construct the general solutions of two families of quad-equations, namely the trapezoidal $H^{4}$ equations and the $H^{6}$ equations. These solutions are obtained exploiting the properties of the first integrals in the Darboux sense, which were derived in [Gubbiotti G., Yamilov R.I., J. Phys. A: Math. Theor. 50 (2017), 345205, 26 pages, arXiv:1608.03506]. These first integrals are used to reduce the problem to the solution of some linear or linearizable non-autonomous ordinary difference equations which can be formally solved.


Introduction
Since its introduction the integrability criterion denoted consistency around the cube (CAC) has been a source of many results in the classification of quad-equations. We define a quad-equation to be a relation of the form: transforms are associated with both linearizable and integrable equations. We point out that to be bona fide a Lax pair has to give rise to a genuine spectral problem [13], otherwise the Lax pair is called fake Lax pair [11,12,14,28,29]. A fake Lax pair is useless in proving (or disproving) the integrability, since it can be equally found for integrable and non-integrable equations. In the linearizable case Lax pairs must be then fake ones, even though proving this statement is usually a nontrivial task [22]. For a complete, pedagogical explanation of the CAC method we refer to [6,32,33].
x x 1 x 2 x 3 x 12 x 13 x 23 x 123 Being algorithmically applicable the CAC criterion proved to be a well suited method to find and classify integrable quad-equations. The first attempt to classify, with some additional assumptions, all the quad-equations possessing CAC was carried out in [1]. The result was the existence of three classes of discrete autonomous equations with this property: the H equations, the Q equations and the A. The A equations can be transformed in particular cases of the Q equations through non-autonomous Möbius transformation, therefore they are usually removed from the general classification. Releasing one of the technical hypothesis of [1], i.e., that face of the cube (Fig. 2) carries the same equation, the same authors in [2] presented some new equations without classification purposes. A complete classification in this extended setting was then accomplished by R. Boll in a series of papers culminating in his Ph.D. Thesis [7,8,9]. In these papers the classification of all the consistent sextuples of quad-equations. The only technical assumption used in [7,8,9] is the tetrahedron property, i.e., the requirement that x 123 is independent from x. The obtained equations may fall into three disjoint families depending on their bi-quadratics h ij = ∂Q ∂y k ∂Q ∂y l − Q ∂ 2 Q ∂y k ∂y l , Q = Q(y 1 , y 2 , y 3 , y 4 ), where we use a special notation for variables of Q, and the pair {k, l} is the complement of the pair {i, j} in {1, 2, 3, 4}. A bi-quadratic is called degenerate if it contains linear factors of the form y i − c, where c is a constant, otherwise a bi-quadratic is called non-degenerate. The three families are classified depending on how many bi-quadratics are degenerate: • Q-type equations: all the bi-quadratics are non-degenerate, • H 4 -type equations: four bi-quadratics are degenerate, • H 6 -type equations: all of the six bi-quadratics are degenerate.
Let us notice that the Q family is the same as the one introduced in [1]. The H 4 equations are divided into two subclasses: rhombic and trapezoidal, depending on their discrete symmetries. We remark that all classification results hold locally in the sense that they relate to a single quadrilateral cell or a single cube. The extension on the whole lattice Z 2 is obtained through reflection considering an elementary cell of size 2×2. This implies that the H 4 and H 6 equations as lattice equations are non-autonomous equations with two-periodic coefficients. For more details on the construction of equations on the lattice from the single cell equations, we refer to [7,8,9,44] and to the Appendix in [21]. A detailed study of all the lattice equations derived from the rhombic H 4 family, including the construction of their three-leg forms, Lax pairs, Bäcklund transformations and infinite hierarchies of generalized symmetries, has been presented in [44]. So there was plenty of results about the Q and the rhombic H 4 equations. On the contrary, besides the CAC property little was known about the integrability features of the trapezoidal H 4 equations and of the H 6 equations. Therefore these equations where thoroughly studied in a series of papers [21,22,23,24,25] with some unexpected results. First in [21] was presented their explicit non-autonomous form.
Remark 1.1. In [35] it was shown that sometimes it is possible to construct different consistent embedding in the Z 2 and in Z 3 lattices. However, in the same paper it was shown that these different embedding need not to be integrable. In this paper we will consider equations (1.2) and (1.3) which are given by the embedding procedure of [2,7,8,9]. As we underlined above this procedure gives equations which, in the sense of the algebraic entropy, are only integrable or linearizable [21,41]. Clearly, it may exist a different embedding in the Z 2 for which the results presented in this paper do not hold. For an example where two different embedding give both rise to linearizable equations, but with different properties, see [20].
In [22] the t H ε 1 equation (1.2a) was studied and it was found that it possessed three-point generalized symmetries depending on arbitrary functions. This property was then linked in [23] to the fact that the t H ε 1 is Darboux integrable [3]. We say that a quad-equation on the Z 2 lattice, possibly non-autonomous: Q n,m (u n,m , u n+1,m , u n,m+1 , u n+1,m+1 ) = 0, (1.5) is Darboux integrable if there exist two independent first integrals, one containing only shifts in the first direction and the other containing only shifts in the second direction. This means that there exist two functions where l 1 < k 1 and l 2 < k 2 are integers, such that the relations hold true identically on the solutions of (1.5). By T n , T m we denote the shift operators in the first and second directions, i.e., T n h n,m = h n+1,m , T m h n,m = h n,m+1 , and by Id we denote the identity operator Id h n,m = h n,m . The number k i − l i , where i = 1, 2, is called the order of the first integral W i . In addition to this result concerning the t H ε 1 equation in [23] it was proved that other quadequations consistent around the cube, which were known to be linearizable [30,31], were in fact Darboux integrable. These facts provide some evidence of an intimate connection between linearizable equations possessing CAC and Darboux integrability. Following these ideas in [26] it was shown that all the trapezoidal H 4 equations and all the H 6 equations are Darboux integrable. This result was proved by explicitly constructing the first integrals with a new algorithm based on those proposed in [17,18,27]. This new algorithm relies on the fact that in the case of non-autonomous quad-equations (1.5) with two-periodic coefficients we can, in general, represent the first integrals in the form In the final section of [26] it was shown, in the case of the t H ε 1 equation, how it is possible to find the general solution using the first integrals, applying a modification of the procedure presented in [17]. In particular we showed how it is possible to obtain a general solution using first integrals of order greater than one. We note that equations with first integrals of first order one are trivial, since possessing a first integral of order one means that the equation itself is a first integral.
In this paper we show that from the knowledge of the first integrals and from the properties of the equations it is possible to construct, maybe after some complicate algebra, the general solutions of all the remaining trapezoidal H 4 equations (1.2) and of the H 6 equations (1.3). By general solution we mean a representation of the solution of any of the equations in (1.2) and (1.3) in terms of the right number of arbitrary functions of one lattice variable n or m. Since the trapezoidal H 4 equations (1.2) and the H 6 equations (1.3) are quad-equations, i.e., the discrete analogue of second-order hyperbolic partial differential equations, the general solution must contain an arbitrary function in the n direction and another one in the m direction, i.e., a general solution is an expression of the form u n,m = F n,m (a n , α m ), (1.9) where a n and α m are arbitrary functions of their variable. Initial conditions are then imposed through substitution in equation (1.9). Nonlinear equations usually possesses also other kinds of solutions, namely the singular solutions which satisfy only specific set of initial values. In this work we outlined when the existence of singular solutions is possible. Moreover we remark that general solutions, in the range of validity of their parameters, enclose also periodic solutions.
Periodic initial values will reflect into periodic solution which will arise by fixing properly the arbitrary functions. Indeed let us consider as an example the (N, −M ) reduction of a quadequation (1.5), with N, M ∈ N * coprime [39,40]. This implies to make the following requirement The existence of the associated periodic solution is subject to the ability to invert formula (1.11). When the integers N and M are not coprime a similar reasoning can be done: taking K = gcd(N, M ) we have just to decompose the reduction condition into K superimposed staircases and convert the scalar condition (1.10) to a vector condition for K fields. The associated reduction will be possible if the associated system possesses a solution.
To obtain the desired solution we will need only the W 1 integrals derived in [26] and the fact that the relation (1.6b) implies W 1 = ξ n with ξ n an arbitrary function of n. The equation W 1 = ξ n can be interpreted as an ordinary difference equation in the n direction depending parametrically on m. Then from every W 1 integral we can derive two different ordinary difference equations, one corresponding to m even and one corresponding to m odd. In both the resulting equations we can get rid of the two-periodic terms by considering the cases n even and n odd and defining This transformation brings both equations to a system of coupled difference equations. This reduction to a system is the key ingredient in the construction of the general solutions for the trapezoidal H 4 equations (1.2) and for the H 6 equations (1.3).
We note that the transformation (1.12) can be applied to the trapezoidal H 4 equations 1 and H 6 equations themselves. This casts these non-autonomous equations with two-periodic coefficients into autonomous systems of four equations. We recall that in this way some examples of direct linearization (i.e., without the knowledge of the first integrals) were produced in [21]. Finally we note that if we apply the even/odd splitting of the lattice variables given by equation (1.12) to describe a general solution we will need two arbitrary functions in both directions, i.e., we will need a total of four arbitrary functions.
In practice to construct these general solutions, we need to solve Riccati equations and nonautonomous linear equations which, in general, cannot be solved in closed form. Using the fact that these equations contain arbitrary functions we introduce new arbitrary functions so that we can solve these equations. This is usually done reducing to total difference, i.e., to ordinary difference equations which can be trivially solved. Let us assume we are given the difference equation depending parametrically on another discrete index m. Then if we can express the function f n as a discrete derivative then the solution of equation (1.13) is simply where γ m is an arbitrary function of the discrete variable m. This is the simplest possible example of reduction to total difference. The general solutions will then be expressed in terms of these new arbitrary functions obtained reducing to total differences and in terms of a finite number of discrete integrations, i.e., the solutions of the simple ordinary difference equation where u n is the unknown and f n is an assigned function. We note that the discrete integration (1.14) is the discrete analogue of the differential equation u (x) = f (x). To give a very simple example of the method of solution we consider how it applies to the prototypical Darboux integrable equation: the discrete wave equation It is easy to check that the discrete wave equation (1.15) is Darboux integrable with two firstorder first integrals From the first integrals (1.16) it is possible to construct the well known discrete d'Alembert solution as follows. From the W 1 first integral we can write W 1 = ξ n with ξ n arbitrary function of its argument. Then we have u n+1,m − u n,m = ξ n . (1.17) This means that choosing the arbitrary function as ξ n = a n+1 − a n , with a n arbitrary function of its argument, we transform (1.17) into the total difference u n+1,m + a n+1 = u n,m + a n , which readily implies where α m is an arbitrary function of its argument. This is of course the discrete analog of the d'Alembert solution of the wave equation and it is the simplest example of solution through the first integrals of a Darboux integrable equation. Now to summarize, in this paper we prove the following result:  [26]. In particular in Section 2.1 we treat the 1 D 2 , 2 D 2 and 3 D 2 equations. In Section 2.2 we treat the D 3 , 1 D 4 and 2 D 4 equations. In Section 2.3 we treat the t H ε 2 and the t H ε 3 equations. The partition in subsection is dictated by the procedure used to obtain the general solution, as we will explain below. Due to the technical nature of the procedures we will present only one example per type. The interested reader will find the remaining procedures of solution in Appendix A. In Section 3 we give some conclusions. Remark 1.3. We remark that the H equations of the ABS classification [1] and their rhombic deformations [2,7,44] should not be Darboux integrable. This can be confirmed directly excluding the existence of integrals up to a certain order as it was done in [17] for some other equations. Moreover it was proved rigorously in [41] using the gcd-factorization method that all the equations of the ABS list [1] possess quadratic growth of the degrees. At heuristic level a similar result was presented in [21] for the rhombic H 4 equations. According to the Algebraic entropy conjecture these result means that the ABS equations and the rhombic H 4 equations are integrable, but not linearizable. Since we have recalled the fact that Darboux integrability for lattice equations implies linearizability we expect that these equations will not possess first integrals of any order. So the results obtained in [26] and in this paper about the trapezoidal H 4 and H 6 equations do not imply anything for H equations and their rhombic deformations.

General solutions of the H 4 and H 6 equations
In this section we present the general solutions of the H 4 equations (1.2) and of the H 6 equations (1.3). We choose to divide this section in three subsections since we have three main different kinds of procedures leading to three different representations of the solution. First in Section 2.1 we present the general solutions of the 1 D 2 , 2 D 2 and 3 D 2 equations (1.3a)-(1.3c). In this case the construction of the general solution is carried out from the sole knowledge of the first integral and the equation acts only as a compatibility condition for the arbitrary functions obtained by solving the equations defined by the first integral. The solution therein obtained is completely explicit and no discrete integration is required.
The in Section 2.2 we present the general solution of the D 3 , 1 D 4 and 2 D 4 equations (1.3d)-(1.3f). In this case the construction of the general solution is carried out through a series of manipulations in the equation itself and from the knowledge of the first integral. The key point will be that the equations defined by the first integrals can be reduced to a single linear equation.
The solution is no longer completely explicit since it is obtained up to two discrete integrations, one in every direction.
Finally in Section 2.3 we present the general solutions of the t H ε 2 equation (1.2b) and of the t H ε 3 equation (1.2c). In this case the construction of the general solution is carried out reducing the equation to a partial difference equation defined on six points and then using the equations defined by the first integrals. The equations defined by the first integrals are reduced to discrete Riccati equations. The solution is then given in terms of four discrete integrations.
Summing up these results and the fact that the general solution of the t H ε 1 equation (1.2a) was presented in [26] we prove Theorem 1.2.

The
We have that the following propositions hold true: the general solution is given by where b k , c k , α l and β l are arbitrary functions of their arguments. If δ = 0 its general solution is given by v k,l = a k + β l , (2.3a) where a k , b k , α l and β l are arbitrary functions of their arguments. If δ 1 = 0 then its general solution is given by where b k , c k , α l and β l are arbitrary functions of their arguments.
Proof . From [26] we know that the 1 D 2 equation (1.3a) is Darboux integrable, and that the form of the first integral depends on the value of the parameter δ 1 . We will begin with the general case when δ 1 = 0 and δ = 0 and then consider the particular cases. Case δ = 0 and δ 1 = 0. In this case the W 1 first integral of the 1 D 2 equation (1.3a) is given by [26] As stated in the introduction, from the relation W 1 = ξ n this first integral defines a three-point, second-order ordinary difference equation in the n direction which depends parametrically on m.
From the parametric dependence we find two different three-point non-autonomous ordinary difference equations corresponding to m even and m odd. We treat them separately. Case m = 2l. If m = 2l we have the following non-autonomous ordinary difference equation where without loss of generality we have chosen α = 1 and β = 1. We can easily see, that once solved for u n+1,2l the equation is linear Tackling this equation directly is very difficult, but we can separate again the cases when n is even and odd and convert (2.6) into a system using the standard transformation (1.12a) where δ is given by (2.1). Now we have two first-order ordinary difference equations. Equation (2.7b) is uncoupled from equation (2.7a). Furthermore, since ξ 2k and ξ 2k+1 are independent functions we can write ξ 2k+1 = a k+1 − a k . So the second equation possesses the trivial solution 2 v k,l = α l + a k .
We define ξ 2k = b k /b k−1 and perform the change of dependent variable: w k,l = b k W k,l . Then W k,l solves the equation The solution of this difference equation is given by where c k is such that Equation (2.9) is not a total difference, but it can be used to define a k in terms of the arbitrary functions b k and c k This means that we have the following solution for the system (2.7)

11b)
Case m = 2l + 1. If m = 2l + 1 we have the following non-autonomous ordinary difference equation We can easily see that the equation is genuinely nonlinear. However we can separate the cases when n is even and odd and convert (2.12) into a system using the standard transformation (1.12b) where we used the definition (2.1). This is a system of two first-order difference equation, and equation (2.13a) is linear and uncoupled from (2.13b). As Hence the solution of (2.14) is given by Inserting (2.15) into (2.13b) and using the definition of ξ 2k+1 in terms of a k , i.e., ξ 2k+1 = a k+1 − a k we obtain We can then represent the solution of (2.16) as where d k satisfies the first-order linear difference equation Inserting the value of a k given by (2.10) inside (2.17) we obtain that this equation is a total difference. Then d k is given by This means that finally we have the following solutions for the fields z k,l and y k,l Equations (2.11), (2.18) provide the value of the four fields, but we have too many arbitrary functions in the l direction, namely α l , β l , γ l and ζ l . Introducing (2.11), (2.18) into (1.3a) and separating the terms even and odd in n and m we obtain two independent equations which allow us to reduce by two the number of independent functions in the l direction. Solving equations (2.19) with respect to γ l and ζ l we obtain Inserting (2.20) into (2.11), (2.18) we obtain that the general solution of 1 D 2 equation (1.3a) is given by (2.2), provided that δ 1 = 0 and δ = 0. Indeed the solution (2.2) is ill-defined if δ 1 = 0 or δ = 0 and we proceed to treat the relevant cases separately.
Case δ = 0. If δ = 0 we can solve (2.1) with respect to δ 1 The first integral (2.5) is not singular for δ 1 given by (2.21). The procedure of solution becomes different only when we arrive to the systems of ordinary difference equations (2.7) and (2.13). So we will present the solution of the systems in this case. Case m = 2l. If δ 1 is given by equation (2.21) the system (2.7) becomes The system (2.22) is uncoupled and imposing

23b)
Case m = 2l + 1. If δ 1 is given by equation (2.21) the system (2.13) becomes where we used the fact that ξ 2k = b k /b k−1 and ξ 2k+1 = a k+1 − a k . The solution to this system is immediate and it is given by As in the general case we obtained the expressions of the four fields, but we have too many arbitrary functions in the l direction, namely α l , β l , γ l and ζ l . Substituting the obtained expressions (2.23), (2.24) in the equation 1 D 2 (1.3a) with δ 1 given by equation (2.21) separating the even and odd terms we obtain two compatibility conditions We can solve this equation with respect to γ l and ζ l and we obtain Case δ 1 = 0. If δ 1 = 0 the first integral (2.5) is singular. Then following [26] the 1 D 2 equation (1.3a) with δ 1 = 0 possesses in the direction n the following three-point, second-order integral In order to solve the 1 D 2 equation (1.3a) in this case we use the first integral (2.26). We start separating the cases even and odd in m. Case m = 2l. If m = 2l we obtain from the first integral (2.26) where we have chosen without loss of generality α = β = 1. This equation is nonlinear. Applying the transformation (1.12a) we transform equation (2.27) into the system The system (2.28) is linear and equation (2.28b) is uncoupled from equation (2.28a). If we put ξ 2k+1 = a k+1 − a k then equation (2.28b) has the solution v k,l = a k + α l .
Substituting into (2.28a) we obtain Equation (2.29) becomes a total difference if we set and then the solution of the system (2.22) is given by Case m = 2l + 1. If m = 2l + 1 we obtain from the first integral (2.26) where we have chosen without loss of generality α = β = 1. This equation is nonlinear. Applying the transformation (1.12b) we transform equation (2.32) into the system where we used the values of ξ 2k and ξ 2k+1 . The system is now linear and equation (2.33a) is solved by Substituting into (2.33b) we obtain Then we have that y k,l is given by where d k solves the ordinary difference equation In (2.34) we inserted the value of a k according to (2.30). Equation (2.34) is a total difference and then d k is given by Then the solution of the system (2.33) is As in the general case we obtained the expressions of the four fields, but we have too many arbitrary functions in the l direction, namely α l , β l , γ l and ζ l . Substituting the obtained expressions (2.31), (2.35) in the equation 1 D 2 (1.3a) with δ 1 = 0 separating the even and odd terms we obtain two compatibility conditions We can solve equation (2.36) with respect to γ l and ζ l and to obtain Inserting (2.37) into (2.31), (2.35) we obtain that the general solution of 1 D 2 equation (1.3a) when δ = 0 is given by (2.4).
This discussion exhausts the possible cases. For any value of the parameters we have the general solution of the 1 D 2 equation (1.3a) and this concludes the proof.
where b k , c k , α l and β l are arbitrary functions of their arguments. If δ 1 = 0 then its general solution is given by where a k , b k , α l and β l are arbitrary functions of their arguments.
Proof . The proof of the two solution (2.38) and (2.39) proceeds as the one outlined in Proposition 2.1. The interested reader can find it in Appendix A.
where b k , c k , α l and β l are arbitrary functions of their arguments. If δ = 0 then its general solution is given by where b k , c k , α l and β l are arbitrary functions of their arguments. If δ 1 = 0 then its general solution is given by where b k , c k , α l and γ l are arbitrary functions of their arguments and ζ 0 is a constant.
Proof . The proof of the three solution (2.40), (2.41) and (2.42) proceeds as the one outlined in Proposition 2.1. The interested reader can find it in Appendix A.

2.2
The D 3 and the i D 4 equations, i = 1, 2 We have the following propositions: We have that the expression of the fields y k,l and v k,l is given by where the function functions c k , d k , α l and β l are arbitrary functions of their arguments, whereas the function e k and γ l are given through the discrete integrations The fields z k,l and w k,l are then given in terms of y k,l and v k,l as z k,l = (y k+1,l y k,l−1 − y k+1,l−1 y k,l ) y k,l − (y k+1,l + y k,l−1 − y k+1,l−1 − y k,l ) v k,l (y k+1,l + y k,l−1 − y k+1,l−1 − y k,l ) y k,l − y k+1,l y k,l−1 + y k+1,l−1 y k,l (2.45a) Remark 2.5. We remark that we can say that equation (2.44b) defines a discrete integration for the function γ l since it can be expressed in the form (1.14). Defining a the new function ζ l by we have that ζ l satisfies the following difference equation i.e., it is given by a discrete integration.
Proof . To find general of the D 3 equation (1.3d) we start from the equation itself. Applying the general transformation (1.12) to the D 3 equation (1.3d) we obtain the following system of four equations: From the system (2.47) we have four different way for calculating z k,l . This means that we have some compatibility conditions. Indeed from (2.47a) and (2.47c) we obtain the following while from (2.47b) and (2.47d) we obtain the following equation for v k+1,l+1 : Equations (2.48) and (2.49) give rise to a compatibility condition between v k+1,l and its shift in the l direction v k+1,l+1 which is given by Discarding the trivial solution v k,l = w 2 k,l we obtain the following value for the field w k,l v k+1,l = y k+1,l − y k+1,l−1 y k,l − y k,l−1 v k,l + (y k+1,l−1 y k,l − y k+1,l y k,l−1 ) 2 (y k+1,l−1 + y k,l − y k+1,l − y k,l−1 )(y k,l − y k,l−1 ) .
we can simplify (2.51) to the equation To go further we need to specify the form of the field y k,l . This can be obtained from the Darboux integrability of the D 3 equation (1.3d). From [26] we know that the D 3 equation (1.3d) possesses the following W 1 four-point, third-order integral Consider now the equation W 1 = ξ n with W 1 given as in (2.54). This relation defines a thirdorder, four-point ordinary difference equation in the n direction depending parametrically on m.
In particular if we choose the case when m = 2l + 1 we have the equation where we have taken without loss of generality α = β = 1. Using the transformation (1.12b) equation (2.55) is converted into the system This system is nonlinear, but if we solve (2.56b) with respect to z k+1,l and we substitute it along with its shift in the k direction into (2.56a) we obtain a linear second-order ordinary difference equation involving only the field y k,l We can lower the order of equation (2.57) by one using the potential transformation Then Y k,l solves the equation Imposing that we obtain that Y k,l can be expressed as From (2.58) we have that Setting we have that the solution of is given by This means that we can write down the following solution for V k,l i.e., up to the condition (2.44a). Substituting the value of V k,l from (2.60) and of y k,l from (2.59) into equation (2.52) we have that v k,l is given by equation (2.43b). Plugging the obtained value of v k,l we can compute w k,l from (2.50). Finally we can compute z k,l from the original system (2.47) and we obtain a single compatibility condition given by i.e., just by (2.44b). Given this conditions all the equations in (2.47) are compatible and z k,l is indifferently given by solving one of the equation. E.g., solving (2.47a) we can say that z k,l is given by equation (2.45a). This ends the procedure of solution of the D 3 equation (1.3d).
Proposition 2.6. The 1 D 4 equation (1.3e) is exactly solvable. We have that the expression of the fields y k,l and v k,l is given by where the function functions c k , d k , α l and β l are arbitrary functions of their arguments, whereas the function e k and γ l are given through the discrete integrations , The fields z k,l and w k,l are then given in terms of y k,l and v k,l as Remark 2.7. We remark that we can say that equation (2.63b) defines a discrete integration for the function γ l since it can be expressed in the form (1.14). Defining a new function ζ l by we have that ζ l is given by the following difference equation i.e., it is given by a discrete integration. Note that (2.65) is exactly the same as (2.46).
Proof . The proof of the solution of the 1 D 2 equation (1.3e) proceeds as the one outlined in Proposition 2.4. The interested reader can find the details in Appendix A.
Proposition 2.8. The 2 D 4 equation (1.3f) is exactly solvable. We have that the expression of the fields y k,l and v k,l is given by where the function functions c k , d k , α l and β l are arbitrary functions of their arguments, whereas the function e k and γ l are given through the discrete integrations The fields z k,l and w k,l are then given in terms of y k,l and v k,l as z k,l = − 1 δ 1 δ 1 δ 3 (y k,l−1 + y k+1,l − y k,l − y k+1,l−1 ) +(v k,l + δ 2 y k,l )(y k+1,l−1 y k,l − y k+1,l y k,l−1 ) y k+1,l−1 y k,l − y k+1,l y k,l−1 +y k,l (y k,l−1 + y k+1,l − y k,l − y k+1,l−1 ) , y k+1,l−1 y k,l − y k+1,l y k,l−1 y k,l−1 + y k+1,l − y k,l − y k+1,l−1 .
(2.68b) Remark 2.9. We remark that we can say that equation (2.63b) defines a discrete integration for the function γ l since it can be expressed in the form (1.14). Defining a new function ζ l by γ l = ζ l δ 2 β l−1 α l − β l α l−1 we have that ζ l is given by the following difference equation i.e., it is given by a discrete integration.
Proof . The proof of the solution of the 2 D 2 equation (1.3f) proceeds as the one outlined in Proposition 2.4. The interested reader can find the details in Appendix A.

2.3
The t H ε 2 and the t H ε 3 equations In this subsection we construct a general solution of the t H ε 2 and the t H ε 3 equations. As we recalled in the introduction, the solution of the t H ε 1 through the first integrals was already presented in [26], so we will not discuss it again. Moreover we also recall that the general solution of the t H ε 1 equation was first found in [21,22] without the knowledge of the first integrals. The first integrals of the t H ε 1 equation were first presented in [23]. The procedure we will follow will make use of the first integrals, in a similar way than in the cases presented in Section 2.2. The main difference is in the fact that the H 4 are nonautonomous only in the direction m, i.e., they depend only on the non-autonomous factors F (±) m as given by (1.4). Therefore instead of the general transformation (1.12) we can use the simplified transformation u n,2l = p n,l , u n,2l+1 = q n,l .
− (α 2 + α 3 )(q n,l − q n,l−1 ) − εα 2 3 (q n,l − q n,l−1 ) + ε α 3 2 + 2α 3 (q n,l−1 + α 2 ) − (q n,l − q n,l−1 )q n+1,l−1 q n+1,l + ε(α 2 + q n,l )(α 2 + q n,l−1 )q n+1,l −εq n+1,l−1 α 3 2 − 2(q n,l + α 2 )α 3 + (α 2 + q n,l )(α 2 + q n,l−1 ) where c n , ζ l and β l are arbitrary functions of their arguments and e n is a solution of the equation while f n and γ l are given by the discrete integrations f n − f n−1 = e n c n−1 − c n e n−1 , If ε = 0, but the field q n,l do not satisfy the discrete wave equation (2.83) then the solution of the t H ε 2 equation (1.2b) is given by where c n , β l and γ l are arbitrary functions of their arguments, ζ 0 is a constant and e n is a solution of (2.72) and f n is a solution of (2.73a). If the field q n,l satisfies the discrete wave equation (2.70) regardless of the value of the parameter ε the solution of the t H ε 2 equation (1.2b) is given by q n,l = a n + ζ 0 , (2.75a) where b n and β l are arbitrary functions of their arguments, ζ 0 is a constant and a n and c n are given by the discrete integration a n+1 − a n + α 2 a n+1 − a n − α Remark 2.11. We remark that the function e n can be obtained from (2.72) as the result of two discrete integrations. Indeed defining and substituting in (2.72) we obtain that E n must solve the equation Note that the right-hand side of (2.78) is not a total difference. So the function e n can be obtained by integrating (2.78) and subsequently integrating (2.77). This provides the value of e n . The obtained value can be plugged in (2.73a) to give f n after discrete integration. This reasoning shows that we can obtain the non-arbitrary functions e n and f n as result of a finite number of discrete integrations. Therefore we can conclude that the solution of the t H ε 2 equation (1.2b) in the general case is given in terms of four discrete integrations. If ε = 0 then the general solution is given in terms of three discrete integrations and finally in the singular case, when q n,l solves the discrete wave equation (2.70), we need only two discrete integrations.
We can introduce a new function b n through discrete integration a n+1 − a n + α 2 a n+1 − a n − α which is just formula (2.76a). Then we have that p n,l must solve the equation The solution of equation (2.84) is given by where c n is given by the discrete integration i.e., through formula (2.76b). This yields the solution of the t H ε 2 equation (1.2b) when q n,l satisfy the discrete wave equation (2.70).
(n − 1, l) (n − 1, l + 1) (n, l) (n, l + 1) (n + 1, l + 1) (n + 1, l) In the general case we have proved that the t H ε 2 equation (1.2b) is equivalent to the system (2.79) which in turn is equivalent to the solution of equations (2.81a) and (2.85). However (2.81a) merely defines p n,l+1 in terms of q n,l and its shifts. Therefore if we find the general solution of equation (2.85) the value of p n,l will follow. Applying T −1 l to equation (2.80a) we obtain then the value of p n,l as displayed in (2.71b). To find the solution for q n,l solution we turn to the first integrals. Like in the case of the H 6 equations (1.3) we will find an expression for q n,l using the first integrals, and then we will insert it into (2.85) to reduce the number of arbitrary functions to the right one. From [26] we know that t H ε 2 equation (1.2b) possesses the following four-point, third-order integral in the n direction . (2.86) We consider the equation W 1 = ξ n , where W 1 is given by (2.86), with m = 2l + 1 Using the substitutions (2.69) we have (q n−1,l − q n+1,l ) (q n+2,l − q n,l ) (q n,l − q n−1,l + α 2 ) (q n+1,l − q n+2,l + α 2 ) = ξ n .

(2.87)
This equation contains only q n,l and its shifts. From equation (2.87) it is very simple to obtain a discrete Riccati equation. Indeed the transformation Q n,l = q n,l − q n−1,l + α 2 q n+1,l − q n−1,l (2.88) brings (2.87) into: which is a discrete Riccati equation. Let us assume a n to be a particular solution of (2.89), then we express ξ n as ξ n = 1 a n (1 − a n+1 ) . (2.90) Using the standard linearization of the discrete Riccati equation Q n,l = a n + 1 Z n,l (2.91) from (2.90) we obtain the following equation for Z n,l Z n+1,l = a n Z n,l + 1 1 − a n+1 . Introducing If we assume that (2.93) can be written as a total difference, i.e., (2.94) So b n must be a total difference and therefore we can represent Z n,l as From (2.91) and (2.92) we obtain the form of Q n,l Q n,l = (c n − c n−1 )(c n+1 + ζ l ) (c n + ζ l )(c n+1 − c n−1 ) . (2.95) Introducing the value of Q n,l from (2.95) into (2.88) we obtain the following equation for q n,l q n+1,l − q n−1,l q n,l − q n−1,l + α 2 = (c n + ζ l )(c n+1 − c n−1 ) (c n − c n−1 )(c n+1 + ζ l ) .
Performing the transformation R n,l = (c n + ζ l )q n,l (2.96) we obtain the following second-order ordinary difference equation for the field R n,l Then we can represent the solutions of the equation (2.97) as R n,l = P n,l + ζ l e n + f n , where e n and f n are particular solutions of (2.99b) P n,l will be then solve the following equation The equations (2.99a) and (2.99b) are not independent. Indeed defining f n − f n−1 = e n c n−1 − c n e n−1 . (2.104b) The system (2.104) is just gives the constraints expressed in formulas (2.72) and (2.73a). Now we turn to the solution of the homogeneous equation (2.100). We can reduce (2.100) to a total difference using the potential substitution T n,l = P n,l − P n−1,l T n+1,l c n+1 − c n − T n,l c n − c n−1 = 0.
This clearly implies where β l is an arbitrary function. The solution to this equation is given by 3 P n,l = (c n + ζ l )β l + γ l , (2.105) where γ l is an arbitrary function. Using (2.96), (2.98), (2.105) we obtain then the following expression for q n,l q n,l = β l + γ l + ζ l e n + f n c n + ζ l , (2.106) where e n and f n are solutions of (2.104). Equation (2.106) is formally the solution presented in (2.71a), but depends on three arbitrary functions in the l direction, namely ζ l , β l and γ l . This means that there is a constraint between these functions, which can be recovered by plugging (2.106) into (2.85). At this point we have a second bifurcation, depending on the value of ε.
Case ε = 0. If ε = 0 inserting (2.106) into (2.85) and factorizing out the n dependent part away we are left with This equation tells us that the function γ l can be expressed after a discrete integration in terms of the two arbitrary functions ζ l and β l . This is just the final constraint expressed in formula (2.73b). This yields the general solution of the t H ε 2 equation (1.2b). Case ε = 0. If ε = 0 inserting (2.106) into (2.85) we obtain the compatibility condition ζ l+1 − ζ l = 0, i.e., ζ l = ζ 0 = const. It is easy to check that the obtained value of q n,l through formula (2.106) is consistent with the substitution of ε = 0 in (2.79). This means that in the case ε = 0 the value of q n,l is given by q n,l = β l + γ l + ζ 0 e n + f n c n + ζ 0 , where the functions e n and f n are defined implicitly and can be found by discrete integration from (2.104), i.e., from formula (2.74a). Since formula (2.71b) is not singular with respect to ε the value of p n,l can be recovered just by substituting ε = 0 and the form of q n,l found in (2.107). This yields equation (2.74b). This concludes the procedure of solution of the t H ε 2 equation (1.2b) in the case when ε = 0. p n,l = α 2 (q n+1,l − q n+1,l−1 ) ε 2 q n,l−1 q n,l + δ 2 α 3 2 + δ 2 α 2 2 α 2 3 (q n,l−1 − q n,l ) + ε 2 q n+1,l q n+1,l−1 (q n,l−1 − q n,l ) (q n+1,l q n,l−1 − q n+1,l−1 q n,l )α 3 α 2 , where e n , β l and γ l are arbitrary functions of their arguments and c n is a solution of the equation while f n and ζ l are given by the discrete integrations p n,l = ε 2 α 2 α 3 α 2 (q n+1,l − q n+1,l−1 )q n,l−1 q n,l + q n+1,l q n+1,l−1 (q n,l−1 − q n,l ) q n+1,l q n,l−1 − q n+1,l−1 q n,l , where e n , ζ l and γ l are arbitrary functions of their arguments β 0 is a constant and c n is a solution of (2.110) and f n is a solution of (2.111a). If the field q n,l satisfies the equation (2.108) regardless of the value of the parameter ε the solution of the t H ε 3 equation (1.2c) is given by q n,l = ζ 0 a n , (2.113a) where b n and β l are arbitrary functions of their arguments, ζ 0 is a constant and a n and c n are given by the discrete integration .
(2.114b) Remark 2.13. We remark that the function c n can be obtained from (2.110) as the result of two discrete integrations. Indeed defining z n = c n+1 − c n , (2.115) and substituting in (2.110) we obtain that z n must solve the equation Note that the right-hand side of (2.116) is not a total difference. So the function c n can be obtained by integrating (2.116) and subsequently integrating (2.115). This provides the value of c n . The obtained value can be plugged in (2.111a) to give f n after discrete integration. This reasoning shows that we can obtain the non-arbitrary functions c n and f n as result of a finite number of discrete integrations. Therefore we can conclude that the solution of the t H ε 3 equation (1.2c) in the general case is given in terms of four discrete integrations. If ε = 0 the general solution is given in terms of three discrete integration and finally in the singular case, when q n,l solves the equation (2.108), we need only two discrete integrations.
Proof . The proof of the solution of the t H ε 3 equation (1.2c) proceeds as the one outlined in Proposition 2.10. The interested reader can find the details in Appendix A.

Conclusions
In this paper we presented a detailed procedure to construct the general solutions for all the H 4 and H 6 equations. As stated in the introduction, these general solutions were obtained in three different ways, but the common feature is that they can be found through some linear or linearizable (discrete Riccati) equations. This is the great advantage of the first integral approach over the direct one which was pursued in [21]. The Darboux integrability therefore yields extra information that it is useful to get the final result, i.e., the general solutions. Moreover linearization arises very naturally from first integrals also in the most complicated cases, whereas in the direct approach can be quite tricky, see, e.g., the examples in [21]. The linearization of the first integrals is another proof of the deep linear nature of the H 4 and H 6 equations. This result is even stronger than Darboux integrability alone, since a priori the first integrals do not need to define linearizable equations.
We also note that our procedure of construction of the general solution, based on the ideas from [18], is likely to be the discrete version of the procedure of linearization and solutions for continuous Darboux integrable equations presented in [46]. The preeminent rôle of the discrete Riccati equation in the solutions is reminiscent of the importance of the usual Riccati equation in the continuous case. Recall, e.g., that the first integrals of the Liouville equation [36] u xt = e u , which is the most famous Darboux integrable system, are Riccati equations. Many other examples of solutions presented in [46] use the reduction of higher-order differential equations to Riccati-like equation in order to obtain the solution, as we have done in the discrete case.
In this paper we constructed the general solutions of the trapezoidal H 4 and of the H 6 equations. Therefore we possess an almost complete theory about these equations ranging from the geometrical background to their analytic properties. For a discussion of the open problems in this field we refer to our previous paper [26].
A Procedure to f ind the general solution in the remaining cases From [26] we know that the first integrals of the 2 D 2 equation (1.3b) can be different depending on the value of the parameter δ 1 . For this reason we treat separately the various cases.
Case δ 1 = 0. If δ 1 = 0 the W 1 first integrals of the 2 D 2 equation (1.3b) is given by [26] As stated in the introduction, from the relation W 1 = ξ n this first integral defines a three-point, second-order ordinary difference equation in the n direction which depends parametrically on m.
From this parametric dependence we find two different three-point non-autonomous ordinary difference equations corresponding to m even and m odd. We treat them separately. Case m = 2l. If m = 2l we have the following non-autonomous nonlinear ordinary difference equation Without loss of generality we set α = 1 and β = δ 1 . Then making the transformation u n,2l = U n,2l − δ 2 (A.2) and putting From the definition (1.12a) applied to U n,2l instead of u n,2l 4 we can separate again the even and the odd part in (A.4). We obtain the following system of two coupled first-order ordinary difference equations Putting ξ 2k = a k /a k−1 the solution to (A.5a) is given by Inserting the value of W k,l from (A.6) into (A.5b) we obtain If we define then (A.7) becomes a total difference. So we obtain the following solutions for the W k,l and the V k,l fields Inverting the transformation (A.2) we obtain for the fields w k,l and v k,l Case m = 2l + 1. If m = 2l + 1 we have the following non-autonomous ordinary difference equation n δδ 1 u n,2l+1 + u n+1,2l+1 δδ 1 u n,2l+1 + u n−1,2l+1 where we already substituted δ as defined in (A.3). Using the standard transformation (1.12b) to get rid of the two-periodic factors we obtain δδ 1 y k,l + z k,l δδ 1 y k,l + z k−1,l = ξ 2k , (A.10a) Both equations in (A.10) are linear in z k,l , y k,l and their shifts. As ξ 2k+1 = b k+1 − b k we have that the solution of equation (A.10b) is given by As ξ 2k = a k /a k−1 and y k,l is given by (A.11) we obtain Recalling the definition of a k in (A.8) we represent z k,l as Equations (A.9), (A.11), (A.12) provide the value of the four fields, but we have too many arbitrary functions in the l direction, namely α l , β l , γ l and ζ l . Inserting (A.9), (A.11), (A.12) into (1.3b) and separating the terms even and odd in n and m we obtain two independent equations δ 1 ζ l + α l δ 2 1 γ l − δ 1 2 λα l + β l α l δ 1 + δα l δ 1 − α l + α l δ 1 = 0, (A.13a) which allow us to reduce by two the number of independent functions in the l direction. Solving (A.13) with respect to γ l and ζ l we find Inserting (A.14) into equations (A.9), (A.11), (A.12) we have the general solution (2.38) of the 2 D 2 equation (1.3b) provided that δ 1 = 0. Indeed the solution of the 2 D 2 equation (1.3b) given by (2.38) is ill-defined if δ 1 = 0. Therefore we now discuss this case separately.
Case δ 1 = 0. Following [26] we have the 2 D 2 equation (1.3b) with δ 1 = 0 possesses the following two-point, first-order first integral in the direction n where we have chosen without loss of generality α = β = 1. Applying the transformation (1.12a) equation (A.16) becomes the system v k,l (δ 2 + w k,l ) = ξ 2k , In this case the system (A.17) do not consist of purely difference equations. Indeed from (A.17a) we can derive immediately the value of the field w k,l we have that (A.19) becomes a total difference. So we have that the system (A.17) is solved by Case m = 2l + 1. If m = 2l + 1 we obtain from the first integral (A.15) ) becomes a total difference. Therefore we can write the solution of the system (A.23) as From [26] we know that the 3 D 2 equation (1.3c) is Darboux integrable, and that the form of the first integral depends on the value of the parameter δ. We will begin with the general case when δ 1 = 0 and δ = 0 and then consider the particular cases.
Case δ 1 = 0 and δ = 0. In this case we know that the W 1 first integrals of the 3 D 2 equation (1.3c) is given by [26] As stated in the introduction, from the relation W 1 = ξ n this first integral defines a three-point, second-order ordinary difference equation in the n direction which depends parametrically on m.
From this parametric dependence we find two different three-point non-autonomous ordinary difference equations corresponding to m even and m odd. We treat them separately. Case m = 2l. If m = 2l we have the following non-autonomous nonlinear ordinary difference equation where we have chosen without loss of generality α = β = 1. We can apply the usual transformation (1.12a) in order to separate the even and odd part in (A.27) .
where δ is given by equation (2.1). Putting ξ 2k = a k /a k−1 we have the following solution for (A.30a) Plugging (A.31) into equation (A.30b) and defining we have that equation (A.30b) becomes a total difference. Then the solution of (A.30b) can be written as So using (A.29) we obtain the following solution for the original system (A.28): Case m = 2l + 1. If m = 2l + 1 we have the following non-autonomous ordinary difference equation where without loss of generality α = β = 1 and δ is given by (2.1). Solving with respect to u n+1,2l+1 it is immediate to see that the resulting equation is linear. Then separating the even and the odd part using the transformation (1.12b) we obtain the following system of linear, first-order ordinary difference equations Then, in the usual way, we can represent the solution as where we have used the explicit definition of a k given in (A.32). So we have the explicit expression for both fields y k,l and z k,l . Equations (A.33), (A.35), (A.36) provide the value of the four fields, but we have too many arbitrary functions in the l direction, namely α l , β l , γ l and ζ l . Inserting (A.33), (A.35), (A.36) into (1.3c) and separating the terms even and odd in n and m we obtain we obtain two equations which allow us to reduce by two the number of independent functions in the l direction. Indeed solving (A.37) with respect to γ l and ζ l we find Inserting (A.38) into (A.33), (A.35), (A.36) we obtain the general solution (2.40) of the 3 D 2 equation (1.3c) provided that δ 1 = 0 and δ = 0. It is easy to see that the solution (2.40) is ill-defined if δ 1 = 0 and if δ = 0. We will treat these two particular cases separately.
Case δ = 0. If δ = 0 we have that δ 1 is given by equation (2.21). In this case the first integral (A.26) is singular since the coefficient of α goes to a constant. Following [26] we have that the 3 D 2 equation with δ 1 given by (2.21) possesses the following first integral in the direction n This first integral is a three-point, second-order first integral. As in the general case we consider separately the m even and odd cases.
Case m = 2l. If m = 2l then the first integral (A.39) becomes the following nonlinear three-point, second-order difference equation where without loss of generality α = β = 1. If we separate the even and the odd part using the general transformation given by (1.12a) we obtain the system This is a system of first-order nonlinear difference equations. However (A.40a) is uncoupled from (A.28b), and it is a discrete Riccati equation which can be linearized through the Möbius transformation (A.29). This linearize the system (A.40) to This yields the following solution of the system (A.40) Case m = 2l + 1. If m = 2l + 1 the first integral (A.39) becomes the following nonlinear, three-point, second-order difference equation where without loss of generality α = β = 1. As usual we can separate the even and odd part in n using the transformation (1.12b). This transformation brings equation (A.46) into the following linear system where we used (A.44) and the definition ξ 2k+1 = a k+1 − a k . Equation (A.47b) is readily solved and gives We can then write for z k,l the following expression where d k solves the equation Equation (A.49) is a total difference with d k given by Therefore we have the following solution to the system (A.47) Equations (A.45), (A.50) provide the value of the four fields, but we have too many arbitrary functions in the l direction, namely α l , β l , γ l and ζ l . Inserting (A.45), (A.50) into (1.3c) with δ 1 given by (2.21) and separating the terms even and odd in n and m we obtain we obtain two equations γ l (δ 2 + 1)α l − λ + ζ l + β l (δ 2 + 1) = 0, Solving this compatibility condition with respect to γ l and ζ l we obtain Case m = 2l. If δ 1 = 0 the system (A.28) becomes . (A.52b) The Case m = 2l + 1. If δ 1 = 0 the system (2.13) becomes where we used (A.55) and ξ 2k = a k−1 /a k . The system is linear and equation (A.57b) is uncoupled from (A.57a). If we put then equation (A.57a) becomes a total difference whose solution is Substituting y k,l given by (A.58) into equation (A.57a) we obtain We can therefore represent the solution as where d k solves the equation Equation (A.59) is a total difference and d k is given by Therefore we have that the solution of the system (A.57) is given by From equations (A.56), (A.60) we have the value of the four fields, but we have too many arbitrary functions in the l direction, namely α l , β l , γ l and ζ l . Inserting (A.56), (A.60) into (1.3c) with δ 1 = 0 and separating the terms even and odd in n and m we obtain we obtain two equations We can solve this compatibility conditions with respect to β l and ζ l we obtain where ζ 0 is a constant. To find the general solution we start from the 1 D 4 equation (1.3e) itself. Applying the general transformation (1.12) we transform the 1 D 4 into the following system v k,l z k,l + w k,l y k,l + δ 1 w k,l z k,l + δ 2 y k,l z k,l + δ 3 = 0, (A.62a) y k,l w k,l+1 + z k,l v k,l+1 + δ 1 z k,l w k,l+1 + δ 2 y k,l z k,l + δ 3 = 0, (A.62b) w k,l y k+1,l + v k+1,l z k,l + δ 1 w k,l z k,l + δ 2 z k,l y k+1,l + δ 3 = 0, (A.62c) From the system (A.62) we have four different way for calculating z k,l . This means that we have some compatibility conditions. Indeed from (A.62a) and (A.62c) we obtain the following equation for v k+1,l v k+1,l = δ 3 + w k,l y k+1,l δ 3 + w k,l y k,l v k,l + (y k+1,l − y k,l ) δ 1 w 2 k,l − δ 2 δ 3 δ 3 + w k,l y k,l , (A.63) while from (A.62b) and (A.62d) we obtain the following equation for v k+1,l+1 v k+1,l+1 = δ 3 + y k+1,l w k,l+1 Equations (A.63) and (A.64) give rise to a compatibility condition between v k+1,l and its shift in the l direction v k+1,l+1 which is given by (y k+1,l+1 y k,l − y k+1,l y k,l+1 )w k,l+1 Discarding the trivial solution v k,l = −δ 1 w k,l + δ 2 δ 3 w k,l we obtain the following value for the field w k,l w k,l = δ 3 y k+1,l−1 − y k+1,l − y k,l−1 + y k,l y k+1,l y k,l−1 − y k+1,l−1 y k,l , (A.65) which makes (A.63) and (A.64) compatible. This gives us the first piece of the solution in (2.64b).
To go further we need to specify the form of the field y k,l . This can be obtained from the Darboux integrability of the 1 D 4 equation (1.3e). From [26] we know that the 1 D 4 equation (1.3e) we have the following four-point, third-order W 1 integral n F (+) m α u 2 n+1,m δ 1 + u n+1,m u n+2,m + u n−1,m (u n,m − u n+2,m ) − δ 2 δ 3 u n+1,m (δ 1 + u n,m ) − δ 2 δ 3 This first integral defines as always the relation W 1 = ξ n which is a third-order, four-point ordinary difference equation in the n direction depending parametrically on m. In particular when m = 2l + 1 we have the equation n (u n,2l+1 − u n+2,2l+1 + δ 1 u n+1,2l+1 )u n−1,2l+1 + u n+1,2l+1 u n+2,2l+1 (u n,2l+1 + δ 1 u n−1,2l+1 )u n+1,2l+1 where we have chosen without loss of generality α = β = 1. Using the transformation (1.12b) then (A.68) is converted into the system (y k,l − y k+1,l + δ 1 z k,l )z k−1,l + y k+1,l z k,l = ξ 2k (y k,l + δ 1 z k−1,l )z k,l , (A.69a) (y k,l − y k+1,l )(z k,l − z k+1,l ) = ξ 2k+1 z k,l (z k+1,l δ 1 + y k+1,l ). (A.69b) If we solve (A.69b) with respect to z k+1,l and then substitute into (A.69a) we obtain a linear, second-order ordinary difference equation for y k,l We can find the solution to this equation in a similar manner than in the case of the D 3 equation. First of all let us consider Y k,l = a k y k,l + b k y k−1,l such that Y k+1,l − Y k,l is equal to the left-hand side of (A.70). To this end we define Therefore y k,l must solve the following equation Then the solution of (A.71) is then given by . This means that the solution of V k,l can be represented as where e k is defined by the discrete integration , which is just equation (2.63a). Inserting the value of V k,l from (A.73) and the value of y k,l (A.72) into equation (A.66) we obtain equation (2.62b). From the obtained value of v k,l we can compute w k,l using (A.65). So finally we can compute z k,l from the original system (A.62), and we obtain a single compatibility condition given by which is just equation (2.63b). Since now the system (A.62) is compatible we can use any of its equations to compute z k,l . E.g., using equation (A.62a) and the value of w k,l from equation (A.65) we obtain equation (2.64b). This concludes the procedure of solution of the 1 D 4 equation (1.3e).
General case: q n,l do not solve (2.108). If the field q n,l do not solve (2.108) we have that we can define p n,l+1 and p n+1,l+1 as in (A.89a) and (A.89b) respectively. Furthermore these two equations must be compatible. The compatibility condition is obtained applying T −1 l to (A.89b) and imposing to the obtained expression to be equal to (A.89a). We then find that q n,l must solve the following equation δ 2 α 2 2 α 2 3 q n+1,l+1 q n,l − q n,l q n−1,l+1 + q n,l+1 (q n−1,l − q n+1,l ) − α 2 ε 2 q n,l q n,l+1 + δ 2 α 2 3 (q n+1,l+1 q n−1,l − q n−1,l+1 q n+1,l ) + ε 2 q n,l q n+1,l+1 (q n−1,l − q n+1,l ) − q n,l+1 q n−1,l q n+1,l q n−1,l+1 + ε 2 q n+1,l+1 q n+1,l q n,l+1 q n−1,l = 0. (A.96) As in the case of the t H ε 2 equation (1.2b) the partial difference equation for q n,l is not defined on the square quad graph of Fig. 1, but it is defined on the six-point lattice shown in Fig. 3. Moreover once equation (A.96) is solved we can use indifferently (A.89a) or (A.89b) to obtain the value of the field p n,l since these two merely defines p n,l+1 in terms of q n,l and its shifts. Therefore if we find the general solution of (A.96) the value of p n,l will follow. E.g., if we solve (A.96) applying T −1 l to (A.89a) we will obtain (2.109b) which is then the first part of the general solution. To find the solution of equation (A.96) we turn to the first integrals. Like in the case of the t H ε 2 equation (1.2b) we will find an expression for q n,l using the first integrals, and then we will insert it into (A.96) to reduce the number of arbitrary functions to the right one. From [26] we know that the t H ε 3 equation (1.2c) possesses a four-point, third-order integral in the n direction We consider the equation W 1 = ξ n /α 2 6 , where W 1 is given by (A.97), with m = 2l + 1 (u n+1,2l+1 − u n−1,2l+1 )(u n+2,2l+1 − u n,2l+1 ) (α 2 u n,2l+1 − u n−1,2l+1 )(u n+2,2l+1 − α 2 u n+1,2l+1 ) = ξ n .
Inserting the definition of R n,l (A.103) into (A.112) we obtain P n,l e n (f n + β l ) = P n−1,l e n−1 (f n−1 + β l ) , i.e., P n,l = γ l e n (f n + β l ). and e n is an arbitrary function, i.e., we have that c n must solve equation (2.110). Formally equation (A.114) has the form of the solution presented in formula (2.109a), but it depends on three arbitrary functions in the l direction, namely ζ l , β l and γ l . Therefore there must be a constraint between these functions. This constraint can be obtained by plugging (A.114) into (A.96), but here we have another bifurcation depending on the value of parameter δ. Indeed it is easy to see that we must distinguish the cases when δ = 0 and when δ = 0.
This equation tells us that the function ζ l can be expressed after a discrete integration in terms of the two arbitrary functions β l and γ l . This condition is just (2.111b). This yields the general solution of the t H ε 3 equation (1.2c) when δ = 0 and the field q n,l do not satisfy equation (2.108). Case δ = 0. Inserting (A.114) into (A.96) if δ = 0 factorizing the n dependent part the compatibility condition is just β l+1 − β l = 0, i.e., β l = β 0 = const. It is easy to check that the obtained value of q n,l through formula (A.114) is consistent with the substitution of δ = 0 in (A.87) and therefore that in the case δ = 0 the value of q n,l is given by q n,l = γ l e n (f n + β 0 ) c n + ζ l , where the functions c n and f n are defined implicitly and can be found by discrete integration from (A.110b) and (A.115) respectively. This is just equation (2.112a). Since equation (2.109b) is not singular if δ = 0 we obtain that in this case the general solution for the field p n,l is given by substituting δ = 0 into (2.109b), i.e., just by equation (2.112b) where q n,l is simply given by (2.112a). This yields the general solution of the t H ε 3 equation (1.2c) in the case when δ = 0.