arxiv:nlin/ v1 [nlin.si] 16 Nov 2003

Size: px
Start display at page:

Download "arxiv:nlin/ v1 [nlin.si] 16 Nov 2003"

Transcription

1 Centre de Recherches Mathématiques CRM Proceedings and Lecture Notes arxiv:nlin/ v1 [nlinsi] 16 Nov 2003 Symbolic Computation of Conserved Densities, Generalized Symmetries, and Recursion Operators for Nonlinear Differential-Difference Equations Willy Hereman, Jan A Sanders, Jack Sayers, and Jing Ping Wang This paper is dedicated to Ryan Sayers ( ) Abstract Algorithms for the symbolic computation of conserved densities, fluxes, generalized symmetries, and recursion operators for systems of nonlinear differential-difference equations are presented In the algorithms we use discrete versions of the Fréchet and variational derivatives, as well as discrete Euler and homotopy operators The algorithms are illustrated for prototypical nonlinear lattices, including the Kac-van Moerbeke (Volterra) and Toda lattices Results are shown for the modified Volterra and Ablowitz-Ladik lattices 1 Introduction The study of complete integrability of nonlinear differential-difference equations (DDEs) largely parallels that of PDEs [AC91, ASY00, W98] Indeed, as in the continuous case, the existence of sufficiently many conserved densities and generalized symmetries is a predictor for complete integrability Based on the first few densities and symmetries a recursion operator (which maps symmetries to symmetries) can be constructed The existence of a recursion operator, which allows one to generate an infinite sequence of symmetries, confirms complete integrability There is a vast body of work on the complete integrability of DDEs Consult eg [ASY00, HH03] for various approaches and references In this article we describe algorithms to symbolically compute polynomial conservation laws, fluxes, generalized symmetries, and recursion operators for DDEs The design of these algorithms heavily relies on related work for PDEs [GH97a, HG99, HGCM98] and work by Oevel et al [OZF93] There exists a close analogy between the continuous and discrete (in space) cases As shown in [HM02, MH02], this analogy can be completely formalized 1991 Mathematics Subject Classification Primary 37J35, 37K10; Secondary 35Q58, 37K05 Key words and phrases Homotopy operators, conserved densities, fluxes, recursion operators, integrability, differential-difference, DDEs, semi-discrete lattices This material is based upon work supported by the National Science Foundation under Grant No CCR c 2003 Willy Hereman

2 2 WILLY HEREMAN, JAN A SANDERS, JACK SAYERS, AND JING PING WANG andboththeoriescanbeformulatedintermsofcomplexes Thesameappliesforthe formulation in terms of Lie algebra complexes in[w98] This allows one to translate by analogy the existing algorithms immediately (although complications arise when there is explicit dependence on the space variable in the discrete case) One of the more useful tools following from the abstract theory is the homotopy operator, which is based on scaling vectorfields, and goes back to Poincaré in the continuous case This operator allows one to directly integrate differential forms and can be straightforwardly implemented in computer algebra packages, since it reduces to integration over one scaling parameter The discrete analogue as described in [HM02, MH02] does the corresponding job in the discrete case and we use it to compute fluxes In this paper we do not explicitly use the abstract framework, yet, it has been a motivating force for the development of our algorithms The algorithms in this paper can be implemented in any computer algebra system Our Mathematica package InvariantsSymmetriesm [GH97b] computes densities and symmetries, and therefore aids in automated testing of complete integrability of semi-discrete lattices Mathematica code to automatically compute recursion operators is still under development The paper is organized as follows In Section 2, we give key definitions and introduce the Kac-van Moerbeke(KvM)[KvM75], Toda[T81] and Ablowitz-Ladik (AL) lattices [AL75], which will be used as prototypical examples throughout the paper The discrete higher Euler operators(variational derivatives) and the discrete homotopy operator are introduced in Section 3 These operators are applied in the construction of densities and fluxes in Section 5 The algorithm for higher-order symmetries is outlined in Section 6 The algorithms for scalar and matrix recursion operators are given in Section 7 and Section 8 The paper concludes with two more examples in Section 9 2 Key Definitions and Prototypical Examples Definition 21 A nonlinear (autonomous) DDE is an equation of the form (21) = F(, 1,,+1,), where and F are vector-valued functions with m components The integer n corresponds to discretization in space; the dot denotes one derivative with respect to the continuous time variable (t) For simplicity, we will denote the components of by (,,w n, ) For brevity, we write F( ), although F typically depends on and a finite number of its forward and backward shifts We assume that F is polynomial with constant coefficients If present, parameters are denoted by lower-case Greek letters No restrictions are imposed on the level of the forward or backward shifts or the degree of nonlinearity in F (22) Example 22 The AL lattice [AL75], = ( )+ ( ), = ( ) ( ), is a completely integrable discretization of the nonlinear Schrödinger equation Definition 23 A DDE is said to be dilation invariant if it is invariant under a scaling (dilation) symmetry

3 DENSITIES, SYMMETRIES, AND RECURSION OPERATORS FOR NONLINEAR DDES 3 Example 24 The KvM lattice [KvM75], (23) = (+1 1 ), is invariant under (t, ) (λ 1 t,λ ), where λ is an arbitrary scaling parameter Example 25 The Toda lattice [T81] in polynomial form [GH98], (24) = 1, = ( +1 ), is invariant under the scaling symmetry (25) (t,, ) (λ 1 t,λ,λ 2 ) Thus, and correspond to one, respectively two derivatives with respect to t, (26) d dt, d2 dt 2 Definition 26 We define the weight, w, of a variable as the number of t- derivatives the variable corresponds to Since t is replaced by t λ, we set w( d dt ) = 1 In view of (26), we have w() = 1, andw( ) = 2fortheTodalattice Weightsofdependentvariablesarenonnegative, rational, and independent of n Definition 27 The rank of a monomial is defined as the total weight of the monomial An expression is uniform in rank if all its terms have the same rank Example 28 In the first equation of (24), all the monomials have rank 2; in the second equation all the monomials have rank 3 Conversely, requiring uniformity inrankforeachequationin(24)allowsonetocomputetheweightsofthedependent variables with simple linear algebra Indeed, (27) w( )+1=w( ), w( )+1=w( )+w( ), yields (28) w( ) = 1, w( ) = 2, which is consistent with (25) Dilation symmetries, which are special Lie-point symmetries, are common to many lattice equations Lattices that do not admit a dilation symmetry can be made scaling invariant by extending the set of dependent variables using auxiliary parameters with scaling Example 29 The AL lattice is not dilation invariant Introducing an auxiliary parameter α, hence replacing (22) by (29) = α( )+ ( ), = α( ) ( ), and requiring uniformity in rank, gives (210) w( )+1 = w(α)+w( ) = 2w( )+w( ), w( )+1 = w(α)+w( ) = 2w( )+w( ) Obviously, (211) w( )+w( ) = w(α) = 1

4 4 WILLY HEREMAN, JAN A SANDERS, JACK SAYERS, AND JING PING WANG Several scales are possible The choice w( ) = w( ) = 1 2,w(α) = 1, corresponds to the scaling symmetry (212) (t,,,α) (λ 1 t,λ 1 2 un,λ 1 2 vn,λα) Definition 210 A scalar function ρ n ( ) is a conserved density of (21) if there exists a scalar function J n ( ) called the associated flux, such that (213) D t ρ n + J n = 0 is satisfied on the solutions of (21) [O93] In (213), we used the (forward) difference operator, = D I, defined by (214) J n = (D I)J n = J n+1 J n, where D denotes the up-shift (forward or right-shift) operator, DJ n = J n+1, and I is the identity operator Operator takes the role of a spatial derivative on the shifted variables as many DDEs arise from discretization of a PDE in (1 +1) variables Most, but not all, densities are polynomial in (215) (216) Example 211 The first three conservation laws for (23) are D t (ln( ))+(+1 + ) ( + 1 ) = 0, D t ( )+(+1 ) ( 1 ) = 0, (217) D t ( 1 2 u2 n ++1 )+ +1 ( ) 1 ( ++1 ) = 0 (218) (219) (220) (221) Example 212 For the Toda lattice (24) the first four density-flux pairs are ρ (0) n = ln( ), J n (0) n, ρ (1) n =, J n (1) n 1, ρ (2) n = 1 2 u2 n +, J (2) n = 1, ρ (3) n = ( 1 + ), J (3) n = The above densities are uniform of ranks 0 through 3 The corresponding fluxes are also uniform in rank with ranks 1 through 4 In general, if in (213) rankρ n = R then rankj n = R+1, since w(d t ) = 1 This comes as no surprise since the conservation law (213) holds on solutions of (21), hence it inherits the dilation symmetry of (21) In Section 5 we will give an algorithm to compute polynomial conserved densities and fluxes and use(24) to illustrate the steps Non-polynomial densities(which are easy to find by hand) can be computed with the method given in [HH03] Definition 213 Compositions of D and D 1 define an equivalence relation ( ) on monomials in the components of All shifted monomials are equivalent Example 214 For example, Definition 215 The main representative of an equivalence class is the monomial of the class with n as lowest label on any component of u (or v if u is missing) Example 216 For example, +2 is the main representative of the class {, 2, 1 +1, +2,+1 +3, } Use lexicographical ordering to resolve conflicts For example, +2 (not 2 ) is the main representative of the class {, 3 1,,+2 +4, }

5 DENSITIES, SYMMETRIES, AND RECURSION OPERATORS FOR NONLINEAR DDES 5 Definition 217 A vector function G( ) is called a generalized symmetry of (21) if the infinitesimal transformation +ǫg leaves (21) invariant up to order ǫ Consequently, G must satisfy [O93] (222) D t G = F ( )[G], on solutions of (21) F ( )[G] is the Fréchet derivative of F in the direction ofg For the scalar case (N = 1), the Fréchet derivative in the direction of G is (223) F ( )[G] = ǫ F( +ǫg) ǫ=0 = k which defines the Fréchet derivative operator (224) F ( ) = k F +k D k F +k D k G, In the vector case with two components, and, the Fréchet derivative operator is (225) F ( ) = Applied to G = (G 1 k F 1 +k D k k k F2 +k D k F 1 +k D k k F2 +k D k G 2 ) T, where T is transpose, one obtains (226) F i ( )[G] = k F i D k G 1 + +k k F i +k D k G 2, i = 1,2 In (223) and (226) summation is over all positive and negative shifts (including k = 0) For k > 0, D k = D D D (k times) Similarly, for k < 0 the down-shift operator D 1 is applied repeatedly The generalization of (225) to a system with N components should be obvious (227) (228) Example 218 The first two symmetries of (23) are G (1) = (+1 1 ), G (2) = +1 ( ) 1 ( ) These symmetries are uniform in rank (rank 2 and 3, respectively) The symmetries of ranks 0 and 1 are both zero Example 219 The first two non-trivial symmetries of (24) are ( ) (229) G (1) v = n 1, ( +1 ) and (230) G (2) = ( vn ( ++1 ) 1 ( 1 + ) (u 2 n+1 u 2 n ) The above symmetries are uniform in rank For example, rankg (2) 1 = 3 and rankg (2) 2 = 4 Symmetries of lower ranks are trivial An algorithm to compute generalized symmetries will be outlined in Section 6 and applied to (24) )

6 6 WILLY HEREMAN, JAN A SANDERS, JACK SAYERS, AND JING PING WANG Definition 220 A recursion operator R connects symmetries (231) G (j+s) = RG (j), where j = 1,2,, and s is the seed In most cases the symmetries are consecutively linked (s = 1) For N-component systems, R is an N N matrix operator The defining equation for R [O93, W98] is (232) D t R+[R,F ( )] = R t +R [F]+R F ( ) F ( ) R = 0, where [, ] denotes the commutator and the composition of operators The operator F ( ) was defined in (225) R [F] is the Fréchet derivative of R in the direction of F For the scalar case, the operator R is often of the form (233) R = U( )O((D I) 1,D 1,I,D)V( ), and then (234) R [F] = k (D k U F) OV + +k k UO(D k F) V +k For the vector case, the elements of the N N operator matrix R are often of the form R ij = U ij ( )O ij ((D I) 1,D 1,I,D)V ij ( ) For the two-component case R [F] ij = k (D k F 1 ) U ij +k O ij V ij + k (D k F 2 ) U ij +k O ij V ij (235) + k U ij O ij (D k F 1 ) V ij +k + k U ij O ij (D k F 2 ) V ij +k (236) Example 221 The KvM lattice (23) has recursion operator R = (I+D)( D D 1 )(D I) 1 1 I = D 1 +( ++1 )I+ D+ (+1 1 )(D I) 1 1 I Example 222 The Toda lattice (24) has recursion operator ( I D 1 I+( 1 )(D I) 1 1 (237) R = I D +1 I+ ( +1 )(D I) 1 1 I I In Section 7 we will give an algorithm for the computation of scalar recursion operators like (236) In Section 8 we cover the matrix case and show how (237) is computed The algorithms complement those for recursion operators of PDEs presented in [HG99] and elsewhere in these proceedings [BHS04] We now introduce two powerful tools which will be used in the computation of densities and fluxes 3 The Discrete Variational Derivative (Euler Operator) Definition 31 A function E n ( ) is a total difference if there exists another function J n ( ), such that E n = J n = (D I)J n Theorem 32 A necessary and sufficient condition for a function E n with positive shifts up to level p 0, to be a total difference is that (31) L (0) (E n ) = 0, )

7 DENSITIES, SYMMETRIES, AND RECURSION OPERATORS FOR NONLINEAR DDES 7 where L (0) is the discrete variational derivative (Euler operator) [ASY00] defined by (32) L (0) = p 0 ( D k ) = (I+D 1 +D 2 + +D p0 ) k=0 A proof is given in eg [HH03] Remark 33 ToverifythatanexpressionE( q,,,,+p )involving negative shifts is a total difference, one must first remove the negative shifts by replacing E n by Ẽn = D q E n Applied to Ẽn, (32) terminates at p 0 = p+q We now introduce a tool to invert the total difference operator = D I 4 The Discrete Homotopy Operator Given is an expression E n (free of negative shifts) We assume that one has verified that E n KerL (0), ie L (0) (E n ) = 0 Consequently, E n Im So, J n such that E n = J n To compute J n = 1 (E n ) one must invert the operator = D I Working with the formal inverse, (41) 1 = D 1 +D 2 +D 3 +, is impractical, perhaps impossible We therefore present the (discrete) homotopy operator which circumvents the above infinite series In analogy to the continuous case [O93], we first introduce the discrete higher Euler operators Definition 41 The discrete higher Euler operators are defined by (42) L u (i) n = p i ( ) k ( D k ) i (43) (44) (45) The higher Euler operators all terminate at some maximal shifts p i k=i Example 42 For scalar component, the first higher Euler operators are: L (0) L (1) L (2) = (I+D 1 +D 2 +D 3 + +D p0 ), = (D 1 +2D 2 +3D 3 +4D 4 + +p 1 D p1 ), = (D 2 +3D 3 +6D 4 +10D p 2(p 2 1)D p2 ) Note that (43) coincides with (32) Similar formulae hold for L (i) Next, we introduce the homotopy operators For notational simplicity we show the formulae for the two component case = (, ) Definition 43 The total homotopy operator is defined as (46) H = 1 where the homotopy operators are given by 0 (j 1,n ( )[λ ]+j 2,n ( )[λ ]) dλ λ, p 1 1 p 2 1 (47) j 1,n ( ) = (D I) i ( L u (i+1) n ), j 2,n ( ) = (D I) i ( L u (i+1) n ) i=0 i=0

8 8 WILLY HEREMAN, JAN A SANDERS, JACK SAYERS, AND JING PING WANG Theorem 44 J n = 1 (E n ) can be computed as J n = H(E n ) A similar theorem and proof for continuous homotopy operators is given in [O93] By constructing a similar variational bicomplex the theorem still holds in the discrete case See [HM02, MH02] and elsewhere in these proceedings [M04] Remark 45 Since the theory for the recursion operator can be formulated in the language of Lie algebra complexes, the results in [SW01a] are immediately applicable This gives one explicit conditions from which the existence of infinite hierarchies of local symmetries and cosymmetries can be concluded Also, the problems with the definition of a recursion operator as treated in [SW01b] also play a role in the discrete case Remark 46 Note that p 1 and p 2 in (47) are the highest shifts for and in E n Furthermore, j r,n ( )[λ ] means that in j r,n ( ) one replaces λ, and, of course, +1 λ+1, +2 λ+2, etc Remark 47 In practice, one does not compute the definite integral(46) since the evaluation at boundary λ = 0 may cause problems Instead, one computes the indefinite integral and evaluates the result at λ = 1 5 Algorithm to Compute Densities and Fluxes As an example, we compute the density-flux pair (221) for (24) Assuming that the weights (28) are computed and the rank of the density is selected (R = 3 here), the algorithm has three steps Step 1: Construct the form of the density List all monomials in and of rank 3 or less: G={ 3, 2,,, } Next, for each monomial in G, introduce the correct number of t-derivatives so that each term has weight 3 Thus, using (24), d 0 dt 0( 3 ) = u 3 d 0 n, dt 0( ) =, d dt ( 2 d ) = 2 = 2 1 2, dt () = = +1, d 2 dt 2() = d dt ( ) = d dt ( 1 ) = Gather all terms in H = { 3, 1,, 1 1,+1 } Identify members belonging to the same equivalence classes and replace them by their main representatives For example, 1 +1, so both are replaced by 1 So, H is replaced by I = { 3, 1, }, which has the building blocks of the density Linear combination of the monomials in I with constant coefficients c i gives (51) ρ n = c 1 3 +c 2 1 +c 3 Step 2: Determine the coefficients Require that (213) holds Compute D t ρ n Use (24) to remove and and their shifts Thus, (52) E n = D t ρ n = (3c 1 c 2 )u 2 n 1 +(c 3 3c 1 )u 2 n +(c 3 c 2 ) 1 +c c 2 v 2 n 1 c 3 +1 c 3 v 2 n

9 DENSITIES, SYMMETRIES, AND RECURSION OPERATORS FOR NONLINEAR DDES 9 Compute Ẽn = DE n First, apply (32) for component to Ẽn : (53) L (0) (Ẽn) = (I+D 1 +D 2 )(Ẽn) Second, apply (32) for component to Ẽn : (54) = 2(3c 1 c 2 ) 1 +2(c 3 3c 1 ) +(c 2 c 3 ) 1 1 +(c 2 c 3 )+1 L (0) (Ẽn) = (I+D 1 )(Ẽn) = (3c 1 c 2 )u 2 n+1 +(c 3 c 2 )+1 +(c 2 c 3 ) +1 +2(c 2 c 3 ) +(c 3 3c 1 )u 2 n +(c 3 c 2 ) 1 Both (53) and (54) must vanish identically Solve the linear system (55) S = {3c 1 c 2 = 0,c 3 3c 1 = 0,c 2 c 3 = 0} The solution is 3c 1 = c 2 = c 3 Substituting c 1 = 1 3,c 2 = c 3 = 1, into (51) (56) ρ n = ( 1 + ) Step 3: Compute the flux In view of (213), we will compute J n = 1 ( E n ) with the homotopy operator introduced in Section 4 Alternative methods are described in [GH98, HH03] Insert c 1 = 1 3,c 2 = c 3 = 1 into (52) and compute (57) Ẽ n = DE n = +1 +v 2 n v 2 n+1 We apply formulae (47) to Ẽn The pieces are listed in Tables 1 and 2 i ( Ẽn) (D I) i ( L (i+1) ( Ẽn)) L (i+1) Table 1 Computation of j 1,n ( )( Ẽn) i L (i+1) ( Ẽn) (D I) i ( L v (i+1) n ( Ẽn)) v 2 n Table 2 Computation of j 2,n ( )( Ẽn) Adding the terms in the right columns in Table 1 and Table 2, (58) j 1,n ( )( Ẽn) = 2 +1, j 2,n ( )( Ẽn) = +1 +2vn 2 Thus, the homotopy operator (46) gives J n = = (59) = +1 +vn 2 0 (j 1,n ( )( Ẽn)[λ ]+j 2,n ( )( Ẽn)[λ ]) dλ λ (3λ λv 2 n ) dλ

10 10 WILLY HEREMAN, JAN A SANDERS, JACK SAYERS, AND JING PING WANG After a backward shift, J n = D 1 Jn, we obtain the final result: (510) ρ n = 1 3 u3 n + ( 1 + ), J n = 1 1 +v 2 n 1 6 Algorithm to Compute Generalized Symmetries As an example, we compute the symmetry (230) of rank (3,4) for (24) The two steps of the algorithm [GH99] are similar to those in Section 5 Step 1: Construct the form of the symmetry Start by listing all monomials in and of ranks 3 and 4, or less: K 1 = {u 3 n,u2 n,,, }, K 2 = {u 4 n,u3 n,u2 n,u 2 n,,,v 2 n,} As in Step 1 in Section 5, for each monomial in K 1 and K 2, introduce the necessary t-derivatives so that each term has rank 3 and 4, respectively At the same time, use (24) to remove all t derivatives Doing so, based on K 1 we obtain (61) L 1 = {u 3 n, 1 1, 1,,+1 } Similarly, based on K 2, we get (62) L 2 = {u 4 n,u2 n 1 1, 1 1,u 2 n 1, 2 1,v 2 n 1,u2 n, +1,u 2 n+1, 1,v 2 n, +1 } In contrast to the strategy for densities, we do not introduce the main representatives, but linearly combine the monomials in L 1 and L 2 to get the form of the symmetry: (63) G 1 = c 1 u 3 n +c c 3 1 +c 4 +c 5 +1, G 2 = c 6 u 4 n +c 7 u 2 n 1 1 +c c 9 u 2 n 1 +c c 11 v 2 n 1 +c 12 u 2 n +c c 14 u 2 n+1 +c c 16 v 2 n +c 17 +1, with constant coefficients c i Step 2: Determine the coefficients To compute the coefficients c i we require that (222) holds on solutions of (21) Compute D t G 1 and D t G 2 and use (24) to remove,, and their shifts This gives the left hand sides of (222) Use (226) to get the right hand sides of (222): (64) F 1 ( )[G] = D 1 G 2 IG 2, F 2 ( )[G] = IG 1 DG 1 +( +1 )IG 2 Substitute (63) into (64) and equate the corresponding left and right hand sides Since all monomials in, and their shifts are independent, one obtains the linear system that determines the coefficients c i The solution is (65) c 1 = c 6 = c 7 = c 8 = c 9 = c 10 = c 11 = c 13 = c 16 = 0, c 2 = c 3 = c 4 = c 5 = c 12 = c 14 = c 15 = c 17 With the choice c 17 = 1, the symmetry (63) finally becomes (66) G 1 = ( ++1 ) 1 ( 1 ), G 2 = (u 2 n+1 u2 n )

11 DENSITIES, SYMMETRIES, AND RECURSION OPERATORS FOR NONLINEAR DDES 11 7 Algorithm to Compute Scalar Recursion Operators In this section we show how to compute the recursion operator (236) of (24) Again, we will use the concept of rank invariance to construct a candidate recursion operator The defining equation (232) is then used to determine the coefficients We observe that (236) in expanded form naturally splits into two pieces: (71) R = R 0 +R 1, where R 0 contains only terms with shift operators D 1,I, and D, and R 1 has terms involving (D I) 1 Step 1: Determine the rank of the recursion operator Inviewof(231)andassumingthat thesymmetriesarelinkedconsecutively(s = 1), the recursion operator R has rank (72) R = rank R = rank G (2) rank G (1) = 3 2 = 1, where we used (227) and (228) If the assumption turns out to be correct, then the recursion operator has rank 1 If the assumption were falls because symmetries are not linked consecutively, then R must be adjusted and the subsequent steps must be repeated See [BHS04] for examples of PDEs for which that happens Step 2: Determine the form of the recursion operator We split this into two sub-steps (i) Determine the form of R 0 The candidate R 0 is a sum of terms involving D 1,I, and D The coefficients of these terms are linear combinations of 1,, and +1, so that all terms have the correct rank (73) R 0 = (c 1 1 +c 2 +c 3 +1 )D 1 +(c 4 1 +c 5 +c 6 +1 )I +(c 7 1 +c 8 +c 9 +1 )D, where the c i are constant coefficients A few remarks are in place First, in R 0 we moved the operators D 1,I, and D all the way to the right Second, the maximum up-shift and down-shift operator that should be included can be determined by comparing two consecutive symmetries Indeed, if the maximum up-shift in the first symmetry is +p, and the maximum up-shift in the next symmetry is +p+r, then R 0 must have D,D 2,,D r The same line of reasoning determines the minimum down-shift operator to be included In our example, there is no need to include terms in D 2,D 2, etc Third, the coefficients of the operators can be restricted to linear combinations of the terms appearing in F Hence, no terms in ±2,±3 and so on occur in (73) (ii) Determine the form of R 1 As in the continuous case [HG99], R 1 is a linear combination (with constant coefficients c jk ) of sums of all suitable products of symmetries and covariants (Fréchet derivatives of densities) sandwiching (D I) 1, ie (74) R 1 = j k c jk G (j) (D I) 1 ρ (k) n For (23), G (1) in (227) and ρ (0) n = ln( ) in (215) are the only suitable pair Indeed, using (224) we have ρ (0) n = ln( ) = 1 I, which has rank -1 Combined

12 12 WILLY HEREMAN, JAN A SANDERS, JACK SAYERS, AND JING PING WANG with G (1) of rank 2, we have a term of rank 1 Other combinations of symmetries and covariants would exceed rank 1 Therefore, (75) R 1 = c 10 (+1 1 )(D I) 1 ( 1 )I where c 10 a constant coefficient Using (71) and renaming c 10 to c 10, R = (c 1 1 +c 2 +c 3 +1 )D 1 +(c 4 1 +c 5 +c 6 +1 )I (76) +(c 7 1 +c 8 +c 9 +1 )D+c 10 (+1 1 )(D I) 1 ( 1 )I is the candidate recursion operator for (23) Step 3: Determine the coefficients Starting from (76), we use (234) with F = (+1 1 ) to compute R [F] The partial results will not be shown due to length Using (223), we compute (77) F = D 1 +(+1 1 )I+ D Then we compute R F and F R After substituting the pieces into (232) we simplify the resulting expression using rules such as (78) (D I) 1 D = D(D I) 1 = I+(D I) 1, (D I) 1 D 1 = D 1 (D I) 1 = D 1 +(D I) 1 We further simplify by recursively using formulas like (79) DU( )(D I) 1 V( )I = U(+1 )V( )I+U(+1 )(D I) 1 V( )I, (D I) 1 U( )V( )D = U( 1 )V( 1 )I+(D I) 1 U( 1 )V( 1 )I Finally, we equate like terms to obtain a linear system for the c i Substituting the solution (710) c 1 = c 3 = c 4 = c 7 = c 9 = 0, c 2 = c 5 = c 6 = c 8 = c 10 = 1 into (76) we obtain the final result (711) R = D 1 +( ++1 )I+ D+ (+1 1 )(D I) 1 1 I A straightforward computation confirms that RG (1) = G (2) with G (1) in (227) and G (2) in (228) 8 Algorithm to Compute Matrix Recursion Operators We construct the recursion operator(237) for(24) Now all the terms in (232) are 2 2 matrix operators Step 1: Determine the rank of the recursion operator The difference in rank of symmetries is again used to compute the rank of the elements of the recursion operator Using (28), (229) and (230), ( (81) rankg (1) 2 = 3 ), rankg (2) = Assuming that RG (1) = G (2), we use the formula (82) rankr ij = rankg (k+1) i ( 3 4 rankg (k) j, )

13 DENSITIES, SYMMETRIES, AND RECURSION OPERATORS FOR NONLINEAR DDES 13 to compute a rank matrix associated to the operator ( ) 1 0 (83) rank R = 2 1 Step 2: Determine the form of the recursion operator As in the scalar case, we build a candidate R 0 : ( ) (R0 ) R 0 = 11 (R 0 ) (84) 12, (R 0 ) 21 (R 0 ) 22 with (R 0 ) 11 = (c 1 +c 2 +1 )I, (R 0 ) 12 = c 3 D 1 +c 4 I, (R 0 ) 21 = (c 5 u 2 n +c c 7 u 2 n+1 +c 8 1 +c 9 )I +(c 10 u 2 n +c c 12 u 2 n+1 +c c 14 )D, (R 0 ) 22 = (c 15 +c )I Analogous to the scalar case, the elements of matrix R 1 are linear combinations with constant coefficients of all suitable products of symmetries and covariants sandwiching (D I) 1 Hence, (85) j k c jk G (j) (D I) 1 ρ (k) n, where denotes the matrix outer product, defined as ( ) ( G (j) ) 1 G (j) (D I) 1 (ρ (k) n,1 ρ (k) n,2 = 2 G (j) 1 (D I) 1 ρ (k) n,1 G (j) 1 (D I) 1 ρ (k) n,2 G (j) 2 (D I) 1 ρ (k) n,1 G (j) 2 (D I) 1 ρ (k) n,2 Only the pair (G (1),ρ (0) n ) can be used, otherwise the ranks in (83) would be exceeded Using (226) and (218) we compute (86) ρ (0) n = ( 1 0 I ), Therefore, using (85) and renaming c 10 to c 17, ( 0 c17 ( 1 )(D I) 1 1 (87) R 1 = I 0 c 17 ( +1 )(D I) 1 1 I Adding (84) and (87) we obtain ( ) R11 R (88) R = 12, R 21 R 22 with R 11 = (c 1 +c 2 +1 )I ) R 12 = c 3 D 1 +c 4 I+c 17 ( 1 )(D I) 1 1 I, R 21 = (c 5 u 2 n +c c 7 u 2 n+1 +c 8 1 +c 9 )I +(c 10 u 2 n +c c 12 u 2 n+1 +c c 14 )D, R 22 = (c 15 +c )I+c 17 ( +1 )(D I) 1 1 I )

14 14 WILLY HEREMAN, JAN A SANDERS, JACK SAYERS, AND JING PING WANG Step 3: Determine the coefficients All the terms in (232) need to be computed The strategy is similar to the scalar case, yet the computations are much more cumbersome Omitting the details, the result is: c 2 = c 5 = c 6 = c 7 = c 8 = c 10 = c 11 = c 12 = c 13 = c 15 = 0, c 1 = c 3 = c 4 = c 9 = c 14 = c 16 = 1, and c 17 = 1 Substitution the constants into (88) gives ( I D 1 I+( 1 )(D I) 1 1 ) (89) R = I I D +1 I+ ( +1 )(D I) 1 1 It is straightforward to verify that RG (1) = G (2) with G (1) in (229) and G (2) in (230) 9 More Examples Example 91 The modified Volterra lattice [ASY00, HH03], (91) = u 2 n (+1 1 ), has two non-polynomial densities ρ (0) n = 1 and ρ (1) n = ln( ), and infinitely many polynomial densities The first two symmetries, (92) (93) G (1) = u 2 n(+1 1 ), G (2) = u 2 nu 2 n+1( ++2 ) u 2 n 1u 2 n( 2 + ), are linked by the recursion operator (94) R = u 2 n D I+u 2 n D+2u2 n (+1 1 )(D I) 1 1 I Example 92 The AL lattice (22) has infinitely many densities [GH98] and symmetries [GH99] The recursion operator is of the form (88) with R 11 = P n D I 1 P n 1 P n I, R 12 = 1 I 1 1 I 1 P n 1 P n I, R 21 = +1 I I++1 P n 1 P n I, R 22 = ( )I+P n D+ 1 1 I++1 P n 1 (95) I, P n where P n = 1+ and = D I This recursion operator has an inverse, which is quite exceptional The investigation of the recursion operator structure of the AL lattice is work in progress Acknowledgements M Hickman and B Deconinck are gratefully acknowledged for valuable discussions D Baldwin is thanked for proof reading the manuscript References [AC91] M J Ablowitz and P A Clarkson, Solitons, nonlinear evolution equations and inverse scattering, London Math Soc Lect Note Ser 149, Cambridge University Press, Cambridge, U K, 1991 [AL75] M J Ablowitz and J F Ladik, Nonlinear differential-difference equations, J Math Phys 16 (1975), [ASY00] V E Adler, A B Shabat and R I Yamilov, Symmetry approach to the integrability problem, Theor Math Phys 125 (2000), I

15 DENSITIES, SYMMETRIES, AND RECURSION OPERATORS FOR NONLINEAR DDES 15 [BHS04] D Baldwin, W Hereman and J Sayers, Symbolic algorithms for the Painlevé test, special solutions and recursion operators for nonlinear PDEs, present proceedings [GH97a] Ü Göktaş and W Hereman, Symbolic computation of conserved densities for systems of nonlinear evolution equations, J Symb Comp 24 (1997), [GH97b], Mathematica package InvariantsSymmetriesm is available since 1997 from [GH98], Computation of conserved densities for nonlinear lattices, Physica D 123 (1998), [GH99], Algorithmic computation of higher-order symmetries for nonlinear evolution and lattice equations, Adv Comp Math 11 (1999), [HG99] W Hereman and Ü Göktaş, Integrability Tests for Nonlinear Evolution Equations Computer Algebra Systems: A Practical Guide, (M Wester, ed), Wiley, New York, 1999, pp [HGCM98] W Hereman, Ü Göktaş, M D Colagrosso and A J Miller, Algorithmic integrability tests for nonlinear differential and lattice equations, Comp Phys Comm 115 (1998), [HH03] M S Hickman and W A Hereman, Computation of densities and fluxes of nonlinear differential-difference equations, Proc Roy Soc Lond A (2003), in press [HM02] P E Hydon and E L Mansfield, A variational complex for difference equations, Found Comp Math (2003), in press [KvM75] M Kac and P van Moerbeke, On an explicitly soluble system of nonlinear differential equations related to certain Toda lattices, Adv Math 16 (1975), [MH02] E L Mansfield and P E Hydon, On a variational complex for difference equations The Geometrical Study of Differential Equations Proc NSF-CBMS Conf (J A Leslie and T P Robart, eds), American Mathematical Society, Providence, Rhode Island, 2002, pp [M04] E L Mansfield and R Quispel, Towards a variational complex for the finite element method, present proceedings [OZF93] W Oevel, H Zhang and B Fuchssteiner, Mastersymmetries and multi-hamiltonian formulations for some integrable lattice systems, Prog Theor Phys 81 (1989), [O93] P J Olver, Applications of lie groups to differential equations, 2nd ed, Graduate Texts in Mathematics 107, Springer, New York, 1993 [SW01a] J A Sanders and J P Wang, Integrable systems and their recursion operators, Nonlinear Anal 47 (2001), [SW01b], On recursion operators, Physica D 149 (2001), 1 10 [T81] M Toda, Theory of nonlinear lattices, Springer, Berlin, 1981 [W98] J P Wang, Symmetries and conservation laws of evolution equations, PhD Thesis, Thomas Stieltjes Institute for Mathematics, Amsterdam, 1998 Department of Mathematical and Computer Sciences, Colorado School of Mines, Golden, Colorado address: whereman@minesedu Department of Mathematics, Faculty of Sciences, Free University, De Boelelaan 1081a, 1081 HV Amsterdam, The Netherlands address: jansa@csvunl Physics Department , California Institute of Technology, Pasadena, California address: jack@itscaltechedu Department of Mathematics, Brock University, St Catharines, Ontario, Canada L2S 3A1 address: jwang@spartanacbrockuca

Symbolic Computation of Recursion Operators of Nonlinear Differential-Difference Equations

Symbolic Computation of Recursion Operators of Nonlinear Differential-Difference Equations Symbolic Computation of Recursion Operators of Nonlinear Differential-Difference Equations Ünal Göktaş and Willy Hereman Department of Computer Engineering Turgut Özal University, Keçiören, Ankara, Turkey

More information

Symbolic Computation of Conservation Laws, Generalized Symmetries, and Recursion Operators for Nonlinear Differential-Difference Equations

Symbolic Computation of Conservation Laws, Generalized Symmetries, and Recursion Operators for Nonlinear Differential-Difference Equations Symbolic Computation of Conservation Laws, Generalized Symmetries, and Recursion Operators for Nonlinear Differential-Difference Equations Ünal Götaş and Willy Hereman Abstract Algorithms for the symbolic

More information

Continuous and Discrete Homotopy Operators with Applications in Integrability Testing. Willy Hereman

Continuous and Discrete Homotopy Operators with Applications in Integrability Testing. Willy Hereman Continuous and Discrete Homotopy Operators with Applications in Integrability Testing Willy Hereman Department of Mathematical and Computer Sciences Colorado School of Mines Golden, Colorado, USA http://www.mines.edu/fs

More information

The first three (of infinitely many) conservation laws for (1) are (3) (4) D t (u) =D x (3u 2 + u 2x ); D t (u 2 )=D x (4u 3 u 2 x +2uu 2x ); D t (u 3

The first three (of infinitely many) conservation laws for (1) are (3) (4) D t (u) =D x (3u 2 + u 2x ); D t (u 2 )=D x (4u 3 u 2 x +2uu 2x ); D t (u 3 Invariants and Symmetries for Partial Differential Equations and Lattices Λ Ünal Göktaοs y Willy Hereman y Abstract Methods for the computation of invariants and symmetries of nonlinear evolution, wave,

More information

Symbolic Computation of Conserved Densities and Symmetries of Nonlinear Evolution and Differential-Difference Equations WILLY HEREMAN

Symbolic Computation of Conserved Densities and Symmetries of Nonlinear Evolution and Differential-Difference Equations WILLY HEREMAN . Symbolic Computation of Conserved Densities and Symmetries of Nonlinear Evolution and Differential-Difference Equations WILLY HEREMAN Joint work with ÜNAL GÖKTAŞ and Grant Erdmann Graduate Seminar Department

More information

Symbolic Computation of. Conserved Densities and Fluxes for. Nonlinear Systems of Differential-Difference. Equations

Symbolic Computation of. Conserved Densities and Fluxes for. Nonlinear Systems of Differential-Difference. Equations Symbolic Computation of Conserved Densities and Fluxes for Nonlinear Systems of Differential-Difference Equations by Holly Eklund Copyright c 003 Holly Eklund, All rights reserved. ii A thesis submitted

More information

arxiv:nlin/ v1 [nlin.si] 16 Nov 2003

arxiv:nlin/ v1 [nlin.si] 16 Nov 2003 Centre de Recherches Mathématiques CRM Proceedings and Lecture Notes arxiv:nlin/0311028v1 [nlin.si] 16 Nov 2003 Symbolic algorithms for the Painlevé test, special solutions, and recursion operators for

More information

ON THE SYMMETRIES OF INTEGRABLE PARTIAL DIFFERENCE EQUATIONS

ON THE SYMMETRIES OF INTEGRABLE PARTIAL DIFFERENCE EQUATIONS Proceedings of the International Conference on Difference Equations, Special Functions and Orthogonal Polynomials, World Scientific (2007 ON THE SYMMETRIES OF INTEGRABLE PARTIAL DIFFERENCE EQUATIONS ANASTASIOS

More information

RESEARCH ARTICLE. A symbolic algorithm for computing recursion operators of nonlinear PDEs

RESEARCH ARTICLE. A symbolic algorithm for computing recursion operators of nonlinear PDEs International Journal of Computer Mathematics Vol. 00, No. 00, Month 200, 1 27 RESEARCH ARTICLE A symbolic algorithm for computing recursion operators of nonlinear PDEs D.E. Baldwin and W. Hereman Department

More information

Application of Euler and Homotopy Operators to Integrability Testing

Application of Euler and Homotopy Operators to Integrability Testing Application of Euler and Homotopy Operators to Integrability Testing Willy Hereman Department of Applied Mathematics and Statistics Colorado School of Mines Golden, Colorado, U.S.A. http://www.mines.edu/fs

More information

Algorithmic Computation of Symmetries, Invariants and Recursion Operators for Systems of. Nonlinear Evolution and Differential-Difference.

Algorithmic Computation of Symmetries, Invariants and Recursion Operators for Systems of. Nonlinear Evolution and Differential-Difference. Algorithmic Computation of Symmetries, Invariants and Recursion Operators for Systems of Nonlinear Evolution and Differential-Difference Equations by Ünal GÖKTAŞ A thesis submitted to the Faculty and the

More information

Symbolic algorithms for the Painlevé test, special solutions, and recursion operators for nonlinear PDEs.

Symbolic algorithms for the Painlevé test, special solutions, and recursion operators for nonlinear PDEs. Centre de Recherches Mathématiques CRM Proceedings and Lecture Notes Symbolic algorithms for the Painlevé test, special solutions, and recursion operators for nonlinear PDEs. Douglas Baldwin, Willy Hereman,

More information

About Integrable Non-Abelian Laurent ODEs

About Integrable Non-Abelian Laurent ODEs About Integrable Non-Abelian Laurent ODEs T. Wolf, Brock University September 12, 2013 Outline Non-commutative ODEs First Integrals and Lax Pairs Symmetries Pre-Hamiltonian Operators Recursion Operators

More information

Symmetries and Group Invariant Reductions of Integrable Partial Difference Equations

Symmetries and Group Invariant Reductions of Integrable Partial Difference Equations Proceedings of 0th International Conference in MOdern GRoup ANalysis 2005, 222 230 Symmetries and Group Invariant Reductions of Integrable Partial Difference Equations A. TONGAS, D. TSOUBELIS and V. PAPAGEORGIOU

More information

INTEGRABLE DISCRETIZATION OF COUPLED ABLOWITZ-LADIK EQUATIONS WITH BRANCHED DISPERSION

INTEGRABLE DISCRETIZATION OF COUPLED ABLOWITZ-LADIK EQUATIONS WITH BRANCHED DISPERSION v..1r0180507 *018.11.15#5f9cb4 INTEGRABLE DISCRETIZATION OF COUPLED ABLOWITZ-LADIK EQUATIONS WITH BRANCHED DISPERSION CORINA N. BABALIC University of Craiova 13 A.I. Cuza, 00585, Craiova, Romania E-mail:

More information

Solving Nonlinear Wave Equations and Lattices with Mathematica. Willy Hereman

Solving Nonlinear Wave Equations and Lattices with Mathematica. Willy Hereman Solving Nonlinear Wave Equations and Lattices with Mathematica Willy Hereman Department of Mathematical and Computer Sciences Colorado School of Mines Golden, Colorado, USA http://www.mines.edu/fs home/whereman/

More information

2. Examples of Integrable Equations

2. Examples of Integrable Equations Integrable equations A.V.Mikhailov and V.V.Sokolov 1. Introduction 2. Examples of Integrable Equations 3. Examples of Lax pairs 4. Structure of Lax pairs 5. Local Symmetries, conservation laws and the

More information

Fibonacci tan-sec method for construction solitary wave solution to differential-difference equations

Fibonacci tan-sec method for construction solitary wave solution to differential-difference equations ISSN 1 746-7233, England, UK World Journal of Modelling and Simulation Vol. 7 (2011) No. 1, pp. 52-57 Fibonacci tan-sec method for construction solitary wave solution to differential-difference equations

More information

The Complete Set of Generalized Symmetries for the Calogero Degasperis Ibragimov Shabat Equation

The Complete Set of Generalized Symmetries for the Calogero Degasperis Ibragimov Shabat Equation Proceedings of Institute of Mathematics of NAS of Ukraine 2002, Vol. 43, Part 1, 209 214 The Complete Set of Generalized Symmetries for the Calogero Degasperis Ibragimov Shabat Equation Artur SERGYEYEV

More information

A MULTI-COMPONENT LAX INTEGRABLE HIERARCHY WITH HAMILTONIAN STRUCTURE

A MULTI-COMPONENT LAX INTEGRABLE HIERARCHY WITH HAMILTONIAN STRUCTURE Pacific Journal of Applied Mathematics Volume 1, Number 2, pp. 69 75 ISSN PJAM c 2008 Nova Science Publishers, Inc. A MULTI-COMPONENT LAX INTEGRABLE HIERARCHY WITH HAMILTONIAN STRUCTURE Wen-Xiu Ma Department

More information

arxiv: v2 [math-ph] 24 Feb 2016

arxiv: v2 [math-ph] 24 Feb 2016 ON THE CLASSIFICATION OF MULTIDIMENSIONALLY CONSISTENT 3D MAPS MATTEO PETRERA AND YURI B. SURIS Institut für Mathemat MA 7-2 Technische Universität Berlin Str. des 17. Juni 136 10623 Berlin Germany arxiv:1509.03129v2

More information

arxiv: v1 [nlin.si] 25 Mar 2009

arxiv: v1 [nlin.si] 25 Mar 2009 Linear quadrilateral lattice equations and multidimensional consistency arxiv:0903.4428v1 [nlin.si] 25 Mar 2009 1. Introduction James Atkinson Department of Mathematics and Statistics, La Trobe University,

More information

B.7 Lie Groups and Differential Equations

B.7 Lie Groups and Differential Equations 96 B.7. LIE GROUPS AND DIFFERENTIAL EQUATIONS B.7 Lie Groups and Differential Equations Peter J. Olver in Minneapolis, MN (U.S.A.) mailto:olver@ima.umn.edu The applications of Lie groups to solve differential

More information

On Local Time-Dependent Symmetries of Integrable Evolution Equations

On Local Time-Dependent Symmetries of Integrable Evolution Equations Proceedings of Institute of Mathematics of NAS of Ukraine 2000, Vol. 30, Part 1, 196 203. On Local Time-Dependent Symmetries of Integrable Evolution Equations A. SERGYEYEV Institute of Mathematics of the

More information

Symbolic Computation of Conservation Laws of Nonlinear Partial Differential Equations in Multi-dimensions

Symbolic Computation of Conservation Laws of Nonlinear Partial Differential Equations in Multi-dimensions . Symbolic Computation of Conservation Laws of Nonlinear Partial Differential Equations in Multi-dimensions Willy Hereman Department of Mathematical and Computer Sciences Colorado School of Mines Golden,

More information

Breaking soliton equations and negative-order breaking soliton equations of typical and higher orders

Breaking soliton equations and negative-order breaking soliton equations of typical and higher orders Pramana J. Phys. (2016) 87: 68 DOI 10.1007/s12043-016-1273-z c Indian Academy of Sciences Breaking soliton equations and negative-order breaking soliton equations of typical and higher orders ABDUL-MAJID

More information

Symmetry Methods for Differential and Difference Equations. Peter Hydon

Symmetry Methods for Differential and Difference Equations. Peter Hydon Lecture 2: How to find Lie symmetries Symmetry Methods for Differential and Difference Equations Peter Hydon University of Kent Outline 1 Reduction of order for ODEs and O Es 2 The infinitesimal generator

More information

Canonical Forms for BiHamiltonian Systems

Canonical Forms for BiHamiltonian Systems Canonical Forms for BiHamiltonian Systems Peter J. Olver Dedicated to the Memory of Jean-Louis Verdier BiHamiltonian systems were first defined in the fundamental paper of Magri, [5], which deduced the

More information

MATH 590: Meshfree Methods

MATH 590: Meshfree Methods MATH 590: Meshfree Methods Chapter 6: Scattered Data Interpolation with Polynomial Precision Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu

More information

The Higher Dimensional Bateman Equation and Painlevé Analysis of Nonintegrable Wave Equations

The Higher Dimensional Bateman Equation and Painlevé Analysis of Nonintegrable Wave Equations Symmetry in Nonlinear Mathematical Physics 1997, V.1, 185 192. The Higher Dimensional Bateman Equation and Painlevé Analysis of Nonintegrable Wave Equations Norbert EULER, Ove LINDBLOM, Marianna EULER

More information

The Modified (G /G)-Expansion Method for Nonlinear Evolution Equations

The Modified (G /G)-Expansion Method for Nonlinear Evolution Equations The Modified ( /-Expansion Method for Nonlinear Evolution Equations Sheng Zhang, Ying-Na Sun, Jin-Mei Ba, and Ling Dong Department of Mathematics, Bohai University, Jinzhou 11000, P. R. China Reprint requests

More information

Symbolic Computation of Conservation Laws of Nonlinear Partial Differential Equations in Multiple Space Dimensions

Symbolic Computation of Conservation Laws of Nonlinear Partial Differential Equations in Multiple Space Dimensions Symbolic Computation of Conservation Laws of Nonlinear Partial Differential Equations in Multiple Space Dimensions Willy Hereman and Douglas Poole Department of Mathematical and Computer Sciences Colorado

More information

Euler characteristic of the truncated order complex of generalized noncrossing partitions

Euler characteristic of the truncated order complex of generalized noncrossing partitions Euler characteristic of the truncated order complex of generalized noncrossing partitions D. Armstrong and C. Krattenthaler Department of Mathematics, University of Miami, Coral Gables, Florida 33146,

More information

arxiv:nlin/ v2 [nlin.si] 9 Oct 2002

arxiv:nlin/ v2 [nlin.si] 9 Oct 2002 Journal of Nonlinear Mathematical Physics Volume 9, Number 1 2002), 21 25 Letter On Integrability of Differential Constraints Arising from the Singularity Analysis S Yu SAKOVICH Institute of Physics, National

More information

arxiv:nlin/ v2 [nlin.si] 15 Sep 2004

arxiv:nlin/ v2 [nlin.si] 15 Sep 2004 Integrable Mappings Related to the Extended Discrete KP Hierarchy ANDREI K. SVININ Institute of System Dynamics and Control Theory, Siberian Branch of Russian Academy of Sciences, P.O. Box 1233, 664033

More information

Symmetry Classification of KdV-Type Nonlinear Evolution Equations

Symmetry Classification of KdV-Type Nonlinear Evolution Equations Proceedings of Institute of Mathematics of NAS of Ukraine 2004, Vol. 50, Part 1, 125 130 Symmetry Classification of KdV-Type Nonlinear Evolution Equations Faruk GÜNGÖR, Victor LAHNO and Renat ZHDANOV Department

More information

TWELVE LIMIT CYCLES IN A CUBIC ORDER PLANAR SYSTEM WITH Z 2 -SYMMETRY. P. Yu 1,2 and M. Han 1

TWELVE LIMIT CYCLES IN A CUBIC ORDER PLANAR SYSTEM WITH Z 2 -SYMMETRY. P. Yu 1,2 and M. Han 1 COMMUNICATIONS ON Website: http://aimsciences.org PURE AND APPLIED ANALYSIS Volume 3, Number 3, September 2004 pp. 515 526 TWELVE LIMIT CYCLES IN A CUBIC ORDER PLANAR SYSTEM WITH Z 2 -SYMMETRY P. Yu 1,2

More information

Compatible Hamiltonian Operators for the Krichever-Novikov Equation

Compatible Hamiltonian Operators for the Krichever-Novikov Equation arxiv:705.04834v [math.ap] 3 May 207 Compatible Hamiltonian Operators for the Krichever-Novikov Equation Sylvain Carpentier* Abstract It has been proved by Sokolov that Krichever-Novikov equation s hierarchy

More information

CLASSIFICATION OF NON-ABELIAN CHERN-SIMONS VORTICES

CLASSIFICATION OF NON-ABELIAN CHERN-SIMONS VORTICES CLASSIFICATION OF NON-ABELIAN CHERN-SIMONS VORTICES arxiv:hep-th/9310182v1 27 Oct 1993 Gerald V. Dunne Department of Physics University of Connecticut 2152 Hillside Road Storrs, CT 06269 USA dunne@hep.phys.uconn.edu

More information

Complete integrability of a modified vector derivative nonlinear Schrödinger equation

Complete integrability of a modified vector derivative nonlinear Schrödinger equation Complete integrability of a modified vector derivative nonlinear Schrödinger equation arxiv:solv-int/9407002v1 14 Jul 1994 Ralph Willox 1, Willy Hereman 2 and Frank Verheest 3 1 Dienst Theoretische Natuurkunde,

More information

Evolutionary Hirota Type (2+1)-Dimensional Equations: Lax Pairs, Recursion Operators and Bi-Hamiltonian Structures

Evolutionary Hirota Type (2+1)-Dimensional Equations: Lax Pairs, Recursion Operators and Bi-Hamiltonian Structures Symmetry, Integrability and Geometry: Methods and Applications Evolutionary Hirota Type +-Dimensional Equations: Lax Pairs, Recursion Operators and Bi-Hamiltonian Structures Mikhail B. SHEFTEL and Devrim

More information

Lax Representations and Zero Curvature Representations by Kronecker Product

Lax Representations and Zero Curvature Representations by Kronecker Product Lax Representations and Zero Curvature Representations by Kronecker Product arxiv:solv-int/9610008v1 18 Oct 1996 Wen-Xiu Ma and Fu-Kui Guo Abstract It is showed that Kronecker product can be applied to

More information

THE NUMBER OF ORTHOGONAL CONJUGATIONS

THE NUMBER OF ORTHOGONAL CONJUGATIONS THE NUMBER OF ORTHOGONAL CONJUGATIONS Armin Uhlmann University of Leipzig, Institute for Theoretical Physics After a short introduction to anti-linearity, bounds for the number of orthogonal (skew) conjugations

More information

arxiv:nlin/ v1 [nlin.si] 13 Apr 2005

arxiv:nlin/ v1 [nlin.si] 13 Apr 2005 Diophantine Integrability R.G. Halburd Department of Mathematical Sciences, Loughborough University Loughborough, Leicestershire, LE11 3TU, UK (Dated: received January 1, revised March 7, 005) The heights

More information

MA22S3 Summary Sheet: Ordinary Differential Equations

MA22S3 Summary Sheet: Ordinary Differential Equations MA22S3 Summary Sheet: Ordinary Differential Equations December 14, 2017 Kreyszig s textbook is a suitable guide for this part of the module. Contents 1 Terminology 1 2 First order separable 2 2.1 Separable

More information

On bosonic limits of two recent supersymmetric extensions of the Harry Dym hierarchy

On bosonic limits of two recent supersymmetric extensions of the Harry Dym hierarchy arxiv:nlin/3139v2 [nlin.si] 14 Jan 24 On bosonic limits of two recent supersymmetric extensions of the Harry Dym hierarchy S. Yu. Sakovich Mathematical Institute, Silesian University, 7461 Opava, Czech

More information

Exact solutions of nonlinear partial differential equations by the method of group foliation reduction

Exact solutions of nonlinear partial differential equations by the method of group foliation reduction Konrad-Zuse-Zentrum für Informationstechnik Berlin Takustraße 7 D-14195 Berlin-Dahlem Germany STEPHEN S. ANCO 1, S. ALI, THOMAS WOLF 3 Exact solutions of nonlinear partial differential equations by the

More information

A new integrable system: The interacting soliton of the BO

A new integrable system: The interacting soliton of the BO Phys. Lett., A 204, p.336-342, 1995 A new integrable system: The interacting soliton of the BO Benno Fuchssteiner and Thorsten Schulze Automath Institute University of Paderborn Paderborn & Germany Abstract

More information

MATH 320, WEEK 7: Matrices, Matrix Operations

MATH 320, WEEK 7: Matrices, Matrix Operations MATH 320, WEEK 7: Matrices, Matrix Operations 1 Matrices We have introduced ourselves to the notion of the grid-like coefficient matrix as a short-hand coefficient place-keeper for performing Gaussian

More information

On a WKB-theoretic approach to. Ablowitz-Segur's connection problem for. the second Painleve equation. Yoshitsugu TAKEI

On a WKB-theoretic approach to. Ablowitz-Segur's connection problem for. the second Painleve equation. Yoshitsugu TAKEI On a WKB-theoretic approach to Ablowitz-Segur's connection problem for the second Painleve equation Yoshitsugu TAKEI Research Institute for Mathematical Sciences Kyoto University Kyoto, 606-8502, JAPAN

More information

arxiv: v1 [math.co] 13 May 2016

arxiv: v1 [math.co] 13 May 2016 GENERALISED RAMSEY NUMBERS FOR TWO SETS OF CYCLES MIKAEL HANSSON arxiv:1605.04301v1 [math.co] 13 May 2016 Abstract. We determine several generalised Ramsey numbers for two sets Γ 1 and Γ 2 of cycles, in

More information

Discriminative Direction for Kernel Classifiers

Discriminative Direction for Kernel Classifiers Discriminative Direction for Kernel Classifiers Polina Golland Artificial Intelligence Lab Massachusetts Institute of Technology Cambridge, MA 02139 polina@ai.mit.edu Abstract In many scientific and engineering

More information

On recursion operators for elliptic models

On recursion operators for elliptic models On recursion operators for elliptic models D K Demskoi and V V Sokolov 2 School of Mathematics and Statistics, University of New South Wales, Sydney, NSW 2052, Australia E-mail: demskoi@maths.unsw.edu.au

More information

A Lie-Group Approach for Nonlinear Dynamic Systems Described by Implicit Ordinary Differential Equations

A Lie-Group Approach for Nonlinear Dynamic Systems Described by Implicit Ordinary Differential Equations A Lie-Group Approach for Nonlinear Dynamic Systems Described by Implicit Ordinary Differential Equations Kurt Schlacher, Andreas Kugi and Kurt Zehetleitner kurt.schlacher@jku.at kurt.zehetleitner@jku.at,

More information

arxiv: v1 [math.ap] 11 Jan 2014

arxiv: v1 [math.ap] 11 Jan 2014 THE UNIFIED TRANSFORM FOR THE MODIFIED HELMHOLTZ EQUATION IN THE EXTERIOR OF A SQUARE A. S. FOKAS AND J. LENELLS arxiv:4.252v [math.ap] Jan 24 Abstract. The Unified Transform provides a novel method for

More information

arxiv:math-ph/ v1 25 Feb 2002

arxiv:math-ph/ v1 25 Feb 2002 FROM THE TODA LATTICE TO THE VOLTERRA LATTICE AND BACK arxiv:math-ph/0202037v1 25 Feb 2002 (1) PANTELIS A DAMIANOU AND RUI LOJA FERNANDES Abstract We discuss the relationship between the multiple Hamiltonian

More information

Symmetry Reductions of Integrable Lattice Equations

Symmetry Reductions of Integrable Lattice Equations Isaac Newton Institute for Mathematical Sciences Discrete Integrable Systems Symmetry Reductions of Integrable Lattice Equations Pavlos Xenitidis University of Patras Greece March 11, 2009 Pavlos Xenitidis

More information

Bäcklund transformation and special solutions for Drinfeld Sokolov Satsuma Hirota system of coupled equations

Bäcklund transformation and special solutions for Drinfeld Sokolov Satsuma Hirota system of coupled equations arxiv:nlin/0102001v1 [nlin.si] 1 Feb 2001 Bäcklund transformation and special solutions for Drinfeld Sokolov Satsuma Hirota system of coupled equations Ayşe Karasu (Kalkanli) and S Yu Sakovich Department

More information

ON POINCARE POLYNOMIALS OF HYPERBOLIC LIE ALGEBRAS

ON POINCARE POLYNOMIALS OF HYPERBOLIC LIE ALGEBRAS ON POINCARE POLYNOMIALS OF HYPERBOLIC LIE ALGEBRAS Meltem Gungormez Dept. Physics, Fac. Science, Istanbul Tech. Univ. 34469, Maslak, Istanbul, Turkey e-mail: gungorm@itu.edu.tr arxiv:0706.2563v2 [math-ph]

More information

Data representation and classification

Data representation and classification Representation and classification of data January 25, 2016 Outline Lecture 1: Data Representation 1 Lecture 1: Data Representation Data representation The phrase data representation usually refers to the

More information

The elliptic sinh-gordon equation in the half plane

The elliptic sinh-gordon equation in the half plane Available online at www.tjnsa.com J. Nonlinear Sci. Appl. 8 25), 63 73 Research Article The elliptic sinh-gordon equation in the half plane Guenbo Hwang Department of Mathematics, Daegu University, Gyeongsan

More information

Symmetry reductions and travelling wave solutions for a new integrable equation

Symmetry reductions and travelling wave solutions for a new integrable equation Symmetry reductions and travelling wave solutions for a new integrable equation MARIA LUZ GANDARIAS University of Cádiz Department of Mathematics PO.BOX 0, 50 Puerto Real, Cádiz SPAIN marialuz.gandarias@uca.es

More information

SMOOTHNESS PROPERTIES OF SOLUTIONS OF CAPUTO- TYPE FRACTIONAL DIFFERENTIAL EQUATIONS. Kai Diethelm. Abstract

SMOOTHNESS PROPERTIES OF SOLUTIONS OF CAPUTO- TYPE FRACTIONAL DIFFERENTIAL EQUATIONS. Kai Diethelm. Abstract SMOOTHNESS PROPERTIES OF SOLUTIONS OF CAPUTO- TYPE FRACTIONAL DIFFERENTIAL EQUATIONS Kai Diethelm Abstract Dedicated to Prof. Michele Caputo on the occasion of his 8th birthday We consider ordinary fractional

More information

Moving Boundary Problems for the Harry Dym Equation & Reciprocal Associates

Moving Boundary Problems for the Harry Dym Equation & Reciprocal Associates Moving Boundary Problems for the Harry Dym Equation & Reciprocal Associates Colin Rogers Australian Research Council Centre of Excellence for Mathematics & Statistics of Complex Systems & The University

More information

Automata on linear orderings

Automata on linear orderings Automata on linear orderings Véronique Bruyère Institut d Informatique Université de Mons-Hainaut Olivier Carton LIAFA Université Paris 7 September 25, 2006 Abstract We consider words indexed by linear

More information

WRONSKIANS AND LINEAR INDEPENDENCE... f (n 1)

WRONSKIANS AND LINEAR INDEPENDENCE... f (n 1) WRONSKIANS AND LINEAR INDEPENDENCE ALIN BOSTAN AND PHILIPPE DUMAS Abstract We give a new and simple proof of the fact that a finite family of analytic functions has a zero Wronskian only if it is linearly

More information

Rational Form Solitary Wave Solutions and Doubly Periodic Wave Solutions to (1+1)-Dimensional Dispersive Long Wave Equation

Rational Form Solitary Wave Solutions and Doubly Periodic Wave Solutions to (1+1)-Dimensional Dispersive Long Wave Equation Commun. Theor. Phys. (Beijing, China) 43 (005) pp. 975 98 c International Academic Publishers Vol. 43, No. 6, June 15, 005 Rational Form Solitary Wave Solutions and Doubly Periodic Wave Solutions to (1+1)-Dimensional

More information

Exact Solutions of Nonlinear Partial Dif ferential Equations by the Method of Group Foliation Reduction

Exact Solutions of Nonlinear Partial Dif ferential Equations by the Method of Group Foliation Reduction Symmetry, Integrability and Geometry: Methods and Applications Exact Solutions of Nonlinear Partial Dif ferential Equations by the Method of Group Foliation Reduction SIGMA 7 (2011), 066, 10 pages Stephen

More information

SYMBOLIC SOFTWARE FOR SOLITON THEORY: INTEGRABILITY, SYMMETRIES CONSERVATION LAWS AND EXACT SOLUTIONS. Willy Hereman

SYMBOLIC SOFTWARE FOR SOLITON THEORY: INTEGRABILITY, SYMMETRIES CONSERVATION LAWS AND EXACT SOLUTIONS. Willy Hereman . SYMBOLIC SOFTWARE FOR SOLITON THEORY: INTEGRABILITY, SYMMETRIES CONSERVATION LAWS AND EXACT SOLUTIONS Willy Hereman Dept. of Mathematical and Computer Sciences Colorado School of Mines Golden, Colorado

More information

Gauge Fixing and Constrained Dynamics in Numerical Relativity

Gauge Fixing and Constrained Dynamics in Numerical Relativity Gauge Fixing and Constrained Dynamics in Numerical Relativity Jon Allen The Dirac formalism for dealing with constraints in a canonical Hamiltonian formulation is reviewed. Gauge freedom is discussed and

More information

Symmetries and solutions of field equations of axion electrodynamics

Symmetries and solutions of field equations of axion electrodynamics Symmetries and solutions of field equations of axion electrodynamics Oksana Kuriksha Petro Mohyla Black Sea State University, 10, 68 Desantnukiv Street, 54003 Mukolaiv, UKRAINE Abstract The group classification

More information

LINEAR EQUATIONS WITH UNKNOWNS FROM A MULTIPLICATIVE GROUP IN A FUNCTION FIELD. To Professor Wolfgang Schmidt on his 75th birthday

LINEAR EQUATIONS WITH UNKNOWNS FROM A MULTIPLICATIVE GROUP IN A FUNCTION FIELD. To Professor Wolfgang Schmidt on his 75th birthday LINEAR EQUATIONS WITH UNKNOWNS FROM A MULTIPLICATIVE GROUP IN A FUNCTION FIELD JAN-HENDRIK EVERTSE AND UMBERTO ZANNIER To Professor Wolfgang Schmidt on his 75th birthday 1. Introduction Let K be a field

More information

Representations of Sp(6,R) and SU(3) carried by homogeneous polynomials

Representations of Sp(6,R) and SU(3) carried by homogeneous polynomials Representations of Sp(6,R) and SU(3) carried by homogeneous polynomials Govindan Rangarajan a) Department of Mathematics and Centre for Theoretical Studies, Indian Institute of Science, Bangalore 560 012,

More information

Conditional symmetries of the equations of mathematical physics

Conditional symmetries of the equations of mathematical physics W.I. Fushchych, Scientific Works 2003, Vol. 5, 9 16. Conditional symmetries of the equations of mathematical physics W.I. FUSHCHYCH We briefly present the results of research in conditional symmetries

More information

(Kac Moody) Chevalley groups and Lie algebras with built in structure constants Lecture 1. Lisa Carbone Rutgers University

(Kac Moody) Chevalley groups and Lie algebras with built in structure constants Lecture 1. Lisa Carbone Rutgers University (Kac Moody) Chevalley groups and Lie algebras with built in structure constants Lecture 1 Lisa Carbone Rutgers University Slides will be posted at: http://sites.math.rutgers.edu/ carbonel/ Video will be

More information

Classification of cubics up to affine transformations

Classification of cubics up to affine transformations Classification of cubics up to affine transformations Mehdi Nadjafikah and Ahmad-Reza Forough Abstract. Classification of cubics (that is, third order planar curves in the R 2 ) up to certain transformations

More information

-state problems and an application to the free particle

-state problems and an application to the free particle -state problems and an application to the free particle Sourendu Gupta TIFR, Mumbai, India Quantum Mechanics 1 2013 3 September, 2013 Outline 1 Outline 2 The Hilbert space 3 A free particle 4 Keywords

More information

JOINING OF CIRCUITS IN P SL(2, Z)-SPACE

JOINING OF CIRCUITS IN P SL(2, Z)-SPACE Bull. Korean Math. Soc. 5 (015), No. 6, pp. 047 069 http://dx.doi.org/10.4134/bkms.015.5.6.047 JOINING OF CIRCUITS IN P SL(, Z)-SPACE Qaiser Mushtaq and Abdul Razaq Abstract. The coset diagrams are composed

More information

Toward Analytic Solution of Nonlinear Differential Difference Equations via Extended Sensitivity Approach

Toward Analytic Solution of Nonlinear Differential Difference Equations via Extended Sensitivity Approach Commun. Theor. Phys. 57 (2012) 5 9 Vol. 57, No. 1, January 15, 2012 Toward Analytic Solution of Nonlinear Differential Difference Equations via Extended Sensitivity Approach G. Darmani, 1, S. Setayeshi,

More information

Exact solutions through symmetry reductions for a new integrable equation

Exact solutions through symmetry reductions for a new integrable equation Exact solutions through symmetry reductions for a new integrable equation MARIA LUZ GANDARIAS University of Cádiz Department of Mathematics PO.BOX, 1151 Puerto Real, Cádiz SPAIN marialuz.gandarias@uca.es

More information

arxiv:hep-th/ v1 16 Jul 1992

arxiv:hep-th/ v1 16 Jul 1992 IC-92-145 hep-th/9207058 Remarks on the Additional Symmetries and W-constraints in the Generalized KdV Hierarchy arxiv:hep-th/9207058v1 16 Jul 1992 Sudhakar Panda and Shibaji Roy International Centre for

More information

EXISTENCE OF NONTRIVIAL SOLUTIONS FOR A QUASILINEAR SCHRÖDINGER EQUATIONS WITH SIGN-CHANGING POTENTIAL

EXISTENCE OF NONTRIVIAL SOLUTIONS FOR A QUASILINEAR SCHRÖDINGER EQUATIONS WITH SIGN-CHANGING POTENTIAL Electronic Journal of Differential Equations, Vol. 2014 (2014), No. 05, pp. 1 8. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu EXISTENCE OF NONTRIVIAL

More information

Periodic and Solitary Wave Solutions of the Davey-Stewartson Equation

Periodic and Solitary Wave Solutions of the Davey-Stewartson Equation Applied Mathematics & Information Sciences 4(2) (2010), 253 260 An International Journal c 2010 Dixie W Publishing Corporation, U. S. A. Periodic and Solitary Wave Solutions of the Davey-Stewartson Equation

More information

BLOCKS IN THE CATEGORY OF FINITE-DIMENSIONAL REPRESENTATIONS OF gl(m n)

BLOCKS IN THE CATEGORY OF FINITE-DIMENSIONAL REPRESENTATIONS OF gl(m n) BLOCKS IN THE CATEGORY OF FINITE-DIMENSIONAL REPRESENTATIONS OF gl(m n) VERA SERGANOVA Abstract. We decompose the category of finite-dimensional gl (m n)- modules into the direct sum of blocks, show that

More information

Gravitational Waves: Just Plane Symmetry

Gravitational Waves: Just Plane Symmetry Utah State University DigitalCommons@USU All Physics Faculty Publications Physics 2006 Gravitational Waves: Just Plane Symmetry Charles G. Torre Utah State University Follow this and additional works at:

More information

New approach for tanh and extended-tanh methods with applications on Hirota-Satsuma equations

New approach for tanh and extended-tanh methods with applications on Hirota-Satsuma equations Volume 28, N. 1, pp. 1 14, 2009 Copyright 2009 SBMAC ISSN 0101-8205 www.scielo.br/cam New approach for tanh and extended-tanh methods with applications on Hirota-Satsuma equations HASSAN A. ZEDAN Mathematics

More information

Symmetry Reductions of (2+1) dimensional Equal Width. Wave Equation

Symmetry Reductions of (2+1) dimensional Equal Width. Wave Equation Authors: Symmetry Reductions of (2+1) dimensional Equal Width 1. Dr. S. Padmasekaran Wave Equation Asst. Professor, Department of Mathematics Periyar University, Salem 2. M.G. RANI Periyar University,

More information

arxiv:math/ v1 [math.co] 6 Dec 2005

arxiv:math/ v1 [math.co] 6 Dec 2005 arxiv:math/05111v1 [math.co] Dec 005 Unimodality and convexity of f-vectors of polytopes Axel Werner December, 005 Abstract We consider unimodality and related properties of f-vectors of polytopes in various

More information

A multiple Riccati equations rational expansion method and novel solutions of the Broer Kaup Kupershmidt system

A multiple Riccati equations rational expansion method and novel solutions of the Broer Kaup Kupershmidt system Chaos, Solitons and Fractals 30 (006) 197 03 www.elsevier.com/locate/chaos A multiple Riccati equations rational expansion method and novel solutions of the Broer Kaup Kupershmidt system Qi Wang a,c, *,

More information

A determinant characterization of moment sequences with finitely many mass-points

A determinant characterization of moment sequences with finitely many mass-points To appear in Linear and Multilinear Algebra Vol 00, No 00, Month 20XX, 1 9 A determinant characterization of moment sequences with finitely many mass-points Christian Berg a and Ryszard Szwarc b a Department

More information

arxiv: v1 [math-ph] 4 May 2016

arxiv: v1 [math-ph] 4 May 2016 arxiv:1605.01173v1 [math-ph] 4 May 2016 On the Classifications of Scalar Evolution Equations with Non-constant Separant Ayşe Hümeyra BİLGE Faculty of Engineering and Natural Sciences, Kadir Has University,

More information

New explicit solitary wave solutions for (2 + 1)-dimensional Boussinesq equation and (3 + 1)-dimensional KP equation

New explicit solitary wave solutions for (2 + 1)-dimensional Boussinesq equation and (3 + 1)-dimensional KP equation Physics Letters A 07 (00) 107 11 www.elsevier.com/locate/pla New explicit solitary wave solutions for ( + 1)-dimensional Boussinesq equation and ( + 1)-dimensional KP equation Yong Chen, Zhenya Yan, Honging

More information

Conservation laws for difference equations. Peter Hydon (Surrey) Elizabeth Mansfield (Kent) Olexandr Rasin (Bar-Ilan) Tim Grant (Surrey)

Conservation laws for difference equations. Peter Hydon (Surrey) Elizabeth Mansfield (Kent) Olexandr Rasin (Bar-Ilan) Tim Grant (Surrey) Conservation laws for difference equations Peter Hydon (Surrey) Elizabeth Mansfield (Kent) Olexandr Rasin (Bar-Ilan) Tim Grant (Surrey) Introduction: What is a conservation law? What is a conservation

More information

A supersymmetric Sawada-Kotera equation

A supersymmetric Sawada-Kotera equation A supersymmetric Sawada-Kotera equation arxiv:0802.4011v2 [nlin.si] 7 Dec 2008 Kai Tian and Q. P. Liu Department of Mathematics, China University of Mining and Technology, Beijing 100083, P.R. China Abstract

More information

THE SEMISIMPLE SUBALGEBRAS OF EXCEPTIONAL LIE ALGEBRAS

THE SEMISIMPLE SUBALGEBRAS OF EXCEPTIONAL LIE ALGEBRAS Trudy Moskov. Matem. Obw. Trans. Moscow Math. Soc. Tom 67 (2006) 2006, Pages 225 259 S 0077-1554(06)00156-7 Article electronically published on December 27, 2006 THE SEMISIMPLE SUBALGEBRAS OF EXCEPTIONAL

More information

Vectors. January 13, 2013

Vectors. January 13, 2013 Vectors January 13, 2013 The simplest tensors are scalars, which are the measurable quantities of a theory, left invariant by symmetry transformations. By far the most common non-scalars are the vectors,

More information

Exact Solutions to the Focusing Discrete Nonlinear Schrödinger Equation

Exact Solutions to the Focusing Discrete Nonlinear Schrödinger Equation Exact Solutions to the Focusing Discrete Nonlinear Schrödinger Equation Francesco Demontis (based on a joint work with C. van der Mee) Università degli Studi di Cagliari Dipartimento di Matematica e Informatica

More information

NILPOTENCY INDEX OF NIL-ALGEBRA OF NIL-INDEX 3

NILPOTENCY INDEX OF NIL-ALGEBRA OF NIL-INDEX 3 J. Appl. Math. & Computing Vol. 20(2006), No. 1-2, pp. 569-573 NILPOTENCY INDEX OF NIL-ALGEBRA OF NIL-INDEX 3 WOO LEE Abstract. Nagata and Higman proved that any nil-algebra of finite nilindex is nilpotent

More information

On the Linearization of Second-Order Dif ferential and Dif ference Equations

On the Linearization of Second-Order Dif ferential and Dif ference Equations Symmetry, Integrability and Geometry: Methods and Applications Vol. (006), Paper 065, 15 pages On the Linearization of Second-Order Dif ferential and Dif ference Equations Vladimir DORODNITSYN Keldysh

More information

Degenerate Perturbation Theory. 1 General framework and strategy

Degenerate Perturbation Theory. 1 General framework and strategy Physics G6037 Professor Christ 12/22/2015 Degenerate Perturbation Theory The treatment of degenerate perturbation theory presented in class is written out here in detail. The appendix presents the underlying

More information