AAA Abstract and Applied Analysis 1687-0409 1085-3375 Hindawi Publishing Corporation 629434 10.1155/2014/629434 629434 Research Article Variational Iteration Method for Solving the Generalized Degasperis-Procesi Equation Lijuan Qian 1 Lixin Tian 2 Kaiping Ma 3 Hong Baojian 1 Faculty of Science Jiangsu University Zhenjiang Jiangsu Province 212013 China ujs.edu.cn 2 Nanjing Normal University Nanjing Jiangsu Province 210097 China njnu.edu.cn 3 College of Engineering Nanjing Agricultural University Nanjing Jiangsu Province 210031 China njau.edu.cn 2014 13 3 2014 2014 31 12 2013 10 02 2014 13 3 2014 2014 Copyright © 2014 Qian Lijuan 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.

We introduce the variational iteration method for solving the generalized Degasperis-Procesi equation. Firstly, according to the variational iteration, the Lagrange multiplier is found after making the correction functional. Furthermore, several approximations of u n + 1 ( x , t ) which is converged to u ( x , t ) are obtained, and the exact solutions of Degasperis-Procesi equation will be obtained by using the traditional variational iteration method with a suitable initial approximation u 0 ( x , t ) . Finally, after giving the perturbation item, the approximate solution for original equation will be expressed specifically.

1. Introduction

The theory of soliton has extensive applications in physics, mechanics, and combustion science. In recent years, many researchers studied the soliton theory in the fields of shock wave [1, 2], light scattering, quantum mechanics, atmospheric physics, neural networks, explosion, and combustion . There are many new methods for searching the soliton solution of nonlinear evolution equations such as hyperbolic tangent function method , the homogeneous balance methods , Jacobi elliptic function expansion method , and pseudo-spectral method .

The variational iteration method (VIM) was developed, in 1999, by He . The VIM gives rapidly convergent successive approximations of the exact solution if such a solution exists; otherwise, a few approximations can be used for numerical purposes. The Adomian decomposition method suffers from the complicated computational work needed for the derivation of Adomian polynomials for nonlinear terms. The VIM has no specific requirements, such as linearization, small parameters for nonlinear operators. Therefore, the VIM can overcome the foregoing restrictions and limitations of perturbation techniques, so that it provides us with a possibility to analyze strongly nonlinear problems. On the other hand, the VIM is capable of greatly reducing the size of calculation while still maintaining high accuracy of the numerical solution . Moreover, the power of the method gives it a wider applicability in handling a huge number of analytical and numerical applications. The VIM was successfully applied to study a variety of differential equations. It is based on Lagrange multiplier, and it has the merits of simplicity and easy execution. As a result, it has been proved by many authors to be a powerful mathematical tool for addressing various kinds of linear and nonlinear problem. For example, this method was used for solving nonlinear wave equations and the Laplace equation by Wazwaz . The VIM for solving linear systems of ODEs with constant coefficients was studied by Khojasteh Salkuyeh . Helmholtz equation was researched by Momani and Abuasad . Geng  introduced the piecewise VIM for solving Riccati differential equation. Fractional vibration equation was researched by Das . Furthermore, higher order boundary value problems were researched by Xu ,Noor, and Mohyud-Din . Noor et al.  applied a modified He’s variational iteration method for solving singular fourth-order parabolic partial differential equations. The proposed modification is made by introducing He’s polynomials in the correction functional. Ghorbani and Saberi-Nadjafi  modified the VIM by constructing an initial trial function without unknown parameters. Sevimlican  constructed approximate Green’s function for a vector equation for the electric field by using VIM.

In this paper, we are concerned with the variational iterations method for solving the generalized Degasperis-Procesi equation. As a review, we will recall the VIM briefly in Section 2.

2. Variational Iteration Method

In this section, the basic concepts of variational iteration method (VIM) are introduced. Here, a description of method  is given to handle the general nonlinear problem. Consider the differential equation of the form (1) L u ( x , t ) + N u ( x , t ) = f ( x , t ) , where L is a linear operator, N is a nonlinear operator, and f ( x , t ) is the inhomogeneous term. According to He’s variational iteration method, we can construct a correction functional for (1) as follows: (2) u n + 1 ( x , t ) = u n ( x , t ) + 0 t λ ( τ ) ( L u n ( x , τ ) + N u ~ n ( x , τ ) - f ( x , τ ) ) d τ , where λ is a general Lagrange multiplier, which can be identified optimally via variational theory [12, 24]. Here u ~ n is considered as a restricted variation [14, 25] which means δ u ~ n = 0 ; the subscript n denotes the n th approximations. The successive approximations u n + 1 ( x , t ) , of the solution u ( x , t ) , can be obtained after using the obtained Lagrange multiplier and the zeroth approximation u 0 ( x , t ) , which are selected from any function that satisfies the initial conditions. With λ determined, several approximations u n + 1 ( x , t ) , n 0 follow. Consequently, the exaction solution may be obtained as (3) u ( x , t ) = lim n + u n + 1 ( x , t ) . In fact, the VIM depends on the suitable selection of the initial approximation u 0 ( x , t ) . Moreover, we use a well-known, powerful tool to prove the convergence of the sequence obtained via the VIM and its rate. It is the Banach’s fixed point theorem that follows.

Theorem 1 (Banach’s fixed point theorem).

Assume that X is a Banach space and (4) A : X X is a nonlinear mapping, and suppose that (5) A [ u ] - A [ v ] α u - v , u , v X for some constant α < 1 . Then A has a unique fixed point. Furthermore, the sequence (6) u n + 1 = A [ u n ] , with an arbitrary choice of u 0 X converges to the fixed point of A .

According to Theorem 1, for the nonlinear mapping (7) A [ u n ( x , t ) ] = u n ( x , t ) + 0 t λ ( τ ) { L u n ( x , τ ) + N ( x , τ ) - f ( x , τ ) } d τ a sufficient condition for the convergence of the variational iteration method is strict contraction of A . Furthermore, the sequence (2) converges to the fixed point of A which is also the solution of problem (1). Some modifications to prove the convergence speed and to lengthen the interval of convergence for VIM series solution are suggested in [17, 2630].

3. The Variational Iteration of Generalized Degasperis-Procesi Equation

Degasperis and Procesi consider the following family of third-order dispersive conservation laws , (8) u t + c 0 u x + γ u x x x - α 2 u x x t = ( c 1 u 2 + c 2 u x 2 + c 3 u u x x ) x , where α , γ , c 0 , c 1 , c 2 , and c 3 are real constants. In this family, only three equations satisfy asymptotic integrability conditions . That is, if c 0 = 1 , c 1 = - 1 / 2 , c 2 = 0 , c 3 = 0 , α 2 = 0 , and γ = 1 , (8) is the KdV equation (9) u t + u x + u u x + u x x x = 0 . If c 0 = 0 , c 1 = - 3 / 2 , c 2 = 1 / 2 , c 3 = 1 , α 2 = 1 , and γ = 0 , (8) is the Camassa-Holm equation (10) u t - u x x t + 3 u u x = 2 u x u x x + u u x x x . If c 0 = 0 , c 1 = - 2 , c 2 = 1 , c 3 = 1 , α 2 = 1 , and γ = 0 , (8) is the Degasperis-Procesi equation (11) u t - u x x t + 4 u u x = 3 u x u x x + u u x x x . It should be mentioned that both C-H and D-P equations are derived as members of a one-parameter family of asymptotic shallow water approximations to the Euler equations. It shows that the two equations are physically relevant; otherwise, the D-P equation would be of purely theoretical interest.

Variational iteration method for KdV-Burgers and Lax’s seventh-order KdV equations has been studied by Soliman . In this paper, we consider the generalized Degasperis-Procesi equation (12) u t - u x x t + 4 u u x - 3 u u x x - u u x x x = f ( u , u x , u t , u x x , u x x x , u x x t ) that was proposed in . f is the generalized perturbation item. We suppose f is a sufficiently smooth function of the variable.

Step 1.

Make the independent variable transformation: (13) ξ = k ( x - ω t ) + ξ 0 . Here, ξ 0 C is an arbitrary complex number. k is wave number; ω is wave velocity. Substituting (13) into (12), we have (14) - ω u + k 2 ω u ′′′ + 4 u u - 3 k u u ′′ - k 2 u u ′′′ = f 1 . Here, u is the derivative of u with respect to ξ ; that is, u = d u / d ξ . f 1 = f 1 ( u , u , u ′′ , u ′′′ ) .

Step 2.

From , we find the special solution, when f is identical to 0 : (15) u 0 ( ξ ) = - 1 2 + 1 2 tanh 2 [ 1 2 ( x - ω t ) + ξ 0 ] .

Remark 2.

Notice that ξ 0 C , i tanh ( i ξ ) = - tan ξ , tanh ( ξ + ( π i / 2 ) ) = coth ξ , i coth ( i ξ ) = cot ξ , and tanh [ ( 1 / 2 ) ( ξ + ( i / 2 ) π ) ] = tanh ξ + i sech ξ , where i = - 1 . These solutions contain the other four types of forms named coth ξ , tan ξ , cot ξ , and tanh ξ + i sech ξ .

Step 3.

Make the correction functional (16) u n + 1 ( ξ )    = u n ( ξ ) - 0 ξ λ ( s ) [ - ω u ~ n + k 2 ω u n ′′′ + 4 u ~ n u ~ n - 3 k u ~ n u ~ n - k 2 u ~ n u ~ n - f ~ 1 ] d s .

Here, u ~ n , u ~ n , u ~ n , and u ~ n are considered as a restricted variation . That is, (17) δ u ~ n = δ u ~ n = δ u ~ n ′′ = δ u ~ n ′′′ = 0 .

Step 4.

Under the above condition, make the correct functional stationary with respect to u n ; noticing that δ u n ( 0 ) = 0 , we have (18) δ u n + 1 ( ξ )    = δ u n ( ξ ) - δ 0 ξ λ ( s ) [ - ω u ~ n + k 2 ω u n ′′′ + 4 u ~ n u ~ n - 3 k 2 u ~ n u ~ n - k 2 u ~ n u ~ n - f ~ ] d s = δ u n ( ξ ) - [ 0 ξ λ δ u ′′ ( s ) | s = ξ - λ ( s ) δ u ( s ) | s = ξ + λ 2 ( s ) δ u ( s ) | s = ξ - 0 ξ λ ′′′ ( s ) δ u d s ] = 0 .

For arbitrary δ u n + 1 , from the above relation, we obtain the Euler-Language equation: (19) 1 - k 2 ω λ ′′ ( s ) | s = ξ = 0 , k 2 ω λ ′′′ ( s ) | s = ξ = 0 , k 2 ω λ ( s ) | s = ξ = 0 , k 2 ω λ ( s ) | s = ξ = 0 .

Solve (19), we derive (20) λ ( s ) = - 1 2 · 1 k 2 ω ( s - ξ ) 2 .

Substituted (20) into (16), we have the integration form: (21) u n + 1 = u n + 0 ξ 1 2 k 2 ω ( s - ξ ) 2 × [ k 2 ω u n - ω u n + 4 u n u n - 3 k u n u n ′′ - k 2 u n u n ′′′ - f 1 ] d s , n = 0,1 , 2 , .

From the above solution procedure, we can see that the approximate solutions converge to its exact solution. That is, u ( x , t ) = lim n u n ( x , t ) , u n ( x , t ) is the approximate solution with arbitrary degree of accurate solitary wave of Degasperis-Procesi equation.

Step 5.

Calculation of the approximate solution.

According to the integration form (21),we can calculate the approximate solution. Firstly, let (15) be the zero-order approximate solution: (22) u 0 ( ξ ) = - 1 2 + 1 2 tanh 2 [ 1 2 ( x - ω t ) + ξ 0 ] . Substitute (15) into (21). We obtain the one-order approximate solution u 1 ( ξ ) (23) u 1 ( ξ ) = u 0 ( ξ ) + 0 ξ 1 2 k 2 ω ( s - ξ ) 2 ( - f 1 ( u 0 ) ) d s = - 1 2 + 1 2 tanh 2 [ 1 2 ( x - ω t ) + ξ 0 ] - 0 ξ 1 2 k 2 ω ( s - ξ ) 2 f 1 ( u 0 ( s ) , u 0 ( s ) , u 0 ′′ ( s ) , u 0 ′′′ ( s ) ) d s = - 1 2 + 1 2 tanh 2 [ 1 2 ( x - ω t ) + ξ 0 ] + v 0 ( ξ ) , in which v 0 ( ξ ) = - 0 ξ ( 1 / 2 k 2 ω ) ( s - ξ ) 2 f 1 ( u 0 ( s ) , u 0 ( s ) , u 0 ′′ ( s ) , u 0 ′′′ ( s ) ) d s .

Then, substitute (23) into (21). We can obtain the second-order approximate solution u 2 ( ξ ) : (24) u 2 ( ξ ) = u 1 ( ξ ) + 0 ξ 1 2 k 2 ω ( s - ξ ) 2 × [ k 2 ω u 1 ′′′ ( s ) - ω u 1 ( s ) + 4 u 1 ( s ) u 1 ( s ) - 3 k 2 u 1 ( s ) u 1 ′′ ( s ) - k 2 u 1 ( s ) u 1 ′′′ ( s ) - f 1 ( u 1 ( s ) , u 1 ( s ) , u 1 ′′ ( s ) , u 1 ′′′ ( s ) ) ] d s = - 1 2 + 1 2 tanh 2 [ 1 2 ( x - ω t ) + ξ 0 ] + v 0 ( ξ ) + 0 ξ 1 2 k 2 ω ( s - ξ ) 2 × [ k 2 ω ( u 0 ( s ) + v 0 ( s ) ) - ω ( u 0 ( s ) + v 0 ( s ) ) + 4 ( u 0 ( s ) + v 0 ( s ) ) × ( u 0 ( s ) + v 0 ( s ) ) - 3 k 2 ( u 0 ( s ) + v 0 ( s ) ) ( u 0 ( s ) + v 0 ( s ) ) - k 2 ( u 0 ( s ) + v 0 ( s ) ) ( u 0 ( s ) + v 0 ( s ) ) ′′′ - f 1 ( u 0 ( s ) + v 0 ( s ) , u 0 ( s ) + v 0 ( s ) , u 0 ′′ ( s ) + v 0 ′′ ( s ) , u 0 ′′′ ( s ) + v 0 ′′′ ( s ) ) ] d s .

Using the same method, we can get the higher order approximate solution.

4. The Optical Soliton Perturbation Solution and the Numerical Example

Specially, we set (25) f = ε g ( u ) = ε u , 0 < ε 1 . Then (12) change to (26) u t - u x x t + 4 u u x - 3 u u x x - u u x x x = ε g ( u ) = ε u . We obtain the zero-order and the one-order approximate solution of (20) (27) u 0 ( ξ ) = - 1 2 + 1 2 tanh 2 [ 1 2 ( x - ω t ) + ξ 0 ] = - 1 2 + 1 2 tanh 2 ( ξ ) , k = 1 2 , u 1 ( ξ ) = u 0 ( ξ ) + 0 ξ 1 2 k 2 ω ( s - ξ ) 2 ( - ε u ( s ) ) d s    = - 1 2 + 1 2 tanh 2 ξ - ε 2 k 2 ω 0 ξ s 2 ( - 1 2 + 1 2 tanh 2 s ) d s + ε ξ 2 k 2 ω 0 ξ s ( - 1 2 + 1 2 tanh 2 s ) d s - ε ξ 2 2 k 2 ω 0 ξ ( - 1 2 + 1 2 tanh 2 s ) d s    = - 1 2 + 1 2 tanh 2 ξ + ε 2 k 2 ω 0 ξ ln cosh s d s = - 1 2 + 1 2 tanh 2 [ 1 2 ( x - ω t ) + ξ 0 ] + 2 ε ω 0 ξ ln cosh s d s , k = 1 2 .

We set (28) ξ 0 = 0 , ω = 1 , t = 0 , ε = 0.1 , ξ 0 = 0 , ω = 1 , t = 0 , ε = 0.05 ; we will obtain the following numerical example. See Tables 1 and 2.

Numerical example for the solution of DP equation, ξ 0 = 0 , ω = 1 , t = 0 , and ε = 0.1 .

x Initial solution u 0 1th approximate solution u 1 Absolute error
0.1 −0.498752 −0.498762 - 4.16770 * 1 0 - 6
0.11 −0.498491 −0.498497 5.54752 * 1 0 - 6
0.12 −0.498204 −0.498211 7.2026 * 1 0 - 6
0.13 −0.497893 −0.497902 9.15804 * 1 0 - 6
0.14 −0.497558 −0.497694 1.53894 * 1 0 - 5
0.15 −0.497198 −0.497212 1.407042 * 1 0 - 5
0.16 −0.496814 −0.496831 1.706776 * 1 0 - 5
0.17 −0.496405 −0.496426 2.04856 * 1 0 - 5
0.18 −0.495972 −0.495996 2.43198 * 1 0 - 5
0.19 −0.495515 −0.495543 2.86050 * 1 0 - 5
0.2 −0.495033 −0.495066 3.33668 * 1 0 - 5

Numerical example for the solution of DP equation, ξ 0 = 0 , ω = 1 , t = 0 , and ε = 0.05 .

x Initial solution u 0 1th approximate solution u 1 Absolute error
0.1 −0.498752 −0.4987524 2.08385 * 1 0 - 6
0.11 −0.498491 −0.498493 2.77376 * 1 0 - 6
0.12 −0.498204 −0.4982076 3.6013 * 1 0 - 6
0.13 −0.497893 −0.4978976 4.57902 * 1 0 - 6
0.14 −0.497558 −0.4975637 5.71902 * 1 0 - 6
0.15 −0.497198 −0.497205 7.03521 * 1 0 - 6
0.16 −0.496814 −0.4968225 8.5388 * 1 0 - 6
0.17 −0.496405 −0.496415 1.02428 * 1 0 - 5
0.18 −0.495972 −0.4959842 1.21599 * 1 0 - 5
0.19 −0.495515 −0.495529 1.43025 * 1 0 - 5
0.2 −0.495033 −0.495049 1.66834 * 1 0 - 5

5. Conclusion

By the analysis of structure on the left side of (8) and the properties of f about the variable and the analytical variational iteration formula, we can prove that the sequence of functions of { u n ( ξ ) } decided by (21) is uniform convergence. So the limit function of { u n ( ξ ) } is the solution of the equation. Moreover, the zero-order approximate solution u 0 ( ξ ) is the soliton of (8) in which f = 0 ; it should be specially pointed out that the more accurate the identification of the multiplier is, the faster the approximations converge to their exact solution.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

The work is supported by the Youth Foundation of National Natural Science Foundation of China (Grant no. 71101072).

McPhaden M. J. Zhang D. Slowdown of the meridional overturning circulation in the upper Pacific Ocean Nature 2002 415 6872 603 608 2-s2.0-0037034026 10.1038/415603a Gu D. Philander S. G. H. Interdecadal climate fluctuations that depend on exchanges between the tropics and extratropics Science 1997 275 5301 805 807 2-s2.0-0031037673 10.1126/science.275.5301.805 Liu S. K. Fu Z. T. Liu S. D. The envelope periodic nonlinear wave equations with Jacobi Elliptic function solutions Acta Physica Sinica 2002 51 1 10 14 MR1894870 ZBL1202.35134 Malfliet W. Solitary wave solutions of nonlinear wave equations American Journal of Physics 1992 60 7 650 654 10.1119/1.17120 MR1174588 ZBL1219.35246 Fan E. Zhang H. A note on the homogeneous balance method Physics Letters A 1998 246 5 403 406 2-s2.0-4344673183 Rosenau P. Hyman J. M. Compactons: solitons with finite wavelength Physical Review Letters 1993 70 5 564 567 2-s2.0-12044253199 10.1103/PhysRevLett.70.564 He J. H. Approximate solution of non linear differential equations with convolution product nonlinearities Computer Methods in Applied Mechanics and Engineering 1998 167 1-2 69 73 10.1016/S0045-7825(98)00109-1 He J.-H. Approximate analytical solution for seepage flow with fractional derivatives in porous media Computer Methods in Applied Mechanics and Engineering 1998 167 1-2 57 68 10.1016/S0045-7825(98)00108-X MR1665221 ZBL0942.76077 He J.-H. Variational iteration method for autonomous ordinary differential systems Applied Mathematics and Computation 2000 114 2-3 115 123 10.1016/S0096-3003(99)00104-6 MR1778658 ZBL1027.34009 He J.-H. Variational principles for some nonlinear partial differential equations with variable coefficients Chaos, Solitons & Fractals 2004 19 4 847 851 10.1016/S0960-0779(03)00265-0 MR2009670 ZBL1135.35303 He J.-H. Some asymptotic methods for strongly nonlinear equations International Journal of Modern Physics B 2006 20 10 1141 1199 10.1142/S0217979206033796 MR2220565 ZBL1102.34039 He J.-H. Wu X.-H. Variational iteration method: new development and applications Computers & Mathematics with Applications 2007 54 7-8 881 894 10.1016/j.camwa.2006.12.083 MR2395625 ZBL1141.65372 He J.-H. Variational iteration method—some recent results and new interpretations Journal of Computational and Applied Mathematics 2007 207 1 3 17 10.1016/j.cam.2006.07.009 MR2332941 ZBL1119.65049 Wazwaz A.-M. The variational iteration method: a reliable analytic tool for solving linear and nonlinear wave equations Computers & Mathematics with Applications 2007 54 7-8 926 932 10.1016/j.camwa.2006.12.038 MR2395630 ZBL1141.65388 Khojasteh Salkuyeh D. Convergence of the variational iteration method for solving linear systems of ODEs with constant coefficients Computers & Mathematics with Applications 2008 56 8 2027 2033 10.1016/j.camwa.2008.03.030 MR2466705 ZBL1165.65376 Momani S. Abuasad S. Application of He's variational iteration method to Helmholtz equation Chaos, Solitons & Fractals 2006 27 5 1119 1123 10.1016/j.chaos.2005.04.113 MR2164845 ZBL1086.65113 Geng F. A modified variational iteration method for solving Riccati differential equations Computers & Mathematics with Applications 2010 60 7 1868 1872 10.1016/j.camwa.2010.07.017 MR2719705 ZBL1205.65229 Das S. Solution of fractional vibration equation by the variational iteration method and modified decomposition method International Journal of Nonlinear Sciences and Numerical Simulation 2008 9 4 361 366 2-s2.0-51749111612 Xu L. The variational iteration method for fourth order boundary value problems Chaos, Solitons and Fractals 2009 39 3 1386 1394 2-s2.0-62949130659 10.1016/j.chaos.2007.06.013 Noor M. A. Mohyud-Din S. T. Variational iteration method for solving higher-order nonlinear boundary value problems using He's polynomials International Journal of Nonlinear Sciences and Numerical Simulation 2007 9 2 141 156 10.1016/j.amc.2006.12.071 MR2332146 ZBL1122.65374 Noor M. A. Noor K. I. Mohyud-Din S. T. Modified variational iteration technique for solving singular fourth-order parabolic partial differential equations Nonlinear Analysis: Theory, Methods & Applications 2009 71 12 e630 e640 10.1016/j.na.2008.11.011 MR2671865 ZBL1238.65103 Ghorbani A. Saberi-Nadjafi J. An effective modification of He's variational iteration method Nonlinear Analysis: Real World Applications 2009 10 5 2828 2833 10.1016/j.nonrwa.2008.08.008 MR2523245 ZBL1168.45301 Sevimlican A. Constructing approximate Green's function for a vector equation for the electric field using the variational iteration method Applied Mathematics Letters 2010 23 5 533 536 10.1016/j.aml.2010.01.006 MR2602404 ZBL1189.35323 Inokuti M. Sekine H. Mura T. Nemat-Nasser S. General use of the Lagrange multiplier in nonlinear mathematical physics Variational Method in the Mechanics of Solids 1978 NewYork, NY, USA Pergamon Press 159 162 He J.-H. Variational iteration method—a kind of non-linear analytical technique: some examples International Journal of Non-Linear Mechanics 1999 34 4 699 708 2-s2.0-0000092673 Wazwaz A.-M. A comparison between the variational iteration method and Adomian decomposition method Journal of Computational and Applied Mathematics 2007 207 1 129 136 10.1016/j.cam.2006.07.018 MR2332954 ZBL1119.65103 Abassy T. A. El-Tawil M. A. El Zoheiry H. Solving nonlinear partial differential equations using the modified variational iteration Padé technique Journal of Computational and Applied Mathematics 2007 207 1 73 91 10.1016/j.cam.2006.07.024 MR2332949 ZBL1119.65095 Odibat Z. M. Reliable approaches of variational iteration method for nonlinear operators Mathematical and Computer Modelling 2008 48 1-2 222 231 10.1016/j.mcm.2007.09.005 MR2431335 ZBL1145.65314 Tatari M. Dehghan M. Improvement of He's variational iteration method for solving systems of differential equations Computers & Mathematics with Applications 2009 58 11-12 2160 2166 10.1016/j.camwa.2009.03.081 MR2557343 ZBL1189.65178 Odibat Z. M. A study on the convergence of variational iteration method Mathematical and Computer Modelling 2010 51 9-10 1181 1192 10.1016/j.mcm.2009.12.034 MR2608904 ZBL1198.65147 Degasperis A. Procesi M. Asymoptotic Integrability, Symmetry and Perturbation Theory 1999 River Edge, NJ, USA World Scientific Soliman A. A. A numerical simulation and explicit solutions of KdV-Burgers' and Lax's seventh-order KdV equations Chaos, Solitons and Fractals 2006 29 2 294 302 10.1016/j.chaos.2005.08.054 MR2211466 ZBL1099.35521 Camassa R. Holm D. D. An integrable shallow water equation with peaked solitons Physical Review Letters 1993 71 11 1661 1664 10.1103/PhysRevLett.71.1661 MR1234453 ZBL0972.35521 Yu L. Q. Tian L. X. A new traveling wave solutions for the Degasperis-Procesi equation Journal of Jiangsu University 2005 26 3 231 234 MR2139822 He J. H. Approximate Nonlinear Analytical Methods in Engineering and Sciences 2002 Zhengzhou, China Henan Science and Technology Press (Chinese)