Next Article in Journal
Improving Deep Interactive Evolution with a Style-Based Generator for Artistic Expression and Creative Exploration
Next Article in Special Issue
Dynamical Invariant and Exact Mechanical Analyses for the Caldirola–Kanai Model of Dissipative Three Coupled Oscillators
Previous Article in Journal
Exergy Analysis of a Bio-System: Soil–Plant Interaction
Previous Article in Special Issue
Dynamical Analysis of a New Chaotic Fractional Discrete-Time System and Its Control
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Simple Equations Method (SEsM): Algorithm, Connection with Hirota Method, Inverse Scattering Transform Method, and Several Other Methods

by
Nikolay K. Vitanov
1,*,
Zlatinka I. Dimitrova
2 and
Kaloyan N. Vitanov
1
1
Institute of Mechanics, Bulgarian Academy of Sciences, Acad. G. Bonchev Str., Block 4, 1113 Sofia, Bulgaria
2
Institute of Solid State Physics, Bulgarian Academy of Sciences, Blvd. Tzarigradsko Chaussee 72, 1784 Sofia, Bulgaria
*
Author to whom correspondence should be addressed.
Entropy 2021, 23(1), 10; https://doi.org/10.3390/e23010010
Submission received: 13 November 2020 / Revised: 19 December 2020 / Accepted: 22 December 2020 / Published: 23 December 2020
(This article belongs to the Special Issue Dynamical Systems, Differential Equations and Applications)

Abstract

:
The goal of this article is to discuss the Simple Equations Method (SEsM) for obtaining exact solutions of nonlinear partial differential equations and to show that several well-known methods for obtaining exact solutions of such equations are connected to SEsM. In more detail, we show that the Hirota method is connected to a particular case of SEsM for a specific form of the function from Step 2 of SEsM and for simple equations of the kinds of differential equations for exponential functions. We illustrate this particular case of SEsM by obtaining the three- soliton solution of the Korteweg-de Vries equation, two-soliton solution of the nonlinear Schrödinger equation, and the soliton solution of the Ishimori equation for the spin dynamics of ferromagnetic materials. Then we show that a particular case of SEsM can be used in order to reproduce the methodology of the inverse scattering transform method for the case of the Burgers equation and Korteweg-de Vries equation. This particular case is connected to use of a specific case of Step 2 of SEsM. This step is connected to: (i) representation of the solution of the solved nonlinear partial differential equation as expansion as power series containing powers of a “small” parameter ϵ ; (ii) solving the differential equations arising from this representation by means of Fourier series, and (iii) transition from the obtained solution for small values of ϵ to solution for arbitrary finite values of ϵ . Finally, we show that the much-used homogeneous balance method, extended homogeneous balance method, auxiliary equation method, Jacobi elliptic function expansion method, F-expansion method, modified simple equation method, trial function method and first integral method are connected to particular cases of SEsM.

1. Introduction

Nature and human society are rich sources of complex systems (examples can be found, e.g., in economics, social sciences, network theory, dynamics of research groups, etc. [1,2,3,4,5,6,7]). Because of this the complex systems attract much research attention in the last decades, see for examples, [8,9,10,11,12,13,14,15,16,17,18,19,20,21,22]. Most of the complex systems are nonlinear—many such examples can be found in the fluid mechanics or solid-state physics [23,24,25,26,27,28,29]. The effects connected to the nonlinearity can be studied, for example, by means of time series analysis or by means of models based on differential or difference equations (additional information about the methodology of the nonlinear time series analysis, some applications of this methodology and basic information about nonlinear differential equations, can be seen, e.g., in [30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45]). Very often the used model equations are nonlinear partial differential equations, and because of this the methodology for solving such equations is interesting to the researchers. Many years ago the methodology for obtaining exact solutions of nonlinear partial differential equations was connected to transformations that transform the solved nonlinear partial differential equation to a linear differential equation. One example is the Hopf-Cole transformation [46,47] which transforms the nonlinear Burgers equation to the linear heat equation. An appropriate transformation of the Korteweg-de Vries equation connected this equation to the famous equation of Schrödinger and led to the development of the Method of inverse scattering transform [48,49,50]. Almost at the same time, Hirota developed a method for obtaining exact solutions of NPDEs—Hirota method [51,52]. This method is based on the bilinearization of the solved nonlinear partial differential equation after appropriate transformation of the nonlinearity of the equation. Truncated Painleve expansions may lead to many of these appropriate transformations [53,54,55,56,57]. We note the work of Kudryashov who formulated the Method of Simplest Equation (MSE) [58] based on determination of singularity order n of the solved NPDE and on searching of a particular solution of this equation as series containing powers of solutions of a simpler equation called the simplest equation. The methodology was extended [59] and applied for obtaining traveling wave solutions of nonlinear partial differential equations (see, e.g., [60,61,62]). We write several more words on the work of Kudryashov as it is of interest for our discussion below in the text and the Method of Simplest Equation leads to interesting results. Kudryashov [63] used various transformations in order to transform the nonlinearity of a generalized evolution equation of the wave dynamics and to obtain exact solutions of this equation. This research was continued in [57,64,65] and recent results connected to the application of the Method of Simplest Equation can be found, for example, in [66,67,68,69,70].
Recently we presented an algorithm for obtaining exact and approximate solutions of nonlinear partial differential equations called Simple Equations Method (SEsM) [71,72,73,74]. We shall discuss this algorithm and its connection to other methods below in this text. We note that some elements of the methodology can be seen in our articles written almost 30 years ago [75,76,77,78,79,80]. More than 10 years ago [81,82] we have used the ordinary differential equation of Bernoulli as simplest equation [83] and applied methodology called Modified Method of Simplest Equation to ecology and population dynamics [84]. In these publications we have used the concept of the balance equation. The Modified Method of Simplest Equation—MMSE [85,86]—is based on the determination of the kind of the simplest equation and truncation of the series of solutions of the simplest equation by means of application of a balance equation and it is equivalent of the Method of Simplest Equation mentioned above. Up to 2018 our contributions to the methodology and its application have been connected to the MMSE [87,88,89,90,91,92,93,94,95]. We note especially the article [94] where we have extended the methodology of the MMSE to simplest equations of the class
d k g d ξ k l = j = 0 m d j g j
where k = 1 , , l = 1 , , and m and d j are parameters. Equation (1) contains as particular cases, for example: (i) trigonometric functions; (ii) hyperbolic functions; (iii) elliptic functions of Jacobi; (iv) elliptic function of Weierstrass.
In the course of time, we extended the algorithm of the Modified Method of Simplest Equation. Here we are going to discuss the last version which is connected to the possibility of use of more than one simple equation. This modification is called SEsM—Simple Equations Method. The reason for the use of this name is that the used simple equations are more simple than the solved nonlinear partial differential equation, but these simple equations in fact can be quite complicated. Thus we have to substitute the word “simplest” by the word “simple”. A variant of SEsM based on two simple equations was applied in [96] and the first description of the algorithm was made in [71] and then in [72,73,74]. For more applications of particular cases of the algorithm see [97,98].
We note that SEsM is not a universal algorithm for obtaining exact analytical solutions of nonlinear differential equations. With very large probability there are many exact solutions of nonlinear differential equation which can not be obtained by SEsM. The presence of such solutions is of large interest for us as it shows possible ways for further development of SEsM. We note also that some steps of SEsM are not rigidly fixed. Examples are the choice of transformations in Step 1 of SEsM and the selection of the simple equation in Step 5. This flexibility increases the number of equations which can be treated by SEsM but the set of these equations is still a subclass of the class of nonlinear differential equations that possess exact solutions.
The organization of the text below is as follows. We describe SEsM in Section 2. In Section 3 we show that the famous method of Hirota is connected to a particular case of SEsM and then we use the connection between SEsM and Hirota method in order to show that SEsM can lead to solutions of integrable differential equations: examples are the three-solution solution of the Korteweg-de Vries equation, two-soliton solution of the nonlinear Schrödinger equation and the soliton solution of the Ishimori equation. In Section 4 we discuss the connection between SEsM and the inverse scattering transform method. This connection is shown in the examples of the Burgers equation and the Korteweg-de Vries equation. In Section 5 we formulate several assumptions which show that numerous methods for obtaining exact particular solutions of nonlinear partial differential equations, namely, homogeneous balance method, extended homogeneous balance method, auxiliary equation method, Jacobi elliptic function expansion method, F-expansion method, modified simple equation method, trial function method, first integral method, are connected to particular cases of SEsM. Generalizations of some of these methods are formulated. Several concluding remarks are summarized in Section 6.

2. Simple Equations Method (SEsM)

Simple Equations Method (SEsM) is an algorithm for obtaining exact and approximate solutions of nonlinear differential equations. In general (Figure 1) the algorithm is designed to obtain solutions of systems of N nonlinear differential equations by the use of solutions of M simple equations. The most applications of the algorithm up to now are for obtaining solutions of 1 nonlinear differential equations by the use of solutions of (a) one simple equation or (b) more than one simple equation (see, e.g., [83,84,85,86,87,88,89,90,91,92,93,94,95] where the corresponding particular case of SEsM called Modified Method of Simplest Equation is applied for obtaining exact solutions of numerous nonlinear partial differential equations). The development of SEsM started with the use of a solution of one simple equation in order to obtain a solution of one nonlinear differential equation. This version of the algorithm was called Modified Method of Simplest Equation. We shall show below that SEsM is capable to lead to particular exact solutions of nonlinear differential equations. In addition, we shall show that SEsM is connected to the method of Hirota and to the inverse scattering transform method and because of this SEsM can lead not only to particular solutions of some nonlinear differential equations.
SEsM has seven steps which are shown in Figure 2.
We consider a system of nonlinear partial differential equations
A i [ u 1 ( x , , t ) , , u n ( x , , t ) ] = 0 , i = 1 , , n ,
where A i [ u 1 ( x , , t ) , , u n ( x , , t ) , ] depend on the functions u 1 ( x , , t ) , , u n ( x , , t ) and some of their derivatives ( u i can be a function of more than 1 spatial coordinates). Then we proceed as follows.
(1.) 
We apply transformations
u i ( x , . . . , t ) = T i [ F i ( x , , t ) , G i ( x , , t ) , ] ,
where T i ( F i , G i , ) is some function of other functions F i , G i , . In general F i ( x , , t ) , G i ( x , , t ) , … are functions of several spatial variables as well as of the time. The transformation has the goal to transform the nonlinearity of the solved differential equations to more treatable kind of nonlinearity or the transformation may even remove the nonlinearity. In the case of one solved equation the transformation T ( F , G , ) can be: the Painleve expansion; u ( x , t ) = 4 tan 1 [ F ( x , t ) ] in the case of the sine–Gordon equation; u ( x , t ) = 4 tanh 1 [ F ( x , t ) ] in the case of sh-Gordon (Poisson–Boltzmann equation) (for applications of the last two transformations, see, e.g., [75,76,77,78]); u ( x , t ) = F ( x , t ) G ( x , t ) ; u ( x , t ) = i = 0 I a i [ F ( x , t ) ] i j = 0 J b j [ G ( x , t ) ] j ; or another transformation.
In numerous particular cases, one may skip this step (then we have just u i ( x , . . . , t ) = F i ( x , . . . , t ) ) but in many cases the step is necessary for obtaining a solution of the studied nonlinear PDE. The application of (3) to (2) leads to a nonlinear PDEs for the functions F i , G i , .
We note that no general form of the transformations T i is known up to now and because of this we cannot write a general relationships for these transformations here. Moreover, some equations can be treated without such transformations. The transformations allow us to extend the class of equations for which exact solutions can be obtained by means of SEsM.
(2.) 
The functions F i ( x , . . . , t ) , G i ( x , , t ) , … are represented as a function of other functions f i 1 , . . . , f i N , g i 1 , , g i M , …. The functions f and g are connected to solutions of some differential equations (these equations can be partial or ordinary differential equations) which are more simple than Equation (2). We note that the possible values of N and M are N = 1, 2, …, M = 1, 2, … (there may be an infinite number of functions f too). The forms of the functions F i ( f 1 , , f N ) , G i ( x , , t ) , … can be different. For an example for the case of a single solved equation the function F can have the form
F = α + i 1 = 1 N β i 1 f i 1 + i 1 = 1 N i 2 = 1 N γ i 1 , i 2 f i 1 f i 2 + i 1 = 1 N i N = 1 N σ i 1 , , i N f i 1 f i N ,
where α , β i 1 , γ i 1 , i 2 , σ i 1 , , i N are parameters. Of course F ( f 1 , , f N ) can have another form too (i.e., the form of F can be different from (4)). SEsM is very flexible with respect to the form of F i , G i , …. We note that the relationship (4) contains, as a particular case, the relationship used by Hirota [51]. The power series i = 0 N μ n f n (where μ is a parameter) used in the previous versions of the methodology based on one simple equation (i.e., the Modified Method of Simplest Equation) are also a particular case of the relationship (4).
(3.) 
In general the functions used in F i , G i , - the functions f i 1 , , f i N , g i 1 , , g i M are solutions of some partial differential equations. These equations are more simple than the solved nonlinear partial differential equation. There are two possibilities: (i) one may use solutions of the simple partial differential equations if such solutions are available, or (ii) one transforms the more simple partial differential equations by means of appropriate ansätze (e.g., traveling-wave ansätze such as ξ = α ^ x + β ^ t ; ζ = μ ^ y + ν ^ t ; ). Then the solved differential equations for f i 1 , …, f i N , g i 1 , , g i M , … may be reduced to differential equations E l , containing derivatives of one or several functions
E l a ( ξ ) , a ξ , a ξ ξ , , b ( ζ ) , b ζ , b ζ ζ , = 0 ; l = 1 , , N + M + .
In many cases (e.g., if the equations for the functions f 1 , are ordinary differential equations) one may skip this step, but the step may be necessary if the equations for f 1 , are complicated partial differential equations.
(4.) 
We assume that the functions a ( ξ ) , b ( ζ ) , etc., are functions of other functions, such as, v ( ξ ) , w ( ζ ) , etc., e.g,
a ( ξ ) = A [ v ( ξ ) ] ; b ( ζ ) = B [ w ( ζ ) ] ; .
Note that SEsM does not prescribe the forms of the functions A, B, …. Thus, different relationships are possible. Often one uses a finite-series relationship, for example,
a ( ξ ) = μ 1 = ν 1 ν 2 q μ 1 [ v ( ξ ) ] μ 1 ; b ( ζ ) = μ 2 = ν 3 ν 4 r μ 2 [ w ( ζ ) ] μ 2 , .
where q μ 1 , r μ 2 , … are parameters. However, other kinds of relationships, and more complicated ones, are also possible.
(5.) 
The functions v ( ξ ) , w ( ζ ) , … are solutions of simple ordinary differential equations. For about 10 years we have used the particular case of the described methodology that was based on the use of just one simple equation. This simple equation was called the simplest equation and the methodology based on one equation was called the Modified Method of Simplest Equation. SEsM contains the Modified Method of Simplest Equation as a particular case.
(6.) 
The application of the steps 1–5 to Equation (2) transforms the left-hand side of these equations. The results of this transformation can be functions which are sums of terms where each term contains some function multiplied by a coefficient. This coefficient contains some of the parameters of the solved equations and some of the parameters of the solutions. In most cases, a balance procedure must be applied in order to ensure that the above-mentioned relationships for the coefficients contain more than one term (e.g., if the result of the transformation is a polynomial, then the balance procedure has to ensure that the coefficient of each term of the polynomial is a relationship that contains at least two terms). This balance procedure may lead to one or more additional relationships among the parameters of the solved equation and parameters of the solution. These relationships are known as balance equations.
(7.) 
We may obtain a nontrivial solution of Equation (2) if all coefficients mentioned in Step 6 are set to 0. This condition usually leads to a system of nonlinear algebraic equations for the coefficients of the solved nonlinear PDE and for the coefficients of the solution. Any nontrivial solution of this algebraic system leads to a solution the studied nonlinear partial differential equation. Usually, the above system of algebraic equations contains many equations and because of this, the support of a computer algebra system is needed.
Figure 3 shows the most frequently used particular cases of SEsM: the cases when one searches for a solution of one nonlinear differential equation (the cases are: (i) solutions of m > 1 simple equations are used, or (ii) the solution of just one simple equation is used). In these cases, we search for solution u of the solved equation and at step 1 of SEsM we transform u by means of a transformation T. The solved equation is transformed into an equation for the function F and at the following steps of SEsM F is represented by functions f i which are constructed by the known solutions of the used simple equations. Thus the solved nonlinear differential equation is reduced to a system of nonlinear algebraic equations and each nontrivial solution of this system leads to a particular solution of the solved differential equation. In the most simple case of SEsM (the Modified Method of Simplest Equation) we skip the Step 1 (the one with the transformation T) and F is represented as a power series of the solution f of a single simple equation. The substitution of these power series in the solved equation reduces this equation to a system of nonlinear algebraic equations and each nontrivial solution of this system leads to an exact particular solution of the solved nonlinear differential equation.
Below we age going to discuss the relation among SEsM and several much-used methods for obtaining exact solutions of nonlinear partial differential equations.

3. Hirota Method and SEsM

3.1. Hirota Method

Hirota method is a very popular method for obtaining soliton solutions of integrable nonlinear partial differential equations. One of the first applications of the method was in the famous article of Hirota [51] and detailed description of the method and its applications can be found in the book of Hirota [52]. Here, we present a very brief summary of this method which is as follows. At Step 1 of the method, one makes a transformation of the solved nonlinear partial differential equation. Hirota often used the transformation
u ( x , t ) = 2 2 x 2 f ( x , t ) = 2 f 2 f x 2 f x 2 f 2 .
The function f is called sometimes the τ -function. At Step 2 of the method, one searches for a solution of the obtained after the transformation equation in the form
f = α + ϵ f 1 + ϵ 2 f 2 + ϵ 3 f 3 +
where α and ϵ are parameters. At Step 3 of the method, (9) is substituted in the solved nonlinear partial differential equation. At Step 4 of the method, the obtained equations for the orders ϵ , ϵ 2 , ϵ 3 , … are solved. In Step 5 of the method, the exact solution is constructed on the basis of the solutions for ϵ , ϵ 2 , ϵ 3 , ….
In order to deal with Step 4 of his method, Hirota introduced the famous bilinear operators,
D t n ( a . b ) = t t n a ( x , t ) b ( x , t ) at   t = t , D x m ( a . b ) = x x m a ( x , t ) b ( x , t ) at   x = x , D x m D t n ( a . b ) = x x m t t n a ( x , t ) b ( x , t ) at   t = t ; x = x ,
which can be written also as follows [99]
D x m ( a . b ) = j = 0 m ( 1 ) ( m j ) m ! j ! ( m j ) ! j a x j m j b x m j , D x m D t n ( a . b ) = j = 0 m i = 0 n ( 1 ) ( m + n j i ) m ! j ! ( m j ) ! n ! i ! ( n i ) ! i + j a t i x j m + n i j b t n i x m j .
The bilinear operators have useful properties such as
D x m ( a . 1 ) = m a x m
D x m ( a . b ) = ( 1 ) m D x m ( b . a )
D x m ( a . a ) = 0 for m odd
D x m D t n ( exp [ k 1 x ω 1 t ] . exp [ k 2 x ω 2 t ] ) = ( k 1 k 2 ) m ( ω 2 ω 1 ) n exp [ ( k 1 + k 2 ) x ( ω 1 + ω 2 ) t ] .
By means of the bilinear operators in many cases, one can find a solution of the sequence of equations for the orders ϵ , ϵ 2 , ϵ 3 , … without much effort.

3.2. Hirota Method and SEsM

Assumption 1.
The method of Hirota is connected to a particular case of SEsM when the transformation in step 1 of SEsM is the same as the transformation used in Hirota’s method, the representation of function f by means of f 1 , f 2 , from Step 2 of SEsM is (9), the differential equations for f 1 , f 2 , from Step 3 of SEsM are the chain of equations obtained for the orders ϵ, ϵ 2 , … within the scope of Hirota’s method and the simple equations which are used in the construction of the solutions for f 1 , f 2 , are differential equations for exponential functions.
We consider SEsM and impose restrictions on its steps in order to reduce this method to a particular case connected to the Hirota method. We proceed as follows. First, we consider the particular case of SEsM where, in Step 1, the transformation (8) is used or any other transformation used within the scope of Hirota’s method. We note that these transformations are a small part of the possible transformations which can be used in SEsM. At Step 2 of SEsM we consider again a particular case: we use the relationship (9) in order to represent the function f by means of the functions f 1 , f 2 , . This is a particular case from the point of view of SEsM as many other kinds of representations of f by f 1 , f 2 , are possible in SEsM. At Step 3 of SEsM, we use one more particular case by considering differential equations for f 1 , f 2 , to be exactly the equations that are obtained within the scope of application of Hirota’s method. This is a particular case as many other kinds of differential equations can be used in SEsM. At Step 4 of SEsM we consider the relationships between f i and the more simple functions to be exactly these ones which arise when Hirota’s method is used. From the point of view of SEsM this is a particular case as much more kinds of relationships can be used in SEsM and one example for this is a relationship of the kind (4). At Step 5 of SEsM we have to determine the simple equations for the functions which participate in the construction of function f i . We take the particular case when these functions are solutions of very simple equations namely differential equations for exponential functions. We note that much more complicated functions can be used in SEsM at this step such as Jacobi elliptic functions, for example. At Step 6 of SEsM we substitute the relationships for f i in the corresponding differential equations and perform a balance procedure if needed. This will lead to a system of nonlinear algebraic equations. We note that in the case of Korteweg-de Vries equation such balance procedure is not needed. Finally, at Step 7 of SEsM one solves the system of nonlinear algebraic equations and obtains the solution of the corresponding nonlinear PDE.
Thus, by means of particular cases of procedures of SEsM we have reduced SEsM to a particular case connected to the Hirota’s method.
Now on the basis of the above, we shall show that SEsM can lead to soliton and multisoliton solutions of several famous equations.

3.3. Example 1: The Three-Soliton Solution of the Korteweg-de Vries Equation

We consider the Korteweg-de Vries equation
u t + 6 u u x + 3 u x 3 = 0 .
At Step 1 of SEsM we transform the nonlinearity in (16). This is made by the transformation (8). The result is
f 2 f x t + f x f t + f 4 f x 4 4 f x 3 f x 3 + 3 2 f x 2 2 = 0 ,
and this can be written by means of the bilinear operators of Hirota as
( D x D t + D x 4 ) ( f · f ) = 0
At Step 2 of SEsM we represent the function f by means of functions f 1 , f 2 , which will be constructed by means of solutions of the simple equations. This representation is (9). At Step 3, of SEsM we substitute (9) in (17) and obtain the sequence of differential equations for f 1 , f 2 , … as follows
2 x t + 3 x 3 f 1 = 0 .
Note that this is a linear equation. The equations for f 2 and f 3 are
2 x t + 3 x 3 f 2 = D x ( D t + D x 3 ) ( f 1 · f 1 ) ,
2 x t + 3 x 3 f 3 = D x ( D t + D x 3 ) ( f 1 · f 2 + f 2 · f 1 ) .
The equations for f 4 , f 5 , … are obtained in a similar way. At steps 4 and 5 of SEsM we have to represent the functions f i , i = 1 , 2 , by means of other functions which are solutions of our simple equations. The form of the corresponding relationships depends on the particular solution of (19) we start with. In order to obtain the single-soliton solution of the Korteweg-de Vries equation, we are starting with a solution for f 1 constructed by just one function which is solution of a very simple differential equation: the differential equation for exponential function
d g 1 d η 1 = g 1 ; g 1 = exp ( η 1 ) ,
where η 1 = λ 1 x + ω 1 t + σ 1 and λ 1 , ω 1 and σ 1 are parameters. We choose f 1 just as
f 1 = exp ( η 1 )
The substitution of (23) in (19) leads to the algebraic equation (the dispersion relation)
ω 1 + λ 1 3 = 0 .
We note that no balance procedure is required here (i.e., Step 6 of SEsM can be skipped). The substitution of f 1 in (20) leads to 0 in the right-hand side of this equation. Then f 2 can be taken to be 0. The same is the situation with f 3 , f 4 , . Thus we obtain the following solution for f
f = 1 + ϵ f 1 .
The parameter ϵ can be absorbed in σ 1 and the solution for f leads to the one-soliton solution for u by use of (8).
In order to obtain the two-soliton solution of the Korteweg-de Vries equation, we take the solution of (19) constructed by solutions of two simple equations
d g 1 d η 1 = g 1 ; g 1 = exp ( η 1 ) ; d g 2 d η 2 = g 2 ; g 2 = exp ( η 2 ) ,
where η i = λ i x + ω i t + σ i and λ i , ω i and σ i , i = 1 , 2 are parameters. Now we construct f 1 as the sum of the solutions of the two simple equations
f 1 = exp ( η 1 ) + exp ( η 2 ) .
The substitution of (27) in (19) leads to a system of two algebraic equations, namely
ω i + λ i 3 = 0 , i = 1 , 2
The substitution of (27) in (20) leads to a nonlinear equation for f 2 . The solution of this equation is
f 2 = a 12 exp ( η 1 + η 2 ) ,
where a 12 = ( λ 1 λ 2 ) 2 ( λ 1 + λ 2 ) 2 . The substitution of f 2 in (21) leads to 0 in the right-hand side of this equation. Thus we can choose f 3 = 0 and then f 4 = f 5 = = 0 . The solution for f becomes
f = 1 + ϵ [ exp ( η 1 ) + exp ( η 2 ) ] + ϵ 2 a 12 exp ( η 1 + η 2 ) .
Again the ϵ and ϵ 2 can be absorbed by σ 1 , 2 and a 1 , 2 and the application of (8) leads to the two-soliton solution of the Korteweg-de Vries equation.
In order to obtain the three-soliton solution of the Korteweg-de Vries equation, we start with the solution for f 1 constructed by means of solutions of three simple equations for exponential functions,
d g i d η i = g i ; g i = exp ( η i ) ; i = 1 , 2 , 3 ,
where η i = λ i x + ω i t + σ i and λ i , ω i and σ i , i = 1 , 2 , 3 are parameters.Then we present f 1 as the sum of the solutions of the three simple equations
f 1 = exp ( η 1 ) + exp ( η 2 ) + exp ( η 3 ) .
The substitution of (32) in (19) leads to the algebraic relationships (dispersion relations)
ω i + λ i 3 = 0 , i = 1 , 2 , 3 ,
and the substitution of (32) in (20) leads to the following relationship for f 2
f 2 = a 12 exp ( η 1 + η 2 ) + a 13 exp ( η 1 + η 3 ) + a 23 exp ( η 2 + η 3 ) ,
where
a i j = ( λ i λ j ) 2 ( λ i + λ j ) 2 , i , j = 1 , 2 , 3 , i < j .
The substitution of the obtained solutions for f 1 and f 2 in (21) leads to the following solution for f 3
f 3 = b 123 exp ( η 1 + η 2 + η 3 ) , b 123 = a 12 a 13 a 23
f 3 is a single exponential function again and then the right-hand side of the equation for f 4 is 0 and we can take f 4 = 0 . Then we can continue with f 5 = f 6 = = 0 and the obtained solution for f is
f = 1 + ϵ [ exp ( η 1 ) + exp ( η 2 ) + exp ( η 3 ) ] + ϵ 2 [ a 12 exp ( η 1 + η 2 ) + a 13 exp ( η 1 + η 3 ) + a 23 exp ( η 2 + η 3 ] + ϵ 3 b 123 exp ( η 1 + η 2 + η 3 )
and after the absorption of epsilons and application of (8) we obtain the three-soliton solution of the Korteweg-de Vries equation.
The procedure can be easily continued and the N-soliton solution of the Korteweg-de Vries equation can be obtained by means of a particular case of the SEsM methodology. We use N simple equations for exponential functions and begin by solution f 1 which is a sum of the solutions of these simple equations. Then step by step we obtain f i up to i = N and for i > N we can set f i = 0 and then we can construct the N-soliton solution of the Korteweg-de Vries equation.

3.4. Example 2: The Two-Soliton Solution of the Nonlinear Schrödinger Equation

Now we show how SEsM leads to the two-soliton solution of the nonlinear Schrödinger equation. Here the transformation at Step 1 of SEsM leads to equations for two functions F and G. The connection of SEsM with the Hirota’s method makes this task of obtaining the two-soliton solution quite easy. The nonlinear Schrödinger equation is
i ψ t + 2 ψ x 2 + q ψ 2 ψ = 0 ,
where q is a parameter. We consider the case q < 0 (and we shall set below q = 2 for convenience) with boundary conditions ψ 2 = ρ 0 2 at x ± . At Step 1 of SEsM we transform the nonlinearity in (37) by the transformation
ψ = G F .
Note that the transformation (38) is different in comparison to the transformation in the case of the Korteweg-de Vries equation. This is an illustration of the feature of SEsM to allow for different kinds of transformations at Step 1 of application of the methodology.
The substitution of (38) to (37) leads to (* means complex conjugated quantity and D x and D t are the operators of Hirota, described above in the text)
i D t G · F F 2 + D x 2 G · F F 2 G F D x 2 F · F F 2 2 G F G G * F 2 = 0 .
This can be written as
{ i [ D t + D x 2 ] G · F } F 2 G F [ D x 2 F · F 2 G G * ] = 0 .
Introducing the constant λ we can write two coupled equations of the basis of (39)
[ i D t + D x 2 ] G · F = λ G F
D x 2 + 2 G G * = λ F 2
λ is a constant which will be determined below. At Step 2 of SEsM we represent F and G by functions which will be then constructed by solutions of simple equations (differential equations for exponential functions in the case discussed here). These representations are made by the following expansions
F = 1 + ϵ f 1 + ϵ 2 f 2 + ,
G = g 0 ( 1 + ϵ g 1 + ϵ 2 g 2 + )
f 1 , f 2 , , g 1 , g 2 , have to go to 0 at x . Then at x g o g 0 * = ρ 0 2 and from (43) and (44) at x
[ i D t + D x 2 λ ] g 0 · 1 = 0 ,
[ D x 2 λ ] 1 · 1 = 2 g 0 g 0 * ,
one obtains
λ = 2 ρ 0 2 ; g 0 = ρ 0 exp ( i η ) ; θ = k x ω t ; ω = k 2 + 2 ρ 0 2 .
Above, k is a real constant. The substitution of (43) and (44) in (41) and (42) leads to relationships for the different powers of ϵ as follows. For the terms of order of ϵ
[ i ( D t + 2 k D x ) + D x 2 ] ( g 1 · 1 + 1 · f 1 ) = 0
[ D x 2 2 ρ 0 2 ] ( f 1 · 1 + 1 · f 1 ) = 2 ρ 0 2 ( g 1 + g 1 * ) .
For the terms of order of ϵ 2
[ i ( D t + 2 k D x ) + D x 2 ] ( g 2 · 1 + g 1 · f 1 + 1 · f 2 ) = 0
[ D x 2 2 ρ 0 2 ] ( f 2 · 1 + f 1 · f 1 + 1 · f 2 ) = 2 ρ 0 2 ( g 2 + g 1 + g 1 * + g 2 * ) .
For the terms of order of ϵ 3
[ i ( D t + 2 k D x ) + D x 2 ] ( g 3 · 1 + g 2 · f 1 + g 1 · f 2 + 1 · f 3 ) = 0
[ D x 2 2 ρ 0 2 ] ( f 3 · 1 + f 2 · f 1 + f 1 · f 2 + 1 · f 3 ) = 2 ρ 0 2 ( g 3 + g 2 + g 1 * + g 1 g 2 * + g 3 * ) .
Next we write relationships if the functions f 1 , , g 1 , by means of functions which are solutions of simple equations. (48) and (49) have many possible solutions. The most simple ones are constructed by means of a simple function which is a solution of a very simple differential equation: the differential equation for exponential function (Steps 3–5 of SEsM). These solutions are
f 1 = exp ( η ) ; g 1 = b exp ( η ) ,
where η = q x ω t and for the satisfaction of (48) and (49) the following relationships must hold (there is no need of balance procedure—Step 6. of SEsM and the relationships follows from the system of algebraic equations obtained from the substitution of (53) in (48) and (49)—Step 7 of SEsM)
ω = q 2 k 4 ρ 0 2 q 2 ; b = q 2 + i ( ω 2 k q ) q 2 i ( ω 2 k q ) .
This choice of f 1 and g 1 leads to f 2 = f 3 = 0 and g 2 = g 3 = = 0 and the solution of the nonlinear Schrödinger equation is
ψ = ρ 0 exp ( i θ ) 1 + b exp ( η ) 1 + exp ( η ) .
In order to obtain the two-soliton solution of the nonlinear Schrödinger equation, we start with solutions of (48) and (49) which depend on solutions of two simple equations for exponential functions
f 1 = exp ( η 1 ) + exp ( η 2 ) ; g 1 = b 1 exp ( η 1 ) + b 2 exp ( η 2 ) ,
where η i = q i x ω i t + σ i , i = 1 , 2 . The substitution of (57) in (48) and (49) (Steps 6 and 7 of SEsM) leads to a system of algebraic equations which solution is
ω i = q i 2 k 4 ρ 0 2 q i 2 ; b i = exp ( 2 i ϕ i ) ; ϕ i = tan 1 q i ω i 2 k q i , i = 1 , 2 .
The substitution of (57), (58) in (50) and (51) leads to solutions containing single exponential function
f 2 = a 12 exp ( η 1 + η 2 ) ; g 2 = b 12 exp ( η 1 + η 2 ) ,
where
a 12 = s i n [ 1 2 ( ϕ 1 ϕ 2 ) ] s i n [ 1 2 ( ϕ 1 + ϕ 2 ) ] 2 ; b 12 = b 1 b 2 a 12 , q i = 2 ρ 0 sin ( ϕ i ) , i = 1 , 2 .
The relationships for f 2 and g 2 leads to f 3 = f 4 = = 0 ; g 3 = g 4 = = 0 and the two-soliton solution of the nonlinear Schrödinger equation becomes
ψ = ρ 0 exp ( i θ ) 1 + exp ( η 1 + 2 i ϕ 1 ) + exp ( η 2 + 2 i ϕ 2 ) + a 12 exp ( η 1 + η 2 + 2 i ϕ i + 2 i ϕ 2 ) 1 + exp ( η 1 ) + exp ( η 2 ) + a 12 exp ( η 1 + η 2 ) .
More complicated soliton solutions of the nonlinear Schrödinger equation can be obtained too. One has just to start by solution constructed by 3,4,… solutions of simple differential equations for the exponential function.

3.5. Example 3: The Soliton Solution of the Ishimori Equation

The last example is connected to the area of solid state physics: dynamics of spin chains. The Ishimori equation [100] is (2+1)-dimensional generalization of the Heisenberg ferromagnetic spin equation and it was introduced in order to explain the dynamics of a classical spin system on a plane. The Ishimori equation is
S t = S 2 S x 2 + σ 2 2 S y 2 + ϕ y S x + ϕ x S y ,
where
2 ϕ x 2 σ 2 2 ϕ y 2 = 2 σ 2 S · S x S y ,
S ( x , y , t ) = ( S 1 , S 2 , S 3 ) is the tree-dimensional spin unit vector, ϕ ( x , y , t ) is a scalar field, σ 2 = ± 1 and ∧ means exterior (wedge) product of corresponding vectors. Step 1 of SEsM is connected with the transformation of the nonlinearity of the studied equation [101,102]. Before Step 1 we introduce the stereographic projection of the spin of the unit sphere on a complex plane, the spin components can be written in terms of the stereographic variable ω ( ω ( r , t ) = S 1 + i S 2 1 + S 3 ) as follows
S + = S 1 + i S 2 = 2 ω 1 + ω , S 3 = 1 ω 2 1 + ω 2 .
The Ishimori equation becomes
i ω t + 2 ω x 2 + σ 2 2 ω y 2 2 ω * 1 + ω 2 ω x 2 + σ 2 ω y 2 i ϕ y ω x i ϕ x ω y = 0 ,
where
2 ϕ x 2 σ 2 2 ϕ y 2 = 4 i σ 2 ( 1 + ω 2 ) 2 ω * x ω y ω x ω * y .
At Step 1 of SEsM, we perform the transformation
ω = g ( x , y , t ) f ( x , y , t ) ,
where f and g are complex functions. This allows us to write the Ishimori equation in the bilinear form by means of the Hirota operators
( i D t D x 2 σ 2 D y 2 ) ( f * · g ) = 0 ( i D t D x 2 σ 2 D y 2 ) ( f * · f g * · g ) = 0 D x [ D x ( f * f + g * g ) · ( f * f + g * g ) ] = σ 2 D y [ D y ( f * f + g * g ) ] · ( f * f + g * g ) ϕ x = 2 i σ 2 D y ( f * · f + g * · g ) f * · f + g * · g ϕ y = 2 i σ 2 D x ( f * · f + g * · g ) f * · f + g * · g
At Step 2 of SEsM, we use two expansions: one expansion for for f and one expansion for g:
f = 1 + n = 1 ϵ 2 n f 2 n ; g = n = 0 ϵ 2 n + 1 g 2 n + 1 .
After substitution of (68) in (67) we arrive at a system of differential equations for the functions f i and g i and we have to solve these systems in a way analogous to the previous examples. In order to obtain the N-soliton solution we start by solution constructed by solutions of simple equations for exponential functions as follows
g 1 = j = 1 N exp ( ξ j ) , ξ j = l j x + m j y + n j t ,
where l j , m j , n j are complex constants. We shall write some details about the single soliton solution of the Ishimori equation. We start from
g 1 = M exp ( ξ 1 ) ; f 2 = e x p [ 2 ( ξ 1 R + ψ ) ] ,
where M is an arbitrary complex constant, ξ 1 = l 1 x + m 1 y + n 1 t , n 1 = i ( l 1 2 + σ 2 m 1 2 ) ; ξ 1 has real and imaginary part: ξ 1 = ξ 1 R + ξ 1 I and exp ( 2 ψ = σ 2 m 1 2 l 1 2 ( l 1 + l 1 * ) 2 σ 2 ( m 1 + m 1 * ) 2 ) . Let us consider further the case σ 2 = 1 . Then we have to solve the system of equations arising by setting to 0 the relationships obtained for the different powers of ϵ . The solution is obtained as in the case of the previous two examples and the obtained results is
S + = 2 E ( l 1 R 2 m 1 R + m 1 R 2 l 1 I + L ) exp ( i ξ 1 I ) sech ( ξ 1 R ) A + 2 B tanh ( ξ 1 R ) + C tanh 2 ( ξ 1 R ) ; S 3 = 1 2 ( l 1 R 2 m 1 R 2 ) 3 sech 2 ( ξ 1 R ) A + 2 B tanh ( ξ 1 R + C tanh 2 ( ξ 1 R ) ) ,
where
A = 2 l 1 R 6 2 m 1 R 2 ( l 1 R 4 + l 1 R 3 m 1 I ) + m 1 R 4 l 1 I 2 + 2 l 1 R 2 l 1 I m 1 R 3 + m 1 I 2 m 1 R 4 + 3 l 1 R 2 m 1 R 4 m 1 R 6 , B = l 1 R 2 m 1 R 2 ( m 1 R 2 + l 1 I 2 + m 1 I 2 + l 1 R 2 ) + l 1 R 3 l 1 I m 1 R + l 1 I ( m 1 R 5 l 1 R 5 ) l 1 R m 1 I m 1 R 4 , C = 2 m 1 R 6 + 2 l 1 I l 1 R 2 m 1 R 3 + l 1 R 4 l 1 I 2 + m 1 I 2 l 1 R 4 + 3 m 1 R 2 l 1 R 4 2 l 1 R 2 m 1 R 4 2 m 1 I l 1 R 3 m 1 R 2 l 1 R 6 , E = l 1 R + i m 1 R l 1 R 2 m 1 R 2 , L = i ( n 1 I m 1 R 2 l 1 R 3 ) + [ m 1 R 3 + l 1 I l 1 R 2 + i ( m i I l 1 R 2 l 1 R m 1 R 2 ) ] .
In a similar way we can obtain two-soliton solution, three-soliton solution, etc. of the Ishimori equation.

4. SEsM and Its Connection with the Inverse Scattering Transform Method

4.1. The Inverse Scattering Transform Method

Below we shall discuss the connection between SEsM and IST method for the case of Korteweg-de Vries equation. First of all we briefly remember application of the IST for the case of the KdV equation:
u t 6 u u x + 3 u x 3 = 0 .
Researchers tried to solve this equation by means of some transformation. Almost 50 years ago Gardner, Greene, Kruskal, and Miura [48] tried the transformation
u = 1 ψ 2 ψ x 2 + λ .
This can be written as
2 ψ x 2 + ( λ u ) ψ = 0 ,
which is the Schrödinger equation where u depends on x and t. However, t can be treated as a parameter and then (75) will be considered as a problem for the scattering of a particle in a potential u ( x , t ) associated with the Korteweg-de Vries equation and described by the linear Schrödinger equation. If the time is a parameter, then for any value of t, there will be a separate scattering problem.
The substitution of (75) in (73) leads to
ψ 2 d λ d t + x ψ P x ψ t P = 0 ,
where
P = ψ t + 3 ψ x 3 3 ( λ + u ) ψ x .
We shall search for solutions of the Korteweg-de Vries equations that decay fast to 0 at x . The associated scattering problem described by (75) has two kinds of eigenvalues for λ : (i) For λ < 0 the values of λ are discrete. We shall write for these eigenvalues λ n = κ n 2 and for the corresponding eigenfunctions ψ n we have ψ n 0 at x 0 and 0 < d x ψ n 2 < . The integration of (76) from to leads to the conclusion that λ n do not depend on t; (ii) For λ > 0 the number of possible values of λ (the spectrum of λ ) is continuous and we shall assume that the values from this spectrum do not depend on t. Then d λ d t = 0 and from (76) and (77) we obtain
P = ψ t + 3 ψ x 3 3 ( λ + u ) ψ x = C ψ ,
where C does not depend on x. Then, (78) describes the evolution of ψ for a fixed value of the parameter λ in the case of continuous spectrum of values λ as well as for the case of discrete spectrum of values of λ .
Next we have to determine C which is done by considering (75) and setting λ = μ 2 . We consider functions which are proportional to exp ( i μ x ) , Im μ 0 at x + and assume that the function of this kind is an asymptotic solution of (78) for any t. This leads to C = 4 i μ 3 and the system of equations becomes
2 ψ x 2 + ( μ 2 u ) ψ = 0 ,
ψ t + 3 ψ x 3 3 ( μ 2 + u ) ψ x + 4 i μ 3 ψ = 0 .
The methodology which is used for calculation of u ( x , t ) comes from the inverse scattering problem from physics: the scattering potential u can be reconstructed on the basis of the knowledge of the scattering coefficient for the waves arriving from x = + and on the basis of knowledge about the spectrum. This methodology has four steps. At step 1, we know u ( x , 0 ) and use this information to solve (79). This leads to discrete eigenvalues μ = i κ n which correspond to the eigenfunctions ψ n and the scattering coefficient β of the incoming waves. The eigenfunctions are chosen as
ψ n ( x ) = χ ( i κ n , x ) .
χ have to satisfy the condition described above, namely, one works with functions that are proportional to exp ( i μ x ) , Im μ 0 at x + . In addition normalization coefficients
γ n = 1 + d x ψ n 2 ,
are introduced. At Step 2 one constructs the solution containing scattering part characterized by the scattering coefficient β . We assume that the solution ψ ( k , x ) of (79) has the following asymptotic behavior
ψ ( k , x ) exp ( i k x ) + β ( k ) exp ( i k x ) , x + , ψ ( k , x ) α ( k ) exp ( i k x ) , x .
k is a positive number, α is the transmission coefficient and β is the reflection coefficient. These coefficients are determined by substitution of the solution (83) in the (79). Thus, one can solve the scattering problem: finding κ n , γ n and β k when the potential u ( x , 0 ) is known. At step 3 one has to determine the time behavior of the scattering parameters. It follows from (79) and (80) (remember that κ n do not change)
d d t + d x χ 2 = 2 χ 2 χ x 2 + 4 χ x 2 + 6 μ 2 χ 2 + 8 i μ 3 + d x χ 2 .
At x ± one has μ = i κ n and ψ n ( x , t ) = χ ( x , t , i κ n ) 0 . The normalization coefficient is
c n ( t ) = γ n exp ( 8 κ n 3 t ) .
The solution of the scattering problem at x has the behavior
ψ ( k , x , t ) f ( k , t ) exp ( i k x ) + g ( k , t ) exp ( i k x ) .
(86) must be an asymptotic solution of (80) when μ = k . This condition leads to the relationships f ( k , t ) = exp ( 8 i k 3 t ) ; g ( k , t ) = β and the reflection coefficient is
b ( k , t ) = g ( k , t ) f ( k , t ) = β ( k ) exp ( 8 i k 3 t ) .
Finally, at step 4, one has to solve the inverse scattering problem. One has to find u ( x , t ) on the basis of known scattering data κ n , c n ( t ) and b ( u , t ) . This happens on the basis of the relationship
u ( x , t ) = 2 d d x K ( x , x , t ) ,
where K ( x , y , t ) is the solution of a linear integral equation known as the equation of Gelfand–Levitan–Marchenko
K ( x , y , t ) + B ( x + y , t ) + x d z K ( x , z , t ) B ( z + y , t ) = 0 ,
where
B ( x + y , t ) = c n ( t ) exp [ κ n ( x + y ) ] + 1 2 π + d k b ( k , t ) exp [ i k ( x + y ) ] = γ n exp [ κ n ( x + y ) + 8 κ n 3 t ] + 1 2 π + d k β ( k ) exp [ i k ( x + y ) + 8 i k 3 t ] .
We note that the parameters κ n , γ n and β ( k ) are determined on the basis of the knowledge of u ( x , 0 ) from step 1 of the schema.

4.2. Connection between SEsM and the Inverse Scattering Transform Method

We shall discuss below the connection between SEsM and IST for the case of Korteweg-de Vries equation. The main points are as follows. We skip Step 1 of SEsM and consider the particular case where no transformation of the solved equation is performed. Then we represent the searched solution as
u ( x , t ) = n = 1 ϵ n u n ( x , t ) .
This is Step 2 of SEsM. One initially treats ϵ as a small parameter, obtains a solution for u and the powers of ϵ are absorbed in the corresponding parameters of the obtained solution. The introduction of (91) in the solved equation leads to a system of equations for u n —Step 3 of SEsM. The obtained system of equations is solved on the basis of solutions constructed from solutions of simple equations for exponential functions (represented by a Fourier series constructed on the basis of these exponential functions. Note that Fourier series method for obtaining solutions of partial differential equations is a particular case of SEsM [74]). Then by means of the appropriate transformation the obtained solution for the case of small values of ϵ is transformed into a solution for finite values of ϵ .
Rosales [103] used the Fourier series in order to obtain the IST methodology for many equations. Below we shall use this excellent work in order to demonstrate the connection between SEsM and IST. We shall use the simple case of Burgers equation and then we shall discuss the methodology for the Korteweg-de Vries equation and in this case we shall arrive at the Gelfand–Levitan–Marchenko equation.

4.3. Example 1: The Burgers Equation

The Burgers equation is
u t + u u x 2 u x 2 = 0 .
We skip Step 1 of SEsM: the transformation of nonlinearity and at Step 2 of SEsM we write the function u ( x , t ) by functions u n ( x , t ) as in (91). For now we shall treat ϵ as a “small” parameter. Our strategy is as follows: First we search for a solution for which u n 0 when x and ϵ is small. Then we rewrite this solution to have a solution valid for all x and we shall no longer require ϵ to be small.
The equations for u n ( x , t ) can be obtained by substitution of (91) in (92) and collection of the terms containing equal powers of ϵ . We obtain
u n t 2 u n x 2 = j = 1 n 1 u n j u j x , n = 2 , 3 ,
Now step by step we can solve the Equation (93) and we can obtain u n . For n = 1 we have
u 1 t 2 u 1 x 2 = 0 .
At Steps 3–5 of SEsM, we are going to represent u n by solutions of more simple equations. As in the case of the method of Hirota above, we shall use as simple equations the equations for exponential functions and we construct u 1 by solutions of these simple equations in the form of a Fourier transform
u 1 = C d λ ( k ) ( 2 i k ) exp [ i k x k 2 t ] ,
where d λ ( k ) is an appropriate measure in the complex plane C and the factor 2 i k is introduced for convenience. The reason for the possibility of representation (95) is that (94) is a linear equation. We note that we can choose the measure above in such a way that u 1 is a superposition of a Fourier transform plus real exponentials
u 1 = m α m exp [ κ m x + κ m 2 ] + d k β ( k ) ( 2 i k ) exp [ i k x k 2 t ] ,
where α m are parameters and β ( k ) is a function.
The choice of solution for u 1 as (95) influences the equations for u 2 , u 3 , …. The equation for u 2 becomes
u 2 t 2 u 2 x 2 = u 1 u 1 x = C 2 d λ ( k 1 ) d λ ( k 2 ) ( 4 i k 1 2 k 2 ) exp [ i ( k 1 + k 2 ) x ( k 1 2 + k 2 2 ) t ] ,
and the equation for u 3 becomes
u 3 t 2 u 3 x 2 = u 2 u 1 x u 1 u 2 x = C 3 d λ ( k 1 ) d λ ( k 2 ) d λ ( k 3 ) 4 i [ k 1 2 k 2 + k 1 k 3 ( k 1 + k 2 ) ] exp { i [ ( k 1 + k 2 + k 3 ) x + i ( k 1 2 + k 2 2 + k 3 2 ) t ] } .
We search for solution of (97) as
u 2 ( x , t ) = C 2 d λ ( k 1 ) d λ ( k 2 ) Φ 2 ( k 1 , k 2 ) exp [ i ( k 1 + k 2 ) x ( k 1 2 + k 2 2 ) t ] ,
and for solution for u 3 we assume
u 3 ( x , t ) = C 3 d λ ( k 1 ) d λ ( k 2 ) d λ ( k 3 ) Φ 3 ( k 1 , k 2 , k 3 ) exp { i [ ( k 1 + k 2 + k 3 ) x + i ( k 1 2 + k 2 2 + k 3 2 ) t ] } .
The substitution of (99) in (97) and of (100) in (98) leads to following relationships for Φ 2 and Φ 3
Φ 2 ( k 1 , k 2 ) = 2 i k 1 ; Φ 3 ( k 1 , k 2 , k 3 ) = 2 i k 1 .
Next we assume the general form of the solution of the equations for u n to be:
u n ( x , t ) = C n [ d λ ( k ) ] n Φ n exp { i Ω n } ,
where [ d λ ( k ) ] n = d λ ( k 1 ) d λ ( k n ) ; Φ n = Φ n ( k 1 , , k n ) and Ω n = ( k 1 + + k n ) x + i ( k 1 2 + + k n 2 ) t . Because we have Φ 2 = Φ 3 = 2 i k 1 we assume further that
Φ n ( k 1 , , k n ) = 2 i k 1 ,
and this is true indeed as it can be seen by direct substitution of (102) in (93). Then we write the solution of Burgers equation at x and small values of ϵ
u = n = 1 ϵ n C n [ d λ ( k ) ] n ( 2 i k 1 ) exp { i Ω n } .
By means of the notation
ζ = ϵ C d λ ( k ) exp { i k x k 2 t } ,
we can write u as
u = n = 1 ϵ n C n [ d λ ( k ) ] n ( 2 i k 1 ) exp { i Ω n } = n = 1 2 ζ n 1 ζ x .
We use the formula for summation of the power series
1 + q + q 2 + q k 1 + = 1 1 q , q < 1 ,
and apply it to (106) with q = ζ . The result is
u = 2 ζ x 1 1 ζ = 2 x ln ( 1 ζ ) .
We remember that ζ and thus 1 ζ where ϵ is a parameter, is the general solution of the heat equation. In addition (107) is exactly the Hopf-Cole transformation [46,47] which linearizes the Burgers equation and reduces it to the heat equation. Thus the obtained solution (106) is a solution of the Burgers equation not only for small ϵ and for x but also for large ϵ and for arbitrary x.
Let us now consider a more complicated equation: the Korteweg-de Vries equation.

4.4. Example 2: The Korteweg-de Vries Equation

We consider the Korteweg-de Vries equation
u t + 6 u u x + 3 u x 3 = 0 .
We skip Step 1 of SEsM (the transformation of the nonlinearity in (108)). At Step 2 of SEsM, we represent u ( x , t ) by means of other functions u 1 , u 2 , which then will be connected to the solutions of the simple equations. We use the relationship (91) where ϵ is considered initially to be a small parameter. The substitution of (91) in (108) leads to equations for u 1 , u 2 , as follows
u n t + 3 u n x 3 = 6 x j = 1 n 1 u j u n j , n = 2 , 3 ,
and the equation for u 1 is
u 1 t + 3 u 1 x 3 = 0
At Step 3 of SEsM, we shall use solutions of the obtained equations for u 1 , u 2 , … which are connected through Fourier series to the solution of a simple equation for exponential functions (Steps 4 and 5 of SEsM). The solution for u 1 is
u 1 = C d λ ( k ) ( k ) exp [ i k x + k 3 t ] ,
where d λ ( k ) is an appropriate measure in the complex plane C and the term ( k ) is introduced for convenience. This choice of u 1 influences the form of equations for u 2 , u 3 , …. For an example the equations for u 2 and u 3 are
u 2 t + 3 u 2 x 3 = 6 i C 2 d λ ( k 1 ) d λ ( k 2 ) k 1 k 2 ( k 1 + k 2 ) exp { i [ ( k 1 + k 2 ) x + ( k 1 3 + k 2 3 ) t ] } ,
u 3 t + 3 u 3 x 3 = 6 x ( u 1 u 2 + u 2 u 3 ) = = 6 i C 3 d λ ( k 1 ) d λ ( k 2 ) d λ ( k 3 ) ( k 1 + k 2 ) ( k 1 + k 2 + k 3 ) exp { i [ ( k 1 + k 2 + k 3 ) x + ( k 1 3 + k 2 3 + k 3 3 ) t ] } .
The solutions for u 2 , u 3 , are searched again in the form
u n ( x , t ) = C n [ d λ ( k ) ] n Φ n ( k 1 , , k n ) exp { i Ω n } ,
and by substitution of this relationship in equations for u 2 , u 3 , we will determine Φ n and Ω n (Steps 6. and 7. of SEsM).
The substitution of (114) in (112) leads to
Φ 2 ( k 1 , k 2 ) = 1 ; Ω 2 = ( k 1 + k 2 ) x + ( k 1 3 + k 2 3 ) t .
The substitution of (114) in (113) leads to
Φ 3 ( K 1 , k 2 , k 3 ) = k 1 + k 2 + k 3 ( k 1 + k 2 ) ( k 2 + k 3 ) ; Ω 3 = ( k 1 + k 2 + k 3 ) x + ( k 1 2 + k 2 2 + k 3 2 ) t .
In addition we can obtain for the parameters of the solution u 4
Φ 4 = k 1 + k 2 + k 3 + k 4 ( k 1 + k 2 ) ( k 2 + k 3 ) ( k 3 + k 4 ) ; Ω 4 = ( k 1 + k 2 + k 3 + k 4 ) x + ( k 1 2 + k 2 2 + k 3 2 + k 4 2 ) t .
The continuation of the calculations leads us to
u n = ( 1 ) n + 1 i x C n [ d λ ( k ) ] n exp ( i Ω n ) j = 1 n 1 ( k j + k j + 1 ) .
The solution of the Korteweg-de Vries equation becomes
u = i x n = 1 ( ϵ ) n C n [ d λ ( k ) ] n exp ( i Ω n ) j = 1 n 1 ( k j + k j + 1 ) .
The straight summation of the infinite series in the sum in (119) is impossible because the different harmonics are coupled by the factors ( k j + k j + 1 ) 1 but the summation still can be done as follows. We write (119) as
u ( x , t ) = x n = 1 n ( ϵ ) n C n d [ λ ( k ) ] n p ^ ( k 1 ) P ^ ( k 1 , k 2 ) P ^ ( k 2 , k 3 ) , P ^ ( k n 1 , k n ) p ^ ( k n ) ,
where
p ^ k = exp i 2 [ k x ω ( k ) t ] ; P ^ ( k , q ) = i p ^ ( k ) p ^ ( q ) k + q ,
and ω ( k ) = k 3 . Now there are two possibilities for the measure d λ ( k ) : to be discrete or to be continuous.
We consider first the case of discrete d λ ( k ) . The integral from (120) can be replaced by a sum as follows
C d λ ( k ) f ( k ) = m a m 2 f ( i k m ) .
Then
u = x n = 1 ( ϵ ) n m j p m 1 P m 1 , m 2 P m 2 , m 3 P m n 1 , m n p m n ,
where
P m q = p m p q k m + k q ; p m = a m p ^ ( i k m ) = a m exp 1 2 ( k m x + k m 3 t ) ,
(123) is a matrix product and if p is the column vector of all p m and P is the square matrix of all P m q we can write p T is the transpose of p
u = ϵ x p T ( I + ϵ P ) 1 p .
For k m > 0 and real a m the matrix P is real symmetric and positive definite. Then u from (125) is nonsingular for ϵ > 0 and ϵ can be absorbed in the coefficients a m . This means that the solution (125) is no longer limited to small values of ϵ and we can set ϵ = 1 (absorbing other values of ϵ in coefficients a m ).
We can write (125) in much more known form. In order to do this we observe that P x = 1 2 p p T and then
u = 2 x Tr [ ( I + P ) 1 P x ] = 2 2 x 2 Tr [ ln ( I + P ) ] = 2 2 x 2 ln det ( I + P ) .
Thus we arrive at
u = 2 2 x 2 ln det ( I + P ) ,
which is the relationship for the multisoliton solution of the Korteweg-de Vries equation.
Up to now we have considered the case of discrete measure. In the case of continuous measure d λ ( k ) we arrive at the Gelfand–Levitan–Marchenko equation from the methodology of the IST as follows. We have to treat P as more general operator as square matrix. We write
( P f ) ( k ) = C d λ ( l ) P ^ ( k , l ) f ( l ) ,
and the solution of KdV equation becomes
u = ϵ x C d λ ( k ) p ^ ( k ) [ ( I + ϵ P ) 1 p ( k ) .
We can write the solution as
u = 2 x K ( x , x ) ,
where
K ( x , y ) = ϵ 2 p T ( x ) [ I + ϵ P ( x ) ] 1 p ( y ) .
After some calculation we obtain
K ( x , y ) = ϵ 2 p T ( x ) p ( y ) ϵ 2 x d z K ( x , z ) p T ( z ) p ( y )
which is the GLM equation with
B ( x + y ) = ϵ 2 p T ( x ) p ( y ) = ϵ 2 d λ ( k ) exp i k 2 ( x + y ) + k 3 t .
Thus we have shown that SEsM is connected to the inverse scattering transform method for the Korteweg-de Vries equation. A similar connection can be proved also for other equations.

5. Several Particular Cases of SEsM and Their Connections with Other Methods

Below we show that many famous methods for obtaining particular exact solutions of nonlinear partial differential equations are connected to particular cases of SEsM. This is important to know in order to avoid errors when applying methods for obtaining exact solutions of nonlinear partial differential equations (for discussion see, e.g., [104,105]).

5.1. Homogeneous Balance Method and SEsM

The homogeneous balance method was discussed by Wang and et al. [106,107,108]. The method is as follows [107]. We consider the partial differential equation
P ( u , u x , u t , u x x , u x t , u t t , ) = 0 ,
where P is in general a polynomial function of its arguments, u = u ( x , t ) , and the subscripts denote the partial derivatives. A function w = w ( x , t ) is called a quasisolution of Equation (133), if there exists a function f = f ( w ) of a single variable so that a suitable linear combination of the following functions
f ( w ) , [ f ( w ) ] x , [ f ( w ) ] t , [ f ( w ) ] x x , [ f ( w ) ] x t , [ f ( w ) ] t t ,
is actually a solution of Equation (133). Four steps are needed to find f ( w ) . At step 1, one chooses a suitable linear combination of the functions from (134). The coefficients in this linear combination must be determined, so that the highest nonlinear terms and the highest order partial derivative terms in the given equation are both transformed into the polynomials with a highest equality degree in partial derivatives of w ( x , t ) in spite of f ( w ) and its various derivatives. These equal highest degrees determine the form of the linear combination. At step 2, after a substitution of the linear combination chosen in the first step into Equation (133), followed by a collection of all terms with the highest degree of derivatives of w(x,t) and setting its coefficient to zero, one obtains an ordinary differential equation for f ( w ) and then one solves it. At step 3, starting from the ODE and its solution obtained above, the nonlinear terms of various derivatives of f ( w ) in the relationship obtained in the second step can be replaced by the corresponding higher order derivatives of f ( w ) . Then by a collection of all terms with the same order derivatives of f ( w ) and by setting the coefficient of each order derivative of f ( w ) to zero respectively, one obtains a set of equations for w ( x , t ) . The left hand sides of these equations are the k degree homogeneous functions in various derivatives of w ( x , t ) , where k is the order of f ( k ) . In view of the homogeneous property of these equations one can expect that w ( x , t ) is an exponential function with some constants to be determined. Substituting the exponential function assumed into each k degree homogeneous equation in partial derivatives of w ( x , t ) , one obtains a set of nonlinear algebraic equations for some constants to be determined. If there exists a solution for these nonlinear algebraic equations, then w ( x , t ) and the coefficients of the linear combination chosen in the first step can be determined. At step 4, by substitution of f ( w ) , w ( x , t ) , and some constants obtained in the second and third steps into the combination chosen in the first step, and after doing some calculations, one obtains an exact solution of Equation (133).
Now we show that the homogeneous balance method is connected to a particular case of SEsM.
Assumption 2.
The homogeneous balance method is connected to the particular case of SEsM where there is no transformation of the nonlinearity of the equation (Step 1 of SEsM is skipped) and the function f ( w ) is the solution of a single simple equation (i.e., we have just one function f ( w ) at Step 2 of SEsM and the form of the function F at Step 2 of SEsM is a linear combination of the functions from (134)). The used simple equation is the equation for the function f ( w ) at Steps 3, 4, and 5 of SEsM. The balance in the homogeneous balance method is a particular case of the balance procedure from Step 6 of SEsM and the algebraic system is the same as the algebraic system from Step 7 of SEsM.
We start from SEsM and then impose restrictions in order to come to the particular case connected to the homogeneous balance method. Let us consider the partial differential equation (133). We skip Step 1 of SEsM (no transformation of the nonlinearity of the equation (133)) and come to the particular class of SEsM methodology without transformation of the nonlinearity of the solved equation. At Step 2 of SEsM, we take the particular case of a single function F which has particular form to be a linear combination of the functions from (134). Note that the functions from (134) depend on a single function w. This is another restriction on the function F from Step 2 of SEsM. At Step 3 of SEsM we make traveling-wave ansatz and thus we consider the particular case of SEsM when the solutions are travelling waves. We skip Step 4 of SEsM. The simple equation from Step 5 of SEsM is a particular case: this equation is chosen to be the differential equation for f ( w ) from the homogeneous balance method. The balance in the homogeneous balance method is a particular case of the balance procedure at Step 6 of SEsM for the case of balance of the highest powers in relationship containing monomials which are made of functions of a single variable. The algebraic system of the homogeneous balance method is the algebraic system from Step 7 of SEsM. Thus we have started from the general SEsM and by applying restrictions on it, we arrived at one particular case of SEsM connected to the homogeneous balance method.
Now let consider the homogeneous balance method in the version of Fan and Zhang [109]. They consider the equation (133) and search for a traveling wave solutions u ( x , t ) = u ( ξ ) = u ( x λ t ) of this equation. The solution is searched in the form
u ( ξ ) = i = 0 m a i v i ( ξ ) ,
where
d v d ξ = k ( 1 v 2 ) ,
and k and λ are parameters. The solution of (136) is
v = tanh ( k ξ ) ; v = coth ( k ξ ) .
Now let us show that this version of the homogeneous balance method is connected to a particular case of SEsM.
Assumption 3.
The homogeneous balance method in the version of Fan and Zhang is connected to a particular case of SEsM where, there is no transformation of the nonlinearity of the equation (Step 1 of SEsM is skipped), just one function u is used and the particular form (135) is used to relate the function u to the solution of the simple equation for the function v-Equation (136).
We start from SEsM and impose restrictions on it in order to reduce SEsM to the particular case connected to the homogeneous balance method in the version of Fan and Zhang. First, we consider the particular case of SEsM where no transformation of the nonlinearity of the solved equation is made. At Step 2 of SEsM, we consider again a particular case where a single function u is used. At Step 3 of SEsM we consider the particular case where we search for a traveling wave solution of the solved equation. At Step 4 of SEsM we assume a particular case of the relationship between the functions u and v, namely (135). At Step 5 of SEsM we again consider the particular case where the simple equation for v is (136). Steps 6 and 7 of SEsM follow and we may obtain an exact traveling wave solution of the solved equation. Thus we have started from SEsM and by considering particular cases of this methodology we have reduced it to a particular case connected to the homogeneous balance method in the version of Fan and Zhang.

5.2. Extended Homogeneous Balance Method and SEsM

The extended homogeneous balance method [110] extends the Fan and Zhang version of the homogeneous balance method by use of two simple equations of the same kind and of two traveling wave variables. Below we consider an assumption for this version of the homogeneous balance method and then we consider a version of the homogeneous balance method, which is based on an arbitrary number of simple equations and for any of these equations there is a separate traveling wave coordinate.
Assumption 4.
The extended homogeneous balance method in the version of El-Wakil et al. [110] is connected to particular case of SEsM for the case when there is no transformation of the nonlinearity of the equation (Step 1 of SEsM is skipped), two functions u and v is used and for these functions particular form (135) is used to relate them to the solution of two simple equations of the same kind
d ϕ d ξ 1 , 2 = a 1 , 2 ϕ 2 + c 1 , 2
where a 1 , 2 and c 1 , 2 are parameters and ξ 1 , 2 = α 1 , 2 x + β 1 , 2 t are two traveling wave coordinates.
We consider SEsM and impose restrictions on it in order to reduce SEsM to the particular case connected to the homogeneous balance method in the version of El-Wakil et al. We consider the particular case of SEsM where no transformation of the nonlinearity of the solved equation is made. At Step 2 of SEsM, we consider another particular case where two functions u and v will be used. At Step 3 of SEsM we consider the particular case where we search for a traveling wave solution of the solved equation. At Step 4 of SEsM we assume a particular case of the relationship between the functions u , v and ϕ , namely (135). At Step 5 of SEsM we again consider the particular case where the simple equations for ϕ are of the kind (138). Steps 6 and 7 of SEsM follow and we may obtain an exact traveling wave solution of the solved equation. Thus we have started from SEsM and by means of considering particular cases of this methodology we have reduced it to a particular case connected to the Extended homogeneous balance method in the version of El Wakil et al.
Now we are going to consider the much more extended homogeneous balance method (MMEHBM) and to show that this method is connected to a particular case of SEsM. In the MMEHBM we do not use the transformation of the nonlinearity of the solved equation. We use arbitrary number of functions u 1 , u 2 , . Any of these functions is represented by a power series of a function of corresponding traveling wave coordinate ξ 1 , ξ 2 , and for the last function there is a separate simple equation (all of these simple equations can be different).
Assumption 5.
The MMEHBM is connected to a particular case of SEsM where there is no transformation of the nonlinearity of the equation (Step 1 of SEsM is skipped), arbitrary number of functions u 1 , u 2 , are used, any of these functions is represented by a power series of a function of corresponding traveling wave coordinate ξ 1 , ξ 2 , and for the last function there is a separate simple equation (all of these simple equations can be different).
In general, MMEHBM will be used for obtaining exact traveling wave solutions of systems of nonlinear partial differential equations. We consider SEsM and impose restrictions on it in order to reduce SEsM to the particular case connected to MMEHBM. First, we consider the particular case of SEsM where no transformation of the nonlinearity of the solved equation is made at Step 1 of SEsM. At Step 2 of SEsM, we consider again a particular case where we use one function u i , ( i = 1 , 2 , ) for any of the i-th solved equations. At Step 3 of SEsM, we consider the particular case where we search for a traveling wave solution of the solved equation and the traveling waves associated with the different functions u i can have different traveling wave coordinates. At Step 4 of SEsM, we assume a particular case of the relationship between the functions u i and ϕ i , which is of the kind (135). We note that this is a large restriction on SEsM. At Step 5 of SEsM, we again consider separate simple equation for ϕ i ( ξ i ) which is of the kind (138). Steps 6 and 7 of SEsM follow and we may obtain an exact traveling wave solution of the solved equation. Thus we have started from SEsM and by means of considering particular cases of this methodology we have reduced it to the particular case connected to MMEHBM.

5.3. Auxiliary Equation Method and SEsM

The summary of the auxiliary equation method is as follows [111]. One considers nonlinear partial differential equation for u ( x , t ) in the form
H ( u , u x , u t , u x x , u x t , u t t , ) = 0 .
Then traveling wave variable ξ = x ω t is introduced and Equation (139) becomes
G ( u , u ξ , u ξ ξ , ) = 0 .
The solution of (140) is searched as
u ( ξ ) = i = 0 n a i z i ( ξ ) ,
where a i are parameters and z ( ξ ) is the solution of the auxiliary equation
d z d ξ 2 = a z ( ξ ) 2 + b z ( ξ ) 3 + c z ( ξ ) 4 .
Equations (141) and (142) are introduced in Equation (140) and then for the obtained relationship, one equates to 0 coefficients of the powers of z ( ξ ) . The solution of the obtained system of algebraic equations leads to solution of solved equation (139).
Now we show that the auxiliary equation method is connected to a particular case of SEsM.
Assumption 6.
The auxiliary equation method is connected to a particular case of SEsM where there is no transformation of the nonlinearity of the equation (Step 1 of SEsM is skipped); Function F at Step 2 of SEsM has particular form—(141)); just one simple equation is used that this simple equation has the form (142) which is a particular case of the form which can be used in SEsM.
We start from SEsM, impose restrictions on it and reduce SEsM to the particular case connected to the auxiliary equation method. At Step 1 of SEsM we do not transform the nonlinearity of the solved equation (we skip this step). At Step 2 of SEsM we use one of the possible forms of the function F—(141). We skip Steps 3 and 4 of SEsM as the function z ( ξ ) in the auxiliary equation method is directly connected to the solution of the used simple equation. Then we consider the particular case of SEsM when only 1 simple equation is used and we restrict further SEsM by the assumption that the simple equation is of the form (142). In such a way we reduce SEsM to a particular case connected to the auxiliary equation method.
Now let us show that a generalized version of the auxiliary equation method is connected to a particular case of SEsM. The generalization is in two directions. First, we can use many relationships for F instead of (141) and second we can use many relationships for the simple equations instead of (142). We call this generalized version of the auxiliary equation method GAEM—general auxiliary equation method. GAEM is formulated as follows. One considers nonlinear partial differential equation for u ( x , t ) in the form
H ( u , u x , u t , u x x , u x t , u t t , ) = 0 .
Then, traveling wave variable ξ = x ω t is introduced and Equation (139) becomes
G ( u , u ξ , u ξ ξ , ) = 0 .
The solution of (140) is searched in any form and not only in the form (141). Equation (141) is just a particular case of the possible forms that can be used in GAEM. z ( ξ ) can be solution of an arbitrary auxiliary equation. (142) is just one possibility for an auxiliary equation which can be used in GAEM. The introduction of the assumed form of the solution and of the assumed form of the auxiliary equation in (143) leads to a relationship consisting of coefficients multiplied by functions. We set to zero these coefficients and obtain a system of nonlinear algebraic equations. Any nontrivial solution of this system leads to an exact solution of (143).
Now we show that GAEM is connected to a particular case of SEsM.
Assumption 7.
The general auxiliary equation method (GAEM) is connected to a particular case of SEsM where there is no transformation of the nonlinearity of the equation (Step 1 of SEsM is skipped); Function F at Step 2 of SEsM has particular form and this is the corresponding form used by GAEM; just one simple equation is used that this simple equation has the particular form which is the corresponding form used by GAEM.
We start from SEsM, impose restrictions on it and reduce SEsM to a particular case connected to GAEM. At Step 1 of SEsM we do not transform the nonlinearity of the solved equation (we skip this step). At Step 2 of SEsM we use a possible form of the function F—the corresponding form used when GAEM is applied. This possible form is just one of the many forms that can be used in SEsM. We skip Steps 3 and 4 of SEsM as the function z ( ξ ) in the auxiliary equation method is directly connected to the solution of the used simple equation. Further, we consider the particular case of SEsM when only 1 simple equation is used and we restrict further SEsM by the assumption that this simple equation is of the form used by GAEM. In such a way we reduce SEsM to a particular case connected to the general auxiliary equation method.
An illustration of a particular case of GAEM can be seen in [112].

5.4. Jacobi Elliptic Function Expansion Method, F-Expansion Method and SESM

We show first that the Jacobi elliptic function expansion method (JEFEM) is connected to a particular case of SEsM. Then we describe general Jacobi elliptic function expansion method (GJEFEM) and show that it is connected to a particular case of SEsM. Finally, we list several methods used in the literature which are particular cases of GJEFEM.
The classic form of JEFEM is as given by [113]. One considers the following nonlinear partial differential equation for u ( x , t )
N ( u , u x , u t , u x x , u x t , u t t , ) = 0 ,
and searches for traveling wave solutions in the form
u = u ( ξ ) : ξ = k ( x c t ) ,
where k and c are parameters. u ( ξ ) is searched in the form of power series of the Jacobi elliptic function sn ( ξ , m ) where m is the modulus of the function sn ,
u ( ξ ) = j = 0 n a j sn ( j , m ) j .
This is a generalization of the tanh-method because for m = 1 sn ( ξ , 1 ) = tanh ( ξ ) . The substitution of (146) and (147) in (145) can lead to a system of nonlinear algebraic equations and any nontrivial solution of this system leads to an exact traveling wave solution of the solved Equation (145).
Assumption 8.
The Jacobi elliptic function expansion method (JEFEM) is connected to a particular case of SEsM where there is no transformation of the nonlinearity of the equation (Step 1 of SEsM is skipped); Function F at Step 2 of SEsM has particular form—(147)); just one simple equation is used that this simple equation is the differential equation for the elliptic function sn .
We start from SEsM, impose restrictions on it and reduce SEsM to the particular case connected to JEFEM. At Step 1 of SEsM we do not transform the nonlinearity of the solved equation (we skip this step). Additional restriction is that we search for a traveling wave solution of the solved equation (145). At Step 2 of SEsM we use a possible form of the function F—(147). This form is just one of the many forms that can be used in SEsM. We skip Steps 3 and 4 of SEsM as the function from (147) in the JEFEM is directly connected to the solution of the used simple equation which is the equation for the Jacobi elliptic function sn . The use of only one simple equation is a further restriction on SEsM. By means of all restrictions above, we reduce SEsM to the particular case connected to JEFEM.
Next we formulate general Jacobi elliptic function expansion method (GJEFEM). By this method we solve in general a system of N nonlinear partial differential equations and search for traveling wave solutions based on different coordinates ξ i = α i x β i t , i = 1 , 2 , , N . The solution is searched as function
u i ( ξ 1 , , ξ n ) = U i [ h 1 ( ξ 1 ) , , h N ( ξ N ) ]
of the functions f 1 , , f N and each of these functions is a solution of a differential equation for the Jacobi elliptic functions
d f i d ξ i 2 = a i f i 4 + b i f i 2 + c i .
We show below that the GJEFEM is connected to a particular case of SEsM.
Assumption 9.
The general Jacobi elliptic function expansion method is connected to a particular case of SEsM where, there is no transformation of the nonlinearity of the equation (Step 1 of SEsM is skipped); Functions u i at Step 2 of SEsM have particular form—(148)); and the simple equations are of the kind of the differential equation for the Jacobi elliptic functions (149).
We start from SEsM, impose restrictions on it and reduce SEsM to JEFEM. At Step 1. of SEsM we do not transform the nonlinearity of the solved equation (we skip this step). An additional restriction is that we search for a traveling wave solution of the solved Equation (145). At Step 2 of SEsM we use a possible form of the functions u i —(148). We skip Steps 3 and 4 of SEsM as the functions from (148) in the JEFEM are directly connected to the solution of the used simple equations which are of the kind of the differential equation for the Jacobi elliptic functions (149). This is an additional restriction on SEsM. By means of all restrictions above, we reduce SEsM to a particular case connected to GJEFEM.
Below, we list several particular cases of GJEFEM.
  • JEFEM is a particular case of GJEFEM for the case of just one solved nonlinear partial differential equation and when the simple equation is the equation for the Jacobi elliptic function sn and in addition the function U is a power series of the function sn .
  • Parks et al. [114] and Fu et al. [115] use expansions based on the elliptic functions cn , dn and cs . This is a particular case of GJEFEM when one simple equation is used and this simple equation is of the kind of (149).
  • Fan and Zhang [116] present an interesting application which is an extension of JEFEM for the case of two functions u 1 , 2 and single simple equation and by means of this extension they obtain solutions of the coupled Schrödinger-KdV system and of two-dimensional Davey–Stewartson equation. This extension of JEFEM is a particular case of GJEFEM when two functions u 1 , 2 are used with the same argument and when the simple equation is the differential equation for the elliptic function sn .
  • Another particular case of GJEFEM is applied by Yan [117] who treated a (2 + 1)-dimensional integrable Davey–Stewartson-type equation for the case of 2 spatial coordinates and travelling wave solutions. We note that SEsM allows treating equations with more that one spatial coordinate and the travelling waves can travel with different velocities which is a more general case than the case discussed by Yan where we have a single traveling wave despite the two spatial coordinates presented. Yan uses the following form of the function u i , i = 1 , 2 , 3
    u i ( ξ ) = a i 0 + j = 1 n f k j 1 ( ξ ) [ a i j f k ( ξ ) + b i j g k ( ξ ) ] ,
    where f k and g k , k = 1 , , 12 are Jacobi elliptic functions (i.e., are functions which satisfy the simple equation of the kind (149)). (150) is a particular form of the function U i from GJEFEM and the simple equations are equations for Jacobi elliptic functions as in GJEFEM.
  • Another particular case of GJEFEM is used in [118]. The simple equations used there are for Jacobi elliptic functions and the particular case of the used single function U is
    U = a 0 + i = 1 N sn 1 ( ξ , m ) [ a i sn ( ξ , m ) + b i cn ( ξ , m ) ] .
  • Liu and Fan [119] apply particular case of GJEFEM for the case of two spatial coordinates and time. These three variables are combined to produce a single traveling wave coordinate which allows the use of single variable simple equations. Wang et al. [120] use also a particular case of GJEFEM for the case of two spatial variables and time and combine all these variables in a single traveling wave variable. The new point in this article is the particular form of the functions U i
    U i = a i 0 + j = 1 m 1 a i j sn j ( ξ , m ) ( μ sn ( ξ , m ) + 1 ) j + b i j sn j 1 ( ξ , m ) cn ( ξ , m ) ( μ sn ( ξ , m ) + 1 ) j .
  • Ye at al. [121] extend (152) and use the following particular case for the functions U i
    U i = a i 0 + j = 1 m 1 a i , 2 j 1 sn j ( ξ , m ) ( μ sn ( ξ , m ) + μ 2 cn ( ξ , m ) + 1 ) j + a i , 2 j sn j 1 ( ξ , m ) cn ( ξ , m ) ( μ sn ( ξ , m ) + μ 2 cn ( ξ , m ) + 1 ) j .
    Other variants for U i are proposed by Wang et al. [122], Chen and Wang [123], Lü [124], Abdou and Elhanbaly [125], El-Sabbagh and Ali [126,127].
  • Another particular case of GJEFEM is the F-expansion method which have the same ideology as JEFEM but only the form of the simple equations for the Jacobi elliptic functions are not specified. In the different variants of the F- expansion method one uses different particular cases for the functions U i from GJEFEM [128,129,130,131].

5.5. Modified Simple Equation Method and SEsM

According to [132] the modified simple equation method is as follows. One considers the nonlinear partial differential equation which can be reduced to an ordinary partial differential equation for the function u ( z )
P ( u , u z , u z z , u z z z , ) = 0
(154) is solved by means of the ansatz
u ( z ) = k = 0 N A k ψ z ψ k .
Above A k are constants and A N 0 . The function ψ is a solution of some ordinary differential equation of a lesser order than (154) (called simplest equation) and solutions of these simplest equations are known. One uses the finite series (155) in order to represent the solution u through the solution of the simplest equation. In doing this, one has to determine the value of N by means of the balance of power of the leading terms of the relationship which is obtained after the substitution of (155) in (154). This relationship is polynomial of ψ z ψ and by setting to 0 of the coefficients to the powers of ψ z ψ one obtains a system of nonlinear algebraic equations which solution leads to an exact solution of (154).
Now we show that the modified simple equation method is connected to a particular case of SEsM.
Assumption 10.
The modified simple equation method is connected to a particular case of SEsM where there is no transformation of the nonlinearity of the equation (Step 1 of SEsM is skipped); Function F at Step 2 of SEsM has particular form—(155)) and just one simple equation is used.
We start from SEsM, impose restrictions on it and reduce SEsM to a particular case connected to the Modified Method of Simple Equation. At Step 1 of SEsM we do not transform the nonlinearity of the solved equation (we skip this step). An additional restriction is that we search for solution of the solved equation which depends on a single coordinate z—(154). At Step 2 of SEsM we use the form (155) of the function F— this form is just one of the many forms that can be used in SEsM. We skip Steps 3 and 4 of SEsM as the function from (155) in the JEFEM is directly connected to the solution of the used simple equation which in this case is called the simplest equation. The use of only one simple equation is a further restriction on SEsM. By means of all restrictions above, we reduce SEsM to a particular case connected to the Modified Method of Simple Equation.

5.6. Trial Function Method and SEsM

The trial function method is described in [133,134] and it is as follows. We consider a nonlinear partial differential equation
N ( u , u x , u t , u x x , u x t , u t t , ) = 0 ,
and take a trial function y ( x , t ) in order to construct a solution u ( y ) of (156). The we substitute u ( y ) in (156) and determine the parameters of the solution.
Assumption 11.
The trial function method is connected to a particular case of SEsM where there is no transformation of the nonlinearity of the equation (Step 1 of SEsM is skipped); Function F at Step 2 of SEsM has particular form— u ( y ) , where y (the trial function) is the solution of a single simple equation.
We start from SEsM, impose restrictions on it and reduce SEsM to a particular case connected to the trial function method. At Step 1 of SEsM we do not transform the nonlinearity of the solved equation (we skip this step). An additional restriction is that we search for a solution of the solved equation which depends on a single coordinate which can be a traveling wave coordinate or another kind of coordinate. At Step 2 of SEsM we use a particular form of the function F which is constructed by means of trial function. In most cases, F is presented by a finite power series of the trial function. The trial function is a solution to one simple equation. Thus, by means of the restrictions above, we reduce SEsM to a particular case connected to the trial function method.

5.7. First Integral Method and SEsM

The first integral method for obtaining exact solutions of nonlinear partial differential equations is as follows [135]. One wants to obtain exact solutions of the nonlinear partial differential equation
P ( u , u x , u t , u x x , u x t , u t t , ) = 0 .
In order to do this one converts (157) to an ordinary differential equation by the traveling wave ansatz u ( x , t ) = U ( z ) = u ( k x ω t ) . Then one introduces X = U and Y = U z and writes (157) as system of equations
Y = X z
Y z = F ( X , Y ) .
The solution is obtained by the assumption that the derivative of the relationship Q ( X , Y ) = i = 0 m a i ( X ) Y i can be represented as
d Q d z = [ g ( X ) + h ( X ) Y ] i = 0 m a i ( X ) Y i
which together with (159) allows computation of the solution.
Assumption 12.
The first integral method is connected to a particular case of SEsM for the case when equations of the kind
X z z = F ( X , X z )
are considered, there is no transformation of the nonlinearity of the equation (Step 1 of SEsM is skipped); one simplest equation is used and this simplest equation is determined by the condition (160).
We stress that the first integral method can be applied to the restricted class of equations (161). This restricted class is obtained from (159) by substitution of (158) there. We start from SEsM, impose restrictions on it and reduce SEsM to the trial function method. At Step 1 of SEsM we do not transform the nonlinearity of the solved equation (we skip this step). An additional restriction is that we search for a solution of the solved equation that depends on a single coordinate which can be a traveling wave coordinate or another kind of coordinate. (160) imposes a further restriction on X and plays the role of the implicit simple equation which together with (161) determine the solution of (157). In this process one has to use a polynomial form of a i ( X ) and to determine the coefficients of these polynomials similar to Steps 6 and 7 of SEsM. Thus first integral method is connected to a particular case of SEsM for obtaining solutions for the limited class of equations (159) under the assumption that (160) holds.

6. Concluding Remarks

In this text, we present the Simple Equations Method (SEsM) for obtaining exact solutions of nonlinear differential equations and discuss its connections with other methods for obtaining exact solutions of nonlinear differential equations. Special importance is given to the connection to two of the most famous method for obtaining exact solutions of integrable nonlinear partial differential equations - method of Hirota and the inverse scattering transform method. We show that one particular case of SEsM is connected to the Hirota method. This means that the SEsM is a useful tool for obtaining exact solutions of nonlinear integrable differential equations and can lead to multisoliton solutions of such equations. Then we show that SEsM can be connected to the famous inverse scattering transform method. This is done in the case of the Korteweg-de Vries equation and much more research on this connection will be presented in the future. In addition to solutions of integrable nonlinear differential equations SEsM can lead also to exact particular solutions of nonintegrable nonlinear differential equations. Actually, the simplest version of SEsM (called Modified Method of Simplest Equation and based on representation of the searched solution as power series of the solution of single simple equation [85,86]) was intended exactly to solve the problem of obtaining particular exact solutions of nonintegrable nonlinear differential equations. Finally, we show that many other methods for obtaining particular exact solutions of nonlinear nonintegrable differential equations are connected to particular cases of SEsM. Many of these particular cases are discussed above in the text. Some more cases are listed in [74]. We note that additional discussion of other particular cases of the SEsM will be presented elsewhere and this is connected to many important questions about the methodology for obtaining exact solutions of nonlinear differential equations raised in [104,105].
We believe that the SEsM algorithm is useful one and it presents a fruitful way for obtaining solutions of nonlinear differential equations. The development of SEsM is far from finished. More research on the possible classes of transformations can be done as well as additional simple equations can be used. We are sure that over the course of the years this algorithm will be further refined. Our intention to go in this direction and we are going to extend and to apply SEsM for obtaining solutions of many nonlinear differential equations of practical significance for natural and social sciences.

Author Contributions

Conceptualization, N.K.V., Z.I.D.; methodology, N.K.V.; software, K.N.V.; validation, Z.I.D., K.N.V.; formal analysis, N.K.V., Z.I.D., K.N.V.; resources, N.K.V.; data curation, Z.I.D., K.N.V.; writing—original draft preparation, N.K.V., Z.I.D. K.N.V.; writing—review and editing, N.K.V.; visualization, K.N.V.; supervision, N.K.V., Z.I.D.; project administration, Z.I.D., K.N.V.; funding acquisition, N.K.V. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially supported by the project BG05 M2OP001-1.001-0008 “National Center for Mechatronics and Clean Technologies”, funded by the Operating Program “Science and Education for Intelligent Growth” of Republic of Bulgaria and by the National Scientific Program “Information and Communication Technologies for a Single Digital Market in Science, Education and Security” (ICTinSES), contract No D01205/23.11.2018, financed by the Ministry of Education and Science in Bulgaria.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Axelrod, R.; Cohen, M. Harnessing Complexity; Basic Books: New York, NY, USA, 2001; ISBN 978-0465005505. [Google Scholar]
  2. Chian, A.C.-L. Complex Systems Approach to Economic Dynamics; Springer: Berlin, Germany, 2007; ISBN 978-3-540-39752-6. [Google Scholar]
  3. Chen, W.-K. Theory of Nets. Flows in Networks; Imperial College Press: London, UK, 2003; ISBN 978-0471851486. [Google Scholar]
  4. Lucas, M.W. Network Flow Analysis; No Starch Press: San Francisco, CA, USA, 2010; ISBN 978-1593272036. [Google Scholar]
  5. Vitanov, N.K. Science Dynamics and Research Production. Indicators, Indexes, Statistical Laws and Mathematical Models; Springer: Cham, Switzerland, 2016; ISBN 978-3-319-41629-8. [Google Scholar]
  6. Treiber, M.; Kesting, A. Traffic Flow Dynamics: Data, Models, and Simulation; Springer: Berlin, Germany, 2013; ISBN 978-3-642-32460-4. [Google Scholar]
  7. Castellani, B.; Hafferty, F. Sociology and Complexity Science; Springer: Berlin, Germany, 2009; ISBN 978-3-540-88462-0. [Google Scholar]
  8. Bahrami, M.; Chinichian, N.; Hosseiny, J.A.G.; Ausloos, M. Optimization of the Post- Crisis Recovery Plans in Scale-Free Networks. Physica A 2020, 540, 123203. [Google Scholar] [CrossRef] [Green Version]
  9. May, R.M.; Levin, S.A.; Sugihara, G. Complex Systems: Ecology for Bankers. Nature 2008, 451, 893–895. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Ivanova, K.; Ausloos, M. Application of the Detrended Fluctuation Analysis (DFA) Method for Describing Cloud Breaking. Physica A 1999, 274, 349–354. [Google Scholar] [CrossRef]
  11. Vitanov, N.K.; Vitanov, K.N. Discrete-time Model for a Motion of Substance in a Channel of a Network with Application to Channels of Human Migration. Physica A 2018, 509, 635–650. [Google Scholar] [CrossRef] [Green Version]
  12. Vitanov, N.K.; Vitanov, K.N. Statistical Distributions Connected to Motion of Substance in a Channel of a Network. Physica A 2019, 527, 121174. [Google Scholar] [CrossRef]
  13. Brehmer, B. Dynamic Decision Making: Human Control of Complex Systems. Acta Psychol. 1992, 81, 211–241. [Google Scholar] [CrossRef]
  14. Blossfeld, H.-P.; Rohwer, G. Techniques of Event History Modeling: New Approaches to Casual Analysis. J. R. Stat. Soc. Ser. D (Stat.) 2003, 52, 236–238. [Google Scholar] [CrossRef]
  15. Vitanov, N.K.; Vitanov, K.N. On the Motion of Substance in a Channel of a Network and Human Migration. Physica A 2018, 490, 1277–1294. [Google Scholar] [CrossRef]
  16. Lambiotte, R.; Ausloos, M. Coexistence of Opposite Opinions in a Network with Communities. J. Stat. Mech. Theory Exp. 2007, P08026. [Google Scholar] [CrossRef] [Green Version]
  17. Ausloos, M. Statistical Physics in Foreign Exchange Currency and Stock Markets. Phys. Stat. Mech. Its Appl. 2000, 285, 48–65. [Google Scholar] [CrossRef]
  18. Kutner, R.; Ausloos, M.; Grech, D.; Di Matteo, T.; Schinckus, C.; Stanley, H.E. Manifesto for a Post-Pandemic Modeling. Physica A 2019, 516, 240–253. [Google Scholar] [CrossRef] [Green Version]
  19. Simon, J.H. The Economic Consequences of Immigration; The University of Michigan Press: Ann Arbor, MI, USA, 1999; ISBN 978-0472086160. [Google Scholar]
  20. Vitanov, N.K.; Ausloos, M.; Rotundo, G. Discrete Model of Ideological Struggle Accounting for Migration. Adv. Complex Syst. 2012, 15 (Suppl. 1), 1250049. [Google Scholar] [CrossRef]
  21. Vitanov, N.K.; Vitanov, K.N. Box Model of Migration Channels. Math. Soc. Sci. 2016, 80, 108–114. [Google Scholar] [CrossRef]
  22. Amaral, L.A.N.; Scala, A.; Barthelemy, M.; Stanley, H.E. Classes of Small-World Networks. Proc. Natl. Acad. Sci. USA 2000, 97, 11149–11152. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Hall, D.A. Nonlinearity in Piezoelectric Ceramics. J. Mater. Sci. 2001, 36, 4575–4601. [Google Scholar] [CrossRef]
  24. Drazin, P.G. Nonlinear Systems; Cambridge University Press: Cambridge, UK, 1992; ISBN 0-521-40489-4. [Google Scholar]
  25. Pikovsky, A.S.; Shepelyansky, D.L. Destruction of Anderson Localization by a Weak Nonlinearity. Phys. Rev. Lett. 2008, 100, 094101. [Google Scholar] [CrossRef] [Green Version]
  26. Dimitrova, Z.I. Numerical Investigation of Nonlinear Waves Connected to Blood Flow in an Elastic Tube with Variable Radius. J. Theor. Appl. Mech. 2015, 45, 79–92. [Google Scholar] [CrossRef] [Green Version]
  27. Boeck, T.; Vitanov, N.K. Low-Dimensional Chaos in Zero-Prandtl-Number Benard–Marangoni Convection. Phys. Rev. E 2002, 65, 037203. [Google Scholar] [CrossRef]
  28. Dimitrova, Z. On Traveling Waves in Lattices: The Case of Riccati Lattices. J. Theor. Appl. Mech. 2012, 42, 3–22. [Google Scholar] [CrossRef]
  29. Kawasaki, K.; Ohta, T. Kink Dynamics in One-Dimensional Nonlinear Systems. Physica A 1982, 116, 573–593. [Google Scholar] [CrossRef]
  30. Kantz, H.; Schreiber, T. Nonlinear Time Series Analysis; Cambridge University Press: Cambridge, UK, 2004; ISBN 978-0511755798. [Google Scholar]
  31. Ashenfelter, K.T.; Boker, S.M.; Waddell, J.R.; Vitanov, N. Spatiotemporal Symmetry and Multifractal Structure of Head Movements During Dyadic Conversation. J. Exp. Psychol. Hum. Percept. Perform. 2009, 35, 1072–1091. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Verhulst, F. Nonlinear Differential Equations and Dynamical Systems; Springer: Berlin, Germany, 2006; ISBN 978-3-540-60934-6. [Google Scholar]
  33. Kantz, H.; Holstein, D.; Ragwitz, M.; Vitanov, N.K. Markov Chain Model for Turbulent Wind Speed Data. Physica A 2004, 342, 315–321. [Google Scholar] [CrossRef]
  34. Ernst, M.H. Nonlinear Model-Boltzmann Equations and Exact Solutions. Phys. Rep. 1981, 78, 1–171. [Google Scholar] [CrossRef]
  35. Benjamin, T.B.; Bona, J.L.; Mahony, J.J. Model Equations for Long Waves in Nonlinear Dispersive Systems. Philos. Trans. R. Soc. Lond. Ser. A 1972, 272, 47–78. [Google Scholar] [CrossRef]
  36. Vitanov, N.K. Upper Bounds on the Heat Transport in a Porous Layer. Physica D 2000, 136, 322–339. [Google Scholar] [CrossRef]
  37. Dimitrova, Z.I.; Vitanov, N.K. Influence of Adaptation on the Nonlinear Dynamics of a System of Competing Populations. Phys. Lett. A 2000, 272, 368–380. [Google Scholar] [CrossRef]
  38. Dimitrova, Z.I.; Vitanov, N.K. Adaptation and its Impact on the Dynamics of a System of Three Competing Populations. Physica A 2001, 300, 91–115. [Google Scholar] [CrossRef]
  39. Struble, R. Nonlinear Differential Equations; Dover: New York, NY, USA, 2018; ISBN 978-0486817545. [Google Scholar]
  40. Dimitrova, Z.I.; Vitanov, N.K. Chaotic Pairwise Competition. Theor. Popul. Biol. 2004, 66, 1–12. [Google Scholar] [CrossRef]
  41. Vitanov, N.K.; Dimitrova, Z.I.; Ausloos, M. Verhulst-Lotka-Volterra Model of Ideological Struggle. Physica A 2010, 389, 4970–4980. [Google Scholar] [CrossRef] [Green Version]
  42. Taulbee, D.B. An Improved Algebraic Reynolds Stress Model and Corresponding Nonlinear Stress Model. Phys. Fluids A 1992, 4, 2555–2561. [Google Scholar] [CrossRef]
  43. Borisov, R.; Dimitrova, Z.I.; Vitanov, N.K. Statistical Characteristics of Stationary Flow of Substance in a Network Channel Containing Arbitrary Number of Arms. Entropy 2020, 22, 553. [Google Scholar] [CrossRef] [PubMed]
  44. Vitanov, N.K.; Vitanov, K.N.; Kantz, H. On the Motion of Substance in a Channel of a Network: Extended Model and New Classes of Probability Distributions. Entropy 2020, 22, 1240. [Google Scholar] [CrossRef]
  45. Grossberg, S. Nonlinear Neural Networks: Principles, Mechanisms, and Architectures. Neural Netw. 1981, 1, 17–61. [Google Scholar] [CrossRef] [Green Version]
  46. Hopf, E. The Partial Differential Equation: ut + uux = ϵuxx. Commun. Pure Appl. Math. 1950, 3, 201–230. [Google Scholar] [CrossRef]
  47. Cole, J.D. On a Quasi-Linear Parabolic Equation Occurring in Aerodynamics. Q. Appl. Math. 1951, 9, 225–236. [Google Scholar] [CrossRef] [Green Version]
  48. Gardner, C.S.; Greene, J.M.; Kruskal, M.D.; Miura, R.R. Method for Solving the Korteweg-de Vries Equation. Phys. Rev. Lett. 1967, 19, 1095–1097. [Google Scholar] [CrossRef]
  49. Ablowitz, M.J.; Kaup, D.J.; Newell, A.C.; Segur, H. The Inverse Scattering Transform -Fourier Analysis for nonlinear problems. Stud. Appl. Math. 1974, 53, 249–315. [Google Scholar] [CrossRef]
  50. Ablowitz, M.J.; Clarkson, P.A. Solitons, Nonlinear Evolution Equations and Inverse Scattering; Cambridge University Press: Cambridge, UK, 1991; ISBN 978-0511623998. [Google Scholar]
  51. Hirota, R. Exact Solution of the Korteweg—De Vries Equation for Multiple Collisions of Solitons. Phys. Rev. Lett. 1971, 27, 1192–1194. [Google Scholar] [CrossRef]
  52. Hirota, R. The Direct Method in Soliton Theory; Cambridge University Press: Cambridge, UK, 2004; ISBN 978-0511543043. [Google Scholar]
  53. Tabor, M. Chaos and Integrability in Dynamical Systems; Wiley: New York, NY, USA, 1989; ISBN 978-0471827283. [Google Scholar]
  54. Carrielo, F.; Tabor, M. Painleve Expansions for Nonintegrable Evolution Equations. Physica D 1989, 39, 77–94. [Google Scholar] [CrossRef]
  55. Carrielo, F.; Tabor, M. Similarity Reductions from Extended Painleve Expansions for Nonintegrable Evolution Equations. Physica D 1991, 53, 59–70. [Google Scholar] [CrossRef]
  56. Weiss, J.; Tabor, M.; Carnevalle, G. The Painleve Property for Partial Differential Equations. J. Math. Phys. 1983, 24, 522–526. [Google Scholar] [CrossRef]
  57. Kudryashov, N.A. On Types of Nonlinear Nonintegrable Equations with Exact Solutions. Phys. Lett. A 1991, 155, 269–275. [Google Scholar] [CrossRef]
  58. Kudryashov, N.A. Simplest Equation Method to Look for Exact Solutions of Nonlinear Differential Equations. Chaos Solitons Fractals 2005, 24, 1217–1231. [Google Scholar] [CrossRef] [Green Version]
  59. Kudryashov, N.A.; Loguinova, N.B. Extended Simplest Equation Method for Nonlinear Differential Equations. Appl. Math. Comput. 2008, 205, 361–365. [Google Scholar] [CrossRef]
  60. Kudryashov, N.A. Partial Differential Equations with Solutions Having Movable First-Order Singularities. Phys. Lett. A 1992, 169, 237–242. [Google Scholar] [CrossRef]
  61. Kudryashov, N.A. Exact Solitary Waves of the Fisher Equation. Phys. Lett. A 2005, 342, 99–106. [Google Scholar] [CrossRef]
  62. Kudryashov, N.A. One Method for Finding Exact Solutions of Nonlinear Differential Equations. Commun. Nonlinear Sci. Numer. Simul. 2012, 17, 2248–2253. [Google Scholar] [CrossRef] [Green Version]
  63. Kudryashov, N.A. Exact Soliton Solutions of the Generalized Evolution Equation of Wave Dynamics. J. Appl. Math. Mech. 1988, 52, 361–365. [Google Scholar] [CrossRef]
  64. Kudryashov, N.A. Exact Solutions of Nonlinear Wave Equations Arising in Mechanics. J. Appl. Math. Mech. 1990, 54, 372–375. [Google Scholar] [CrossRef]
  65. Kudryashov, N.A. Exact Solutions and Integrability of the Duffing—Van der Pol Equation. Regul. Chaotic Dyn. 2018, 23, 471–479. [Google Scholar] [CrossRef]
  66. Kudryashov, N.A. Exact Solutions of the Equation for Surface waves in a Convecting Fluid. Appl. Math. Comput. 2019, 344–345, 97–106. [Google Scholar] [CrossRef]
  67. Kudryashov, N.A. A Generalized Model for Description of Propagation Pulses in Optical Fiber. Optik 2019, 189, 42–52. [Google Scholar] [CrossRef]
  68. Kudryashov, N.A. First Integrals and Solutions of the Traveling Wave Reduction for the Triki–Biswas Equation. Optik 2019, 185, 275–281. [Google Scholar] [CrossRef]
  69. Kudryashov, N.A. Highly Dispersive Optical Solitons of the Generalized Nonlinear Eighth-Order Schrödinger Equation. Optik 2020, 206, 164335. [Google Scholar] [CrossRef]
  70. Kudryashov, N.A. The Generalized Duffing Oscillator. Commun. Nonlinear Sci. Numer. Simul. 2021, 93, 105526. [Google Scholar] [CrossRef]
  71. Vitanov, N.K. Recent Developments of the Methodology of the Modified Method of Simplest Equation with Application. Pliska Stud. Math. Bulg. 2019, 30, 29–42. [Google Scholar]
  72. Vitanov, N.K. Modified Method of Simplest Equation for Obtaining Exact Solutions of Nonlinear Partial Differential Equations: History, recent development and studied classes of equations. J. Theor. Appl. Mech. 2019, 49, 107–122. [Google Scholar] [CrossRef]
  73. Vitanov, N.K. The Simple Equations Method (SEsM) For Obtaining Exact Solutions Of Nonlinear PDEs: Opportunities Connected To The Exponential Functions. AIP Conf. Proc. 2019, 2159, 030038. [Google Scholar] [CrossRef]
  74. Vitanov, N.K.; Dimitrova, Z.I. Simple Equations Method (SEsM) and Other Direct Methods for Obtaining Exact Solutions of Nonlinear PDEs. AIP Conf. Proc. 2019, 2159, 030039. [Google Scholar] [CrossRef]
  75. Martinov, N.; Vitanov, N. On the Correspondence Between the Self-consistent 2D Poisson-Boltzmann Structures and the Sine- Gordon Waves. J. Phys. A Math. Gen. 1992, 25, L51–L56. [Google Scholar] [CrossRef]
  76. Martinov, N.; Vitanov, N. On Some Solutions of the Two-Dimensional Sine-Gordon Equation. J. Phys. A Math. Gen. 1992, 25, L419–L426. [Google Scholar] [CrossRef]
  77. Martinov, N.K.; Vitanov, N.K. New Class of Running-Wave Solutions of the (2 + 1)-Dimensional Sine-Gordon Equation. J. Phys. A Math. Gen. 1994, 27, 4611–4618. [Google Scholar] [CrossRef]
  78. Martinov, N.K.; Vitanov, N.K. On Self-Consistent Thermal Equilibrium Structures in Two-Dimensional Negative-Temperature Systems. Can. J. Phys. 1994, 72, 618–624. [Google Scholar] [CrossRef]
  79. Vitanov, N.K. On Travelling Waves and Double-Periodic Structures in Two-Dimensional Sine—Gordon Systems. J. Phys. A Math. Gen. 1996, 29, 5195–5207. [Google Scholar] [CrossRef]
  80. Vitanov, N.K. Breather and Soliton Wave Families for the Sine-Gordon Equation. Proc. R. Soc. Lond. A 1998, 454, 2409–2423. [Google Scholar] [CrossRef]
  81. Vitanov, N.K.; Jordanov, I.P.; Dimitrova, Z.I. On Nonlinear Dynamics of Interacting Populations: Coupled Kink Waves in a System of Two Populations. Commun. Nonlinear Sci. Numer. Simul. 2009, 14, 2379–2388. [Google Scholar] [CrossRef]
  82. Vitanov, N.K.; Jordanov, I.P.; Dimitrova, Z.I. On Nonlinear Population Waves. Appl. Math. Comput. 2009, 215, 2950–2964. [Google Scholar] [CrossRef]
  83. Vitanov, N.K. Application of Simplest Equations of Bernoulli and Riccati Kind for Obtaining Exact Traveling-Wave Solutions for a Class of PDEs with Polynomial Nonlinearity. Commun. Nonlinear Sci. Numer. Simul. 2010, 15, 2050–2060. [Google Scholar] [CrossRef]
  84. Vitanov, N.K.; Dimitrova, Z.I. Application of The Method of Simplest Equation for Obtaining Exact Traveling-Wave Solutions for Two Classes of Model PDEs from Ecology and Population Dynamics. Commun. Nonlinear Sci. Numer. Simul. 2010, 15, 2836–2845. [Google Scholar] [CrossRef]
  85. Vitanov, N.K.; Dimitrova, Z.I.; Kantz, H. Modified Method of Simplest Equation and its Application to Nonlinear PDEs. Appl. Math. Comput. 2010, 216, 2587–2595. [Google Scholar] [CrossRef]
  86. Vitanov, N.K. Modified Method of Simplest Equation: Powerful Tool for Obtaining Exact and Approximate Traveling-Wave Solutions of Nonlinear PDEs. Commun. Nonlinear Sci. Numer. Simul. 2011, 16, 1176–1185. [Google Scholar] [CrossRef]
  87. Vitanov, N.K.; Dimitrova, Z.I.; Vitanov, K.N. On the Class of Nonlinear PDEs That Can be Treated by the Modified Method of Simplest Equation. Application to Generalized Degasperis–Processi Equation and B-Equation. Commun. Nonlinear Sci. Numer. Simul. 2011, 16, 3033–3044. [Google Scholar] [CrossRef]
  88. Vitanov, N.K. On Modified Method of Simplest Equation for Obtaining Exact and Approximate Solutions of Nonlinear PDEs: The Role of the Simplest Equation. Commun. Nonlinear Sci. Numer. Simul. 2011, 16, 4215–4231. [Google Scholar] [CrossRef]
  89. Vitanov, N.K. On Modified Method of Simplest Equation for Obtaining Exact Solutions of Nonlinear PDEs: Case of Elliptic Simplest Equation. Pliska Stud. Math. Bulg. 2012, 21, 257–266. [Google Scholar]
  90. Vitanov, N.K.; Dimitrova, Z.I.; Kantz, H. Application of the Method of Simplest Equation for Obtaining Exact Traveling-Wave Solutions for the Extended Korteweg–de Vries Equation and Generalized Camassa–Holm Equation. Appl. Math. Comput. 2013, 219, 7480–7492. [Google Scholar] [CrossRef]
  91. Vitanov, N.K.; Dimitrova, Z.I.; Vitanov, K.N. Traveling Waves and Statistical Distributions Connected to Systems of Interacting Populations. Comput. Math. Appl. 2013, 66, 1666–1684. [Google Scholar] [CrossRef]
  92. Vitanov, N.K.; Vitanov, K.N. Population Dynamics in Presence of State Dependent Fluctuations. Comput. Math. Appl. 2013, 68, 962–971. [Google Scholar] [CrossRef]
  93. Vitanov, N.K.; Dimitrova, Z.I. Solitary Wave Solutions for Nonlinear Partial Differential Equations that Contain Monomials of Odd and Even Grades with Respect to Participating Derivatives. Appl. Math. Comput. 2014, 247, 213–217. [Google Scholar] [CrossRef]
  94. Vitanov, N.K.; Dimitrova, Z.I.; Vitanov, K.N. Modified Method of Simplest Equation for Obtaining Exact Analytical Solutions of Nonlinear Partial Differential Equations: Further Development of the Methodology with Applications. Appl. Math. Comput. 2015, 269, 363–378. [Google Scholar] [CrossRef] [Green Version]
  95. Vitanov, N.K.; Dimitrova, Z.I.; Ivanova, T.I. On Solitary Wave Solutions of a Class of Nonlinear Partial Differential Equations Based on the Function 1/cosh(αx + βt)n. Appl. Math. Comput. 2017, 315, 372–380. [Google Scholar] [CrossRef]
  96. Vitanov, N.K.; Dimitrova, Z.I. Modified Method of Simplest Equation Applied to the Nonlinear Schrödinger Equation. J. Theor. Appl. Mech. Sofia 2018, 48, 59–68. [Google Scholar] [CrossRef] [Green Version]
  97. Nikolova, E.V.; Jordanov, I.P.; Dimitrova, Z.I.; Vitanov, N.K. Evolution of Nonlinear Waves in a Blood-Filled Artery with an Aneurysm. AIP Conf. Proc. 2017, 1895, 07002. [Google Scholar] [CrossRef] [Green Version]
  98. Jordanov, I.P.; Vitanov, N.K. On the Exact Traveling Wave Solutions of a Hyperbolic Reaction- Diffusion Equation. Stud. Comput. Intell. 2019, 793, 199–210. [Google Scholar] [CrossRef]
  99. Hereman, W.; Zhuang, W. Symbolic Computation of Solitons via Hirota’s Bilinear Method, Preprint; Department of Mathematical and Computer Sciences Colorado School of Mines: Golden, CO, USA, 1994. [Google Scholar]
  100. Ishimori, Y. Multi-Vortex Solutions of a Two-Dimensional Nonlinear Wave Equation. Prog. Theor. Phys. 1984, 72, 33–37. [Google Scholar] [CrossRef] [Green Version]
  101. Vijayalakshmi, S.; Lakshmanan, M. Localized Coherent Structures of Ishimori Equation I through Hirota’s Bilinearization Method: Time Dependent/Stationary Boundaries. Chaos Solitons Fractals 2007, 33, 203–216. [Google Scholar] [CrossRef] [Green Version]
  102. Lakshmanan, M.; Daniel, M. On the Evolution of Higher Dimensional Heisenberg Continuum Spin Systems. Physica A 1981, 107, 533–552. [Google Scholar] [CrossRef]
  103. Rosales, R. Exact Solutions of Some Nonlinear Evolution Equations. Stud. Appl. Math. 1978, 59, 117–151. [Google Scholar] [CrossRef]
  104. Kudryashov, N.A. Seven Common Errors in Finding Exact Solutions of Nonlinear Differential Equations. Commun. Nonlinear Sci. Numer. Simul. 2009, 14, 3507–3529. [Google Scholar] [CrossRef] [Green Version]
  105. Kudryashov, N.A. A Note on the G’/G-expansion method. Appl. Math. Comput. 2010, 217, 1755–1758. [Google Scholar] [CrossRef]
  106. Wang, M.-L. Solitary Wave Solutions for Variant Boussinesq equations. Phys. Lett. A 1995, 199, 169–172. [Google Scholar] [CrossRef]
  107. Wang, M.-L. Exact Solutions for a Compound KdV-Burgers Equation. Phys. Lett. A 1996, 213, 279–287. [Google Scholar] [CrossRef]
  108. Wang, M.-L.; Zhou, Y.; Li, Z. Application of a Homogeneous Balance Method to Exact Solutions of Nonlinear Equation in Mathematical Physics. Phys. Lett. A 1996, 216, 67–75. [Google Scholar] [CrossRef]
  109. Fan, E.; Zhang, H. A Note on the Homogeneous Balance Method. Phys. Lett. A 1998, 246, 403–406. [Google Scholar] [CrossRef]
  110. El-Wakil, S.A.; Abulwafa, E.M.; Elhanbaly, A.; Abdou, M.A. The Extended Homogeneous Balance Method and its Applications for a Class of Nonlinear Evolution Equations. Chaos Solitons Fractals 2007, 33, 1512–1522. [Google Scholar] [CrossRef]
  111. Sirendaoreji, J.S. Auxiliary Equation Method for Solving Nonlinear Partial Differential Equations. Phys. Lett. A 2003, 309, 387–396. [Google Scholar] [CrossRef]
  112. Kilic, B.; Bulut, H. A New Method with a Different Auxiliary Equation to Obtain Solitary Wave Solutions for Nonlinear Partial Differential Equations. Adv. Math. Phys. 2013, 890784. [Google Scholar] [CrossRef] [Green Version]
  113. Liu, S.; Fu, Z.; Liu, S.; Zhao, Q. Jacobi Elliptic Function Expansion Method and Periodic Wave Solutions of Nonlinear Wave Equations. Phys. Lett. A 2001, 289, 69–74. [Google Scholar] [CrossRef]
  114. Parkes, E.J.; Duffy, B.R.; Abbott, P.C. The Jacobi Elliptic-Function Method for Finding Periodic-Wave Solutions to Nonlinear Evolution Equations. Phys. Lett. A 2002, 295, 280–286. [Google Scholar] [CrossRef]
  115. Fu, Z.; Liu, S.; Liu, S.; Zhao, Q. New Jacobi Elliptic Function Expansion and New Periodic Solutions of Nonlinear Wave Equations. Phys. Lett. A 2001, 290, 72–76. [Google Scholar] [CrossRef]
  116. Fan, E.; Zhang, J. Applications of the Jacobi Elliptic Function Method to Special-Type Nonlinear Equations. Phys. Lett. A 2002, 305, 383–392. [Google Scholar] [CrossRef]
  117. Yan, Z. Abundant Families of Jacobi Elliptic Function Solutions of the (2 + 1)-Dimensional Integrable Davey–Stewartson-Type Equation Via a New Method. Chaos Solitons Fractals 2003, 18, 299–309. [Google Scholar] [CrossRef]
  118. Shen, S.; Pan, Z. A Note on the Jacobi Elliptic Function Expansion Method. Phys. Lett. A 2003, 308, 143–148. [Google Scholar] [CrossRef]
  119. Liu, G.-T.; Fan, T.-Y. New Applications of Developed Jacobi Elliptic Function Expansion Methods. Phys. Lett. A 2005, 345, 161–166. [Google Scholar] [CrossRef]
  120. Wang, Q.; Chen, Y.; Hongqing, Z. A New Jacobi Elliptic Function Rational Expansion Method and its Application to (1+1)-Dimensional Dispersive Long Wave Equation. Chaos Solitons Fractals 2005, 23, 477–483. [Google Scholar] [CrossRef]
  121. Yu, Y.; Wang, Q.; Zhang, H. The Extended Jacobi Elliptic Function Method to Solve a Generalized Hirota–Satsuma Coupled KdV Equations. Chaos Solitons Fractals 2005, 26, 1415–1421. [Google Scholar] [CrossRef]
  122. Wang, Q.; Chen, Y.; Zhang, H. An Extended Jacobi Elliptic Function Rational Expansion Method and its Application to (2+1)-Dimensional Dispersive Long Wave Equation. Phys. Lett. A 2005, 340, 411–426. [Google Scholar] [CrossRef]
  123. Chen, Y.; Wang, Q. Extended Jacobi Elliptic Function Rational Expansion Method and Abundant Families of Jacobi Elliptic Function Solutions to (1 + 1)-Dimensional Dispersive Long Wave Equation. Chaos Solitons Fractals 2005, 24, 745–757. [Google Scholar] [CrossRef]
  124. Lü, D. Jacobi Elliptic Function Solutions for Two Variant Boussinesq Equations. Chaos Solitons Fractals 2005, 24, 1373–1385. [Google Scholar] [CrossRef]
  125. Abdou, M.A.; Elhanbaly, A. Construction of Periodic and Solitary Wave Solutions by the Extended Jacobi Elliptic Function Expansion Method. Commun. Nonlinear Sci. Numer. Simul. 2007, 12, 1229–1241. [Google Scholar] [CrossRef]
  126. El-Sabbagh, M.F.; Ali, A.T. New Generalized Jacobi Elliptic Function Expansion Method. Commun. Nonlinear Sci. Numer. Simul. 2008, 13, 1758–1766. [Google Scholar] [CrossRef]
  127. Ali, A.T. New Generalized Jacobi Elliptic Function Rational Expansion Method. J. Comput. Appl. Math. 2011, 235, 4117–4127. [Google Scholar] [CrossRef] [Green Version]
  128. Zhou, Y.; Wang, M.; Wang, Y. Periodic Wave Solutions to Coupled KdV Equations With Variable Coefficients. Phys. Lett. A 2003, 308, 31–36. [Google Scholar] [CrossRef]
  129. Wang, M.; Zhou, Y. The Periodic Wave Solutions for the Klein–Gordon–Schrödinger Equations. Phys. Lett. A 2003, 318, 84–92. [Google Scholar] [CrossRef]
  130. Wang, M.; Li, X. Applications of F-Expansion to Periodic Wave Solutions for a New Hamiltonian Amplitude Equation. Chaos Solitons Fractals 2005, 24, 1257–1268. [Google Scholar] [CrossRef]
  131. Ren, Y.-J.; Zhang, H.-Q. A Generalized F-Expansion Method to Find Abundant Families of Jacobi Elliptic Function Solutions of the (2+1)-Dimensional Nizhnik-Novikov-Veselov Equation. Chaos Solitons Fractals 2006, 27, 959–979. [Google Scholar] [CrossRef]
  132. Jawad, A.J.M.; Petkovic, M.D.; Biswas, A. Modified Simple Equation Method for Nonlinear Evolution Equations. Appl. Math. Comput. 2010, 217, 869–877. [Google Scholar] [CrossRef]
  133. Liu, S.; Fu, Z.; Liu, S.-D.; Zhao, Q. A Simple Fast Method in Finding Particular Solutions of Some Nonlinear PDE. Appl. Math. Mech. 2001, 22, 326–331. [Google Scholar] [CrossRef]
  134. Xie, Y.; Tang, J. A Unified Trial Function Method in Finding the Explicit and Exact Solutions to Three NPDEs. Phys. Scr. 2006, 74, 197–200. [Google Scholar] [CrossRef]
  135. Feng, Z. The First-Integral Method to Study the Burgers–Korteweg–de Vries Equation. J. Phys. Math. Gen. 2002, 35, 343–349. [Google Scholar] [CrossRef]
Figure 1. The general case of simple equations method (SEsM) and its particular cases. The general case of the SEsM is for a system of N differential equations and the solution is constructed on the basis of solutions of M simple equations (we note that he parameter M may depend on the parameter N). A particular case of the general SEsM is the case when one has to solve one differential equation and the solution is constructed on the basis of solutions of M simple equations. The simplest case of SEsM is when one has to solve one differential equation and the solution is constructed by solutions of one simple equation. This particular case is known as the Modified Method of Simplest Equation.
Figure 1. The general case of simple equations method (SEsM) and its particular cases. The general case of the SEsM is for a system of N differential equations and the solution is constructed on the basis of solutions of M simple equations (we note that he parameter M may depend on the parameter N). A particular case of the general SEsM is the case when one has to solve one differential equation and the solution is constructed on the basis of solutions of M simple equations. The simplest case of SEsM is when one has to solve one differential equation and the solution is constructed by solutions of one simple equation. This particular case is known as the Modified Method of Simplest Equation.
Entropy 23 00010 g001
Figure 2. The seven steps in the general case of SEsM. For a more detailed description, see the text.
Figure 2. The seven steps in the general case of SEsM. For a more detailed description, see the text.
Entropy 23 00010 g002
Figure 3. The steps for the particular case of SEsM for solving one differential equation by means of m simple equations. On the right-hand side of the figure one sees the particular case of the methodology called the Modified Method of Simplest Equation.
Figure 3. The steps for the particular case of SEsM for solving one differential equation by means of m simple equations. On the right-hand side of the figure one sees the particular case of the methodology called the Modified Method of Simplest Equation.
Entropy 23 00010 g003
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Vitanov, N.K.; Dimitrova, Z.I.; Vitanov, K.N. Simple Equations Method (SEsM): Algorithm, Connection with Hirota Method, Inverse Scattering Transform Method, and Several Other Methods. Entropy 2021, 23, 10. https://doi.org/10.3390/e23010010

AMA Style

Vitanov NK, Dimitrova ZI, Vitanov KN. Simple Equations Method (SEsM): Algorithm, Connection with Hirota Method, Inverse Scattering Transform Method, and Several Other Methods. Entropy. 2021; 23(1):10. https://doi.org/10.3390/e23010010

Chicago/Turabian Style

Vitanov, Nikolay K., Zlatinka I. Dimitrova, and Kaloyan N. Vitanov. 2021. "Simple Equations Method (SEsM): Algorithm, Connection with Hirota Method, Inverse Scattering Transform Method, and Several Other Methods" Entropy 23, no. 1: 10. https://doi.org/10.3390/e23010010

APA Style

Vitanov, N. K., Dimitrova, Z. I., & Vitanov, K. N. (2021). Simple Equations Method (SEsM): Algorithm, Connection with Hirota Method, Inverse Scattering Transform Method, and Several Other Methods. Entropy, 23(1), 10. https://doi.org/10.3390/e23010010

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop