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

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.


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 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 solu-tions 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 threesolution 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.

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 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.
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)  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 i1 , ..., f iN , g i1 , . . . , g iM , . . . . 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 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  . . , f iN , g i1 , . . . , g iM 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 i1 , . . . , f iN , g i1 , . . . , g iM , . . . 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, 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, where q µ 1 , r µ 2 , . . . are parameters. However, other kinds of relationships, and more complicated ones, are also possible.  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.

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 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 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, which can be written also as follows [99] D m The bilinear operators have useful properties such as D m x (a.a) = 0 for m odd (14) 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.

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.

Example 1: The Three-Soliton Solution of the Korteweg-de Vries Equation
We consider the Korteweg-de Vries equation At Step 1 of SEsM we transform the nonlinearity in (16). This is made by the transformation (8). The result is and this can be written by means of the bilinear operators of Hirota as 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 Note that this is a linear equation. The equations for f 2 and f 3 are 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 where η 1 = λ 1 x + ω 1 t + σ 1 and λ 1 , ω 1 and σ 1 are parameters. We choose f 1 just as The substitution of (23) in (19) leads to the algebraic equation (the dispersion relation) 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 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 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 The substitution of (27) in (19) leads to a system of two algebraic equations, namely The substitution of (27) in (20) leads to a nonlinear equation for f 2 . The solution of this equation is where 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 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, Then we present f 1 as the sum of the solutions of the three simple equations The substitution of (32) in (19) leads to the algebraic relationships (dispersion relations) and the substitution of (32) in (20) leads to the following relationship for f 2 where The substitution of the obtained solutions for f 1 and f 2 in (21) leads to the following solution for f 3 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 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.
where q is a parameter. We consider the case q < 0 (and we shall set below q = −2 for convenience) with boundary conditions | ψ | 2 = ρ 2 0 at x → ±∞. At Step 1 of SEsM we transform the nonlinearity in (37) by the transformation 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) This can be written as Introducing the constant λ we can write two coupled equations of the basis of (39) λ 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 2 , . . . , g 1 , g 2 , . . . have to go to 0 at x → −∞. Then at x → −∞ g o g * 0 = ρ 2 0 and from (43) and (44) 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 For the terms of order of 2 For the terms of order of 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 where η = qx − ω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) 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 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 where The substitution of (57) in (48) and (49) (Steps 6 and 7 of SEsM) leads to a system of algebraic equations which solution is The substitution of (57), (58) in (50) and (51) leads to solutions containing single exponential function where 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 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.

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 where 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 The Ishimori equation becomes i ∂ω ∂t where At Step 1 of SEsM, we perform the transformation 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 At Step 2 of SEsM, we use two expansions: one expansion for for f and one expansion for g: 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 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 where M is an arbitrary complex constant, ξ 1 = l 1 x + m 1 y + n 1 t, n 1 = i(l 2 1 + σ 2 m 2 1 ); ξ 1 has real and imaginary part: . 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 where In a similar way we can obtain two-soliton solution, three-soliton solution, etc. of the Ishimori equation.

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: Researchers tried to solve this equation by means of some transformation. Almost 50 years ago Gardner, Greene, Kruskal, and Miura [48] tried the transformation This can be written as 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 where P = ∂ψ ∂t 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 = −κ 2 n and for the corresponding eigenfunctions ψ n we have | ψ n |→ 0 at | x |→ 0 and 0 < ∞ −∞ dxψ 2 n < ∞. 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λ dt = 0 and from (76) and (77) we obtain 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 = −4iµ 3 and the system of equations becomes ∂ψ ∂t 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 χ 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 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 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) At x → ±∞ one has µ = iκ n and ψ n ( The solution of the scattering problem at x → ∞ has the behavior (86) must be an asymptotic solution of (80) when µ = k. This condition leads to the relationships f (k, t) = exp(−8ik 3 t); g(k, t) = β and the reflection coefficient is 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 where K(x, y, t) is the solution of a linear integral equation known as the equation of Gelfand-Levitan-Marchenko where 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.

Connection between SEsM and the Inverse Scattering Transform Method
We shall discuss below the connection between SEsM and IST for the case of Kortewegde 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 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.

Example 1: The Burgers Equation
The Burgers equation is ∂u ∂t 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 Now step by step we can solve the Equation (93) and we can obtain u n . For n = 1 we have 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 where dλ(k) is an appropriate measure in the complex plane C and the factor 2ik 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 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 (97) and the equation for u 3 becomes We search for solution of (97) as and for solution for u 3 we assume The substitution of (99) in (97) and of (100) in (98) leads to following relationships for Φ 2 and Φ 3 Φ 2 (k 1 , k 2 ) = 2ik 1 ; Φ 3 (k 1 , k 2 , k 3 ) = 2ik 1 .
Next we assume the general form of the solution of the equations for u n to be: 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 2 1 + · · · + k 2 n )t. Because we have Φ 2 = Φ 3 = 2ik 1 we assume further that 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 By means of the notation we can write u as We use the formula for summation of the power series and apply it to (106) with q = ζ. The result is We remember that ζ and thus 1 − ζ where is a parameter, is the general solution of the heat equation. In addition (107) (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.

Example 2: The Korteweg-de Vries Equation
We consider the Korteweg-de Vries equation 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 and the equation for u 1 is ∂u 1 ∂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 The solutions for u 2 , u 3 , . . . are searched again in the form 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 The substitution of (114) in (113) leads to ; In addition we can obtain for the parameters of the solution u 4 ; The continuation of the calculations leads us to The solution of the Korteweg-de Vries equation becomes 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 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 Then where (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 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 pp T and then 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), (127) and the solution of KdV equation becomes We can write the solution as where After some calculation we obtain which is the GLM equation with 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.

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]).

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 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 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 where and k and λ are parameters. The solution of (136) is Now let us show that this version of the homogeneous balance method is connected to a particular case of SEsM. 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.

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.
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). 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.

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 Then traveling wave variable ξ = x − ωt is introduced and Equation (139) becomes The solution of (140) is searched as where a i are parameters and z(ξ) is the solution of the auxiliary equation Equations (141)  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 Then, traveling wave variable ξ = x − ωt is introduced and Equation (139) becomes The solution of (140)  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].

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 xx , u xt , u tt , . . . ) = 0, and searches for traveling wave solutions in the form 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, 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 of the functions f 1 , . . . , f N and each of these functions is a solution of a differential equation for the Jacobi elliptic functions We show below that the GJEFEM is connected to a particular case of SEsM. 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.

1.
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.

2.
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).

3.
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.

4.
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 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.

5.
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 6. 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 7.
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].

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 zz , u zzz , . . . ) = 0 (154) (154) is solved by means of the ansatz 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 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 Fthis 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.

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 xx , u xt , u tt , . . . ) = 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.

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 xx , u xt , u tt , . . . ) = 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(kx − ωt). Then one introduces X = U and Y = U z and writes (157) as system of equations The solution is obtained by the assumption that the derivative of the relationship Q(X, Y) = m ∑ i=0 a i (X)Y i can be represented as 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 zz = 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.

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.