In this paper, we shall show that the Hamiltonian structure can be defined for any nonlinear evolution equations which describe surfaces of a constant negative curvature, so that the densities of conservation laws can...In this paper, we shall show that the Hamiltonian structure can be defined for any nonlinear evolution equations which describe surfaces of a constant negative curvature, so that the densities of conservation laws can be considered as corresponding Hamiltonians. This paper obtains the soliton solution and conserved quantities of a new fifth-order nonlinear evolution equation by the aid of inverse scattering method.展开更多
With the urgent need to resolve complex behaviors in nonlinear evolution equations,this study makes a contribution by establishing the local existence of solutions for Cauchy problems associated with equations of mixe...With the urgent need to resolve complex behaviors in nonlinear evolution equations,this study makes a contribution by establishing the local existence of solutions for Cauchy problems associated with equations of mixed types.Our primary contribution is the establishment of solution existence,illuminating the dynamics of these complex equations.To tackle this challenging problem,we construct an approximate solution sequence and apply the contraction mapping principle to rigorously prove local solution existence.Our results significantly advance the understanding of nonlinear evolution equations of mixed types.Furthermore,they provide a versatile,powerful approach for tackling analogous challenges across physics,engineering,and applied mathematics,making this work a valuable reference for researchers in these fields.展开更多
An evolution inequality of Sobolev type involving a nonlinear convolution term is considered.By using the nonlinear capacity method and the contradiction argument,the non-existence of the nontrivial local weak solutio...An evolution inequality of Sobolev type involving a nonlinear convolution term is considered.By using the nonlinear capacity method and the contradiction argument,the non-existence of the nontrivial local weak solution is proved.展开更多
The(2 + 1)-dimensional generalized fifth-order Kd V(2GKd V) equation is revisited via combined physical and mathematical methods. By using the Hirota perturbation expansion technique and via setting the nonzero backgr...The(2 + 1)-dimensional generalized fifth-order Kd V(2GKd V) equation is revisited via combined physical and mathematical methods. By using the Hirota perturbation expansion technique and via setting the nonzero background wave on the multiple soliton solution of the 2GKd V equation, breather waves are constructed, for which some transformed wave conditions are considered that yield abundant novel nonlinear waves including X/Y-Shaped(XS/YS),asymmetric M-Shaped(MS), W-Shaped(WS), Space-Curved(SC) and Oscillation M-Shaped(OMS) solitons. Furthermore, distinct nonlinear wave molecules and interactional structures involving the asymmetric MS, WS, XS/YS, SC solitons, and breathers, lumps are constructed after considering the corresponding existence conditions. The dynamical properties of the nonlinear molecular waves and interactional structures are revealed via analyzing the trajectory equations along with the change of the phase shifts.展开更多
This paper studies a generalized nonlinear evolution equation. Using the homotopic mapping method, it constructs a corresponding homotopic mapping transform. Selecting a suitable initial approximation and using homoto...This paper studies a generalized nonlinear evolution equation. Using the homotopic mapping method, it constructs a corresponding homotopic mapping transform. Selecting a suitable initial approximation and using homotopic mapping, it obtains an approximate solution with an arbitrary degree of accuracy for the solitary wave. From the approximate solution obtained by using the homotopic mapping method, it possesses a good accuracy.展开更多
The process of evolution, especially that of nonlinear evolution, of C-type instability of laminar-turbulent flow transition in nonparallel boundary layers are studied by means of a newly developed method called parab...The process of evolution, especially that of nonlinear evolution, of C-type instability of laminar-turbulent flow transition in nonparallel boundary layers are studied by means of a newly developed method called parabolic stability equations (PSE). Initial conditions, which are very important for the nonlinear problem, are investigated by computing initial solution of the harmonic waves, modifying the mean-flow-distortion, and giving initial value of TS wave and its subharmonic waves at initial station by solving linear PSE. A numerical method with high-order accuracy are developed in the text, the key normalization conditions in the PSE are satisfied, and nonlinear PSE are solved efficiently and implemented stably by the spatial marching. It has been shown that the computed process of nonlinear evolution of C-type instability in Blasius boundary layer is in good agreement with the experimental results.展开更多
A Riccati equation involving a parameter and symbolic computation are used to uniformly construct the different forms of travelling wave solutions for nonlinear evolution equations.It is shown that the sign of the pa...A Riccati equation involving a parameter and symbolic computation are used to uniformly construct the different forms of travelling wave solutions for nonlinear evolution equations.It is shown that the sign of the parameter can be applied in judging the existence of various forms of travelling wave solutions.An efficiency of this method is demonstrated on some equations,which include Burgers Huxley equation,Caudrey Dodd Gibbon Kawada equation,generalized Benjamin Bona Mahony equation and generalized Fisher equation.展开更多
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.展开更多
It has still been difficult to solve nonlinear evolution equations analytically.In this paper,we present a deep learning method for recovering the intrinsic nonlinear dynamics from spatiotemporal data directly.Specifi...It has still been difficult to solve nonlinear evolution equations analytically.In this paper,we present a deep learning method for recovering the intrinsic nonlinear dynamics from spatiotemporal data directly.Specifically,the model uses a deep neural network constrained with given governing equations to try to learn all optimal parameters.In particular,numerical experiments on several third-order nonlinear evolution equations,including the Korteweg-de Vries(KdV)equation,modified KdV equation,KdV-Burgers equation and Sharma-Tasso-Olver equation,demonstrate that the presented method is able to uncover the solitons and their interaction behaviors fairly well.展开更多
Solving nonlinear evolution partial differential equations has been a longstanding computational challenge.In this paper,we present a universal paradigm of learning the system and extracting patterns from data generat...Solving nonlinear evolution partial differential equations has been a longstanding computational challenge.In this paper,we present a universal paradigm of learning the system and extracting patterns from data generated from experiments.Specifically,this framework approximates the latent solution with a deep neural network,which is trained with the constraint of underlying physical laws usually expressed by some equations.In particular,we test the effectiveness of the approach for the Burgers'equation used as an example of second-order nonlinear evolution equations under different initial and boundary conditions.The results also indicate that for soliton solutions,the model training costs significantly less time than other initial conditions.展开更多
A trial equation method to nonlinear evolution equation with rank inhomogeneous is given. As appncations, the exact traveling wave solutions to some higher-order nonlinear equations such as generalized Boussinesq equa...A trial equation method to nonlinear evolution equation with rank inhomogeneous is given. As appncations, the exact traveling wave solutions to some higher-order nonlinear equations such as generalized Boussinesq equation, generalized Pochhammer-Chree equation, KdV-Burgers equation, and KS equation and so on, are obtained. Among these, some results are new. The proposed method is based on the idea of reduction of the order of ODE. Some mathematical details of the proposed method are discussed.展开更多
With the use of computer algebra, the method that straightforwardly leads to travelling wave solutions is presented. The compound KdV-Burgers equation and KP-B equation are chosen to illustrate this approach. As a res...With the use of computer algebra, the method that straightforwardly leads to travelling wave solutions is presented. The compound KdV-Burgers equation and KP-B equation are chosen to illustrate this approach. As a result, their abundant new soliton-like solutions and period form solutions are found.展开更多
To construct the infinite sequence new exact solutions of nonlinear evolution equations and study the first kind of elliptic function, new solutions and the corresponding B^cklund transformation of the equation are pr...To construct the infinite sequence new exact solutions of nonlinear evolution equations and study the first kind of elliptic function, new solutions and the corresponding B^cklund transformation of the equation are presented. Based on this, the generalized pentavalent KdV equation and the breaking soliton equation are chosen as applicable examples and infinite sequence smooth soliton solutions, infinite sequence peak solitary wave solutions and infinite sequence compact soliton solutions are obtained with the help of symbolic computation system Mathematica. The method is of significance to search for infinite sequence new exact solutions to other nonlinear evolution equations.展开更多
This paper is concerned with the existence and the nonlinear asymptotic stabil- ity of traveling wave solutions to the Cauchy problem for a system of dissipative evolution equations {θt=vζx+(ζθ)x+aθxx,ζt=-θ...This paper is concerned with the existence and the nonlinear asymptotic stabil- ity of traveling wave solutions to the Cauchy problem for a system of dissipative evolution equations {θt=vζx+(ζθ)x+aθxx,ζt=-θx+βζxx;with initial data and end states (ζθ)(x,0)=(ζ0,θ0)(x)→(ζ±,θ±)as x→∞.We obtain the existence of traveling wave solutions by phase plane analysis and show the asymptotic nonlinear stability of traveling wave solutions without restrictions on the coeffi- cients a and v by the method of energy estimates.展开更多
A unified approach is presented for finding the travelling wave solutions to one kind of nonlinear evolution equation by introducing a concept of 'rank'. The key idea of this method is to make use of the arbit...A unified approach is presented for finding the travelling wave solutions to one kind of nonlinear evolution equation by introducing a concept of 'rank'. The key idea of this method is to make use of the arbitrariness of the manifold in Painlevé analysis. We selected a new expansion variable and thus obtained a rich variety of travelling wave solutions to nonlinear evolution equation, which covered solitary wave solutions, periodic wave solutions, Weierstrass elliptic function solutions, and rational solutions. Three illustrative equations are investigated by this means, and abundant travelling wave solutions are obtained in a systematic way. In addition, some new solutions are firstly reported here.展开更多
A new algorithm for symbolic computation of polynomial-type conserved densities for nonlinear evolution systems is presented. The algorithm is implemented in Maple. The improved algorithm is more efficient not only in...A new algorithm for symbolic computation of polynomial-type conserved densities for nonlinear evolution systems is presented. The algorithm is implemented in Maple. The improved algorithm is more efficient not only in removing the redundant terms of the genera/form of the conserved densities but also in solving the conserved densities with the associated flux synchronously without using Euler operator. Furthermore, the program conslaw.mpl can be used to determine the preferences for a given parameterized nonlinear evolution systems. The code is tested on several well-known nonlinear evolution equations from the soliton theory.展开更多
In this paper, the authors discuss the global existence and blow-up of the solution to an evolution p-Laplace system with nonlinear sources and nonlinear boundary condition. The authors first establish the local exist...In this paper, the authors discuss the global existence and blow-up of the solution to an evolution p-Laplace system with nonlinear sources and nonlinear boundary condition. The authors first establish the local existence of solutions, then give a necessary and sufficient condition on the global existence of the positive solution.展开更多
We give the generalized definitions of variable separable solutions to nonlinear evolution equations, and characterize the relation between the functional separable solution and the derivative-dependent functional sep...We give the generalized definitions of variable separable solutions to nonlinear evolution equations, and characterize the relation between the functional separable solution and the derivative-dependent functional separable solution. The new definitions can unify various kinds of variable separable solutions appearing in references. As application, we classify the generalized nonlinear diffusion equations that admit special functional separable solutions and obtain some exact solutions to the resulting equations.展开更多
For the nonconservative schemes of the nonlinear evolution equations, taking the one-dimensional shallow water wave equation as an example, the necessary conditions of computational stability are given. Based on numer...For the nonconservative schemes of the nonlinear evolution equations, taking the one-dimensional shallow water wave equation as an example, the necessary conditions of computational stability are given. Based on numerical tests, the relationship between the nonlinear computational stability and the construction of difference schemes, as well as the form of initial values, is further discussed. It is proved through both theoretical analysis and numerical tests that if the construction of difference schemes is definite, the computational stability of nonconservative schemes is decided by the form of initial values.展开更多
A new numerical method of integrating the nonlinear evolution equations, namely the Taylor expansion method, was presented. The standard Galerkin method can be viewed as the 0_th order Taylor expansion method; while t...A new numerical method of integrating the nonlinear evolution equations, namely the Taylor expansion method, was presented. The standard Galerkin method can be viewed as the 0_th order Taylor expansion method; while the nonlinear Galerkin method can be viewed as the 1_st order modified Taylor expansion method. Moreover, the existence of the numerical solution and its convergence rate were proven. Finally, a concrete example, namely, the two_dimensional Navier_Stokes equations with a non slip boundary condition,was provided. The result is that the higher order Taylor expansion method is of the higher convergence rate under some assumptions about the regularity of the solution.展开更多
文摘In this paper, we shall show that the Hamiltonian structure can be defined for any nonlinear evolution equations which describe surfaces of a constant negative curvature, so that the densities of conservation laws can be considered as corresponding Hamiltonians. This paper obtains the soliton solution and conserved quantities of a new fifth-order nonlinear evolution equation by the aid of inverse scattering method.
基金Supported by the National Natural Science Foundation of China(12201368,62376252)Key Project of Natural Science Foundation of Zhejiang Province(LZ22F030003)Zhejiang Province Leading Geese Plan(2024C02G1123882,2024C01SA100795).
文摘With the urgent need to resolve complex behaviors in nonlinear evolution equations,this study makes a contribution by establishing the local existence of solutions for Cauchy problems associated with equations of mixed types.Our primary contribution is the establishment of solution existence,illuminating the dynamics of these complex equations.To tackle this challenging problem,we construct an approximate solution sequence and apply the contraction mapping principle to rigorously prove local solution existence.Our results significantly advance the understanding of nonlinear evolution equations of mixed types.Furthermore,they provide a versatile,powerful approach for tackling analogous challenges across physics,engineering,and applied mathematics,making this work a valuable reference for researchers in these fields.
基金Supported by Scientific Research Fund of Hunan Provincial Education Departmen(t23A0361)。
文摘An evolution inequality of Sobolev type involving a nonlinear convolution term is considered.By using the nonlinear capacity method and the contradiction argument,the non-existence of the nontrivial local weak solution is proved.
基金provided by the National Natural Science Foundation of China (Grant No. 12271324)the Natural Science Basic Research Program of Shaanxi Province (Grant No. 2024JC-YBQN-0069)+2 种基金the China Postdoctoral Science Foundation (Grant No. 2024M751921)the 2023 Shaanxi Province Postdoctoral Research Project (Grant No.2023BSHEDZZ186)the Fundamental Research Funds for the Central Universities (Grant No. 1301032598)。
文摘The(2 + 1)-dimensional generalized fifth-order Kd V(2GKd V) equation is revisited via combined physical and mathematical methods. By using the Hirota perturbation expansion technique and via setting the nonzero background wave on the multiple soliton solution of the 2GKd V equation, breather waves are constructed, for which some transformed wave conditions are considered that yield abundant novel nonlinear waves including X/Y-Shaped(XS/YS),asymmetric M-Shaped(MS), W-Shaped(WS), Space-Curved(SC) and Oscillation M-Shaped(OMS) solitons. Furthermore, distinct nonlinear wave molecules and interactional structures involving the asymmetric MS, WS, XS/YS, SC solitons, and breathers, lumps are constructed after considering the corresponding existence conditions. The dynamical properties of the nonlinear molecular waves and interactional structures are revealed via analyzing the trajectory equations along with the change of the phase shifts.
基金supported by the National Natural Science Foundation of China(Grant Nos 40676016 and 40876010)the Knowledge Innovation Project of Chinese Academy of Sciences(Grant No KZCX2-YW-Q03-08)+1 种基金LASG State Key Laboratory Special fundE-Institutes of Shanghai Municipal Education Commission of China(Grant No E03004)
文摘This paper studies a generalized nonlinear evolution equation. Using the homotopic mapping method, it constructs a corresponding homotopic mapping transform. Selecting a suitable initial approximation and using homotopic mapping, it obtains an approximate solution with an arbitrary degree of accuracy for the solitary wave. From the approximate solution obtained by using the homotopic mapping method, it possesses a good accuracy.
文摘The process of evolution, especially that of nonlinear evolution, of C-type instability of laminar-turbulent flow transition in nonparallel boundary layers are studied by means of a newly developed method called parabolic stability equations (PSE). Initial conditions, which are very important for the nonlinear problem, are investigated by computing initial solution of the harmonic waves, modifying the mean-flow-distortion, and giving initial value of TS wave and its subharmonic waves at initial station by solving linear PSE. A numerical method with high-order accuracy are developed in the text, the key normalization conditions in the PSE are satisfied, and nonlinear PSE are solved efficiently and implemented stably by the spatial marching. It has been shown that the computed process of nonlinear evolution of C-type instability in Blasius boundary layer is in good agreement with the experimental results.
基金Supported by the Postdoctoral Science Foundation of ChinaChinese Basic Research Plan"MathematicsMechanization and A Platform
文摘A Riccati equation involving a parameter and symbolic computation are used to uniformly construct the different forms of travelling wave solutions for nonlinear evolution equations.It is shown that the sign of the parameter can be applied in judging the existence of various forms of travelling wave solutions.An efficiency of this method is demonstrated on some equations,which include Burgers Huxley equation,Caudrey Dodd Gibbon Kawada equation,generalized Benjamin Bona Mahony equation and generalized Fisher equation.
文摘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.
基金the support of the National Natural Science Foundation of China(No.11675054)the Shanghai Collaborative Innovation Center of Trustworthy Software for Internet of Things(Grant No.ZF1213)the Science and Technology Commission of Shanghai Municipality(No.18dz2271000)。
文摘It has still been difficult to solve nonlinear evolution equations analytically.In this paper,we present a deep learning method for recovering the intrinsic nonlinear dynamics from spatiotemporal data directly.Specifically,the model uses a deep neural network constrained with given governing equations to try to learn all optimal parameters.In particular,numerical experiments on several third-order nonlinear evolution equations,including the Korteweg-de Vries(KdV)equation,modified KdV equation,KdV-Burgers equation and Sharma-Tasso-Olver equation,demonstrate that the presented method is able to uncover the solitons and their interaction behaviors fairly well.
基金supported by the National Natural Science Foundation of China(No.11675054)Shanghai Collaborative Innovation Center of Trustworthy Software for Internet of Things(Grant No.ZF1213)Science and Technology Commission of Shanghai Municipality(No.18dz2271000)。
文摘Solving nonlinear evolution partial differential equations has been a longstanding computational challenge.In this paper,we present a universal paradigm of learning the system and extracting patterns from data generated from experiments.Specifically,this framework approximates the latent solution with a deep neural network,which is trained with the constraint of underlying physical laws usually expressed by some equations.In particular,we test the effectiveness of the approach for the Burgers'equation used as an example of second-order nonlinear evolution equations under different initial and boundary conditions.The results also indicate that for soliton solutions,the model training costs significantly less time than other initial conditions.
文摘A trial equation method to nonlinear evolution equation with rank inhomogeneous is given. As appncations, the exact traveling wave solutions to some higher-order nonlinear equations such as generalized Boussinesq equation, generalized Pochhammer-Chree equation, KdV-Burgers equation, and KS equation and so on, are obtained. Among these, some results are new. The proposed method is based on the idea of reduction of the order of ODE. Some mathematical details of the proposed method are discussed.
基金The project supported by the National Key Basic Research Development Project Program under Grant No.G1998030600the Foundation of Liaoning Normal University
文摘With the use of computer algebra, the method that straightforwardly leads to travelling wave solutions is presented. The compound KdV-Burgers equation and KP-B equation are chosen to illustrate this approach. As a result, their abundant new soliton-like solutions and period form solutions are found.
基金supported by the National Natural Science Foundation of China(Grant No.10862003)the Science Research Foundation of Institution of Higher Education of Inner Mongolia Autonomous Region,China(Grant No.NJZZ07031)the Natural Science Foundation of Inner Mongolia Autonomous Region,China(Grant No.2010MS0111)
文摘To construct the infinite sequence new exact solutions of nonlinear evolution equations and study the first kind of elliptic function, new solutions and the corresponding B^cklund transformation of the equation are presented. Based on this, the generalized pentavalent KdV equation and the breaking soliton equation are chosen as applicable examples and infinite sequence smooth soliton solutions, infinite sequence peak solitary wave solutions and infinite sequence compact soliton solutions are obtained with the help of symbolic computation system Mathematica. The method is of significance to search for infinite sequence new exact solutions to other nonlinear evolution equations.
基金supported by the Natural Science Foundation of China(11001095)the Ph.D.specialized grant of the Ministry of Education of China(20100144110001)+2 种基金the Special Fund for Basic Scientific Research of Central Colleges(CCNU12C01001)supported by the Fundamental Research Funds for the Central Universities(2015IA009)the Natural Science Foundation of China(61573012)
文摘This paper is concerned with the existence and the nonlinear asymptotic stabil- ity of traveling wave solutions to the Cauchy problem for a system of dissipative evolution equations {θt=vζx+(ζθ)x+aθxx,ζt=-θx+βζxx;with initial data and end states (ζθ)(x,0)=(ζ0,θ0)(x)→(ζ±,θ±)as x→∞.We obtain the existence of traveling wave solutions by phase plane analysis and show the asymptotic nonlinear stability of traveling wave solutions without restrictions on the coeffi- cients a and v by the method of energy estimates.
文摘A unified approach is presented for finding the travelling wave solutions to one kind of nonlinear evolution equation by introducing a concept of 'rank'. The key idea of this method is to make use of the arbitrariness of the manifold in Painlevé analysis. We selected a new expansion variable and thus obtained a rich variety of travelling wave solutions to nonlinear evolution equation, which covered solitary wave solutions, periodic wave solutions, Weierstrass elliptic function solutions, and rational solutions. Three illustrative equations are investigated by this means, and abundant travelling wave solutions are obtained in a systematic way. In addition, some new solutions are firstly reported here.
文摘A new algorithm for symbolic computation of polynomial-type conserved densities for nonlinear evolution systems is presented. The algorithm is implemented in Maple. The improved algorithm is more efficient not only in removing the redundant terms of the genera/form of the conserved densities but also in solving the conserved densities with the associated flux synchronously without using Euler operator. Furthermore, the program conslaw.mpl can be used to determine the preferences for a given parameterized nonlinear evolution systems. The code is tested on several well-known nonlinear evolution equations from the soliton theory.
基金supported by a grant from the National High Technology Researchand and Development Program of China (863 Program) (2009AA044501)by NSFC (10776035+2 种基金10771085)by Key Lab of Symbolic Computation and Knowledge Engineering of Ministry of Educationby the 985 program of Jilin University
文摘In this paper, the authors discuss the global existence and blow-up of the solution to an evolution p-Laplace system with nonlinear sources and nonlinear boundary condition. The authors first establish the local existence of solutions, then give a necessary and sufficient condition on the global existence of the positive solution.
基金The project supported by National Natural Science Foundation of China under Grant Nos.10447007 and 10671156the Natural Science Foundation of Shaanxi Province of China under Grant No.2005A13
文摘We give the generalized definitions of variable separable solutions to nonlinear evolution equations, and characterize the relation between the functional separable solution and the derivative-dependent functional separable solution. The new definitions can unify various kinds of variable separable solutions appearing in references. As application, we classify the generalized nonlinear diffusion equations that admit special functional separable solutions and obtain some exact solutions to the resulting equations.
基金supported by the project"Global Changefor Regional Response"of the Important Study Project of the National Natural Science Foundation of China (Grant No.902110041)the Key Innovation Project of the Chinese Academy of Sciences (KZCX3-SW-213).
文摘For the nonconservative schemes of the nonlinear evolution equations, taking the one-dimensional shallow water wave equation as an example, the necessary conditions of computational stability are given. Based on numerical tests, the relationship between the nonlinear computational stability and the construction of difference schemes, as well as the form of initial values, is further discussed. It is proved through both theoretical analysis and numerical tests that if the construction of difference schemes is definite, the computational stability of nonconservative schemes is decided by the form of initial values.
文摘A new numerical method of integrating the nonlinear evolution equations, namely the Taylor expansion method, was presented. The standard Galerkin method can be viewed as the 0_th order Taylor expansion method; while the nonlinear Galerkin method can be viewed as the 1_st order modified Taylor expansion method. Moreover, the existence of the numerical solution and its convergence rate were proven. Finally, a concrete example, namely, the two_dimensional Navier_Stokes equations with a non slip boundary condition,was provided. The result is that the higher order Taylor expansion method is of the higher convergence rate under some assumptions about the regularity of the solution.