arxiv:nlin/ v1 [nlin.si] 7 Sep 2005

Transcription

1 NONSINGULAR POSITON AND COMPLEXITON SOLUTIONS FOR THE COUPLED KDV SYSTEM arxiv:nlin/5918v1 [nlin.si] 7 Sep 25 H. C. HU 1,2, BIN TONG 1 AND S. Y. LOU 1,3 1 Department of Physics, Shanghai Jiao Tong University, Shanghai, 23, P. R. China 2 College of Science, University of Shanghai for Science and Technology, Shanghai, 293, P. R. China 3 Department of Physics, Ningbo University, Ningbo, , P. R. China Abstract. Taking the coupled KdV system as a simple eample, analytical and nonsingular compleiton solutions are firstly discovered in this letter for integrable systems. Additionally, the analytical and nonsingular positon-negaton interaction solutions are also firstly found for S-integrable model. The new analytical positon, negaton and compleiton solutions of the coupled KdV system are given out both analytically and graphically by means of the iterative Darbou transformations. PACS.2.3.Ik, 2.3.Jr, 5.45.Yv. 1. Introduction On the eact solutions of integrable models, there is a new classification way recently based on the property of the spectral parameter [1]: If the spectral parameter is positive, the related solution of the nonlinear system is called positon which is usually epressed by means of the triangle functions. If the spectral parameter is negtive, the related solution of the nonlinear system is called negaton which is usually epressed with help of the hyperbolic functions. The term positon and negaton may be pursued back to Matveev et al [2, 3] though the singular solutions had been studied earlier [4]. The so-called compleiton, which is epressed by combinations of the triangle functions and the hyperbolic functions, is related to the comple spectral parameters. It is meaningful that the compleiton eists for real integrable systems. The most of known positon solutions are singular. For various important integrable systems such as the KdV equation there is no nonsingular positon found. The negaton can be both singular and nonsingular. Especially nonsingular negatons are just the usual solitons. The interaction solutions among positons and negatons have also been 1

2 studied recently[5, 6]. Unfortunately, ecept for a special C-integrable model, the STO (Sharma-Tass-Olver) model [5], in our knowledge, there is no known nonsingular positonnegaton interaction solution. Though the compleiton is proposed for several years, but up to now, one has not yet found a nonsingular compleiton solution for any (1+1)-dimensional integrable systems. In high dimensions, (2+1)- and (3+1)-dimensions, some special types of analytic nonsingular positon, posilton-negaton interaction solutions and the compleiton solutions can be easily obtained because of the eistence of the arbitrary functions in their epressions of eact solutions[7, 8]. Then an important and interesting question is naturally arised: Are there any nonsingular compleiton and positon-negaton solutions for S-integrable models? As we all know, the Korteweg-de Vries (KdV) equation is a typical and classical model to describe the weakly shallow long waves in the soliton theory. Many important properties of KdV equation, such as infinitely many symmetries, infinitely many conservation laws [9], inverse scattering transform, Bäcklund transformation, Darbou transformation, etc. have been studied etensively for a long time and it can be viewed as a completely integrable Hamiltonian system[1, 11]. Nevertheless, it is known that both the positon and compleiton solutions of the KdV equation are singuler[2, 1]. On the other hand, coupled integrable systems, which come up in many physical fields, have also been studied etensively and a lot of interesting results are given out from both the classification view and application field [12]. The first coupled KdV system was put forward by Hirota and Satsuma (HS) in 1981 [13]. The multi-soliton solutions, conservation laws and Darbou transformation, etc for the HS model are studied in detail [14]. Many other coupled KdV systems such as the Ito s system [15], the Drinfeld and Sokolov [16] model, the Fuchssteiner equation [17], the Nutk-O guz model[18], the Zharkov system [19] and the Foursov model [2] are also constructed. Vladimir gives a classification of the integrable coupled KdV systems by using a symmetry approach [12]. 2

3 In the real application aspect of the coupled KdV systems, a quite general coupled KdV equation, u t 6uu +u +ǫ 1 vv +ǫ 2 (uv) +ǫ 3 v =, v t 6βvv +βv V v +α[ǫ 1 (vu) +ǫ 2 uu +ǫ 3 u ] = (1) where ǫ 1, ǫ 2, ǫ 3, β and V are arbitrary constants, is derived from a two-wave modes in a shallow stratified liquid [21]. Another type of general coupled KdV system u t +α 1 vu +(α 2 v 2 +α 3 uv +α 4 u +α 5 u 2 ) =, v t +δ 1 vu +(δ 2 u 2 +δ 3 uv +δ 4 v +δ 5 v 2 ) =, (2) with arbitrary constants α i, δ i, i = 1, 2,..., 5 is obtained from a two layer model of the atmospheric dynamical system [22]. A more general coupled KdV system τ q j + ζ γj kl q k q l + ζ 3 βjq k k =, j = 1, 2 (3) k,l is derived to describe two-component Bose-Einstein condensates [23]. k In this Letter, we study the eact solutions, especially the positons, negatons and compleitons of a special coupled KdV system u t +6vv 6uu +u =, v t 6uv 6vu +v =. (4) It is clear that the coupled KdV system (4) is a special case for all the physical systems (1), (2) and (3). It is also interesting that, the system (4) can be read off simply from the real and imaginary parts of the single comple KdV equation where U t 6UU +U =, (5) U = u+iv, I 1. From the fact that (4) can be derived from the comple KdV equation, one can reasonable believe that all the integrable properties of the real KdV equation can be preserved. Actually, 3

4 we have really proved its Painlevé property though we haven t write down the proof procedure here because the method is standard and even can be performed simply by pressing the enter key with help of some eistent Maple or Mathematica programmes related to the Painlevé test, say, Ptest by Xu and Li [24]. In the remained part of this paper, we use the Darbou transformation which is one of the most powerful tools in constructing eact solutions for many integrable nonlinear equations to the coupled KdV system (4) to study its positon, negaton, compleiton solutions and the interaction property among these types of ecitations. 2. La pair and Darbou transformations of the coupled KdV system (4). In this section, we apply the Darbou transformation to the coupled KdV system (4) in order to construct eact solutions. Because the coupled KdV system is obtained by imposing comple variables on the KdV equation(5), we can easily obtain the Darbou transformation for the coupled KdV system (4) from those of the KdV equation. It is known that the La pair for the KdV equation (5) reads Φ +UΦ = λφ, (6) Φ t = 4Φ +6UΦ +3U Φ. (7) Substituting Φ = φ 1 iφ 2 into Eqs. (6) and (7), the La pair for Eq. (4) is as follows: φ 1 = vφ 2 +uφ 1 λφ 1, (8) φ 2 = uφ 2 vφ 1 λφ 2, (9) and φ 1t = 4φ 1 +6vφ 2 +6uφ 1 +3v φ 2 +3u φ 1, (1) φ 2t = 4φ 2 +6uφ 2 6vφ 1 +3u φ 2 3v φ 1, (11) then it is easy to prove that the compatibility conditions of four equations (8) (11) are eactly the equation system (4). 4

5 With the help of the Darbou transformation for the KdV equation (5), we can easily obtain the first step Darbou transformation for the coupled KdV system u[1] = u ln(f 2 +g 2 ), (12) v[1] = v 2arctan and the new wave functions are transformed to φ 1 = φ 1 ln(f2 +g 2 ) φ 1 arctan 2 ( ) g, (13) f ( ) f φ 2, (14) g ( ) φ 2 = φ 2 ln(f2 +g 2 ) f φ 2 +arctan φ 1, (15) 2 g where{f,g}isawavefunctionsolutionofthelapair(8) (11)withλ = λ, f = φ 1, g = φ 2. In order to construct the second step Darbou transformation, we need two wave function seed solutions {φ 11,φ 12 } and {φ 21,φ 22 } with two different spectral parameters λ 1 and λ 2 respectively, then the second step Darbou transformation is u[2] = u [ ln(f 2 +G 2 ) ], (16) [ v[2] = v 2 arctan where two functions F and G are given by and ( )] G, (17) F F = W(φ 11,φ 21 ) W(φ 12,φ 22 ), (18) G = W(φ 11,φ 22 ) W(φ 12,φ 21 ), (19) respectively with W(a,b) = ab ba is the usual Wronskian determinant. Similarly, the third and fourth step Darbou transformations are written down as follows simply, u[3] = u [ln(h 2 +L 2 )], (2) 5

6 where [ ( )] L v[3] = v 2 arctan, (21) H H = W(φ 11,φ 21,φ 31 ) W(φ 11,φ 22,φ 32 ) W(φ 12,φ 22,φ 31 ) W(φ 12,φ 21,φ 32 ), and L = W(φ 12,φ 22,φ 32 ) W(φ 12,φ 21,φ 31 ) W(φ 11,φ 22,φ 31 ) W(φ 11,φ 21,φ 32 ) with {φ i1, φ i2 }(i = 1,2,3) are solutions of the La pairs (8) (11) with three different parameters λ 1,λ 2 and λ 3. The fourth step Darbou transformation is more complicated because more wave functions are included in. The final results reads u[4] = u [ln(p 2 +Q 2 )], (22) [ ( )] Q v[4] = v 2 arctan, (23) P where P P(,t) and Q Q(,t) are and P(,t) = W(φ 11,φ 21,φ 31,φ 41 )+W(φ 12,φ 22,φ 32,φ 42 ) W(φ 12,φ 21,φ 32,φ 41 ) W(φ 11,φ 21,φ 32,φ 42 ) W(φ 12,φ 22,φ 31,φ 41 ) W(φ 12,φ 21,φ 31,φ 42 ) W(φ 11,φ 22,φ 31,φ 42 ) W(φ 11,φ 22,φ 32,φ 41 ), Q(,t) = W(φ 12,φ 22,φ 32,φ 41 )+W(φ 12,φ 21,φ 32,φ 42 )+W(φ 11,φ 22,φ 32,φ 42 ) +W(φ 12,φ 22,φ 31,φ 42 ) W(φ 11,φ 21,φ 32,φ 41 ) W(φ 11,φ 21,φ 31,φ 42 ) W(φ 12,φ 21,φ 31,φ 41 ) W(φ 11,φ 22,φ 31,φ 41 ), respectively, while {φ i1,φ i2 }, (i = 1,2,3,4) are four wave function vector of the La pair (8) (11) corresponding to four parameters λ i,(i = 1,2,3,4). by In general, the N-step Darbou transformation for the coupled KdV system (4) is given u[n] = u [ln(w 2 r +W2 I )], (24) 6

7 [ v[n] = v 2 arctan ( Wr with W r and W I being the real and imaginary part of the Wronskian of the N comple wave functions i.e., W I )] Φ i φ i1 +Iφ i2, i = 1, 2,..., N, (25) W = W(Φ 1, Φ 2,..., Φ N ) W r +IW I. (26) 3. Positon, Negaton and compleiton solutions In this section, we systematically study the positon, negaton and compleiton solutions for the coupled KdV system by means of the Darbou transformations given in the last section Positon solution, λ >. Based on the first step Darbou transformation, analytical positon solutions for the coupled KdV system can be constructed directly. Taking the seed solution as {u =,v = }, then solving the La pair (8)-(11) directly, one can find f = φ 1 = C 2 sin( λ+4 λ 3 t)+c 1 cos( λ+4 λ 3 t) = C 1 cos( λ+4 λ 3 t+δ 1 ) C 1 cosξ 1, (27) g = φ 2 = C 4 sin( λ+4 λ 3 t)+c 3 cos( λ+4 λ 3 t) = C 3 cos( λ+4 λ 3 t+δ 2 ) C 3 cosξ 2. (28) Substituting the result equations (27) and (28) into the first step Darbou transformation yields the single general positon solution for the coupled KdV system (4): u = 2λ[(C2 1 C 2 3)(C 2 1 cos 2 ξ 1 C 2 3 cos 2 ξ 2 )+4C 2 1C 2 3 cos(δ 1 δ 2 )cosξ 1 cosξ 2 ] (C 2 1 cos2 ξ 1 +C 2 3 cos2 ξ 2 ) 2, (29) v = 4C 1C 3 λ[(c 2 1 cos2 ξ 1 C 2 3 cos2 ξ 2 )cos(δ 1 δ 2 ) (C 2 1 C2 3 )cosξ 1cosξ 2 ] (C 2 1 cos2 ξ 1 +C 2 3 cos2 ξ 2 ) 2. (3) 7

8 2 u Fig.1a Fig.1b 3 v Figure 1. Structur of the Positon solution (a) for the field u epressed by (29) and (b) for the quantity v given by (3) with the parameter selections (31) at time t =. It is easy to find that the positon solution (29) (3) is always nonsingular for v. Actually, from (29) (3) we know that the positon is singular only for δ 2 = δ 1 +nπ, n =, ±1, ±2,..., while this constant condition leads to v =. This result coincides with the fact that the positon of the real KdV is singular. Fig. 1 shows the nonsingular positon structure (29) (3) with the constant parameter selections at time t =. C 1 = 2, C 3 = 4, δ 1 =, δ 2 = 1, λ = 4 (31) From Fig. 1, and the epressions (29) and (3), we know that a nonsingular positon is just a special nonsingular periodic wave solution. 8

9 3.2. Negaton solution (λ = k 2 < ). Inordertoobtainthenegatonsolutionsforcoupled KdV system (4), we only need to substitute λ = k 2 into Eqs. (27) and (28), then the wave functions are f = C 2 e k 4k3t +C 1 e k+4k3 t = c 1 sinh(k 4k 3 t+ 1 ) c 1 sinhη 1, (32) g = C 4 e k 4k3t +C 3 e k+4k3 t = c 2 sinh(k 4k 3 t+ 2 ) c 2 sinhη 2. (33) Substituting the above wave functions into the epressions of the first step Darbou transformation, we have the single negaton solution for the coupled KdV system (4) u = 2k2[ (c 2 1 c 2 2)(c 2 1sinh 2 η 1 c 2 2sinh 2 η 2 )+4c 2 1c 2 2cosh( 1 2 )sinhη 1 sinhη 2 ] (c 2 1 sinh2 η 1 +c 2 2 sinh2 η 2 ) 2,(34) v = 2c 1c 2 k 2 [(c 2 1 sinh(2η 1)+c 2 2 sinh(2η 2)]sinh( 1 2 ) (c 2 1sinh 2 η 1 +c 2 2sinh 2 η 2 ) 2. (35) Fig. 2 shows a special negaton structure under the parameter selections c 1 = 1, c 2 = 5, 1 =, 2 = 2, k = 2 (36) at time t =. From the negaton epression (34) (35), one can see that the negaton solution of the model is nonsingular ecept that 2 = 1 +nπ 1 (n =, ±1, ±2,...) which is corresponding to the real KdV case, v = Negaton interaction solutions. With the N-step Darbou transformations, one may obtain much richer solution structures including the interactions among different types of localized ecitations. A type of soliton solution including two negatons can be constructed easily from the two-step Darbou transformation. Let two wave function vectors, {φ 11,φ 12 } and {φ 21,φ 22 }, are solutions of the La pair (8) (11) with the seed {u =, v = } and two spectral parameters {λ 1 = k 2 1,λ 2 = k 2 2 } respectively, then we have the two negaton 9

10 Fig.2a u 3 Fig.2b 2 1 v Figure 2. A typical negaton structure of the coupled KdV system (4) (a) for u epressed by (34) and (b) for v given by (35) under the parameter selections (36) at time t =. solution epressed by (16) (19) with φ 11 = c 1 cosh(k 1 4k 3 1 t+δ 1), (37) φ 12 = c 2 cosh(k 1 4k 3 1 t+δ 2), (38) φ 21 = c 3 cosh(k 2 4k 3 2 t+δ 3), (39) φ 22 = c 4 cosh(k 2 4k 3 2t+δ 4 ). (4) Fig. 3 and Fig.4 display the special two-negaton interaction procedure for the fields u and v described by (16) (19) with zero seed, wave function selections (37) (4) and the parameter selections δ 1 =, c 1 = k 1 = δ 2 = 1, k 2 = 2, c 2 = δ 3 = 3, δ 4 = 4, c 3 = c 4 = 5 (41) at times t =.25,.5,.1, and.35 respectively. 1

11 Fig.3a u Fig.3b u Fig.3c u Fig.3d u 3 4 Figure 3. Negaton interaction solution epressed by (16) (19) with zero seed, wave function selections (37) (4) and the parameter selections (41) for the field u at times (a)t =.25; (b).5; (c).1 and (d).35 respectively. 11

12 Fig.4a 2 1 v Fig.4b 3 v Fig.4c 3 v Fig.4d v Figure 4. Negaton interaction solution for the quantity v with the same wave function and parameter selections as Fig. 3 12

13 Fig.5 1 u t.4 Figure 5. Positon-negaton interaction solution epressed by (16) and (18)- (19), (42)-(45) for the coupled KdV system (4) Positon-negaton interaction solution. In the same way, by selecting the spectral parameters appropriately in the N-step Darbou transformations one can obtain the interaction solutions among positons and negatons. For instance, based on the second step Darbou transformation by choosing two different spectral parameters with opposite sign, λ 1 = κ 2, λ 2 = k 2, then we can obtain the positon-negaton interaction solutions for the coupled KdV system (4). To show a concrete eample, we further fi the spectral parameters as λ 1 = 1, λ 2 = 4 and select the corresponding spectral functions in the forms φ 11 = 4sin(+4t)+3cos(+4t), (42) φ 12 = 2sin(+4t) 6cos(+4t). (43) φ 21 = 4e 2+32t +e 2 32t, (44) φ 22 = 4e 2+32t +.1e 2 32t. (45) Substituting (42)-(45) and (18)-(19) into (16)-(17) with the initial seed solution {u =, v = }, a special positon-negaton solution follows immediately. Fig. 6 and Fig. 7 show the interaction procedure epressed by (16)-(19), (42)-(45) and the initial solution{u =, v = } for the fields u and v respectively. 13

14 Fig.6 v t.4 Figure 6. Positon-negaton solution epressed by (17)-(19), (42)-(45) for the coupled KdV system (4) Analytical compleiton solution. For many integrable system, compleiton solutions have been found [1]. But it should be pointed out that the known compleiton solutions have singularities. It is lucky that the analytical compleiton solutions are found in this Letter for the first time. Let two spectral parameters λ and λ be conjugated comple number, that is λ = λ 1 +Iλ 2 and λ = λ 1 Iλ 2, then we arrive at two coupled wave solutions with the initial solution {u =,v = }, φ 11 = IC 1 e (4Iλ3 1 12λ2 1 λ 2 12Iλ 1 λ λ3 2 )t+(λ 1I λ 2 ) +IC 2 e (12λ2 1 λ 2 4Iλ Iλ 1λ 2 2 4λ3 2 )t+(λ 2 Iλ 1 ), (46) φ 12 = IC 3 e (4Iλ3 1 12λ2 1 λ 2 12Iλ 1 λ λ3 2 )t+(λ 1I λ 2 ) +IC 4 e (12λ2 1 λ 2 4Iλ Iλ 1λ 2 2 4λ3 2 )t+(λ 2 Iλ 1 ), (47) φ 21 = c 1 e (4Iλ λ2 1 λ 2 12Iλ 1 λ 2 2 4λ3 2 )t+(λ 1I+λ 2 ) +c 2 e ( Iλ 1 λ 2 )+(12Iλ 1 λ λ3 2 4Iλ3 1 12λ2 1 λ 2)t, (48) φ 22 = c 3 e (4Iλ λ2 1 λ 2 12Iλ 1 λ 2 2 4λ3 2 )t+(λ 1I+λ 2 ) +c 4 e ( Iλ 1 λ 2 )+(12Iλ 1 λ λ3 2 4Iλ3 1 12λ2 1 λ 2)t. (49) Substituting (46)-(49) into (18)-(19), we can obtain F = 4λ 2 (c 1 C 1 c 3 C 3 )sin[2λ 1 (4tλ tλ 2 2)] +4λ 1 (c 4 C 3 c 2 C 1 )cosh[2λ 2 (+12tλ 2 1 4tλ 2 2)], (5) G = 4λ 2 (c 3 C 1 +c 1 C 3 )sin[2λ 1 (4tλ tλ2 2 )] +4λ 1 (c 2 C 3 +c 4 C 1 )cosh[2λ 2 (+12tλ 2 1 4tλ2 2 )], (51) 14

15 where {c 1,c 2,c 3,c 4,C 1,C 2,C 3,C 4,λ 1,λ 2 } are all constants. To show a detailed structure of the compleiton, we fi the constant parameters further. The substitution of (5) and (51) into (16)-(17) with the fiing constants {c 1 = 5,c 2 =,c 3 =,c 4 = 5,C 1 = 1,C 2 = 2,C 3 = 2,C 4 = 1,λ 1 = 2,λ 2 = 1 }, yields the following analytical nonsingular compleiton solution 2 U = [ ln(f 2 +G 2 ) ] = A(,t) B(, t) (52) and where [ V = 2 arctan ( )] G = C(,t) F B(,t), (53) A(,t) = [272cos(14t+8)+24+24cos(14t+8)cosh(2+94t) +272cosh(2+94t) 128sin(14t+8)sinh(2+94t)], B(,t) = cos(14t+8)+ 1 8 cos(28t+16)+136cosh(2+94t) 8cos(14t+8)cosh(2+94t)+32cosh(4+188t), and C(,t) = [1862sin(52t+4)cosh(+47t)+3cosh(+47t)sin(156t+12) +24cos(52t+4)sinh(+47t)+16sinh(+47t)cos(156t+12) +48sin(52t+4)cosh(3+141t)+256cos(52t+4)sinh(3+141t)]. The compleiton solution displayed above is analytical without any singularity and it is the first time to find this type of analytical compleiton solutions for the integrable systems to our knowledge. Fig. 7 and Fig. 8 show the detailed structures of the fields u and v epressed by (52) and (53) respectively. 4. Summary and discussion In this Letter, eact solutions of a coupled KdV system which is obtained as a special case of the general coupled KdV systems derived from the two layer fluid dynamical systems and two-component Bose-Einstein condensates are studied in detail. The model can also be 15

16 Fig.7a u Fig.7b 4 -u t Figure 7. Analytical compleiton solution epressed by (52) for the coupled KdV system (4). (a) The structure plot at the fied time t =.2; (b) The evolution plot. simply read off from the well known single KdV equation by assuming its field is comple. Thanks to the model can be obtained from the comple KdV equation, the La pair of the coupled KdV system follows immediately and so does the Darbou transformation. To study the detailed structures of the localized and periodic ecitations of the model, the first four step Darbou transformations are eplicitly given in terms of the real Wronskian determinant forms while the general N-step Darbou transformation is written down with help of the comple Wronskian determinant form. Starting from the trivial seed solution, the first step of the Darbou transformation leads to two types of eact solutions, the negaton and the positon which are related to negative and positive spectral parameters respectively. Though the single negaton solution for other integrable models, say, the single real KdV system, can be both singular and nonsingular, 16

17 8 Fig.8a v Fig.8b v t Figure 8. Nonsingular compleiton solution given by (53) for the coupled KdV system (4). (a) The structure plot at the fied time t =.2; (b) The evolution plot. it can be only analytic for the coupled KdV system (4). Similarly, the positon solution can also be only analytic while for the single real KdV equation there is no analytic positon. The solutions obtained from the second step Darbou transformation and the trivial zero seed can be two negaton interaction solution, two positon solution, negaton-positon solution and a single compleiton dependent on the selections of the spectral parameters. The interaction among negatons is completely elastic. The positon-negaton solution can be considered as a single soliton solution with a periodic wave background. The single real compleiton solution yields if two spectral parameters are comple conjugates. The known compleiton solutions for other integrable systems are singular[1]. It is interesting that the compleiton solution obtained here for the coupled KdV system is nonsingular! This type of nonsingular compleiton solution may eist for other types of coupled integrable systems. 17

18 Though the coupled KdV system studied here can be read off from the comple KdV equation, the solution properties are quite different. For the real single component KdV equation, there is no nonsingular positons and compleitons. However, for the coupled KdV system, negatons, positons and compleitons can be nonsingular. For the single KdV system, various other types of eact solutions, such as the rational solutions, cnoidal wave solutions, algebraic geometry solutions and τ function solutions have been obtained by many authors under different approaches. The similar solutions for the coupled KdV system (4) can be easily obtained in some similar ways. However, the detailed properties of the solutions, especially there analytic behavior, must be quite different. The more about the model and its possible applications in fluids, atmospherical dynamics and the Bose-Einstein condensates are worthy of studying further. The authors are indebt to thanks Dr. X. Y. Tang and Prof. Y. Chen for their helpful discussions. The work was supported by the National Natural Science Foundation of China (No and No ). References [1] W. X. Ma, Phys. Lett. A 31 (22) 35; W. X. Ma and K. Maruno, Physica A 343 (24) 219. [2] V. B. Matveev and M. A. Salle Darbou transformations and Solitons, Springer-Verlag 1991; V. B. Matveev, Phys. Lett. A 166 (1992) 25; 29; 168 (1992) 463; J. Math. Phys. 35 (1994) [3] C. Rasinariu, U. Sukhatme and A. Khare, J. Phys. A: Math. Gen, 29 (1996) 183. [4] M. D. Kruskal, Lect. Appl. Math. 15 (1974) 61; H. Airault, H. P. Mckean and J. Moser, Commun. Pure Appl. Math. 3 (1977) 95; M. J. Ablowitz and H. Cornille, Phys. lett. A 72 (1979) 277. [5] Z. J. Lian and S. Y. Lou, Nonlinear Analysis, (25) in press. [6] Y. Chen and Q. Wang, Phys. Lett. A. (25) in press. [7] S. Y. Lou, H. C. Hu and X. Y. Tang, Phys. Rev. E., 71 (25) 3664; H. C. Hu and S. Y. Lou, Phys. Lett. A 341 (25) 422. [8] S. Y. Lou, Phys. Lett. A 277 (2) 94; X. Y. Tang, S. Y. Lou and Y. Zhang, Phys. Rev. E., 66 (22) 4661; H. C. Hu, X. Y. Tang, S. Y. Lou and Q. P. Liu, Chaos, Solitons and Fractals, 22 (24) 327. [9] R. M. Miura, C. S. Gardner and M. D. Kruskal J. Math. Phys. 9 (1968) 124; S. Y. Lou, J. Math. Phys. 35 (1994) 239. [1] F. Magri J. Math. Phys. 19 (1978) [11] V. E. Zakharov and L. D. Faddeev, Functional Anal. Appl. 5 (1971) 28. [12] V. V. Sokolov and T. Wolf, J. Phys. A: Math. Gen. 34 (21) V. V. Sokolov and S. I. Svinolupov, Theor. Math. Phys. 1 (1994) 959. [13] R. Hirota and J. Satsuma, Phys. Lett. A 85 (1981) 47. [14] R. Hirota and J. Satsuma, J. Phys. Soc. Japan 51 (1982) 339. R. Dodd and A. P. Fordy Phys. Lett. A 89 (1982) 168. H. C. Hu and Q. P. Liu, Chaos, Solitons & Fractals 17 (23) 921. [15] M. Ito, Phys. Lett. A 91 (1982)

19 [16] V. G. Drinfeld and V. V. Sokolov, Sov. J. Math. 3 (1985) [17] B. Fuchssteiner, Prog. Theor. Phys. 68 (1982) 182. [18] Y. Nutku and O. O guz, Nuovo Cimento Soc. Ital. Fis. B 15 (199) [19] A. Y. Zharkov, J. Symb. Comput. 15 (1993) 85. [2] M. V. Foursov, Inverse Problems 16 (2) 259. [21] J. A. Gear and R. Grimshaw, Stud. Appl. Math. 7 (1984) 235; J. A. Gear, Stud. Appl. Math. 72 (1985) 95. [22] S. Y. Lou, B. Tong, H. C. Hu and X. Y. Tang, Coupled KdV equations derived from atmospherical dynamics, preprint ArXiv: nlin.si/5829. [23] V. A. Brazhnyi and V. V. Konotop, Phys. Rev. E 72 (25) [24] G. Q. Xu and Z. B. Li, Comp. Phys. Commun. 161 (24)