Abstract and Applied Analysis

Volume 2013 (2013), Article ID 656297, 15 pages

http://dx.doi.org/10.1155/2013/656297

Research Article

A Generalized KdV Equation of Neglecting the Highest-Order Infinitesimal Term and Its Exact Traveling Wave Solutions

Xianbin Wu,¹ Weiguo Rui,² and Xiaochun Hong³

¹Junior College, Zhejiang Wanli University, Ningbo 315100, China

²College of Mathematics, Honghe University, Mengzi, Yunnan 661100, China

³College of Statistics and Mathematics, Yunnan University of Finance and Economics, Kunming, Yunnan 650221, China

Received 3 November 2012; Accepted 15 January 2013

Academic Editor: Julian López-Gómez

Copyright © 2013 Xianbin Wu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

We study a generalized KdV equation of neglecting the highest order infinitesimal term, which is an important water wave model. Some exact traveling wave solutions such as singular solitary wave solutions, semiloop soliton solutions, dark soliton solutions, dark peakon solutions, dark loop-soliton solutions, broken loop-soliton solutions, broken wave solutions of U-form and C-form, periodic wave solutions of singular type, and broken wave solution of semiparabola form are obtained. By using mathematical software Maple, we show their profiles and discuss their dynamic properties. Investigating these properties, we find that the waveforms of some traveling wave solutions vary with changes of certain parameters.

1. Introduction

In 1995, based on the physical and asymptotic considerations, Fokas [1] derived the following generalized KdV equation:

which is an important water wave model, where , , , , , , , , . Regarding the as free parameters and using the to replace the (i.e., ), (1) becomes the following PDE:

which is given by Tzirtzilakis et al. in [2]; they called it high-order wave equation of KdV type. Just as Tzirtzilakis [2] said these two equations are both water wave equations of KdV type, which are more physically and practically meaningful. The motion described by the model (1) or (2) is a 2-dimensional, inviscid, incompressible fluid (water) lying above a horizontal fiat bottom located at ( is a constant), and letting the air above the water. It turns out that, for such a system if the vorticity is zero initially, it remains zero. The fluids (water) analyzed by Fokas are only irrotational flows. In addition, this system is characterized by two parameters , with , , where and are two typical values of the amplitude and of the wavelength of the waves. The parameters , , satisfy the condition because the system is a model of short amplitude and long wavelength. Therefore the parameters and satisfy the condition , .

Assuming that the waves are unidirectional and neglecting terms of , (1) can be reduced to the classical KdV equation:

In [1], Fokas assumed that is less than ; this implies . According to this assumption, we easily know that and . Neglecting two high-order infinitesimal terms of , (1) can be reduced to another high-order wave equation of KdV type [2–5] as follows:

Equation (4) is a special case of (1) for . In [1], it was observed that (4) can be reduced by the local transformation of coordinates

to a completely integrable PDE as follows:

Equation (6) was first derived in [6] by using the method of bi-Hamiltonian systems, which based on the physical considerations, its Lax Pair was also given in [7].

If only neglecting the highest order infinitesimal term of , then (1) can be reduced to a new generalized KdV equation as follows:

We call it a generalized KdV equation of neglecting the highest order infinitesimal term. In fact, (7) is another special case of (1) for ; it is also third-order approximate equation of KdV type. Of course, on describing dynamical behaviors of water waves, (4) is only a rough approximative model of (1) compared with (7); that is, the precision of model (7) is better than that of model (4) on describing dynamical behaviors of water waves. In other words, model (7) exhibits much richer phenomenology than model (4). Therefore, the investigation of exact traveling wave solutions for (7) are more practically meaningful than that of (4).

On the other hand, under the local transformation

(1) can be reduced to the following generalized Gardner equation [1]:

where are certain expressions of . Clearly, (9) is equivalent to the following generalization form of the modified KdV equation [1, 8] under the transformation ,

From the above references and the references cited therein, we can know that (1) and (2) are very important water wave models. However, (1) and (2) are too complex to obtain their exact solution under universal conditions. Only under some special parametric conditions, their exact solutions were obtained in the existing literature. Next, let us briefly review the research backgrounds for the above equations.

In [2], Tzirtzilakis et al. only obtained two soliton-like solutions of (2) in the two groups of special conditions , , , , and , , , , . In [9], by using the planar bifurcation method of dynamical systems, under the four groups of special conditions , , , ; , , ; , , ; , , , , Li et al. studied (2); the existence of all kinds of traveling wave solutions were discussed completely, but its exact solutions were not investigated although some results of numerical simulation were obtained in this literature. In [10], in order to answer what is the dynamical behavior of one-loop soliton solution, Li studied the special case of (1) for . In [11], under different kinds of parametric conditions, Marinakis discussed two integrable cases for the third-order approximation model (1). In [12], Marinakis proved that (1) and some its of special cases are integrable. In [13], Gandarias and Bruzon proved that (1) is self-adjoint if and only if . In [14], by using the method as in [9], Li et al. studied (10); the existence of solitary wave, kink, and antikink wave solutions and uncountably infinite many smooth and nonsmooth periodic wave solutions was discussed except exact solutions. In [8], Bi also obtained some results of numerical simulation of (10); it is a pity that the exact traveling wave solutions of (10) were not obtained yet. In [4, 5], we obtained some exact traveling wave solutions of (4) under the parametric conditions or . In [15], under a new ansätze, Khuri studied (4); some exact solitary wave solutions and periodic wave solutions were obtained. In [16], by using method of planar dynamical system, Long et al. studied (6); the existence of smooth solitary wave and uncountably infinite many smooth and nonsmooth periodic wave solutions was proved in this literature.

From the above research backgrounds of (1) and (2), we can see that their exact solutions in universal conditions are hard to obtain because they are highly nonlinear equations and most probably they are not integrable equations in general. Thus, large numbers of research results are still concentrated in the classical KdV equation and some other high-order equations with KdV type, such as KdV-Burgers equation [17, 18] and KdV-Burgers-Kuramoto equation [19], at present. Therefore, the investigation of the more exact solutions for (1) is very important and necessary. However, by using the current methods, we can not obtain exact solutions of (1) in universal conditions; the next best thing is the investigation of exact solutions of (7). In this paper, still regarding the as free parameters and by using the integral bifurcation method [20, 21], we will investigate exact traveling wave solutions and their properties of (7).

The rest of this paper is organized as follows. In Section 2, we will derive two-dimensional planar system which is equivalent to (7) and give its first integral equation. In Section 3, by using the integral bifurcation method, we will obtain some new traveling wave solutions and discuss their dynamic properties.

2. Two-Dimensional Planar Dynamical System of (7) and Its First Integral and Conservation of Energy

Making a transformation with , (7) can be reduced to the following ODE:

where is wave velocity which moves along the direction of -axis and . Integrating (11) once and setting the integral constant as zero yields

Let . Thus (12) can be reduced to a planar system

and a linear equation

Obviously, the solutions of (12) cover the solutions of (13) and (14). We notice that the second equation in (13) is not continuous when ; that is, the function is not defined by . So (13) is a singular system; the line is called the singular line. Thus, we make the following transformation:

where is a free parameter. Under the transformation (15), (13) can be rewritten as the following system

Clearly, (16) is equivalent to (12). Except for the singular line , (13) and (16) have the same first integral.

First, from the parametric conditions in the model (1), we easily know that . Second, from [9–11, 13], we know that (1) contains many good properties including loop soliton, integrability, and self-adjoint property when . In addition, in the condition , it becomes easy when we solve (7). Therefore, we only consider the special case of (7) for in this paper. Thus, by using this assumption, we obtain the first integral of (13) and (16) as follows:

where is an integral constant. For the convenience of discussion, taking the integral constant in (17) yields

where denotes the particle's mass of system. Let , , . Thus, (18) can be rewritten as the following equation of conservation of energy:

where denotes kinetic energy, denotes external energy, and denotes potential energy. In (19), the kinetic energy and potential energy are not conserved because the external energy , but the global energy () are subject to the energy conservation. if only if . Obviously, once ; under this case, (7) becomes the following equation:

Similarly, (20) can be reduced to the following planar system:

Obviously, the system (21) is a regular system; it has no singular line. Taking the integral constant as zero, we obtain (21)'s first integral as follows:

Equation (22) can be rewritten as the following equation of conservation of energy:

where are given above. The kinetic energy and potential energy are conserved in (23); in other words, the particle motion satisfies the conservation of kinetic energy and potential energy, which converts the kinetic energy from the potential energy then it converts the potential energy from the kinetic energy and go round and round. Therefore, (20) only has nonsingular traveling wave solutions, and all solutions are smooth; this is very different from (7).

In order to discuss singular or nonsingular traveling wave solutions of (7), we will consider (17). When , (17) can be reduced to

Equation (24) can be further simplified under the following parametric conditions.

Case 1. One has .

Case 2. One has .

Case 3. One has .

Case 4. One has .

Case 5. One has or or .

Under the parametric conditions of Case 1, (24) can be reduced to

Under the parametric conditions of Case 2, (24) can be reduced to

Under the parametric conditions of Case 3, (24) can be reduced to

Under the parametric conditions of Case 4, (24) can be reduced to

Under the the parametric conditions of Case 5, (24) can be reduced to

3. Exact Traveling Wave Solutions of (7) and Their Dynamic Properties

In this section, by using (25), (26), (27), (28), and (29), we will derive many exact traveling wave solutions and discuss their dynamic properties.

3.1. The Exact Solutions under the Parametric Conditions of Case 1

Substituting (25) into the left expression (i.e., ) of (16) yields

For the sake of convenience, in all the following discussions, we only discuss the case where the right of equation is “+” sign; the case of “−” sign can be similarly discussed. By the way, the solutions obtained by taking “+” sign and the solutions obtained by taking “−” sign are same when the is even function.

Taking “+” in (30), it can be reduced to

Let . Equation (31) can be rewritten as

Taking the as initial values, respectively, integrating (32) we obtain

where are two roots of equation .

(i) When , completing the integrals (33) yields

where , and with .

By using (34) and the transformation , we obtain four exact solutions of (31) as follows:

where the constants are given above.

Substituting (35) into (15) yields

where is a normal elliptic integral of the first kind and .

Substituting (36) into (15) yields

where , and .

Combining (35) with (37), we obtain (7)'s two exact solutions of parametric type as follows:

Combining (36) with (38), we obtain (7)'s another two exact solutions of parametric type as follows:

(ii) When , completing the integrals (33) yields

Similarly, by using (41) and the transformation , we obtain four periodic solutions of (31) as follows:

where ; that is, . Respectively, substituting (42) and (43) into (15) yields

where .

Combining (42) with (44), we obtain (7)'s two exact solutions of parametric type as follows:

where .

Combining (43) with (45), we obtain (7)'s another two exact solutions of parametric type as follows:

Among the above traveling wave solutions, it is very worthy to mention solutions (40), (46), and (47); they have some peculiar dynamical properties and their traveling wave phenomena are very interesting. Solution (40) denotes a singular solitary wave; its shape is very similar to bright soliton, but it is not a normal soliton; the left parts of waveform degenerate into a crook. Solution (46) denotes a bounded wave with level asymptote; its shape is a half of a whole loop soliton; we call it semiloop soliton. In order to describe the dynamic properties of these two traveling wave solutions intuitively, we, respectively, plot profile of solutions (40) and (46) by using software Maple, which are shown in Figure 1.

(a) Singular solitary wave 

(b) Semiloop soliton

Figure 1: Profiles of solutions (40) and (46) for the fixed parameters.

Figure 1(a) shows a shape of singular solitary wave under the fixed parameters , , , . Figure 1(b) shows a shape of semiloop soliton under the fixed parameters , , , .

It is very interesting that the solution (47) has six kinds of waveforms. Solution (47) with “+”, respectively, denotes a dark peakon, a dark loop soliton, and a broken loop soliton when the parameter decreases from to , which are shown in Figures 2(a)–2(d). Solution (47) with “−”, respectively, denotes a bright soliton and a singular compacton when the parameter decreases from to , which are shown in Figures 2(e) and 2(f). In other words, its waveforms depend on parameter extremity. Similar traveling wave phenomena that one solution contains multiwaveform first appeared in the investigation of the Degasperis-Procesi equation [22].

(a) Dark soliton

(b) Dark peakon

(c) Dark loop soliton

(d) Broken loop soliton (broken wave of C-form)

(e) Bright soliton

(f) Singular compacton

Figure 2: Different kinds of waveforms on solution (47) under different parametric conditions.

Figure 2(a) shows a shape of dark soliton under the fixed parameters , , , and . Figure 2(b) shows a shape of dark peakon under the fixed parameters , , , and . Figure 2(c) shows a shape of dark loop soliton under the fixed parameters , , , and . Similar traveling wave phenomena also appear in [5], in “Concluding remarks” of this literature; by using these phenomena, we successfully explain the movement of water waves. Figure 2(d) shows a shape of broken wave of C-form (upward hatch C) under the fixed parameters , , , and ; we also call it broken loop soliton or open-orbicular wave. Figure 2(e) shows a shape of bright soliton under the fixed parameters , , , and . Figure 2(f) shows a shape of singular compacton under the fixed parameters , , , and ; we call it singular compacton because it is not a normal compacton. Specifically, both sides of compacton wave become upward crook; it very much likes solitary waves but it is not a solitary wave after all because the left and right parts of waveform do not extend as .

3.2. The Exact Solutions under the Parametric Conditions of Case 2

Notice that the characters of (26)–(29) are quite different from (25). Therefore, by using the expression and (26)–(29), we will investigate the implicit solutions of (7) next.

Substituting (26) into the expression of (13) yields

Taking “+”, (48) can be reduced to

(i) When , (49) can be reduced to

By the Ferrari method or the Descartes method, it is easy to know that the equation

has four roots as follows:

where with . Especially, when , these four roots can be reduced to . For the convenience of discussion, we denote the above two real roots by the signs and always assume that the biggest root is denoted by . We also denote the above two complex roots by the signs . Thus (50) can be rewritten as

where . Taking as the initial values and integrating (53), we obtain implicit solution of (7) as follows:

where , , , ,

(ii) When , (49) can be reduced to

Similarly, (56) can be rewritten as

where are four roots of (51) and . Taking as the initial values then integrating (56), we obtain another implicit solution of (7) as follows:

where the constants are given above, the , and

3.3. The Exact Solutions under the Parametric Conditions of Case 3

As in Section 3.2, substituting (27) into the expression of (13) yields

Taking “+”, (60) can be rewritten as

When , (61) can be reduced to

When , the equation

has four roots as follows:

where with . When , (63) has another four roots .

(i) Particularly, when and , (62) can be reduced to

According to the above cases, we find that (63) still has two real roots and two complex roots; this case is very similar to that in (51) except their especial cases are different. Therefore, the types of their solutions are the same when and . Thus, we omit these parts of discussions. The especial case of (62) (i.e., (65)) can be reduced to

where . Taking the as initial values and integrating (66), we obtain (7)'s periodic wave solution of implicit function as follows:

(ii) Especially, when and , (61) can be reduced to

where . Taking the as initial values and integrating (68), we obtain (7)'s implicit solitary wave solution as follows:

3.4. The Exact Solutions under the Parametric Conditions of Case 4

Substituting (28) into the expression of (13) yields

Taking “+”, (70) can be rewritten as

When , (71) can be reduced to

Next, we discuss the roots of the following equation:

When , (73) has four real roots as follows:

When , (73) has still four real roots as follows:

where is a double root.

When , (73) has two real roots and two complex roots as follows:

Thus, by using the above parametric conditions, we can obtain three kinds of exact solutions of implicit function type; see the following discussions.

Under the conditions and , (72) can be reduced to

where , , , and . Taking as initial values and integrating (77), we obtain (7)'s periodic wave solution of implicit function type as follows:

where .

Under the conditions and , (72) can be reduced to

Equation (79) can be rewritten as

Taking as initial values and integrating (80), we obtain (7)'s solitary wave solution of implicit function type as follows:

Under the conditions and , (72) can be reduced to

where and . Taking as initial values and integrating (82), we obtain (7)'s periodic wave solution of implicit function type as follows:

where , , , , .

(ii) When , (72) can be rewritten as

Under the conditions and , (84) can be reduced to

where , and . Taking as initial values and integrating (85), we obtain (7)'s periodic wave solution of implicit function type as follows:

where .

Under the conditions and , (84) can be reduced to

Respectively, taking and (0,0) as initial values, integrating (87), we obtain (7)'s periodic wave solution of implicit function type as follows:

Under the conditions and , (84) can be reduced to

where and . Taking as initial values and integrating (89), we obtain (7)'s periodic wave solution of implicit function type as follows:

where constants are given above, the , and

3.5. The Exact Solutions under the Parametric Conditions of Case 5

Under different parameters, (29) can be reduced different types. As an example, here we only discuss one reduced to type under the parametric conditions for (29); the others can be similarly discussed.

Especially, when and , (29) can be reduced to

Also, we only consider one case where the equation has three real roots here. The other two cases (i.e., two real roots and one real root) can be similarly discussed. Because the expressions of these three roots are very complex, we use the symbols instead of them. In fact, once the parameters are fixed, these three roots can be solved by the above equation. For example, if , then . In addition, if , then ; if , then . Thus, under the above parametric conditions, (92) can be rewritten as

or

where .

Respectively, substituting (93) and (94) into the equation of (13) yields

Taking as initial values and integrating (95), respectively, we obtain (7)'s two groups of exact solutions of implicit function type as follows:

where , , , , , , , .

Among the above traveling wave solutions, it is also worthy to mention solutions (81), (88), (96), and (97); they have some peculiar dynamical properties, too. Solution (81) denotes a bounded wave with level asymptote under the fixed parameters , , . Solution (88) denotes a periodic wave of singular type under the fixed parameters , , . Solution (96) denotes a broken wave of U-form (leftward fallen U) under the fixed parameters , , , , , , , , , . Solution (97) denotes a broken wave of semi-parabola form under the fixed parameters , , , , , , , , , . In order to describe the dynamic properties of the above traveling wave solutions intuitively, similarly, we, respectively, draw profile of solutions (81), (88), (96), and (97) by using software Maple, which are shown in Figure 3.

(a) Bounded wave

(b) Periodic wave of singular type

(c) Broken wave of U-form

(d) Broken wave of semiparabola form

Figure 3: The profiles of solutions (81), (88), (96), and (97) for the fixed parameters.

4. Conclusion

In this work, by using the integral bifurcation method, we studied a generalized KdV equation of neglecting the highest order infinitesimal term. Some singular and nonsingular traveling wave solutions including singular solitary wave solutions, semiloop soliton solutions, dark soliton solutions, dark peakon solutions, dark loop-soliton solutions, broken wave solutions of C-form and U-form, periodic wave solutions of singular type, and broken wave solution of semi-parabola form were obtained. Among these traveling wave solutions, the solutions such as (40), (46), (47), (81), (88), (96), and (97) have some peculiar dynamical properties and very interesting traveling wave phenomena. It is very worthy to mention solution (47); it denotes a dark peakon, a dark loop soliton, a broken loop soliton, a bright soliton, and a singular compacton when the parameter varies. This course of changes of waveform shows that solution (47) is a continuous solution. Indeed, the loop soliton which contains (47) is also continuous; this case has proved again that the loop soliton solution of some nonlinear traveling wave equation is one continuous solution, not three breaking solutions; see the arguments in [23, 24]. However, the phenomenon of half loop soliton exists in troth. For example, the semiloop soliton and broken loop soliton show that their waveforms are only a part of a whole loop soliton. Maybe, a loop soliton can be divided into two parts or three parts like the semiloop soliton or three breaking waves which appeared in [24], but every part of them is still continuous, and then joining them together, they come into being a continuous loop soliton again. Indeed, in 2006, Vakhnenko and Parkes's work [25] successfully explained similar phenomena. In [25], the graphical interpretation of the solution for gDPE is presented. In this analysis, the 3D-spiral (whether one loop from a spiral or a half loop of a spiral) has the different projections that is the essence of the possible solutions.

The results obtained in this work are very different from those in other literatures because this special case of (1) for is not studied in other references. Especially, the results of (7) are different from those of (4) because of . In addition, we hope more and more researchers pay attention to investigations of exact solutions for (7) under general conditions in the future work.

Acknowledgments

The authors thank reviewers very much for their useful comments and helpful suggestions. This work was financially supported by the Natural Science Foundation of China (Grant no. 11161038). It was also supported by the Natural Science Foundations of Zhejiang Province (Grant no. Y2111160) and Yunnan Province (Grant no. 2011FZ193).

References

A. S. Fokas, “On a class of physically important integrable equations,” Physica D. Nonlinear Phenomena, vol. 87, no. 1–4, pp. 145–150, 1995.

E. Tzirtzilakis, V. Marinakis, C. Apokis, and T. Bountis, “Soliton-like solutions of higher order wave equations of the Korteweg-de Vries type,” Journal of Mathematical Physics, vol. 43, no. 12, pp. 6151–6165, 2002.

E. Tzirtzilakis, M. Xenos, V. Marinakis, and T. C. Bountis, “Interactions and stability of solitary waves in shallow water,” Chaos, Solitons & Fractals, vol. 14, no. 1, pp. 87–95, 2002.

J. Li, W. Rui, Y. Long, and B. He, “Travelling wave solutions for higher-order wave equations of KdV type. III,” Mathematical Biosciences and Engineering, vol. 3, no. 1, pp. 125–135, 2006.

W. Rui, Y. Long, and B. He, “Some new travelling wave solutions with singular or nonsingular character for the higher order wave equation of KdV type (III),” Nonlinear Analysis: Theory, Methods & Applications, vol. 70, no. 11, pp. 3816–3828, 2009.

B. Fuchssteiner and A. S. Fokas, “Symplectic structures, their Bäcklund transformations and hereditary symmetries,” Physica D. Nonlinear Phenomena, vol. 4, no. 1, pp. 47–66, 1981.

R. Camassa and D. D. Holm, “An integrable shallow water equation with peaked solitons,” Physical Review Letters, vol. 71, no. 11, pp. 1661–1664, 1993.

Q. Bi, “Peaked singular wave solutions associated with singular curves,” Chaos, Solitons & Fractals, vol. 31, no. 2, pp. 417–423, 2007.

J. Li, J. Wu, and H. Zhu, “Traveling waves for an integrable higher order KdV type wave equations,” International Journal of Bifurcation and Chaos in Applied Sciences and Engineering, vol. 16, no. 8, pp. 2235–2260, 2006.

J. Li, “Dynamical understanding of loop soliton solution for several nonlinear wave equations,” Science in China. Series A. Mathematics, vol. 50, no. 6, pp. 773–785, 2007.

V. Marinakis, “Integrable third order equations of the KdV type,” Journal of Mathematical Sciences. Advances and Applications, vol. 5, no. 2, pp. 317–332, 2010.

V. Marinakis, “Higher-order equations of the KdV type are integrable,” Advances in Mathematical Physics, vol. 2010, Article ID 329586, 5 pages, 2010.

M. Gandarias and M. S. Bruzon, “Conservation laws for a third order equation by using freesoftware MAXIMA,” in Recent Researches in Engineering Education and Software Engineering, pp. 48–51, WSEAS Press, 2012.

J. Li and J. Zhang, “Bifurcations of travelling wave solutions for the generalization form of the modified KdV equation,” Chaos, Solitons & Fractals, vol. 21, no. 4, pp. 889–913, 2004.

S. A. Khuri, “Soliton and periodic solutions for higher order wave equations of KdV type. I,” Chaos, Solitons & Fractals, vol. 26, no. 1, pp. 25–32, 2005.

Y. Long, W. Rui, and B. He, “Travelling wave solutions for a higher order wave equations of KdV type. I,” Chaos, Solitons and Fractals, vol. 23, no. 2, pp. 469–475, 2005.

R. Kumar, K. Chong, and J. Park, “Exact null controllability of KdV-burgers equation with memory effect systems,” Abstract and Applied Analysis, vol. 2012, Article ID 531659, 20 pages, 2012.

W. Zhang and X. Li, “Approximate damped oscillatory solutions for generalized KdV-Burgers equation and their error estimates,” Abstract and Applied Analysis, vol. 2011, Article ID 807860, 26 pages, 2011.

J.-M. Kim and C. Chun, “New exact solutions to the KdV-Burgers-Kuramoto equation with the Exp-function method,” Abstract and Applied Analysis, vol. 2012, Article ID 892420, 10 pages, 2012.

W. Rui, B. He, Y. Long, and C. Chen, “The integral bifurcation method and its application for solving a family of third-order dispersive PDEs,” Nonlinear Analysis: Theory, Methods & Applications, vol. 69, no. 4, pp. 1256–1267, 2008.

R. Weiguo, L. Yao, H. Bin, and L. Zhenyang, “Integral bifurcation method combined with computer for solving a higher order wave equation of KdV type,” International Journal of Computer Mathematics, vol. 87, no. 1–3, pp. 119–128, 2010.

V. O. Vakhnenko and E. J. Parkes, “Periodic and solitary-wave solutions of the Degasperis-Procesi equation,” Chaos, Solitons & Fractals, vol. 20, no. 5, pp. 1059–1073, 2004.

K. K. Victor, B. B. Thomas, and C. K. Timoleon, “Reply to "comment on "on two-loop soliton solution of the Schäfer-Wayne short-pulse equation using Hirota's method and Hodnett-Moloney approach"",” Journal of the Physical Society of Japan, vol. 76, no. 11, Article ID 116002, pp. 116002-1–116002-2, 2007.

Y. Zhang and J. Li, “Comment on “on two-loop soliton solution of the Schäfer-Wayne shortpulse equation using hirota’s method and hodnett-molony approach,” Journal of the Physical Society of Japan, vol. 76, no. 11, pp. 116001-1–116001-2, 2007.

V. O. Vakhnenko and E. J. Parkes, “The solutions of a generalized Degasperis-Procesi equation,” Reports of the National Academy of Sciences of Ukraine, no. 8, pp. 88–94, 2006.