On the Use of Composite Functions in the Simple Equations Method to Obtain Exact Solutions of Nonlinear Differential Equations

: We discuss the Simple Equations Method (SEsM) for obtaining exact solutions of a class of nonlinear differential equations containing polynomial nonlinearities. We present an amended version of the methodology, which is based on the use of composite functions. The number of steps of the SEsM was reduced from seven to four in the amended version of the methodology. For the case of nonlinear differential equations with polynomial nonlinearities, SEsM can reduce the solved equations to a system of nonlinear algebraic equations. Each nontrivial solution of this algebraic system leads to an exact solution of the solved nonlinear differential equations. We prove the theorems and present examples for the use of composite functions in the methodology of the SEsM for the following three kinds of composite functions: (i) a composite function of one function of one independent variable; (ii) a composite function of two functions of two independent variables; (iii) a composite function of three functions of two independent variables.


Introduction
We discuss in this article the mathematical problem for obtaining exact analytical solutions of nonlinear differential equations. The discussion is based on an amended version of a methodology called the SEsM. The emphasis in this amended version is on the use of the composite functions and their derivatives. The motivation of the study is as follows: 1.
Often, the effects connected to the nonlinearity are studied by means of time series analysis or by means of models based on differential or difference equations [14][15][16][17][18][19].
The corresponding model equations are nonlinear differential equations.
The following points from the history of the methodology for obtaining exact solutions of nonlinear differential equations are relevant for our study: 1.
Initially, the efforts were directed toward removing the nonlinearity of the solved equation by means of an appropriate transformation. An example is the Hopf-Cole transformation [20,21]. It transforms the nonlinear Burgers equation to the linear heat equation; 2.
Another transformation connects the Korteweg-de Vries equation to the famous equation of Schrödinger. Thus, the method of inverse scattering transform was born [22][23][24]; 3.
Almost at the same time, Hirota developed a method for obtaining exact solutions of nonlinear partial differential equations [25,26]. The Hirota method is connected also 1.
We note the work of Kudryashov. He formulated the Method of the Simplest Equation (MSE) [32]. The method is based on the determination of the singularity order n of the solved equation. Then, a particular solution of this equation is searched as a series containing powers of the solution of a simpler equation. This simpler equation is called the simplest equation. The methodology was extended [33] and applied to obtain traveling wave solutions of nonlinear partial differential equations (see, e.g., [34][35][36]). Kudryashov [37] used various transformations in order to transform the nonlinearity of a generalized evolution equation of the wave dynamics. Then, he obtained exact solutions of this equation. This research was continued in [31,38,39]. For recent results connected to the application of the method of the simplest equation, see [40][41][42][43][44][45][46]; 2.
We developed a methodology for obtaining the exact and approximate solutions of nonlinear partial differential equations. The methodology is called the Simple Equations Method (SEsM) [47][48][49][50][51]. Some elements of the methodology can be seen in our publications written a long time ago [52][53][54][55]. At the beginning [56,57], we used the ordinary differential equation of Bernoulli as the simplest equation [58]. This version of the methodology was called the Modified Method of the Simplest Equation (MMSE). It was used to obtain exact solutions of model nonlinear partial differential equations from ecology and population dynamics [59]; 3.
In these early publications, we used the concept of the balance equation. This helped us determine the kind of the simplest equations, as well as the form of the solution as a series of the solution of the simplest equation [60,61]. We note that the MMSE leads to results that are equivalent to the methodology of Kudryashov mentioned above. Our contributions to the methodology and its application till 2018 were connected to the MMSE [62][63][64][65][66][67]. We note especially the article [66]. It is connected to the part of the topics discussed below in the text; 4.
In the course of the years, the MMSE was extended to the SEsM [47]. The SEsM is connected to the possibility of the use of more than one simple equation. Thus, the solution of the solved nonlinear differential equation can be constructed on the basis of many simple equations. A version of the SEsM based on two simple equations was applied in [68]. The first description of the methodology was made in [48] and then in [47,[49][50][51]69]. For more applications of specific cases of the SEsM, see [70][71][72].
Our idea in this study is as follows. We used the SEsM to study the mathematical problem for obtaining exact solutions of certain classes of nonlinear differential equations. Below, we were not interested in imposing boundary conditions on the solved differential equations. An important part of the SEsM is the construction of a solution of the solved equation. This solution is presented as a composite function of the solutions of simpler differential equations.
The text is organized as follows. We briefly describe the previous version of the SEsM from [47] in Section 2. Our interest below was to search for exact traveling wave solutions of one nonlinear differential equation. In Section 3, we discuss the amended version of the SEsM. It has a simplified schema because of the use of the composite function and its derivatives. One new theorem is discussed. In addition, a consequence of another theorem (proved by us) is considered. The theorem and the consequence are connected to the application of the SEsM to a specific case. This is the case of a nonlinear differential equation with polynomial nonlinearities and simple equations containing polynomial nonlinearities as well. Illustrative examples for the application of the amended version of the SEsM are discussed. Finally, several concluding remarks are given in Section 4.

Materials and Methods
Below, we discuss the simple equations method. It is an algorithm for obtaining exact and approximate solutions of nonlinear differential equations. The method was designed for the search of the solutions of systems of n nonlinear differential equations. The solution was constructed by the solutions of m simpler differential equations ( Figure 1). The most applications of the algorithm up to now have been to obtain solutions of one nonlinear differential equation. This solution is searched as a function of the solution of one simple equation.  Figure 1. The general case of the simple equations method and two of its specific cases. The general case of the SEsM is to search for the solution of a system of N differential equations. The solution is constructed on the basis of the solutions of M simpler differential equations (the parameter M may depend on the parameter N). A specific case of the SEsM is the case when one has to solve one differential equation and the solution is constructed on the basis of the solutions of M simpler differential equations. The simplest case of the SEsM is as follows. One has to solve one differential equation. The solution is constructed by the solution of one simple equation. This specific case is known as the modified method of simplest equation.

USE OF M SIMPLE EQUATIONS
The version of the SEsM discussed in [47] has 7 steps. They are shown in Figure 2. We observed that this version of the methodology can be amended by the use of composite functions and their derivatives. The result was that the number of steps of the SEsM decreased from 7 to 4. This amended version is discussed below.  [47]. For details, see [47].

The Amended Version of the SEsM
The amended version of the method of simple equations consists of four steps (see Figure 3). We unified Steps 2-4 from the previous version into Step 2 of the amended version. In addition, Steps 5 and 6 from the previous version were unified into Step 3 of the amended version. The description of the amended version of the methodology in is presented in more detail as follows.  We considered a system of nonlinear partial differential equations: In Equation (1), B 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 several spatial coordinates). We performed four steps: (1) We applied transformations: (2) They can transform the nonlinearity of the solved differential equations to a more treatable kind of nonlinearity (e.g., to polynomial nonlinearity).
In the case of one solved equation, the transformation T(F, G, . . . ) can be: the in the case of sh-Gordon (Poisson-Boltzmann equation) (for applications of the last two transformations, see, e.g., [52][53][54] ; or another transformation. In numerous cases, one may skip this step (then, we have u i (x, . . . , t) = F i (x, . . . , t)). In many other cases, the transformation is needed to obtain a solution of the studied nonlinear PDE. The application of Equation (2) to Equation (1) leads to a nonlinear differential equations for the functions F i , G i , . . . . No general form is known for the transformations T i . The reason is that the nonlinearities in the solved equations can be of different kinds. The most studied cases of transformations are transformations that result in differential equations containing polynomial nonlinearities; (2) This step is based on the use of composite functions. It unifies Steps 2-4 from the version of the SEsM from [47]. In this step, the functions F i (x, . . . , t), G i (x, . . . , t), . . . are chosen as composite functions of the functions f i1 , . . . , f iN , g i1 , . . . , g iM , . . . , which are solutions of simpler differential equations (Step 2 in Figure 3). There are two possibilities: (a) The construction relationship for the composite function is not fixed. Then, we have to use the Faa di Bruno relationship for the derivatives of a composite function; The construction relationship for the composite function is fixed. For example, for the case of one solved equation and one function F, the construction relationship can be chosen to be: Then, one can directly calculate the corresponding derivatives from the solved differential equation; (3) In this step, we have to select the simple equations for the functions f i1 , . . . , f iN , g i1 , . . . , g iM , . . . . In addition, if we are in the hypothesis of Point (a) of Step 2, we have to fix the relationship between the composite functions F i (x, . . . , t), G i (x, . . . , t), . . . , and the functions f i1 , . . . , f iN , g i1 , . . . , g iM , . . . . We note that the fixation of the simple equations and the fixation of the relationships for the composite functions are connected. The fixations transform the left-hand sides of Equation (1). The result of this transformation can be functions that are the sum of terms. Each of these terms contains some function multiplied by a coefficient. This coefficient is a relationship containing some of the parameters of the solved equations and some of the parameters of the solutions of the simple equations used. The fixation mentioned above is performed by a balance procedure that ensures that the relationships for the coefficients contain more than one term. This balance procedure leads to one or more additional relationships among the parameters of the solved equation and parameters of the solutions of the simple equations used. These additional relationships are known as balance equations; (4) We may obtain a nontrivial solution of Equation (1) if all coefficients mentioned in Step 3 are set to zero. This condition usually leads to a system of nonlinear algebraic equations. The unknown variables in these equations are the coefficients of the solved nonlinear differential equation and the coefficients of the solutions of the simple equations. Any nontrivial solution of this algebraic system leads to a solution of the studied nonlinear partial differential equation.

Faa di Bruno Relationship for Derivatives of a Composite Function
The composite functions play an important role in the amended version of the SEsM. Below, we provide information about the derivatives of the composite functions. These derivatives were used in the process of the application of the SEsM.

The General Case: Composite Function of Many Functions of Many Independent Variables
Let us consider the function h(x 1 , . . . , x d ). It is a function of d independent variables x 1 , . . . , x d . We assumed that the function h is a composite function of m other functions We note that in this notation, D 0 x is the identity operator; 7.
|| z ||= max 1≤i≤d | z i | is the maximum value component of the multicomponent vari- Then, we define: Ordering of vector indexes: For two vector indexes µ = (µ 1 , . . . , µ d ) and ν = (ν 1 , . . . , ν d ), we have µ ≺ ν when one of the following holds: Below, we also use the notation: Then, the Faa di Bruno relationship for the composite derivative of a function containing functions of many variables is [73]: In Equation (4), n =| ν |. In addition, Finally, Equation (4) can be simplified by a change of the notation [73]. We introduce: In addition: Finally, 0 ≺ l n−s+1 · · · ≺ l n are such that Equation (4) can be written as:

Several Specific Cases of the General Relationship
We discuss below the specific case when the composite function h is a function of two independent variables x 1 and x 2 . In this case: Then, from Equation (4), the Faa di Bruno formula for the composite function containing functions of two variables is: The version of the formula that arises from Equation (5) is: Let us consider the case of the composite function containing three functions of two independent variables. In this case: Then, from Equation (4), the Faa di Bruno formula for a composite function containing three functions of two independent variables is: The version of Equation (8) arising from Equation (5) is: Next, we considered the case of the composite function containing two functions of two independent variables. In this case, the composite function is a function of the functions g (1) (x 1 , x 2 ) and g (2) (x 1 , x 2 ). The Faa di Bruno formula for the composite function containing two functions that are functions of two variables is: The version of Equation (10) arising from Equation (5) is: Finally, we considered the case of a composite function containing one function of one variable. For the case of one function of one variable h = f [g(x)], the Faa di Bruno formula is: In Equation (12): • h (n) = d n h dx n is the n-th derivative of the function h; • f (k) = d k f dg k is the k-th derivative of the function f ; • g (i) = d i g dx i is the i-th derivative of the function g; • p(n, k) = {λ 1 , λ 2 , . . . , λ n }: the set of numbers such that:

Several Results Relevant for Applications of the SEsM
Here, we present two theorems. They are for some features of the SEsM connected to the use of composite functions by this methodology. Theorem 1 is for the case of a differential equation containing polynomial nonlinearities where the unknown function h depends on two independent variables. The theorem states that for the case when the unknown function is a composite function, constructed by exponential functions, the solved equation can be reduced to a system of nonlinear algebraic equations. The theorem ensures that the SEsM is applicable in this case. The theorem is connected to the solutions of (1 + 1)dimensional equations constructed on the basis of exponential functions. Many equations that have multisoliton solutions are of this class. The applicability of the discussed version of the SEsM to such equations is shown on the basis of a classic example-the Korteweg-de Vries equation. Another application is connected to obtaining standing wave solutions of differential equations for the case when these solutions are constructed by trigonometric and hyperbolic functions.
Theorem 2 shows that the SEsM can be applied also when the composite function is constructed by a function of a single variable and this function of a single variable satisfies the differential equation containing polynomial nonlinearities. We discuss a consequence of this theorem for the case when the simple equation used contains the equations of Bernoulli and Riccati as specific cases. This case ensures that the SEsM can be applied for the solution of differential equations containing polynomial nonlinearities for the cases when the simple equations used are the equations of Bernoulli, Riccati, Jacobi elliptic functions, etc. In such a way, we show why a large number of methods for obtaining exact solutions of nonlinear differential equations can be treated as specific cases of the SEsM. Theorem 1. Let us consider a nonlinear partial differential equation that contains a polynomial P of the function h(x 1 , x 2 ) and its derivatives. The relationship for this equation is: Above, N can be any natural number. We searched for the solution of the above equation in the form: h(x 1 , where h is the polynomial of the functions g (1) (x 1 , x 2 ), . . . , g (m) (x 1 , x 2 ). Let each function g (i) (x 1 , x 2 ) satisfy the simple equation: where α i,j is a constant parameter. Then, the solved nonlinear partial differential equation is reduced to a polynomial of the functions g (1) Proof. Let us consider Equation (7). We substituted Equation (2) in (7). The result is: We looked at Equation (15). The term: is a polynomial of g (1) , . . . , g (m) . It is multiplied by another polynomial of g (1) , . . . , g (m) , namely by . Then, any of h ν 1 , . . . , h ν N is a polynomial of g (1) , . . . , g (m) . This means that P(h, h ν 1 , . . . , h ν N ) is a polynomial of g (1) , . . . , g (m) .
We note that the resulting polynomial P(h, h ν 1 , . . . , h ν N ) contains monomials multiplied by coefficients that are nonlinear algebraic relationships between the parameters of the solved nonlinear partial differential equation and parameters α i,j . We can set these coefficients to zero and obtain a system of nonlinear algebraic equations. Any nontrivial solution of this algebraic system (if such a solution exists) leads to a solution of the solved nonlinear partial differential Equation (13). Thus, Theorem 1 ensures that we can reduce to a nonlinear algebraic system any nonlinear differential equation that depends on two independent variables and has polynomial nonlinearities. In order to do this, we have to search for the solution as a composite function containing exponential functions. Theorem 1 is an indication of how to apply the SEsM in such cases. The class of equations concerned is a large one. Many important equations are in this class (e.g., many equations that have multisoliton solutions). For the convenience of the reader, we write the form of the derivatives of the composite function constructed by the exponential function for two cases. For the case of a composite function containing three exponential functions that are functions of two independent variables, we considered Equation (9) and substituted Equation (2) in (9). The result is: For the case of a composite function containing two exponential functions that are functions of two independent variables, we considered Equation (11). We substituted Equation (2) in (11). The result is: (k j,1 !k j,2 !)(l j,1 ! + l j,2 !) k j,1 +k j,2 α l j,1 k j,1 1,1 α l j,1 k j,2 2,1 × α l j,2 k j,1 1,2 α l j,2 k j,2 2,2 g (1) k j,1 g (2) k j,2 .
We can use simple equations that are more complicated than the equation for the exponential function (Equation (2)) in the SEsM. Let us consider the case of the composite function of one function, which depends on a single independent variable. Let us consider a nonlinear partial differential equation with nonlinearities that are polynomials of the unknown function h(x, t) and its derivatives. We searched for a solution of the kind: where µ and ν are parameters. The basis of our search is a solution g(ξ) of a certain simplest equation. Hence: h from Equation (16) is a composite function. For the n-th derivative of h, we have the Faa di Bruno formula (12). Let us assume that f is a polynomial of g. Then: We used the following simple equation: In Equation (18), k, l, m are integers. The nonlinearity in Equation (18) is of the polynomial kind. As mentioned above, this simple equation is more complicated than Equation (2). The solution of Equation (18)  This function V has interesting properties. Its specific cases are trigonometric, hyperbolic, elliptic functions of Jacobi, etc. Below, we used the function V a 0 ,a 1 ,...,a m (ξ; 1, 2, m), which is the solution of the simple equation: We have proven the following theorem [66].

Theorem 2. If g 2
(1) is given by Equation (19) and f is a polynomial of g given by Equation (17), then for h[ f (g)], the following relationship holds: where K n (q, m)(g) and Z n (q, m)(g) are polynomials of the function g(ξ).
The polynomials K n (q, m)(g) and Z n (q, m)(g) can be calculated as follows.
Then, starting from Equation (20), we obtain: Several of the polynomials K n and Z n from Equation (21) are calculated in Appendix A. Theorem 2 can explain why many methods for obtaining exact traveling wave solutions of nonlinear partial differential equations that contain polynomial nonlinearities are specific cases of the SEsM. Many of these methods use simple equations, which are the equation of Bernoulli or the equation of Riccati. These equations are specific cases of Equation (19). This can be proven as follows. The equations of Bernoulli and Riccati are specific cases of the simple equation: In Equation (22), n and c j are constant parameters. However, Equation (22) is a specific case of (19). We prove this as follows. The idea is that Equation (19) contains all cases of (22). From Equation (22): In Equation (23), a k are appropriate combinations of the coefficients c i . Equation (23) is of the kind (19). Then, Equation (19) contains all possible relationships of the kind (22). However, Equation (19) contains more than this. For example, Equation (19) contains the case: Equation (24) cannot be reduced to a relationship of the kind (22). Then, Equation (22) is a specific case of (19).
Thus, the methods based on the use of the equations of Bernoulli and Riccati for the construction of solutions of nonlinear differential equations with polynomial nonlinearities are specific cases of the SEsM where the composite function can be represented according to Theorem 2. The discussed consequence of Theorem 2 shows, for example, that the homogeneous balance method [74] and tanh method [75] are specific cases of the SEsM.
The direct application of Theorem 2 is convenient when, for example, the simple equation is an equation for a Jacobi elliptic function. Theorem 2 shows that, for example, the Jacobi elliptic function expansion method [76] and the F-expansion method [77] are specific cases of the SEsM.
For the case when the simple equation has the specific form of Equation (22), we have a simpler situation. Instead of the two kinds of polynomials Z n and K n , we have a single kind of polynomial L n . In other words, for the case when the simple equation is of the kind (Equation (22)), h (n) is a polynomial of g: h (n) = L n (g). These polynomials can be calculated as follows. We start from: Then, we use the recurrence relationship: Several of the polynomials L i are calculated in Appendix B. In addition, we present in Appendix C several derivatives of a composite function that is function of two functions of two independent variables. In Appendix D, we present several derivatives of a composite function that is function of three functions of two independent variables. The information from these appendices are used in the examples below.

Illustrative Examples
The goal of the discussion below is to show how the methodology of the SEsM works. For an example, multisoliton solutions can be obtained. We show this for the Korteweg-de Vries equation. In addition, we show that the methodology can lead to specific solutions of nonlinear partial differential equations that are not integrable.
Below, we used the composite functions in the SEsM as follows. We used Equation (3) to fix the relationship for the composite function. On the basis of Equation (3), h can be written as: N can be any natural number. Note that the value of N can also be equal to ∞. The relationships for T i are: . . .
We note that the functions g (i n ) are the solutions of the corresponding simple equations. For the cases considered below, we used the following specific forms of Equations (26) and (27). For the case of the composite function of a function of one variable h = f [g (1) In Equation (28), N is a parameter, which is fixed by means of balance equation. For the case of the composite function of a function of two variables h = f [g (1) (x, t), g (2) (x, t)], we used the simple relationship: For the case of the composite function of a function of three variables h = f [g (1) (x, t), g (2) (x, t), g (3) (x, t)], we used the simple relationship: Equations (28) and (29) were obtained by selecting only several terms from the large number of possible terms in Equation (26).
Next, we considered illustrative examples for the case of the composite function that is a function of one function of one variable. First of all, we considered the Korteweg-de Vries equation: Following the steps of the SEsM, we applied the transformation: This is Step 1 of the SEsM. The Korteweg-de Vries Equation (31) was transformed to: Next, we considered h as a composite function of one function of one variable ξ = λx + ωt + σ. λ, ω, and σ are parameters (Step 2 of the SEsM). Equation (33) was transformed to: In Step 3 of the SEsM, we used Equation (28) as the relationship for the composite function h. We set N = 1 in (28). In addition, the simple function g (1) (ξ) was assumed to satisfy the simple Equation (22) for the case n = 1 and c 0 = 0. Then, the derivatives h (n) of the composite function h are given by the polynomials L i discussed above in the text (the relationships for L i are given in Appendix D).
In such a way, the Korteweg-de Vries Equation (34) was reduced to a polynomial that has to be equal to zero. This was achieved by setting the coefficients of the terms of the polynomial to zero. The result is a system of algebraic equations (Step 4 of the SEsM). This system of algebraic equations contains one equation: in addition to β 0 = 1. Equation (35) can be solved for ω, for example. The solution is nontrivial. This, nontrivial solution leads to the one-soliton solution of the Korteweg-de Vries equation. Another example is for the equation of the class: In Equation (36), p 1 , p 2 , p 3 , p 4 , p 5 , p 6 , p 7 , p 8 are parameters. We skipped Step 1 of the SEsM (no transformation of the nonlinearity). In Step 2 of the SEsM, we used Equation (28) as the relationship for the composite function h. In Step 3 of the SEsM, we fixed the value of N. The form of the composite function h[g (1) (ξ)] (ξ = α 1 x + α 2 t) becomes: At the same time, the form of the simple equation for g (1) was fixed by means of a balance equation: By means of Equations (37) and (38), Equation (36) was reduced to a system of seven nonlinear algebraic equations (Step 4 of the SEsM). One nontrivial solution of this system is: .
The solution of the simple Equation (38) is: Then, a particular exact solution of the equation: is: where: Next, we present examples for the case of the composite function that is a function of two functions of two variables. First, we considered again the Korteweg-de Vries equation. In this case, we applied the SEsM as follows. We used the transformation (32) and obtained Equation (33) from (31) (Step 1 of the SEsM). Then, we considered h as a composite functions of two functions of two variables. We used the relationship (29) (Step 2 of the SEsM). In addition, the simple functions g (1) (η 1 ) and g (2) (η 2 ) satisfy simple equations of the kind (22) for the case n = 1 and c 0 = 0 (Step 3 of the SEsM) η i = λ i x + ω i t + σ i , i = 1, 2.
The further application of the methodology (Step 4 of the SEsM) led to γ 1 = γ 2 = 0 and reduced the Korteweg-de Vries equation to the system of two algebraic equations: in addition to α = 1; β 1 = β 2 = β, γ 3 = β 2 µ, where µ is a parameter. Equation (39) can be solved for ω i , i = 1, 2, for example. The solutions are nontrivial. These nontrivial solutions lead to the two-soliton solution of the Korteweg-de Vries equation. Another example is connected to the equation: We applied the SEsM and skipped Step 1 (no transformation of the nonlinearity). The composite function h is of the kind (29) where α = 0, β 1 = β 2 = 0, γ 1 = γ 2 = 0 (Step 2 of the SEsM). In addition, g (1) does not depend on t and g (2) does not depend on x. Let γ 3 = A. The composite function becomes: h(x, t) = Ag (1) (αx)g (2) (δγt), δ = ±1. (41) We note that the composite function (41) allows for complicated simple equations for g (1) and g (2) . These equations can be of the kind of equations for the elliptic functions of Jacobi (Step 3 of the SEsM): In Step 4 of the SEsM, Equation (40) is reduced to a system of algebraic equations: Equation (43) has various nontrivial solutions. For example, one of these solutions is when α 2 − γ 2 < 1. We can consider A as a free parameter. Then, In Equation (44), cn(αx; k 1 ) and cn(γt; k 2 ) are the corresponding Jacobi elliptic functions of modulus 0 ≤ k 1 ≤ 1 and 0 ≤ k 2 ≤ 1, respectively. Equation (43) has an interesting specific case when k 1 = 1 and k 1 = 0. In this case, cn(αx; k 1 ) = sech(αx) and cn(δγt) = cos(δγt). Then, Equation (45) can also be obtained straightforwardly on the basis of the composite function (41) if one takes for g (1) and g (2) the corresponding simple equations for the trigonometric and hyperbolic functions, respectively.
Finally, we discuss an example use case of the composite function that is a function of three functions of two variables. This led us to the three-soliton solution of the Korteweg-de Vries equation. In this case, we applied the SEsM as follows. We used the transformation (32) and obtained Equation (33) from (31). Then, we considered h as the composite function of three functions of two variables. We used the relationship (30). In addition, the simple functions g (1) (η 1 ), g (2) (η 2 ), and g (3) (η 3 ) satisfy simple equations of the kind of Equation (22) for the case n = 1 and c 0 = 0. η i = λ i x + ω i t + σ i , i = 1, 2, 3. The application of the SEsM led to γ 1 = γ 2 = 0 and reduced the Korteweg-de Vries equation to the system of three algebraic equations: in addition to α = 1; β 1 = β 2 = β 3 = β, γ i,i = 0, δ i,j,k = 0 except for δ 1,2,3 . Moreover, Finally, δ 1,2,3 = βa 1,2 a 1,3 a 2,3 . Equation (46) can be solved for ω i , i = 1, 2, 3, for example. The solutions are nontrivial. These nontrivial solutions led to the three-soliton solution of the Korteweg-de Vries equation.

Concluding Remarks
This article was devoted to a discussion of the methodology called the Simple Equations Method (SEsM) and several results connected to this methodology. The new results were as follows:

1.
We presented an amended version of the SEsM in comparison to the version from [47]. We discussed a theorem that states that under certain conditions, a nonlinear differential equation with polynomials nonlinearities can be reduced to a polynomial containing monomials consisting of exponential functions. This theorem justified the application of the SEsM as the setting of the coefficients of the obtained polynomial to zero led to a system of nonlinear algebraic equation, which led exactly to Step 4 of the SEsM. We note that in such a way, the SEsM can lead to multisoliton solutions of a large class of equations. An illustrative example for the Korteweg-de Vries equation was given; 3.
A consequence of a theorem proven in [66] was used in order to show that the simple equation of the SEsM can contain polynomial nonlinearities of large power for the case when the composite function used in the SEsM is a function of one independent variable. This consequence showed that many methods that search for exact traveling wave solutions of nonlinear differential equations on the basis, for example, of the equations of Riccati and Bernoulli are specific cases of the SEsM; 4.
We presented many illustrative examples for the application of the amended version of the SEsM.
Finally, we want to stress that the SEsM is an interesting methodology for obtaining exact solutions of nonlinear differential equations. We showed already that the SEsM is connected to the inverse scattering transform method and to the method of Hirota [47]. The SEsM can lead to specific solutions also of nonintegrable differential equations. The use of composite functions in the methodology of the SEsM opens possibilities for obtaining additional results on the methodology, as well as specific solutions of many nonlinear differential equations. In addition, the kinds of nonlinearities different from polynomial ones can be treated. This will be one of the goals of our future research.

Conflicts of Interest:
The authors declare no conflict of interest.