arxiv: v1 [math-ph] 9 Dec 2009

Transcription

1 Darboux transformation of the generalized coupled dispersionless integrable system arxiv:092.67v [math-ph] 9 Dec 2009 M. Hassan Department of Mathematics University of Glasgow Glasgow G2 8QW Scotland and Department of Physics University of the Punjab Quaid-e-Azam Campus Lahore Pakistan. Abstract The Darboux transformation on matrix solutions to the generalized coupled dispersionless integrable system based on some non-abelian Lie group is studied and the solutions are shown to be expressed in terms of quasideterminants. As an explicit example the Darboux transformation on scalar solutions to the system based on the Lie group SU(2) is discussed in detail and the solutions are shown to be expressed as ratios of determinants. PACS: 02.30Zz Ik Keywords: Integrable systems dispersionless hierarchies Darboux transformation mhassan@physics.pu.edu.pk mhassan@maths.gla.ac.uk

2 Introduction There has been a considerable interest in dispersionless integrable hierarchies for the last couple of decades partly because of their emergence in diverse areas of mathematical and theoretical physics such as quantum field theory conformal field theory string theory soliton theory etc. (see e.g. []-[7]). These integrable systems arise as semi-classical limit of the ordinary integrable systems. In other words the quasi-classical limit corresponds to the solutions which slowly depend on the independent variables resulting in the elimination of dispersion term of the original integrable equation. In recent past the coupled dispersionless integrable system and its generalization has also attracted a great deal of interest because of its nice integrability structure and its soliton solutions [8]-[5]. However it should be mentioned here that the coupled dispersionless integrable system introduced in Refs. [8]-[3] is referred to asdispersionless in the sense that it does not contain dispersion term andnot in the sense of quasi-classical limit of some ordinary integrable system where the dispersion term is eliminated by taking an appropriate limit. The coupled dispersionless integrable system and its generalization is solvable by the inverse scattering method possesses infinite number of conservation laws and has Painleve property [8]-[5]. For the case of Lie group SU(2) the multi-soliton solutions have been investigated using Backlund and Darboux transformations [6]-[7]. In this paper we study the Darboux transformation of the generalized coupled dispersionless integrable system based on some Lie group G and obtain its solutions expressed in terms of quasideterminants of Gelfand and Retakh [8]-[20]. For the special case of G = SU(2) we show that the scalar solutions to the system are expressed as ratios of determinants showing the transition from the matrix case to the scalar case. We compare our results with those of [7] and also with the previously known results of the sine-gordon system. For the construction of multi-solitons we use the method of Darboux matrix and by iteration of Darboux transformation we obtain the quasideterminant solutions to the generalized coupled dispersionless integrable system. The layout of the paper is as follows. In section 2 we give a review of the generalized coupled dispersionless integrable system by writing its lagrangian lagrange equations of motion and the Lax pair. In section 2 we also introduce the notion of quasideterminants which shall be used in section 3 to express the solutions. In section 3 we define a Darboux transformation via a Darboux matrix on matrix solutions to the generalized coupled dispersionless

3 2 integrable system and express the solutions to the system in a closed form as quasideterminants. In section 4 we take an explicit example of the system when the underlying Lie group is SU(2) and write its multisoliton solutions as ratios of determinants by using the Darboux transformation already discussed in section 3. At the end we compare our results with the already known results and make concluding remarks in the last section. 2 The generalized coupled dispersionless integrable system The generalized coupled dispersionless integrable system based on some non-abelian Lie group G is a classical lagrangian field theory defined by the following action functional [2] I = dtdxl(s x S t S) (2.) with the lagrangian density defined by ( L = Tr 2 xs t S ) 3 G[S[ xss]] (2.2) where S is some matrix field and G is a constant matrix; taking values in non-abelian Lie algebra g of the Lie group G. Let us introduce a basis of anti-hermitian generators T a for g with real totally antisymmetric structure constants f abc and normalization given by [T a T b ] = f abc T c Tr(T a T b ) = δ ab. (2.3) For any X g we write X = T a X a X a = Tr(T a X). (2.4) The matrix fields S and G of the generalized coupled dispersionless system (2.2 ) take values in the Lie algebra g therefore we write S = φ a T a G = κ a T a where φ = φ(xt) is a vector field with components {φ a a = 2 dimg} and κ is a constant vector with components {κ a a = 2 dimg}. The lagrange equation of motion resulting from the action is given by t x S [[SG] x S] = 0. (2.5)

4 3 The field equation (2.5) is known as the generalized coupled dispersionless integrable system based on a Lie group G. The equation (2.5) has an n n AKS (Ablowitz-Kaup-ewell- Segur) representation and the nonlinearity of the equation is due to the non-abelian character of the underlying Lie group. The equation (2.5) is a generalization of the coupled dispersionless equation originally introduced in [8] for G = SL(2R) and SU(2). For G = SL(2R) the equation (2.5) reduces to t x q + x (rs) = 0 t x r 2 x qr = 0 t x s 2 x qs = 0 where qr and s are real valued functions of x and t. For G = SU(2) the system (2.5) is equivalent to t x q +r x r = 0 t x r 2 x qr = 0 t x r 2 x q r = 0 (2.6) where r is a complex valued function and r denotes complex conjugation of r. The system (2.6) is in fact equivalent to Pohlmeyer-Lund-Regge system [0]. If r is real then the system (2.6) is the one originally proposed in [8] and is equivalent to the sine-gordon equation [9]. With reference to Darboux transformation we shall discuss it in detail in section 4. The generalized coupled dispersionless integrable system (2.5) can be expressed as the compatibility condition of the following Lax pair [2] x ψ = U(xtλ)ψ (2.7) t ψ = V(xtλ)ψ (2.8) where ψ G λ is the spectral parameter and the matrix fields U(xtλ) and V(xtλ) are given by U(xtλ) = λ x S V(xtλ) = [SG]+λ G. The compatibility condition of (2.7)-(2.8) is the zero-curvature condition for the matrices U(xtλ) and V(xtλ) t U(xtλ) x V(xtλ)+[U(xtλ)V(xtλ)] = 0 (2.9)

5 4 which is equivalent to the equation of motion (2.5). In the next section we define the Darboux transformation via a Darboux matrix on matrix solutions ψ of the Lax pair (2.7)- (2.8). To write down the explicit expressions for matrix solutions to the generalized coupled dispersionless integrable system we will use the notion of quasideterminant introduced by Gelfand and Retakh [8]-[20]. LetX beann nmatrixover aringr(noncommutative ingeneral). Forany i j n let r i be the ith row and c j be the jth column of X. The matrix X has n 2 quasideterminants denoted by X ij for ij =...n and are defined by X ij = x ij r j i ( X ij ) c i j (2.0) where x ij is the ijth entry of X r j i represents the ith row of X without the jth entry cj i represents the jth column of X without the ith entry and X ij is the submatrix of X obtained by removing from X the ith row and the jth column. The quasideterminats are also denoted by the following notation X ij = X ij If the entries of the matrix X all commute then r j i cj i x ij. (2.) i+j detx X ij = ( ) detxij. (2.2) For a detailed account of quasideterminants and their properties see e.g. [8]-[20]. In this paper we will consider only quasideterminants that are expanded about an n n matrix over a commutative ring. In our case the ring R is the (noncommutative) ring of n n matrices over another commutative ring. Let ( A B C D be a block decomposition of any matrix where the matrix D is n n and A is invertible. The quasideterminant expanded about the matrix D is defined by A B C D = D CA B. (2.3) )

6 5 3 Darboux transformation The Darboux transformation on the matrix solutions of the Lax pair (2.7)-(2.8) is defined in termsofann nmatrixd(xtλ) calledthedarbouxmatrix. TheDarbouxmatrixconnects the two matrix solutions of the Lax pair (2.7)-(2.8) such that the Lax pair is covariant under the Darboux transformation (see e.g. [22]-[28]). Let us denote the new matrix solution to the Lax pair (2.7)-(2.8) by ψ[] so that the Darboux transformation is defined by ψ[] = D(xtλ)ψ. (3.) For the present case we make the following ansatz for the Darboux matrix D(xtλ) = λ I M(xt) (3.2) where M(xt) is some n n matrix field to be determined and I is an n n identity matrix. For the covariance of the Lax pair we require that x ψ[] = λ x S[]ψ[] (3.3) t ψ[] = [S[]G[]]ψ[]+λ G[]ψ[] (3.4) where S[] and G[] are the Darboux transformed matrix fields defined by S[] = S M G[] = G (constant matrix) (3.5) such that the matrix M is required to satisfy the following conditions x MM = [ x SM] (3.6) t M = [[SG]M]+[GM]M. (3.7) The next step is to determine the matrix M in terms of a particular solution of the Lax pair such the M satisfies the conditions (3.6)-(3.7). For this pupose we proceed as follows. Let λ λ n be non-zero distinct real or complex constant parameters and e e n be constant column vectors such that Θ = (ψ(λ )e ψ(λ n )e n ) = (θ θ n ) (3.8)

7 6 be an invertible n n matrix. Each column θ i = ψ(λ i )e i of the matrix Θ is a column solution of the Lax pair (2.7)-(2.8) at λ = λ i i.e. the columns θ i i = n satisfy x θ i = λ i x Sθ i (3.9) t θ i = [SG]θ i +λ i Gθ i. (3.0) Let us define an invertible constant diagonal matrix with entries being the eigenvalues λ i corresponding to the eigenvectors θ i Λ = diag(λ λ n ). (3.) The Lax pair (2.7)-(2.8) can now be written in the matrix form as x Θ = x SΘΛ (3.2) t Θ = [SG]Θ+GΘΛ (3.3) where Θ is a particular matrix solution of the Lax pair with a matrix Λ of particular eigenvalues. ow we show that the matrix M = ΘΛ Θ expressed in terms of particular matrix solution Θ of the Lax pair (2.7)-(2.8) satisfies the conditions (3.6)-(3.7) on the matrix M imposed by the covariance of the Lax pair under the Darboux transformation(3.) and(3.5). For this purpose let us take the x derivative of the matrix M = ΘΛ Θ x M = x ΘΛ Θ +ΘΛ x Θ = x S ΘΛ Θ x SΘΛΘ = x S M x SM (3.4) which is equation (3.6). ow take the t derivative of M = ΘΛ Θ to get t M = t ΘΛ Θ +ΘΛ t Θ = [SG]ΘΛ Θ +GΘΛ 2 Θ ΘΛ Θ [SG] ΘΛ Θ GΘΛ Θ = [SG]M +GM 2 M[SG] MGM = [[SG]M]+[GM]M (3.5) which is equation (3.7). This shows that the matrix M = ΘΛ Θ satisfies the conditions (3.6)-(3.7). So we have established that the transformation ψ[] = D(xtλ)ψ = (λ I ΘΛ Θ )ψ (3.6)

8 7 constitutes the Darboux transformation on the matrix solution ψ of the Lax pair (2.7)- (2.8) of the generalized coupled dispersionless integrable system (2.5). The corresponding Darboux transformations on the matrix fields S and G are S[] = S ΘΛ Θ G[] = G. (3.7) By introducing the notation ψ () = λ ψ and Θ () = ΘΛ the one-fold Darboux transformations (3.6)-(3.7) can also be expressed in terms of quasideterminants as Θ ψ ψ[] = Θ () ψ () (3.8) Θ I S[] = S + O (3.9) where O is an n n null matrix. The iteration of Darboux transformation times gives the quasideterminant matrix solution to the generalized coupled dispersionless integrable system (2.5). For each k = 2 let Θ k be an invertible matrix solution of the Lax pair (2.7)-(2.8) at Λ = Λ k. ow using the notation Θ (k) = ΘΛ k Θ[0] = Θ ψ[0] = ψ and for we write Θ () ψ[] = ψ[ ] Θ () [ ](Θ[ ]) ψ () [ ] Θ Θ ψ = Θ ( ) Θ ( ) ψ ( ). (3.20) Θ () Θ () ψ () Similarly the expression for S[] is S[] = S = S + j= Θ[j] () (Θ[j]) Θ Θ O Θ ( 2) Θ ( 2) O Θ ( ) Θ ( ) I Θ () Θ () O. (3.2) These results can be proved by induction and the proof is identical to the one for noncommutative KP equation (see e.g. [20]-[2] ).

9 8 4 The SU(2) system In order to get explicit solutions and compare our results with the already known soliton solutions of the coupled dispersionless integrable system we proceed as follows. Let us make general remarks about the case when G = SU(n) and then focus on the specific case of G = SU(2). The Lie group G = SU(n) consists of unimodular anti-hermitian n n matrices. Since the matrix fields S and G are valued in the Lie algebra su(n) of the Lie group SU(n) therefore S = S G = G TrS = 0 TrG = 0. (4.) In order to have S[] and G[] to be valued in su(n) we need to have M = M and TrM = 0. For the particular solutions θ i at λ = λ i let us compute x (θ i θ j) = λ i θ i xs θ j +λ j θ i xsθ j x (θ i θ ( j) = θ i [SG] +[SG] ) θ j + λ i θ i G θ j +λ j θ i Gθ j. (4.2) Since S and G are anti-hermitian therefore we get ) x (θ i θ j = t θ i θ j = 0 (4.3) when λ i λ j (i.e. λ i = λ j ). Also from the definition of the the matrix M we can calculate ( θ i M +M ) θ j = ( λ ) i +λ j θ i θ j (4.4) which implies θ i θ j = 0 (4.5) when λ i λ j. Given this weconclude that thecolumnvectors θ i arealllinearly independent and the condition (4.5) holds at every point of space and time. For the system with G = SU(2) the corresponding equation of motion consists of the hierarchy of localized induction equation of a thin vortex filament which is equivalent to the Heisenberg spin equation and the hierarchy of dispersionless equations which are equivalent to Pohlmeyer-Lund-Regge system and the sine-gordon equation. Here we shall consider the coupled dispersionless integrable system which is equivalent to the sine-gordon equation [8]. Let us introduce a vector φ = (φ φ 2 φ 3 ) in such a way that the matrix field S is defined by φ 3 φ iφ 2 S = i (4.6) φ +iφ 2 φ 3

10 9 which is traceless and anti-hermitian. The matrices U and V are then given by ( ) x φ 3 x φ i x φ 2 U = iλ x φ +i x φ 2 x φ 3 0 φ iφ 2 V = i 0. (4.7) φ iφ 2 0 2λ 0 By writing φ = r φ 2 = 0 φ 3 = q we get the coupled dispersionless integrable system of ref. [8] x t q +2 x rr = 0 x t r 2 x qr = 0. (4.8) The soliton solutions to this system have been computed in [7] by using Backlund and Darboux transformations. In ref. [7] the Darboux matrix D(λ) has been obtained using the method described in the section 3 and is given by ( ) λ λ cosω λ sinω D(λ) = (4.9) λ sinω λ +λ cosω where tan ω 2 = β α and (αβ)t is a particular column solution of the Lax pair at λ = λ. The corresponding matrix M is M = λ ( cosω sinω sinω cosω ) (4.0) which is traceless and is anti-hermitian when λ = λ. The one-fold Darboux transformation (3.7) gives the Darboux transformation on the fields q and r q[] = q +iλ cos ω r[] = r +iλ sin ω. (4.) By assuming the seed solutions to be x q = and r = 0 the one-soliton solutions turn out to be x q[] = +2sech 2 ( 2iλ x i r[] = i ( sech 2iλ x i t λ λ ) t λ ) (4.2) where x q[] is called a soliton of dark-type and r[] is called a soliton of bright-type [8].

11 0 To get simpler expressions of Darboux transformation on scalar solutions (eigenfunctions) of the Lax pair and the scalar fields q and r we make a gauge transformation on matrix fields S U and V. Let us introduce a matrix Ω = (4.3) 2 i i so that the gauge equivalent matrix field S is given by S = Ω 0 q +ir SΩ = i. (4.4) q ir 0 The Lax pair is then expressed as X 0 x (q +ir) X x = iλ Y x (q ir) 0 Y ( X ir i ) X t = 2λ Y i. (4.5) ir Y 2λ Let X and Y be the particular scalar solutions of the Lax pair (4.5) at λ = λ then the one-fold Darboux matrix for the Lax pair (4.5) is given by λ λ X Y D(λ) = λ Y X λ. (4.6) In other words the Lax pair (4.5) is covariant under the following Darboux transformations X[] = X () X() Y = (X Y XY)[2] Y 2 (X Y )[] Y[] = Y () Y () q[] = q i 2 r[] = r 2 X () = 2(X Y XY)[2] X (X Y )[] X () + Y () = q + t log( (X Y )[] 2 (X Y )[]) Y X X () Y () = r +i t log Y X (X Y )[] (4.7) 2 (X Y )[] where we have used the notation X () = λ X and the determinants and 2 are defined by Y Y (X Y XY)[2] = X () X () X X 2(X Y XY)[2] = Y () Y () (X Y )[] = X 2 (X Y )[] = Y.

12 Similarly the two-fold Darboux transformations on scalar solutions X Y of the lax pair (4.5) and q r of the coupled dispersionless integrable system (4.8) are given by X[2] = X () [] X() [] Y [] Y[] = (X k Y k XY)[3] 2 (X k Y k )[2] Y[2] = Y () Y () X () = 2(X Y XY)[3] X (X Y )[2] q[2] = q[] i X () [] + Y () [] = q + t log( (X k Y k )[2] 2 (X k Y k )[2]) 2 Y [] X [] r[2] = r[] X () [] Y () [] (X k Y k )[2] = r +i t log (4.8) 2 Y [] X [] 2 (X k Y k )[2] where for each k = 2 we have assumed X k Y k to be particular solutions of the Lax pair (4.5) at λ = λ k and the determinants and 2 are defined by X X 2 X (X k Y k XY)[3] = Y () Y () 2 Y () X (2) X (2) 2 X (2) Y Y 2 Y 2 (X k Y k XY)[3] = X () X () 2 X () Y (2) Y (2) 2 Y (2) Y Y 2 (X k Y k )[2] = X () X () 2 X X 2 2 (X k Y k )[2] = Y () Y (). (4.9) 2 We can iterate the Darboux transformation times to arrive at the -soliton solutions to the coupled dispersionless integrable system (4.8) expressed in closed form as ratios of determinants. For each k = let X k and Y k be particular scalar solutions of the Lax pair at λ = λ k respectively then by using the notation X (k) j = λ k j X j and X (k) = λ k X we get the following result X[] = (X k Y k XY)[ +] 2 (X k Y k )[] Y[] = 2(X k Y k XY)[ +] (4.20) (X k Y k )[] where the determinants are defined by the following expressions. For odd the determinants

13 2 and 2 are given by (X k Y k XY)[ +] = 2 (X k Y k XY)[ +] = Y Y Y Y ( ) Y ( ) X () X () Y ( ) X () X X X X ( ) X ( ) Y () Y () X ( ) For even we have the following determinants X X X. (X k Y k XY)[ +] =..... Y ( ) Y ( ) X () X () Y Y Y. 2 (X k Y k XY)[ +] =..... X ( ) X ( ) Y () Y () Y () Y ( ) X () X ( ) Y (). (4.2). (4.22) Similarly after performing Darboux transformations times we get the following expressions of the -soliton solutions to the couple dispersionless integrable system (4.8) q[] = q + t log( (X k Y k )[] 2 (X k Y k )[]) r[] = ( ) (X k Y k )[] r +i t log. (4.23) 2 (X k Y k )[] Since the system (4.8) is equivalent to the sine-gordon system therefore the solutions to both the system are related with each other. The natural boundary conditions for the system (4.8) are x q (const.) r 0 x. With these boundary conditions the solution ϕ to the sine-gordon equation x t ϕ = 2sinϕ and the solutions q and r to the system (4.8) are related by x q = cosϕ r = ± 2 tϕ.

14 3 The -fold Darboux transformation on ϕ is given by (X k Y k )[] ϕ[] = ϕ+2ilog 2 (X k Y k )[] 2 (X k Y k )[] exp( iϕ[]) = exp( iϕ) (4.24) 2 (X k Y k )[] which are the well-known expressions of -soliton solutions to the sine-gordon equation [23]. 5 Concluding remarks In this paper we have studied the Darboux transformation on matrix solutions to the generalized coupled dispersionless integrable system based on some non-abelian Lie group and expressed the matrix solutions in terms of quasideterminants. For the particular case of the coupled dispersionless system which is equivalent to the sine-gordon equation we have defined a Darboux transformation on scalar solutions of the linear problem and expressed them as ratios of determinants as is the case with the sine-gordon equation. There are a number of directions where the present work could be extended and the immediate one is to study the effect of Moyal noncommutativity of space coordinates on the soliton solutions of the system. The other interesting direction to pusue is to study the supersymmetric generalization of the generalized coupled dispersionless integrable system. We shall address these and related problems in some later work. Acknowledgements I thank Jonathan immo for helpful discussions and the Department of Mathematics University of Glasgow for warm hospitality and provision of facilities while this work was done. I gratefully acknowledge financial support of the Higher Education Commission of Pakistan. References [] I. Krichever The dispersionless Lax equations and topological minimal models Commun. Math. Phys. 43 (992) 45.

15 4 [2] S. Aoyama and Y. Kodama Topological conformal field theory with a rational W potential and the dispersionless KP hierarchy Mod. Phys. Lett. A 9(994) 248. [3] R. Carroll and Y. Kodama Solution of the dispersionless Hirota equations J. Phys. A 28 (995) 6373 [arxiv:hep-th/ ]. [4] K. Takasaki and T. Takebe Integrable hierarchies and dispersionless limit Rev. Math. Phys. 7 (995) 743 [arxiv:hep-th/ ]. [5] K. Takasaki Dispersionless Toda hierarchy and two-dimensional string theory Commun. Math. Phys. 70 (995) 0 [arxiv:hep-th/940390]. [6] P. B. Wiegmann and A. Zabrodin Conformal maps and dispersionless integrable hierarchies Commun. Math. Phys. 23 (2000) 523 [arxiv:hep-th/990947]. [7] M. Dunajski Interpolating dispersionless integrable system J. Phys. A 4 (2008) [arxiv: [nlin.si]]. [8] K. Konno and H. Oono ew coupled integrable dispersionless equations J. Phys. Soc. Jpn. 63 (994) 377. [9] R. Hirota and S. Tsujimoto ote on ew coupled integrable dispersionless equations J. Phys. Soc. Jpn. 63 (994) [0] V. P. Kotlyarov On equations gauge equivalent to the sine-gordon and Pohlmeyer- Lund-Regge equations J. Phys. Soc. Jpn. 63 (994) [] H. Kakuhata and K. Konno Lagrangian hamiltonian and conserved quantities for coupled integrable dispersionless equations J. Phys. Soc. Jpn. 65 (995). [2] H. Kakuhata and K. Konno A generalization of coupled integrable dispersionless system J. Phys. Soc. Jpn. 65 (996) 340. [3] H. Kakuhata and K. Konno Canonical formulation of generalized coupled dispersionless system J. Phys. A 30 (997) L40. [4] T. Alagesan and K. Porsezian Painleve analysis and the integrability properties of coupled integrable dispersionless equations Chaos Solitons and Fractals 7 (996) 209.

16 5 [5] T. Alagesan and K. Porsezian Singularity structure analysis and Hirota s bilinearisation of the coupled integrable dispersionless equations Chaos Solitons and Fractals 8 (997) 645. [6] T. Alagesan Y. Chung and K. akkeeran Backlund transformation and soliton solutions for the coupled dispersionless equations Chaos Solitons and Fractals 2 (2004) 63. [7] A. Chen and X. Li Soliton solutions of the coupled dispersionless equation Phys. Lett. A 370 (2007) 28. [8] I. Gelfand and V. Retakh Determinants of matrices over noncommutative rings Funct. Anal. Appl. 25 (99) no [9] I. Gelfand and V. Retakh A theory of noncommutative determinants and characteristic functions of graphs Funct. Anal. Appl. 26 (992) no [20] I. Gelfand S. Gelfand V. Retakh and R. L. Wilson Quasideterminants Adv. Math. 93 (2005) 56. [2] C. R. Gilson and J. J. C. immo On a direct approach to quasideterminant solutions of a noncommutative KP equation J. Phys. A: Math. Theor. 40 (2007) [22] D. Levi O. Ragnisco A. Sym Dressing method vs. classical Darboux transformation Il uovo Cim. 83B() (984) 34. [23] V. B. Matveev and M. A. Salle Darboux Transformations and Soliton Springer-Verlag (99). [24] A. L. Sakhnovich Dressing procedure for solutions of non-linear equations and the method of operator indentities Inv. Problems 0 (994) 699. [25] J. Cieslihski An algebraic method to construct the Darboux matrix J. Math. Phys. 36 (995) [26] M. Manas Darboux transformations for the nonlinear Schrodinger equations J. Phys. A: Math. Gen. 29 (996) 772. [27] C. Rogers and W. K. Schief Backlund and Darboux Transformations: Geometry and Modern Applications in Soliton Theory Cambridge (2003).

17 6 [28] C. Gu H. Hu and Z. Zhou Darboux transformations in integrable systems theory and their applications to geometry Springer (2005).