Making use of a new generalized ans?tze and a proper transformation, we generalized the extended tanh-function method. Applying the generalized method with the aid of Maple, we consider some nonlinear evolution equati...Making use of a new generalized ans?tze and a proper transformation, we generalized the extended tanh-function method. Applying the generalized method with the aid of Maple, we consider some nonlinear evolution equations. As a result, we can successfully recover the previously known solitary wave solutions that had been found by the extended tanh-function method and other more sophisticated methods. More importantly, for some equations, we also obtain other new and more general solutions at the same time. The results include kink-profile solitary-wave solutions, bell-profile solitary-wave solutions, periodic wave solutions, rational solutions, singular solutions and new formal solutions.展开更多
To seek new infinite sequence soliton-like exact solutions to nonlinear evolution equations (NEE(s)), by developing two characteristics of construction and mechanization on auxiliary equation method, the second ki...To seek new infinite sequence soliton-like exact solutions to nonlinear evolution equations (NEE(s)), by developing two characteristics of construction and mechanization on auxiliary equation method, the second kind of elliptie equation is highly studied and new type solutions and Backlund transformation are obtained. Then (2+ l )-dimensional breaking soliton equation is chosen as an example and its infinite sequence soliton-like exact solutions are constructed with the help of symbolic computation system Mathematica, which include infinite sequence smooth soliton-like solutions of Jacobi elliptic type, infinite sequence compact soliton solutions of Jacobi elliptic type and infinite sequence peak soliton solutions of exponential function type and triangular function type.展开更多
Applying the generalized method, which is a direct and unified algebraic method for constructing multipletravelling wave solutions of nonlinear partial differential equations (PDEs), and implementing in a computer alg...Applying the generalized method, which is a direct and unified algebraic method for constructing multipletravelling wave solutions of nonlinear partial differential equations (PDEs), and implementing in a computer algebraicsystem, we consider the generalized Zakharov-Kuzentsov equation with nonlinear terms of any order. As a result, wecan not only successfully recover the previously known travelling wave solutions found by existing various tanh methodsand other sophisticated methods, but also obtain some new formal solutions. The solutions obtained include kink-shapedsolitons, bell-shaped solitons, singular solitons, and periodic solutions.展开更多
The impacts of initial perturbations on the computational stability of nonlinear evolution equations for non-conservative difference schemes and non-periodic boundary conditions are studied through theoretical analysi...The impacts of initial perturbations on the computational stability of nonlinear evolution equations for non-conservative difference schemes and non-periodic boundary conditions are studied through theoretical analysis and numerical experiments for the case of onedimensional equations.The sensitivity of the difference scheme to initial values is further analyzed.The results show that the computational stability primarily depends on the form of the initial values if the difference scheme and boundary conditions are determined.Thus,the computational stability is sensitive to the initial perturbations.展开更多
With the aid of symbolic computation system Mathematica, several explicit solutions for Fisher's equation and CKdV equation are constructed by utilizing an auxiliary equation method, the so called G′/G-expansion met...With the aid of symbolic computation system Mathematica, several explicit solutions for Fisher's equation and CKdV equation are constructed by utilizing an auxiliary equation method, the so called G′/G-expansion method, where the new and more general forms of solutions are also constructed. When the parameters are taken as special values, the previously known solutions are recovered.展开更多
With symbolic computation, the Hirota method and Riemann theta function are employed to directly construct the periodic wave solutions for the Hirota-Satsuma equation for shallow water waves and Boiti-Leon-Manna- Pemp...With symbolic computation, the Hirota method and Riemann theta function are employed to directly construct the periodic wave solutions for the Hirota-Satsuma equation for shallow water waves and Boiti-Leon-Manna- Pempinelli equation. Then, the corresponding figures of the periodic wave solutions are given. Fhrthermore, it is shown that the known soliton solutions can be reduced from the periodic wave solutions.展开更多
In this paper, based on a new more general ansitz, a new algebraic method, named generalized Riccati equation rational expansion method, is devised for constructing travelling wave solutions for nonlinear evolution eq...In this paper, based on a new more general ansitz, a new algebraic method, named generalized Riccati equation rational expansion method, is devised for constructing travelling wave solutions for nonlinear evolution equations with nonlinear terms of any order. Compared with most existing tanh methods for finding travelling wave solutions, the proposed method not only recovers the results by most known algebraic methods, but also provides new and more general solutions. We choose the generalized Burgers-Fisher equation with nonlinear terms of any order to illustrate our method. As a result, we obtain several new kinds of exact solutions for the equation. This approach can also be applied to other nonlinear evolution equations with nonlinear terms of any order.展开更多
In the paper, a generalized sub-equation method is presented to construct some exact analytical solutions of nonlinear partial differential equations. Making use of the method, we present rich exact analytical solutio...In the paper, a generalized sub-equation method is presented to construct some exact analytical solutions of nonlinear partial differential equations. Making use of the method, we present rich exact analytical solutions of the onedimensional nonlinear Schrfdinger equation which describes the dynamics of solitons in Bose-Einstein condensates with time-dependent atomic scattering length in an expulsive parabolic potential. The solutions obtained include not only non-traveling wave and coefficient function's soliton solutions, but also Jacobi elliptic function solutions and Weierstra.ss elliptic function solutions. Some plots are given to demonstrate the properties of some exact solutions under the Feshbachmanaged nonlinear coefficient and the hyperbolic secant function coefficient.展开更多
In this paper, the Adomian decomposition method is developed for the numerical solutions of a class of nonlinear evolution equations with nonlinear term of any order, utt+auxx + bu + cu^p+ du^2p-1=0, which contain...In this paper, the Adomian decomposition method is developed for the numerical solutions of a class of nonlinear evolution equations with nonlinear term of any order, utt+auxx + bu + cu^p+ du^2p-1=0, which contains some important famous equations. When setting the initial conditions in different forms, some new generalized numerical solutions: numerical hyperbolic solutions, numerical doubly periodic solutions are obtained. The numerical solutions are compared with exact solutions. The scheme is tested by choosing different values of p, positive and negative, integer and fraction, to illustrate the efficiency of the ADM method and the generalization of the solutions.展开更多
Based on the Reynolds equation with Reynolds boundary conditions, the Castelli method was employed to solve the Reynolds equation for oil lubrication upon bearings. By doing so, a profile of nonlinear oil film force o...Based on the Reynolds equation with Reynolds boundary conditions, the Castelli method was employed to solve the Reynolds equation for oil lubrication upon bearings. By doing so, a profile of nonlinear oil film force of single-pad journal bearings is established. According to the structure of combination journal bearings, nonlinear oil film force of combination journal bearing is obtained by retrieval, interpolation and assembly techniques. As for symmetrical flexible Jeffcott rotor systems supported by combination journal bearings, the nonlinear motions of the center of the rotor are calculated by the self-adaptive Runge-Kutta method and Poincar6 mapping with different rotational speeds. The numerical results show that the system performance is slightly better when the pivot ratio changes from 0.5 to 0.6, and reveals nonlinear phenomena of periodic, period-doubing, quasi-periodic motion, etc.展开更多
The existence and representation of the exact solution are given for a nonlinear functional equation in the reproducing kernel space. For a numerical computation, we present a large-range convergence iterative method ...The existence and representation of the exact solution are given for a nonlinear functional equation in the reproducing kernel space. For a numerical computation, we present a large-range convergence iterative method for solving the nonlinear functional equation. In the iterative method, the convergent condition is simple and the convergence is irrespective to the choice of the initial function. It is worthy to note that the presented method can be generalized to solve other nonlinear operator equations.展开更多
An improved algorithm for symbolic computation of Hirota bilinear forms of KdV-type equations withlogarithmic transformations is presented.In the algorithm,the general assumption of Hirota bilinear form is successfull...An improved algorithm for symbolic computation of Hirota bilinear forms of KdV-type equations withlogarithmic transformations is presented.In the algorithm,the general assumption of Hirota bilinear form is successfullyreduced based on the property of uniformity in rank.Furthermore,we discard the integral operation in the traditionalalgorithm.The software package HBFTrans is written in Maple and its running effectiveness is tested by a variety solitonequations.展开更多
Combining Adomian decomposition method (ADM) with Pade approximants, we solve two differentiaidifference equations (DDEs): the relativistic Toda lattice equation and the modified Volterra lattice equation. With t...Combining Adomian decomposition method (ADM) with Pade approximants, we solve two differentiaidifference equations (DDEs): the relativistic Toda lattice equation and the modified Volterra lattice equation. With the help of symbolic computation Maple, the results obtained by ADM-Pade technique are compared with those obtained by using ADM alone. The numerical results demonstrate that ADM-Pade technique give the approximate solution with faster convergence rate and higher accuracy and relative in larger domain of convergence than using ADM.展开更多
For a specific combustion problem involving calculations of several species at the equilibrium state, it is simpler to write a general computer program and calculate the combustion concentration. Original work describ...For a specific combustion problem involving calculations of several species at the equilibrium state, it is simpler to write a general computer program and calculate the combustion concentration. Original work describes, an adaptation of Newton-Raphson method was used for solving the highly nonlinear system of equations describing the formation of equilibrium products in reacting of fuel-additive-air mixtures. This study also shows what possible of the results. In this paper, to be present the efficient numerical algorithms for. solving the combustion problem, to be used nonlinear equations based on the iteration method and high order of the Taylor series. The modified Adomian decomposition method was applied to construct the numerical algorithms. Some numerical illustrations are given to show the efficiency of algorithms. Comparisons of results by the new Matlab routines and previous routines, the result data indicate that the new Matlab routines are reliable, typical deviations from previous results are less than 0.05%.展开更多
Using the extended homogeneous balance method, we obtained abundant exact solution structures ofthe (3 + 1)-dimensional Nizhnik-Novikov-Veselov (NNV) equation. By means of the leading order term analysis, thenonlinear...Using the extended homogeneous balance method, we obtained abundant exact solution structures ofthe (3 + 1)-dimensional Nizhnik-Novikov-Veselov (NNV) equation. By means of the leading order term analysis, thenonlinear transformations of the (3+1)-dimensional NNV equation are given first, and then some special types of singlesolitary wave solution and the multisoliton solutions are constructed.展开更多
Using difference quotient instead of derivative, the paper presents the solution method and procedure of the nonlinear least square estimation containing different classes of measurements. In the meantime, the paper s...Using difference quotient instead of derivative, the paper presents the solution method and procedure of the nonlinear least square estimation containing different classes of measurements. In the meantime, the paper shows several practical cases, which indicate the method is very valid and reliable.展开更多
With the aid of computation, we consider the variable-coefficient coupled nonlinear Schrodinger equations with the effects of group-velocity dispersion, self-phase modulation and cross-phase modulation, which have pot...With the aid of computation, we consider the variable-coefficient coupled nonlinear Schrodinger equations with the effects of group-velocity dispersion, self-phase modulation and cross-phase modulation, which have potential applications in the long-distance communication of two-pulse propagation in inhomogeneous optical fibers. Based on the obtained nonisospectral linear eigenvalue problems (i.e. Lax pair), we construct the Darboux transformation for such a model to derive the optical soliton solutions. In addition, through the one- and two-soliton-like solutions, we graphically discuss the features of picosecond solitons in inhomogeneous optical fibers.展开更多
An improved algorithm for symbolic computation of Hirota bilinear form of nonlinear equations by a logarithm transformation is presented. The improved algorithm is more efficient by using the property of Hirota-D oper...An improved algorithm for symbolic computation of Hirota bilinear form of nonlinear equations by a logarithm transformation is presented. The improved algorithm is more efficient by using the property of Hirota-D operator. The software package HBFTrans2 is written in Maple and its running efficiency is tested by a variety of soliton equations.展开更多
Investigated in this paper is the generalized nonlinear Schrodinger equation with radial symmetry. With the help of symbolic computation, the one-, two-, and N-soliton solutions are obtained through the bilinear metho...Investigated in this paper is the generalized nonlinear Schrodinger equation with radial symmetry. With the help of symbolic computation, the one-, two-, and N-soliton solutions are obtained through the bilinear method. B^cklund transformation in the bilinear form is presented, through which a new solution is constructed. Graphically, we have found that the solitons are symmetric about x = O, while the soliton pulse width and amplitude will change along with the distance and time during the propagation.展开更多
文摘Making use of a new generalized ans?tze and a proper transformation, we generalized the extended tanh-function method. Applying the generalized method with the aid of Maple, we consider some nonlinear evolution equations. As a result, we can successfully recover the previously known solitary wave solutions that had been found by the extended tanh-function method and other more sophisticated methods. More importantly, for some equations, we also obtain other new and more general solutions at the same time. The results include kink-profile solitary-wave solutions, bell-profile solitary-wave solutions, periodic wave solutions, rational solutions, singular solutions and new formal solutions.
基金Supported by the Natural Natural Science Foundation of China under Grant No.10461006the Science Research Foundation of Institution of Higher Education of Inner Mongolia Autonomous Region,China under Grant No.NJZZ07031the Natural Science Foundation of Inner Mongolia Autonomous Region,China under Grant No.2010MS0111
文摘To seek new infinite sequence soliton-like exact solutions to nonlinear evolution equations (NEE(s)), by developing two characteristics of construction and mechanization on auxiliary equation method, the second kind of elliptie equation is highly studied and new type solutions and Backlund transformation are obtained. Then (2+ l )-dimensional breaking soliton equation is chosen as an example and its infinite sequence soliton-like exact solutions are constructed with the help of symbolic computation system Mathematica, which include infinite sequence smooth soliton-like solutions of Jacobi elliptic type, infinite sequence compact soliton solutions of Jacobi elliptic type and infinite sequence peak soliton solutions of exponential function type and triangular function type.
基金The project supported by National Natural Science Foundation of China under Grant No.10072013the National Key Basic Research Development Program under Grant No.G1998030600
文摘Applying the generalized method, which is a direct and unified algebraic method for constructing multipletravelling wave solutions of nonlinear partial differential equations (PDEs), and implementing in a computer algebraicsystem, we consider the generalized Zakharov-Kuzentsov equation with nonlinear terms of any order. As a result, wecan not only successfully recover the previously known travelling wave solutions found by existing various tanh methodsand other sophisticated methods, but also obtain some new formal solutions. The solutions obtained include kink-shapedsolitons, bell-shaped solitons, singular solitons, and periodic solutions.
基金supported by the"Strategic Priority Research Program-Climate Change:Carbon Budget and Relevant Issues"of the Chinese Academy of Sciences (Grant No.XDA01020304)
文摘The impacts of initial perturbations on the computational stability of nonlinear evolution equations for non-conservative difference schemes and non-periodic boundary conditions are studied through theoretical analysis and numerical experiments for the case of onedimensional equations.The sensitivity of the difference scheme to initial values is further analyzed.The results show that the computational stability primarily depends on the form of the initial values if the difference scheme and boundary conditions are determined.Thus,the computational stability is sensitive to the initial perturbations.
文摘With the aid of symbolic computation system Mathematica, several explicit solutions for Fisher's equation and CKdV equation are constructed by utilizing an auxiliary equation method, the so called G′/G-expansion method, where the new and more general forms of solutions are also constructed. When the parameters are taken as special values, the previously known solutions are recovered.
基金Supported by the National Natural Science Foundation of China under Grant No.60772023the Open Fund of the State Key Laboratory of Software Development Environment under Grant No.BUAA-SKLSDE-09KF-04+1 种基金Beijing University of Aeronautics and Astronautics,by the National Basic Research Program of China (973 Program) under Grant No.2005CB321901the Specialized Research Fund for the Doctoral Program of Higher Education under Grant Nos.20060006024 and 200800130006,Chinese Ministry of Education
文摘With symbolic computation, the Hirota method and Riemann theta function are employed to directly construct the periodic wave solutions for the Hirota-Satsuma equation for shallow water waves and Boiti-Leon-Manna- Pempinelli equation. Then, the corresponding figures of the periodic wave solutions are given. Fhrthermore, it is shown that the known soliton solutions can be reduced from the periodic wave solutions.
基金The project partially supported by the State Key Basic Research Program of China under Grant No. 2004CB318000
文摘In this paper, based on a new more general ansitz, a new algebraic method, named generalized Riccati equation rational expansion method, is devised for constructing travelling wave solutions for nonlinear evolution equations with nonlinear terms of any order. Compared with most existing tanh methods for finding travelling wave solutions, the proposed method not only recovers the results by most known algebraic methods, but also provides new and more general solutions. We choose the generalized Burgers-Fisher equation with nonlinear terms of any order to illustrate our method. As a result, we obtain several new kinds of exact solutions for the equation. This approach can also be applied to other nonlinear evolution equations with nonlinear terms of any order.
基金The project supported by the Natural Science Foundation of Zhejiang Province of China under Grant Nos. 605408 and Y604056, the Doctoral Foundation of Ningbo City under Grant No. 2005A61030, and the Postdoctoral Science Foundation of China under Grant No. 2005038441
文摘In the paper, a generalized sub-equation method is presented to construct some exact analytical solutions of nonlinear partial differential equations. Making use of the method, we present rich exact analytical solutions of the onedimensional nonlinear Schrfdinger equation which describes the dynamics of solitons in Bose-Einstein condensates with time-dependent atomic scattering length in an expulsive parabolic potential. The solutions obtained include not only non-traveling wave and coefficient function's soliton solutions, but also Jacobi elliptic function solutions and Weierstra.ss elliptic function solutions. Some plots are given to demonstrate the properties of some exact solutions under the Feshbachmanaged nonlinear coefficient and the hyperbolic secant function coefficient.
基金supported by National Natural Science Foundation of China under Grant No.10735030Natural Science Foundation of Zhejiang Province of China under Grant No.Y604056Doctoral Science Foundation of Ningbo City under Grant No.2005A61030
文摘In this paper, the Adomian decomposition method is developed for the numerical solutions of a class of nonlinear evolution equations with nonlinear term of any order, utt+auxx + bu + cu^p+ du^2p-1=0, which contains some important famous equations. When setting the initial conditions in different forms, some new generalized numerical solutions: numerical hyperbolic solutions, numerical doubly periodic solutions are obtained. The numerical solutions are compared with exact solutions. The scheme is tested by choosing different values of p, positive and negative, integer and fraction, to illustrate the efficiency of the ADM method and the generalization of the solutions.
基金Project(2007CB707706) supported by the National Basic Research Program of China Projects(51075327,10972179) supported by the National Natural Science Foundation of China+2 种基金 Project(SKLMT-KFKT-201011) supported by Open Foundation of State Key Laboratory of Mechanical Transmission,China Projects(2009JQ7006,2007E203) supported by the Natural Science Foundation of Shaanxi Province of China Projects(09JK680,07JK340) supported by the Natural Science Foundation of Department of Education of Shaanxi Province of China
文摘Based on the Reynolds equation with Reynolds boundary conditions, the Castelli method was employed to solve the Reynolds equation for oil lubrication upon bearings. By doing so, a profile of nonlinear oil film force of single-pad journal bearings is established. According to the structure of combination journal bearings, nonlinear oil film force of combination journal bearing is obtained by retrieval, interpolation and assembly techniques. As for symmetrical flexible Jeffcott rotor systems supported by combination journal bearings, the nonlinear motions of the center of the rotor are calculated by the self-adaptive Runge-Kutta method and Poincar6 mapping with different rotational speeds. The numerical results show that the system performance is slightly better when the pivot ratio changes from 0.5 to 0.6, and reveals nonlinear phenomena of periodic, period-doubing, quasi-periodic motion, etc.
基金Sponsored by the Education Department Science and Technology Foundation of Heilongjiang Province (Grant No.11531324)
文摘The existence and representation of the exact solution are given for a nonlinear functional equation in the reproducing kernel space. For a numerical computation, we present a large-range convergence iterative method for solving the nonlinear functional equation. In the iterative method, the convergent condition is simple and the convergence is irrespective to the choice of the initial function. It is worthy to note that the presented method can be generalized to solve other nonlinear operator equations.
基金Supported by Scientific Research Fund of Zhejiang Provincial Education Department under Grant No.20070979the National Natural Science Foundations of China under Grant Nos.10675065 and 10735030+1 种基金the Scientific Research Found of Ningbo University under Grant No.XKL09059the K.C.Wong Magana Fund in Ningbo University
文摘An improved algorithm for symbolic computation of Hirota bilinear forms of KdV-type equations withlogarithmic transformations is presented.In the algorithm,the general assumption of Hirota bilinear form is successfullyreduced based on the property of uniformity in rank.Furthermore,we discard the integral operation in the traditionalalgorithm.The software package HBFTrans is written in Maple and its running effectiveness is tested by a variety solitonequations.
基金supported by the National Natural Science Foundation of China under Grant No. 10735030Shanghai Leading Academic Discipline Project under Grant No. B412Program for Changjiang Scholars and Innovative Research Team in University under Grant No. IRT0734
文摘Combining Adomian decomposition method (ADM) with Pade approximants, we solve two differentiaidifference equations (DDEs): the relativistic Toda lattice equation and the modified Volterra lattice equation. With the help of symbolic computation Maple, the results obtained by ADM-Pade technique are compared with those obtained by using ADM alone. The numerical results demonstrate that ADM-Pade technique give the approximate solution with faster convergence rate and higher accuracy and relative in larger domain of convergence than using ADM.
文摘For a specific combustion problem involving calculations of several species at the equilibrium state, it is simpler to write a general computer program and calculate the combustion concentration. Original work describes, an adaptation of Newton-Raphson method was used for solving the highly nonlinear system of equations describing the formation of equilibrium products in reacting of fuel-additive-air mixtures. This study also shows what possible of the results. In this paper, to be present the efficient numerical algorithms for. solving the combustion problem, to be used nonlinear equations based on the iteration method and high order of the Taylor series. The modified Adomian decomposition method was applied to construct the numerical algorithms. Some numerical illustrations are given to show the efficiency of algorithms. Comparisons of results by the new Matlab routines and previous routines, the result data indicate that the new Matlab routines are reliable, typical deviations from previous results are less than 0.05%.
基金The project supported by National Natural Science Foundation of China and the Natural Science Foundation of Shandong Province
文摘Using the extended homogeneous balance method, we obtained abundant exact solution structures ofthe (3 + 1)-dimensional Nizhnik-Novikov-Veselov (NNV) equation. By means of the leading order term analysis, thenonlinear transformations of the (3+1)-dimensional NNV equation are given first, and then some special types of singlesolitary wave solution and the multisoliton solutions are constructed.
文摘Using difference quotient instead of derivative, the paper presents the solution method and procedure of the nonlinear least square estimation containing different classes of measurements. In the meantime, the paper shows several practical cases, which indicate the method is very valid and reliable.
基金Supported by the National Natural Science Foundation of China under Grant No.60772023 the Open Fund of the State Key Laboratory of Software Development Environment under Grant No.BUAA-SKLSDE-09KF-04+2 种基金Beijing University of Aeronautics and Astronautics,by the National Basic Research Program of China (973 Program) under Grant No.2005CB321901the Specialized Research Fund for the Doctoral Program of Higher Education under Grant Nos.20060006024 and 200800130006Chinese Ministry of Education
文摘With the aid of computation, we consider the variable-coefficient coupled nonlinear Schrodinger equations with the effects of group-velocity dispersion, self-phase modulation and cross-phase modulation, which have potential applications in the long-distance communication of two-pulse propagation in inhomogeneous optical fibers. Based on the obtained nonisospectral linear eigenvalue problems (i.e. Lax pair), we construct the Darboux transformation for such a model to derive the optical soliton solutions. In addition, through the one- and two-soliton-like solutions, we graphically discuss the features of picosecond solitons in inhomogeneous optical fibers.
基金Supported by Scientific Research Fund of Zhejiang Provincial Education Department under Grant No.Y201017148the National Natural Science Foundations of China under Grant No.10735030+2 种基金the Natural Science Fund of Ningbo under Grant No.2009B21003the Scientific Research Fund of Ningbo University under Grant No.XKL09059the K.C.Wong Magana Fund in Ningbo University
文摘An improved algorithm for symbolic computation of Hirota bilinear form of nonlinear equations by a logarithm transformation is presented. The improved algorithm is more efficient by using the property of Hirota-D operator. The software package HBFTrans2 is written in Maple and its running efficiency is tested by a variety of soliton equations.
基金Supported by the National Natural Science Foundation of China under Grant No.60772023the Open Fund (No.BUAASKLSDE-09KF-04)+2 种基金Supported Project (No.SKLSDE-2010ZX-07) of the State Key Laboratory of Software Development Environment,Beijing University of Aeronautics and Astronauticsthe National Basic Research Program of China (973 Program) under Grant No.2005CB321901the Specialized Research Fund for the Doctoral Program of Higher Education under Grant No.200800130006,Chinese Ministry of Education
文摘Investigated in this paper is the generalized nonlinear Schrodinger equation with radial symmetry. With the help of symbolic computation, the one-, two-, and N-soliton solutions are obtained through the bilinear method. B^cklund transformation in the bilinear form is presented, through which a new solution is constructed. Graphically, we have found that the solitons are symmetric about x = O, while the soliton pulse width and amplitude will change along with the distance and time during the propagation.