Next Article in Journal
A Preprocessing Perspective for Quantum Machine Learning Classification Advantage in Finance Using NISQ Algorithms
Next Article in Special Issue
Influence of Information Blocking on the Spread of Virus in Multilayer Networks
Previous Article in Journal
Tsallis Entropy for Loss Models and Survival Models Involving Truncated and Censored Random Variables
Previous Article in Special Issue
Flows of Substances in Networks and Network Channels: Selected Results and Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Simple Equations Method (SEsM): An Effective Algorithm for Obtaining Exact Solutions of Nonlinear Differential Equations

by
Nikolay K. Vitanov
1,2
1
Institute of Mechanics, Bulgarian Academy of Sciences, Acad. G. Bonchev Str., Bl. 4, 1113 Sofia, Bulgaria
2
Climate, Atmosphere and Water Research Institute, Bulgarian Academy of Sciences, Blvd. Tzarigradsko Chaussee 66, 1784 Sofia, Bulgaria
Entropy 2022, 24(11), 1653; https://doi.org/10.3390/e24111653
Submission received: 6 October 2022 / Revised: 2 November 2022 / Accepted: 11 November 2022 / Published: 14 November 2022

Abstract

:
Exact solutions of nonlinear differential equations are of great importance to the theory and practice of complex systems. The main point of this review article is to discuss a specific methodology for obtaining such exact solutions. The methodology is called the SEsM, or the Simple Equations Method. The article begins with a short overview of the literature connected to the methodology for obtaining exact solutions of nonlinear differential equations. This overview includes research on nonlinear waves, research on the methodology of the Inverse Scattering Transform method, and the method of Hirota, as well as some of the nonlinear equations studied by these methods. The overview continues with articles devoted to the phenomena described by the exact solutions of the nonlinear differential equations and articles about mathematical results connected to the methodology for obtaining such exact solutions. Several articles devoted to the numerical study of nonlinear waves are mentioned. Then, the approach to the SEsM is described starting from the Hopf–Cole transformation, the research of Kudryashov on the Method of the Simplest Equation, the approach to the Modified Method of the Simplest Equation, and the development of this methodology towards the SEsM. The description of the algorithm of the SEsM begins with the transformations that convert the nonlinearity of the solved complicated equation into a treatable kind of nonlinearity. Next, we discuss the use of composite functions in the steps of the algorithms. Special attention is given to the role of the simple equation in the SEsM. The connection of the methodology with other methods for obtaining exact multisoliton solutions of nonlinear differential equations is discussed. These methods are the Inverse Scattering Transform method and the Hirota method. Numerous examples of the application of the SEsM for obtaining exact solutions of nonlinear differential equations are demonstrated. One of the examples is connected to the exact solution of an equation that occurs in the SIR model of epidemic spreading. The solution of this equation can be used for modeling epidemic waves, for example, COVID-19 epidemic waves. Other examples of the application of the SEsM methodology are connected to the use of the differential equation of Bernoulli and Riccati as simple equations for obtaining exact solutions of more complicated nonlinear differential equations. The SEsM leads to a definition of a specific special function through a simple equation containing polynomial nonlinearities. The special function contains specific cases of numerous well-known functions such as the trigonometric and hyperbolic functions and the elliptic functions of Jacobi, Weierstrass, etc. Among the examples are the solutions of the differential equations of Fisher, equation of Burgers–Huxley, generalized equation of Camassa–Holm, generalized equation of Swift–Hohenberg, generalized Rayleigh equation, etc. Finally, we discuss the connection between the SEsM and the other methods for obtaining exact solutions of nonintegrable nonlinear differential equations. We present a conjecture about the relationship of the SEsM with these methods.

1. Short Introduction and Overview of the Literature

Nature and human societies are full of complex systems. Atomic chains and lattices, biological systems, the dynamics of research groups, networks, social sciences, and economics are some of the many examples [1,2,3,4,5,6,7,8]. These complex systems require sophisticated research instruments and powerful computers for their study. The research on complex systems is complicated as many of these systems are nonlinear. Examples of this nonlinearity can be seen in fluid mechanics, solid-state physics, etc. [9,10,11,12,13]. Time-series analysis and models based on differential or difference equations are often used in the studies of nonlinear complex systems [14,15,16,17,18,19].
This review is focused on one of the aspects of the study of complex systems: the use of nonlinear differential equations for modeling such systems. The research on nonlinear differential equations has a long history [20,21,22,23,24,25,26]. In particular, nonlinear waves have been attracting the attention of researchers for many decades [27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43]. Several examples of the occurrence of such waves are:
  • elastic media, specifically elastic crystals [44,45];
  • fluids [46,47,48,49,50,51,52], especially the interaction of nonlinear waves and offshore structures [53], and plasmas [54];
  • solids [55], for example, nonlinear acoustic waves in solids [56] or nonlinear spin waves in magnetic films [57];
  • optics [58], especially nonlinear waves in waveguides [59];
  • chemical systems [60].
Much attention has been devoted to the research of nonlinear waves in dispersive systems [61,62,63,64,65,66]. These waves often belong to a special class of nonlinear waves: the solitons [67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,82,83,84,85,86,87,88,89,90,91,92,93,94,95,96]. Solitons are predicted and observed in many systems such as conducting polymers [97], plasmas [98,99], and helium films [100], and can be used, for example, in long-distance communications [101].
Below, we describe a methodology for obtaining exact solutions of nonlinear differential equations. The methodology is called the SEsM (Simple Equations Method). The methodology is a recent development in a long line of studies of methods for obtaining exact solutions of nonlinear differential equations. The classic approach is the Inverse Scattering Transform method. The basis for the Inverse Scattering Transform method was described in the paper by Gardner, Greene, Kruskal, and Miura [102]. An excellent introduction to the methodology is [103,104]. Additional results on the methodology can be found in [105]. Lax [106] proposed a general principle for associating the evolution of nonlinear equations with linear operators so that the eigenvalues of the linear operator are integrals of the nonlinear equation. We also note the contribution of Zakharov and Shabat [107,108] and the results of Fokas and Ablowitz [109,110]; Manakov [111,112]; Dod and Bullough [113]; Ablowitz, Kruskal, and Segur [114]; Kaup and Newell [115,116]; Beals [117]; Kaup [118,119]; and other authors [120,121]. This short list of contributors is in no way complete.
The Inverse Scattering Transform method stimulated much research on nonlinear differential equations. We have already mentioned above the nonlinear waves in dispersive systems and solitons. Additional results have been obtained on the integrability of dynamical systems [122,123] and the integrability of evolution equations [124,125,126,127,128,129,130,131], the construction of higher-dimensional integrable systems [132], the spectral transformations method for solving nonlinear evolution equations [133], etc. [134,135]. Other methods for obtaining exact solutions of nonlinear partial differential equations have also appeared. A prominent example of such a method is the method of Hirota [136,137,138,139,140,141].
Numerous equations have been studied using the above-mentioned methods. One of the most popular is the Korteweg-de Vries equation [142,143,144,145,146,147,148,149,150,151,152,153,154,155,156,157,158,159,160,161,162,163,164,165,166,167,168,169,170]. A short list of some other studied equations is as follows:
The obtained exact solutions have been used to describe various phenomena. Several examples are:
  • the description of optical pulses [258,259,260,261,262,263,264,265], the propagation of optical solitons in fibers [266,267,268,269,270,271,272,273,274,275], and discrete solitons in optical systems [276,277];
  • solitons in condensed matter physics [278]; solitons in lattices [279,280,281,282,283,284,285,286,287,288,289,290,291,292,293,294,295,296,297,298,299]; solitons in Heisenberg chains [300,301,302,303,304]; spin-wave solitons in magnetic films [305]; solitons in liquid crystals [306]; solitons in elastic solids [307]; solitons in rods, plates, and shells [308]; and collisions of solitons in spin systems [309];
  • solitons in plasmas [310], for example, ion-acoustic solitons [311,312,313,314,315] and Langmuir solitons [316];
  • solitons in fluid mechanics, for example, water waves [317,318,319,320,321,322,323,324]. This is a large area of research as the model equations used in fluid mechanics can often be solved analytically and the corresponding solutions describe solitons and multisolitons. Several examples are the observation of solitons in a sea [325], the description of internal waves in the atmosphere and atmospheric solitons [326,327,328], solitons in film flows [329], solitary waves and liquid drops [330], nonlinear soliton wave–beach interaction [331], long internal waves in deep fluid [332], solitons in rotating baroclinic fluid [333], ocean waves [334], solitary waves in two-fluid media [335], the propagation of bores [336,337], Peregrine soliton [338], solitary Rossby waves [339], and nonlinear focusing of waves in a fluid of finite depth [340];
  • solitons in molecular structures [341,342,343] such as solitary waves in large blood vessels [344], dissipative solitons [345], and solitary excitations in muscle proteins [346];
  • solitons connected to the Thirring model [347,348], the relationship between quantum mechanics and solitons on an elastic road [349], and the interpretation of solitons as particles and oscillators [350].
In addition, there is a large amount of literature on the various properties of solitons and processes with the participation of solitons. Several examples are breaking solitons [351], solitons and collapses [352]; solitons in non-uniform and disordered media [353,354]; bifurcations of solitons [355]; dynamics of solitons under perturbations [356,357,358]; vortex solitons [359], positons [360]; coupled solitons [361]; multidimensional solitons [362]; collisions of solitary waves [363]; the connection between strings, vortices, and solitons [364]; the interaction between short and long waves [365]; and wave–wave scattering [366].
The methodology of the Inverse Scattering Transform method stimulated the research and obtainment of many interesting mathematical results. Several examples are:
  • Results on integrability: integrable Hamiltonian systems [367,368,369,370], asymptotic integrability [371,372], integrable nonlocal nonlinear equations [373], Lax pair for the one-dimensional Hubbard model [374] and higher-order models of water waves [375], the integrability of higher-dimensional discrete systems [376], geometrical aspects of solitons [377,378], conservation laws for nonlinear evolution equations [379], perturbed integrable equations [380], and the relationships between different methods for obtaining exact solutions of integrable systems [381];
  • Stability of solitary waves [382,383,384];
  • Pseudopotentials connected to nonlinear differential equations [385,386];
  • Algebraic results and results connected to symmetries: the relationship between integrable systems and groups of Lee [387], the results from algebraic theory applicable to the theory of nonlinear equations [388], symmetries and the integrability of nonlinear differential equations [389,390,391], similarity reductions of the Boussinesq equation [392], symmetries, Lax pairs and solutions for AKNS systems [393,394,395], and Katz–Moodey–Lee algebras and solitons [396];
  • Painlevé analysis of nonlinear differential equations [397,398,399,400,401], the connection between nonlinear evolution equations and ordinary differential equations of P-type [402,403], and perturbation theory for solitons [404];
  • The results from rational solutions and other solutions of nonlinear evolution equations [405,406] and connections between the AKNS and KP hierarchies [407];
  • The results from transformations of nonlinear differential equations: the direct linearization of nonlinear difference equations [408], gauge transformations applicable to solitons [409], Bäcklund transformations [410,411], the inverse spectral transform [412], generalized Fourier transforms for soliton equations [413], the singular manifold method for recovering Lax pairs and Darboux transformations [414], and the study of Bäcklund–Darboux transformations [415];
  • The results from solutions of nonlinear differential equations: Wronskian representation of N-soliton solutions [416], asymptotic solutions of nonlinear evolution equations [417], and the dynamics of classical solitons in nonintegrable systems [418].
Numerical analysis of nonlinear waves is of special interest [419]. Several results in this research area are connected to methods for the symplectic integration of Hamiltonian systems [420], the finite difference and finite element methods for the Korteweg-de Vries equation [421,422,423,424], the Galerkin method [425,426,427,428,429], and the collocation and radial basis function methods for the Korteweg-de Vries equation [430]. The numerical methods have been successfully applied to, for example, the study of the evolution of solitary waves over a shelf [431], Lax pairs symbolic computations [432], the nonisospectral scattering problems [433], the numerical solution of the nonlinear Schrödinger equation [434], the numerical study of internal wave solitons [435], and the study of the soliton decay of nonlinear Alvfen waves [436].
This concludes the brief survey of the literature, noting that it should be considered an introduction to the vast area of research on nonlinear differential equations and the methods for obtaining exact analytical solutions of such equations. Next, we present a brief history of the research that led to the appearance of the methodology, which is the focus of this review article.
As we have seen, research on the methodology for obtaining exact solutions of nonlinear differential equations is important. In the early years of these studies, there was an attempt to remove the nonlinearity of the studied equation by an appropriate transformation. An example is the Hopf–Cole transformation [437,438]. This transformation converts the nonlinear Burgers equation to the linear heat equation. Another transformation connected the nonlinear Korteweg-de Vries equation to the linear Schrödinger equation. Efforts to obtain such transformations led to the method of the Inverse Scattering Transform [439], as well as the method of Hirota [136,141]. The truncated Painlevé expansions also supplied suitable transformations of the nonlinearities [440,441,442,443]. The possibilities connected to the Painlevé expansions were also explored by Kydryashov [444], who used a truncated Painlevé expansion. In these expansions, the truncation occurs after the constant whose term is kept in the expansion. Kudryashov [445] formulated the Method of the Simplest Equation (MSE). The method is based on the determination of the singularity order n of the solved NPDE. Then, a particular solution of this equation is searched as a power series of solutions of a simpler equation called the simplest equation. For further results on this methodology and its applications, see [446,447,448,449,450,451,452,453,454,455,456,457,458,459].
About 13 years ago we started the development of a methodology for obtaining exact and approximate solutions of nonlinear differential equations. This methodology is now called the SEsM (Simple Equations Method) [460,461,462,463,464,465,466,467,468,469,470,471]. Elements of this methodology were used in our articles written decades ago [472,473,474,475]. In 2009 and 2010 [476,477], we used the ODE of Bernoulli as the simplest equation [478] in the first version of the method, the Modified Method of the Simplest Equation (MMSE), and applied the methodology to population dynamics and ecology [479]. The MMSE [480,481] uses the concept of balance equations for the fixation of the simplest equation. After this, the searched solution is constructed as a truncated power series of the solution of the simplest equation. This methodology led to equivalent results with respect to the Method of the Simplest Equation mentioned above.
The MMSE was actively used by us until 2018 [482,483,484,485,486,487]. A special role in this period was played by article [486]. There, the methodology of the MMSE was extended by the use of a class of equations d k g d ξ k l = j = 0 m d j g j as the simplest equations. Here, k = 1 , , l = 1 , , and m and d j are parameters. The solution of this equation contains specific cases of many well-known functions such as hyperbolic functions, trigonometric functions, the elliptic functions of Weierstrass and Jacobi, etc. The capacity of the methodology has been extended over the years. The current version of the methodology (SEsM) can use more than one simple equation. The SEsM based on two simple equations was used in [488]. The first discussion of the SEsM was in [460]. Further discussions on the SEsM can be seen in [461,462,463,464,465]. Applications of specific cases of the SEsM can be seen in [489,490,491,492,493,494,495].
Below, we discuss the SEsM in detail. The discussion covers the use of transformations to convert the nonlinearity of the solved equation to a more treatable kind of nonlinearity, the use of composite functions in the SEsM and the application of the Faa di Bruno formula for derivatives of these functions, the role of simple equations and a special class of simple equations that define an interesting special function, as well as a series of polynomials, which appear frequently in the application of the SEsM.
The text is organized as follows. In Section 2, the methodology of the SEsM is described. In Section 3, we discuss the use of the transformations of the nonlinearity of the solved equation using the SEsM. These transformations convert various kinds of nonlinearities to a polynomial kind of nonlinearity which is more easily treatable by the SEsM. In Section 4, the use of the composite functions of the SEsM is discussed with examples of such use. In Section 5, the role of simple equations in the SEsM is described. In Section 6, the connection of the SEsM to the inverse scattering transform method and the method of Hirota is discussed. Section 7 is devoted to the specific application of the SEsM for obtaining the exact solution of a differential equation connected to the SIR model of epidemic spreading, with a possible application for analysis of epidemic waves of a class of epidemics, including COVID-19. In Section 8, additional applications of the SEsM are presented. In Section 9, we discuss the conjecture about the connection of the SEsM with other methods for obtaining specific solutions on nonlinear nonintegrable differential equations. Several concluding remarks are summarized in Section 10.

2. The Simple Equations Method (SEsM)

The SEsM has the goal of obtaining exact and approximate solutions of nonlinear differential equations. The algorithm of the SEsM is designed for solving systems of n nonlinear differential equations. The solution is built by known exact solutions of m simpler differential equations [470]. Here, we describe a specific case of the SEsM for the solution of a single nonlinear differential equation based on the solution of m simpler equations. The MMSE is the specific case. There, m = 1 .
The discussed version of the SEsM has four steps. The solved differential equation is
A [ u 1 ( x , , t ) , , u n ( x , , t ) ] = 0 .
A [ u 1 ( x , , t ) , , u n ( x , , t ) , ] depends on the functions u 1 ( x , , t ) , , u n ( x , , t ) and some of their derivatives. The functions u i can depend on several spatial coordinates. The idea of the SEsM is to transform (1) to
r i ( ) B i = 0 .
Here, B i are some functions. r i are the relationships among the parameters of the solved equation and the parameters of the solution. We set
r i = 0 .
Thus, a system of nonlinear algebraic equations is obtained. Each nontrivial solution of (3) leads to a solution of the solved differential equation.
The steps for the transition from (1) to (3) are as follows. First, the transformation
u ( x , . . . , t ) = T [ F ( x , , t ) , G ( x , , t ) , ] ,
is applied. T ( F , G , ) is a composite function of the other functions F , G , . These functions can be the functions of several spatial variables and time. If possible, the transformation can remove the nonlinearity of the solved equation. An example is the Hopf–Cole transformation, which reduces the nonlinear Burgers equation to the linear heat equation [437,438]. The removal of the nonlinearity is rarely achieved. In most cases, the transformation T transforms the nonlinearity of the solved equation to an easily treatable kind of nonlinearity such as polynomial nonlinearity. Examples of the transformation T are u ( x , t ) = 4 tan 1 [ F ( x , t ) ] for the sine-Gordon equation and u ( x , t ) = 4 tanh 1 [ F ( x , t ) ] for the Poisson–Boltzmann equation [472,473,474,475]. The Painlevé expansion can be another appropriate transformation. u ( x , . . . , t ) = F ( x , . . . , t ) (no transformation) is also a possibility for certain classes of nonlinear differential equations. The application of (4) to (1) leads to a treatable nonlinear differential equation for F i , G i , . The transformation T may remain unfixed in the first step of the SEsM. Thus, the function T will be unknown and we have to determine it in some of the other three steps of the methodology.
Step 2 of the SEsM deals with the selection of the functions F ( x , . . . , t ) , G ( x , , t ) , …. They are the composite functions of the known solutions of simpler differential equations. We can let the form of the composite functions be undetermined. This, together with the form of the transformation T, can be determined in Step 3 of the SEsM. Often, F, G, …, are fixed in Step 2 of the SEsM. A form that is often used is
F = α + i 1 = 1 N β i 1 f i 1 + i 1 = 1 N i 2 = 1 N γ i 1 , i 2 f i 1 f i 2 + i 1 = 1 N i N = 1 N σ i 1 , , i N f i 1 f i N .
α , β i 1 , γ i 1 , i 2 , σ i 1 , , i N are parameters. The relationship used by Hirota [141] is a specific case of (5). Another specific case of (5) is the power series F = i = 0 N μ n f n ( μ is a parameter). This is frequently used in the MMSE.
Step 3 of the SEsM is often the most important in the application of the methodology. Here, one can determine the forms of the more simple equations and the composite functions T, F, G, …. All these are determined from the requirements for the satisfaction of the balance equations. The balance equations occur due to the requirement that any of r i from (3) must have at least two terms in order to ensure that a nontrivial solution of the solved equation is obtained. This determines the shapes of all composite functions and the forms of all used simpler equations.
The satisfaction of the balance equations in Step 3 leads to a system of nonlinear algebraic equations (3). This system contains relationships between the parameters of the solved equation and the parameters of the solution. In Step 4 of the SEsM, one solves the system (3). Each nontrivial solution of (3) leads to a nontrivial solution of the solved Equation (1).

3. The SEsM and the Transformation of the Nonlinearity of the Solved Equation

Step 1 of the SEsM is connected to the transformation of the nonlinearity of the solved equation (if necessary). If the nonlinearity of the solved equation is a polynomial one, there is no need for a transformation. Steps 2–4 of the SEsM can be applied directly in this case. In the case of the solved equation, which contains nonpolynomial nonlinearity, in Step 1 of the SEsM, one can attempt the application of a transformation that can transform the nonpolynomial nonlinearity to a polynomial nonlinearity. The following proposition can be performed, which was proved in [469].
Proposition 1.
Let us consider a differential equation for the function u ( x , , t ) , which contains two kinds of terms:
1. 
Terms containing only derivatives of u;
2. 
Terms containing one or several nonpolynomial nonlinearities of the function u and these nonpolynomial nonlinearities are of the same kind.
Let u = T ( F ) be a transformation with the following properties:
1. 
Property 1: The transformation T transforms any of the nonpolynomial nonlinearities to a function that contains only polynomials of F.
2. 
Property 2: The transformation T transforms the derivatives of u to terms containing only the polynomials of the derivatives of F or the polynomials of the derivatives of F multiplied or divided by the polynomials of F.
Then, the transformation T transforms the studied differential equation into a differential equation containing only the polynomial nonlinearity of F.
The above proposition allows us to transform numerous nonpolynomial nonlinearities into polynomial nonlinearities. A list of some nonpolynomial nonlinearities that can be transformed into polynomial nonlinearities can be found in [469]
Case 1:
N ( u ) = exp ( u ) ; N ( u ) = [ exp ( u ) ] m . The transformation is u = ln ( F ) .
Case 2:
N ( u ) = sin ( u ) ; N ( u ) = [ sin ( u ) ] m . A possible transformation is u = 4 tan 1 ( F ) .
Case 3:
N ( u ) = cos ( u ) ; N ( u ) = [ cos ( u ) ] m . The transformation is u = 4 tan 1 ( F ) .
Case 4:
N ( u ) = tan ( u ) ; N ( u ) = [ tan ( u ) ] m . A possible transformation is u = tan 1 ( F ) .
Case 5:
N ( u ) = cot ( u ) ; N ( u ) = [ cot ( u ) ] m . The transformation in this case is u = cot 1 ( F ) .
Case 6:
N ( u ) = sinh ( u ) ; N ( u ) = [ sinh ( u ) ] m . The transformation is u = 4 tanh 1 ( F ) .
Case 7:
N ( u ) = cosh ( u ) ; N ( u ) = [ cosh ( u ) ] m . The transformation is u = 4 tanh 1 ( F ) .
Case 8:
N ( u ) = tanh ( u ) ; N ( u ) = [ tanh ( u ) ] m . The transformation is u ( F ) = tanh 1 ( F ) . N ( u ) .
Case 9:
N ( u ) = coth ( u ) ; N ( u ) = [ coth ( u ) ] m . The transformation is u ( F ) = coth 1 ( F ) .
Case 10:
N ( u ) = sin ( m u ) ; N ( u ) = cos ( m u ) . The transformation is u = 4 tan 1 F .
One illustrative example of the application of a transformation in Step 1 of the SEsM is for the equation
b u x x 2 + d u t t 2 = l sin 2 ( u ) .
b, d, and l are parameters. In Step 1 of the SEsM, we use the transformation u = 4 tan 1 ( F ) . The result is the conversion of (6) to an equation that contains only polynomial nonlinearities:
4 F 2 ( b F x 4 + d F t 4 ) 4 ( F + F 3 ) ( b F x 2 F x x + d F t 2 F t t ) + ( 1 + 2 F 2 ) ( b F x x 2 + d F t t 2 ) + F 4 ( b F x x + d F t t ) l ( F 6 2 F 4 + F 2 ) = 0 .
Step 2 of the SEsM is connected with the choice of F as a composite function of simpler functions. In our case, the choice is a specific case of the composite function F ( x , t ) = F [ T 1 ( x , t ) , T 2 ( x , t ) ] , namely F ( x , t ) = A T 1 ( μ ) T 2 ( ξ ) . Here, μ = α x , ξ = γ t .
In Step 3 of the SEsM, one has to determine the form of the functions T 1 and T 2 in such a way that the balance equation leads to a balanced system of nonlinear algebraic relationships. The choice is
T 1 μ 2 = i = 0 N 1 δ i T 1 i ; T 2 ξ 2 = j = 0 N 2 ϵ i T 2 j .
Above, δ i and ϵ i are parameters. The balance equation has a specific solution N 1 = N 2 = 4 . Then, T 1 μ 2 = δ 4 T 1 4 + δ 3 T 1 3 + δ 2 T 1 2 + δ 1 T 1 + δ 0 ; T 2 ξ 2 = ϵ 4 T 2 4 + ϵ 3 T 2 3 + ϵ 2 T 2 2 + ϵ 1 T 1 + ϵ 0 . For simplicity, specific cases of the above simple equations are considered, namely
T 1 μ 2 = p T 1 4 + q T 1 2 + r ; T 2 ξ 2 = s T 2 4 + v T 2 2 + w .
The form of the solution and the choice of the simple equations lead to the following balance system of nonlinear relationships in Step 4 of the SEsM:
b α 4 q 2 + d γ 4 v 2 = l , d γ 4 v A 2 w b α 4 p q = 0 , 4 d γ 4 s w b α 4 q 2 4 b α 4 p r + l d γ 4 v 2 = 0 , b α 4 A 4 r 2 + d γ 4 s 2 = 0 , b α 4 q A 2 r + d γ 4 s v = 0 , b α 4 p 2 + d γ 4 A 4 w 2 = 0 , d γ 4 s v + b α 4 q A 2 r = 0 , d γ 4 v A 2 w + b α 4 p q = 0 .
A nontrivial solution of this system is
p = w = r = s = 0 , q = δ [ b ( d γ 2 v 2 l ) ] 1 / 2 α 2 b , δ = ± 1 , b ( d γ 2 v 2 l ) 0 .
v, A, l, b, d, α , γ are free parameters that satisfy the condition b ( d γ 2 v 2 l ) 0 . The corresponding solution of (6) is
u ( x , t ) = 4 tan 1 A exp δ 1 α v x + γ δ [ b ( d γ 2 v 2 l ) ] 1 / 2 α 2 b t , δ 1 = ± 1 .
We note that in (12), one can group the free parameters. These groups are A, α v , and γ δ [ b ( d γ 2 v 2 l ) ] 1 / 2 α 2 b . Thus, the free parameters for Solution (12) are reduced to three. Let us obtain this result by obtaining (12) in another way. Above, we used a composite function of the kind F ( x , t ) = A T 1 ( μ ) T 2 ( ξ ) ; μ = α x , ξ = γ t , in other words, we searched for standing waves. We can obtain (12) based on a search for the traveling wave solutions of (6), where u ( x , t ) = u ( z ) = u ( x c t ) . c is the velocity of the traveling wave. We introduce the parameter a 2 = l / ( b + d c 2 ) and the velocity becomes
c = l b a 2 d a 2 .
Then, (6) can be factorized: u z z 2 a 2 sin 2 ( u ) = ( u z z + sin ( u ) ) ( u z z sin u ) = 0 , and this leads to the following solution of (6):
u = 4 tan 1 1 + k sn a 1 / 2 x l b a 2 d a 2 t + z 0 , k dn a 1 / 2 x l b a 2 d a 2 t + z 0 , k .
In (14), sn and dn are the Jacobi elliptic functions of modulus k, and z 0 is a constant of integration. Now, let k = 1 . The Jacobi elliptic functions reduce to hyperbolic functions and (14) is reduced to its specific case
u = 4 tan 1 2 exp a 1 / 2 x l b a 2 d a 2 t + z 0 .
Here, we have a and A = 2 exp ( a 1 / 2 z 0 ) as free parameters. However, we have another free parameter and it is connected to the group of parameters l b a 2 d a 2 . For this group of parameters, we have the relationship (13). Let us fix a and z 0 . This does not fix c in (13) as l, b, and d are not fixed. This (15) has effectively three free parameters, similar to (12). We are grateful to the anonymous reviewer who suggested discussing the number of free parameters for the solution (12) using the considerations described in (13)–(15).

4. Composite Functions and Their Role in the Algorithm of SEsM

One encounters composite functions in several places of the algorithm of the SEsM. The presence of composite functions leads to preferences in the choice of the simple equations and the occurrence of specific sets of polynomials in the SEsM. Below, we discuss this in detail.
The differential equations contain derivatives of unknown functions. If the unknown functions are composite, then one has to use the Faa di Bruno formula for the derivative of a composite function. The Faa di Bruno formula is written below for the function h ( x 1 , , x d ) of d independent variables x 1 , , x d , which is a composite function of m other functions g 1 ( 1 ) , , g ( m )
h ( x 1 , , x d ) = f [ g ( 1 ) ( x 1 , , x d ) , , g ( m ) ( x 1 , , x d ) ] .
The following notations are used:
  • ν = ( ν 1 , , ν d ) : d-dimensional index containing integer non-negative numbers ν 1 , , ν d .
  • z = ( z 1 , , z d ) : d-dimensional object containing real numbers z 1 , , z d .
  • ν = i = 1 d ν i : sum of the elements of the d-dimensional index ν .
  • ν ! = i = 1 d ν i ! : factorial of the multicomponent index ν .
  • z ν = i = 1 d z i ν i : ν -th power of the multicomponent variable z .
  • D x ν = ν x 1 ν 1 x d ν d , ν > 0 : ν -th derivative with respect to the multicomponent variable x . Then, the identity operator is denoted as D x 0 .
  • z = max z i : maximum value component of the multicomponent variable z for the interval 1 i d .
  • For the d-dimensional index l = ( l 1 , , l d ) ( l 1 , , l d are integers), it follows that l ν when l i ν i , i = 1 , , d . The following notation is used:
    ν l = i = 1 d ν i l i = ν ! l ! ( ν l ) ! .
  • Ordering of vector indexes: for two vector indexes μ = ( μ 1 , , μ d ) and ν = ( ν 1 , , ν d ) , we have μ ν when one of the following holds:
    (a)
    μ < ν .
    (b)
    μ = ν and μ 1 < ν 1 .
    (c)
    μ = ν , μ 1 = ν 1 , … μ k = ν k and μ k + 1 < ν k + 1 for some 1 k < d .
In addition, the following notation is used:
h ( ν ) = D x ν h ; f ( λ ) = D y λ f ; g ( μ ) ( i ) = D x μ g ( i ) ; g ( μ ) = ( g ( μ ) ( 1 ) , , g ( μ ) ( m ) ) .
Based on the above notation, the Faa di Bruno for the composite derivative of a function containing the functions of many variables is [496]
h ( ν ) = 1 λ n f ( λ ) s = 1 n p s ( ν , λ ) ( ν ! ) j = 1 s [ g ( l j ) ] k j ( k j ! ) [ l j ! ] k j .
In (18) n = ν and
p s ( ν , λ ) = { k 1 , , k s ; l 1 , , l s } , k i > 0 ; 0 l 1 l s , i = 1 s k i = λ , i = 1 s k i l i = ν .
For example, the derivatives of the composite function of two independent variables x 1 and x 2 h ( x 1 , x 2 ) = f [ g ( 1 ) ( x 1 , x 2 ) , , g ( m ) ( x 1 , x 2 ) ] are
h ( ν ) = ν 1 + ν 2 h x 1 ν 1 x 2 ν 2 = 1 ( λ 1 + + λ m ) ν 1 + ν 2 λ 1 + + λ m f g ( 1 ) λ 1 g ( m ) λ m { s = 1 ν 1 + ν 2 p s ( ν , λ ) ( ν 1 ! ν 2 ! ) × j = 1 s 1 ( k j , 1 ! k j , m ! ) ( l j , 1 ! + l j , 2 ! ) k j , 1 + + k j , m i = 1 m l j , 1 + l j , 2 x 1 l j , 1 x 2 l j , 2 g ( i ) k j , i } .
The specific case of a composite function containing one function of one variable h = f [ g ( x ) ] is of special interest, as it is often used in practice. For this case, the Faa di Bruno formula is
h ( n ) = k = 1 n f ( k ) p ( k , n ) n ! i = 1 n g ( i ) λ i ( λ i ! ) ( i ! ) λ i .
In (21), h ( n ) = d n h d x n ; f ( k ) = d k f d g k ; and g ( i ) = d i g d x i . In addition, p ( n , k ) = { λ 1 , λ 2 , , λ n } is a set of numbers satisfying i = 1 n λ i = k ; i = 1 n i λ i = n .
Two results of the composite functions are of interest to the methodology of the SEsM. The first is connected to the use of simple equations for exponential functions [470]
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
P ( h , h ν 1 , , h ν N ) = 0 .
Above, N can be any natural number. We search for the solution of the above equation in the form
h ( x 1 , x 2 ) = f [ g ( 1 ) ( x 1 , x 2 ) , , g ( m ) ( x 1 , x 2 ) ] ,
where h is a 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
g ( x j ) ( i ) = α i , j g ( i ) ,
where α i , j is a constant parameter. Then, the solved nonlinear PDE is reduced to a polynomial of the functions g ( 1 ) ( x 1 , x 2 ) , , g ( m ) ( x 1 , x 2 ) containing monomials multiplied by some coefficients. These coefficients are nonlinear algebraic relations between the parameters of the solved equation and the parameters α i , j . We set the above coefficients to 0. The result is a system of nonlinear algebraic equations. Any nontrivial solution of this system leads to a solution of the solved nonlinear PDE.
The second result is connected to the use of the composite function of a function of a single variable. In this case, one can use a more complicated simple equation in the SEsM. In more detail, we consider a nonlinear partial differential equation with nonlinearities that are polynomials of the unknown function h ( x , t ) and its derivatives. We search for a solution of the kind
h ( x , t ) = h ( ξ ) ; ξ = μ x + ν t ,
where μ and ν are parameters and H is a composite function of another function g:
h = f [ g ( ξ ) ] .
The assumption is that f is a polynomial of g:
f = r = 0 q b r g r ,
and the general form of the simple equation is
g ( k ) l = d k g d ξ k l = j = 0 m a j g j .
Equation (27) defines the function V a 0 , a 1 , , a m ( ξ ; k , l , m ). In this notation, k is the order of derivative of g, l is the degree of the derivative, and m is the highest degree of the polynomial of g in the defining ODE. V has as specific cases the trigonometric, hyperbolic, and elliptic functions of Weierstrass and Jacobi, etc. The result we describe is connected to a specific case of the equation for the V-function V a 0 , a 1 , , a m ( ξ ; 1 , 2 , m ) :
g ( 1 ) 2 = d g d ξ 2 = j = 0 m a j g j .
The result is [486].
Theorem 2.
If g ( 1 ) 2 is given by Equation (28) and f is a polynomial of g given by Equation (26), then for h [ f ( g ) ] , the following relationship holds:
h ( n ) = K n ( q , m ) ( g ) + g ( 1 ) Z n ( q , m ) ( g )
where K n ( q , m ) ( g ) and Z n ( q , m ) ( g ) are polynomials of the function g ( ξ ) .
Note that for some values of n, one of the polynomials K n ( q , m ) or Z n ( q , m ) can be equal to 0.
The theorem states that the derivatives of the composite function h are constructed by specific polynomials. These polynomials are calculated by the recurrence relationships [486]
K n + 1 = Z n 2 j = 0 m j a j g j 1 + d Z n d g j = 0 m a j g j ; Z n + 1 = d K n d g ,
with K 0 = r = 0 q b r g r and Z 0 = 0 . In this way, K 1 = 0 ; Z 1 = r = 0 q r b r g r 1 , and, for example,
K 4 = r = 0 q j = 0 m u = 0 m 1 2 j r + r ( r 1 ) ( j + r 2 ) 1 2 u + j + r 3 a j b r a u g j + r + u 4 , Z 4 = 0 .
Another kind of simple equation leads to another set of polynomials. The simple equation (n and c j are constant parameters)
g ( 1 ) = j = 0 n c j g j ,
is frequently used. In this case [470], one obtains the following set of polynomials. One starts from L 0 = r = 0 q b r q r and uses the recurrence relationship
L i + 1 = d L i d g j = 0 m c j g j .
The derivative of the composite function is h ( n ) = L n ( g ) , for example,
L 3 = r = 0 q j = 0 m k = 0 m l = 0 m b r r ( r + j 1 ) ( r + j + k 2 ) c j c k c l g r + j + k + l 3 .
Let us consider two examples of the application of the SEsM with the participation of composite functions. The first example is for the generalized Korteweg-de Vries equation
u t + A u p u x + 3 u x 3 = 0 ,
In (34), p is a positive integer number and A is a parameter. The solution of Equation (34) is searched as a composite function u = h [ f ( g ( ξ ) ) ] , where ξ = α x + β t , g is the solution of the simplest Equation (28) and f is given by Equation (26). The substitution of u = h [ f ( g ( ξ ) ) ] in Equation (34) and the application of the SEsM leads to an equation of the kind
W 0 ( g ) + W 1 ( g ) g ( 1 ) = 0 .
with
W 0 ( g ) = ν K 1 ( g ) + μ A K 0 ( g ) p K 1 ( g ) + μ 3 K 3 ( g ) , W 1 ( g ) = ν Z 1 ( g ) + μ A K 0 ( g ) p Z 1 ( g ) + μ 3 Z 3 ( g ) .
K 1 = K 3 = 0 and the first equation of (35) is satisfied. The second equation needs balance. It is m = 2 + p q . For the case q = 1 , m = 2 + p . Then, from Equations (26) and (28) one obtains h = b 0 + b 1 g ; g ( 1 ) 2 = j = 0 2 + p a j g j . Thus, the second equation from (36) is reduced to the system of nonlinear algebraic relationships among the parameters of Equation (34) and the parameters of the solution (Step 4 of the SEsM):
ν b 1 δ 0 , k + p k μ A b 0 p k b 1 k + 1 + 1 2 μ 3 ( k + 1 ) ( k + 2 ) a k + 2 b 1 = 0 , k = 0 , , p .
δ is the delta-symbol of Kronecker.
The solution of Equation (37) is obtained for b 0 = 0 . Then, one obtains the system of algebraic equations:
k = 0 : ν + μ 3 a 2 = 0 , k = 1 : a 3 = 0 , k = p 1 : a p + 1 = 0 , k = p : A b 1 p + 1 2 μ 2 ( p + 1 ) ( p + 2 ) a p + 2 = 0 .
The solution is a 2 = ν μ 3 ; a p + 2 = 2 A b 1 p μ 2 ( p + 1 ) ( p + 2 ) . This fixes the simple equation, Equation (28), as
g ( 1 ) 2 = ν μ 3 A g 2 2 b 1 p μ 2 A ( p + 1 ) ( p + 2 ) g p + 2 .
The solution of (39) is g ( ξ ) = Ω cosh ω ( ξ ) , with ω = 2 p ; μ 2 = p 2 A Ω p b 1 p 2 ( p + 1 ) ( p + 2 ) ; ν = 4 μ 3 p 2 . The solution of (34) becomes
u ( ξ ) = Ω b 1 cosh 2 / p ( ξ ) ; ξ = μ x + ν t .
For the particular case p = 1 , A = 6 (34) coincides with the classical Korteweg-de Vries equation. The values b 1 = 1 and Ω = 2 of the parameters lead to μ = 1 , ν = 4 . Then, (40) reduces to the single soliton solution of the KdV equation u ( x , t ) = 2 cosh 2 ( x 4 t ) One can show that (40) is the solution of Equation (34) for an arbitrary real nonzero value of p (and not only for positive integer values of p) [486].
In a similar manner, one can obtain numerous solutions of the Olver equation [486]
u t + u x + α 0 u u x + α 1 u x 2 u x 2 + α 2 u 3 u x 3 + α 3 u 2 u x + α 4 3 u x 3 + α 5 5 u x 5 = 0 .

5. The Role of the Simple Equations in the SEsM

In general, the form of the solution of the solved nonlinear differential equation is connected to the form of the solution of the used simple equations by means of the balance equations, which have to be satisfied in Step 3 of the SEsM. We are not obligated to choose the form of the simple equations before Step 3 of the SEsM. However, if we do, additional restrictions are imposed on the form of the searched solution of the solved equation [483]. Below, we discuss this in more detail.
We consider a simple equation of the kind
Q Φ , d Φ d ξ , d 2 Φ d ξ 2 , , d n Φ d ξ n = 0 .
Using the solution Φ of this simple equation, we build a function F = F ( Φ ) , which has to solve the more complicated equation
P F , d F d ξ , d 2 F d ξ 2 , , d n F d ξ n = 0 .
We can choose Q and F. Thus, we have a function and the question is as follows: What is the manifold of equations that have this function as a solution? Below, are the responses to this question for the case of simple equations of the kind
d Φ d ξ ϵ = α = 0 β γ α [ Φ ( ξ ) ] α ,
where ϵ and β are positive integers and the function F ( Φ ) is a polynomial
F ( Φ ) = μ = 0 ν θ μ [ Φ ( ξ ) ] μ .
For the construction of equations of the kind (43), one needs to know the derivatives of F with respect to ξ . Some derivatives are
d F d ξ = d F d Φ d Φ d ξ .
d 2 F d ξ 2 = d 2 F d Φ 2 d Φ d ξ 2 + d F d Φ d 2 Φ d ξ 2 .
d 3 F d ξ 3 = d 3 F d Φ 3 d Φ d ξ 3 + 3 d 2 F d Φ 2 d Φ d ξ d 2 Φ d ξ 2 + d F d Φ d 3 Φ d ξ 3 .
d 4 F d ξ 4 = d 4 F d Φ 4 d Φ d ξ 4 4 + 6 d 3 F d Φ 3 d Φ d ξ 2 d 2 Φ d ξ 2 + 3 d 2 F d Φ 2 d 2 Φ d ξ 2 2 + 4 d 2 F d Φ 2 d Φ d ξ d 3 Φ d ξ 3 + d F d Φ d 4 Φ d ξ 4 .
d 5 F d ξ 5 = d 5 F d Φ 5 d Φ d ξ 4 5 + 10 d 4 F d Φ 4 d Φ d ξ 3 d 3 Φ d ξ 2 + 15 d 3 F d Φ 3 d Φ d ξ d 2 Φ d ξ 2 2 + 10 d 3 F d Φ 3 d Φ d ξ 2 d 3 Φ d ξ 3 + 10 d 2 F d Φ 2 d 2 Φ d ξ 2 d 3 Φ d ξ 3 + 5 d 2 F d Φ 2 d Φ d ξ d 4 Φ d ξ 4 + d F d Φ d 5 Φ d ξ 5 .
Derivatives (46)–(50) have different forms for different forms of (45) and different versions of the simple equation, Equation (44). This results in different hierarchies of nonlinear PDEs, which have as their solutions the corresponding form of F ( Φ ) . For an illustration of this, we consider as a simple equation the extended tanh-equation
d Φ d ξ = b 2 a 2 Φ 2 ,
which has the solution
Φ ( ξ ) = b a tanh [ a b ( ξ + ξ 0 ) ] .
In (52), a 2 Φ ( ξ ) 2 < d 2 and ξ 0 is a constant of integration. Equation (52) can be used for the construction of traveling wave solutions of numerous nonlinear PDEs. Solution (51) leads to specific forms of Derivatives (46)–(50). The simplest specific case of Equation (45) is d F d Φ = p 1 . Then,
F ( Φ ) = p 0 + p 1 Φ = p 0 + p 1 b a tanh [ a b ( ξ + ξ 0 ) ] .
One obtains Derivatives (46)–(50):
d F d ξ = p 1 a 2 Φ 2 + p 1 b 2 .
d 2 F d ξ 2 = 2 p 1 a 4 Φ 3 2 p 1 a 2 b 2 Φ .
d 3 F d ξ 3 = 6 p 1 a 6 Φ 4 + 8 p 1 a 4 b 2 Φ 2 2 p 1 a 2 b 4 .
d 4 F d ξ 4 = 24 p 1 a 8 Φ 5 40 p 1 a 6 b 2 Φ 3 + 16 p 1 a 4 b 4 Φ .
d 5 F d ξ 5 = 120 p 1 a 10 Φ 6 + 240 p 1 a 8 b 2 Φ 4 136 p 1 a 6 b 4 Φ 2 + 16 p 1 a 4 b 6 .
The function (53) is a solution of numerous nonlinear PDEs. Let us demonstrate the construction of some of these equations. The first equation is
q 1 d 2 F d ξ 2 + q 2 F d F d ξ + q 3 F 3 + q 4 F 2 + q 5 F + q 6 = 0 .
For a solution of the kind (53), the derivatives of F are given by Equations (54) and (55). The substitution of (53)–(55) in Equation (59) leads to the proof of the following proposition.
Proposition 2
([483]). Each solution of the system of nonlinear algebraic relationships
2 q 1 a 4 p 1 2 + q 3 p 1 3 q 2 p 1 3 a 2 = 0 , 2 q 1 a 4 p 1 p 0 3 q 2 p 0 a 2 p 1 2 + 3 q 3 p 0 p 1 2 + q 4 p 1 2 = 0 , q 5 p 1 2 q 2 p 0 2 a 2 p 1 + q 2 p 1 ( b 2 a 2 p 0 2 ) 2 q 1 a 2 p 1 2 b 2 + 3 q 3 p 0 2 p 1 + 2 q 4 p 0 p 1 = 0 , q 2 p 0 ( b 2 a 2 p 0 2 ) + q 6 + q 5 p 0 2 q 1 a 2 p 0 p 1 b 2 + q 3 p 0 3 + q 4 p 0 2 = 0 .
leads to an exact solution of the kind (53) of Equation (59).
Let us set, for example, q 3 = q 4 = q 6 = 0 . Then, we obtain the conditions of a function of the kind (53) to be the solution of the equation
q 1 d 2 F d ξ 2 + q 2 F d F d ξ + q 5 F = 0 .
For this, the system of the following nonlinear algebraic equations must be satisfied
q 2 p 1 3 a 2 + 2 q 1 a 4 p 1 2 = 0 , 3 q 2 p 0 a 2 p 1 2 + 2 q 1 a 4 p 1 p 0 = 0 , 2 q 2 p 0 2 a 2 p 1 + q 2 p 1 ( b 2 a 2 p 0 2 ) 2 q 1 a 2 p 1 2 b 2 + q 5 p 1 = 0 , q 2 p 0 ( b 2 a 2 p 0 2 ) 2 q 1 a 2 p 0 p 1 b 2 + q 5 p 0 = 0 .
Let us now set q 6 = 0 . Thus, we search for a solution of
q 1 d 2 F d ξ 2 + q 2 F d F d ξ + q 3 F 3 + q 4 F 2 + q 5 F = 0 .
One solution of (60) is
p 0 = q 4 2 ( a 2 q 2 q 3 ) ; p 1 = 2 a 4 q 1 a 2 q 2 q 3 ; b = 1 2 ( q 4 2 4 q 3 q 5 + 4 a 2 q 2 q 5 ) ( q 2 q 3 a 2 q 2 2 + 4 q 1 2 a 6 ) .
The substitution of (64) in (53) leads to the following solution of Equation (63)
F ( ξ ) = q 4 2 ( a 2 q 2 q 3 ) + a 3 q 1 a 2 q 2 q 3 ( q 4 2 4 q 3 q 5 + 4 a 2 q 2 q 5 ) ( q 2 q 3 a 2 q 2 2 + 4 q 1 2 a 6 ) tanh a 2 ( q 4 2 4 q 3 q 5 + 4 a 2 q 2 q 5 ) ( q 2 q 3 a 2 q 2 2 + 4 q 1 2 a 6 ) ( ξ + ξ 0 ) .
The solution (65) is an exact solution of each nonlinear partial differential equation that can be reduced to Equation (63).
One can consider more complicated differential equations such as
q 1 d 2 F d ξ 2 + q 2 F d F d ξ + q 3 d F d ξ + q 4 F 3 + q 5 F 2 + q 6 F + q 7 = 0 ,
and to search for conditions for the existence of a solution of the kind (53). One can proceed by using more complicated forms of F ( Φ ) . One such form is d F d Φ = p 1 + 2 p 2 Φ . Then,
F = p 0 + p 1 Φ + p 2 Φ 2 .
Φ ( ξ ) is (52). The substitution of (52) and (67) in (66) leads to a system of seven nonlinear algebraic equations that have to be solved in order to construct a solution of the kind (67) and (52) of (66). A nontrivial solution of this algebraic system is
a = q 4 q 2 ; q 1 = 0 ; q 3 = q 5 q 2 q 4 ; q 6 = q 4 q 7 q 5 ; b = q 4 q 7 q 2 q 5 .
Thus, the equation
q 2 F d F d ξ + q 5 q 2 q 4 d F d ξ + q 4 F 3 + q 5 F 2 + q 4 q 7 q 5 F + q 7 = 0 ,
has the exact solution
F = p 0 + p 1 a q 4 q 7 q 2 q 5 tanh a q 4 q 7 q 2 q 5 ( ξ + ξ 0 ) p 2 q 4 q 7 a 2 q 2 q 5 tanh 2 a q 4 q 7 q 2 q 5 ( ξ + ξ 0 ) .
All nonlinear partial differential equations that can be reduced to (69) have the solution (70).
We can continue the process by taking, for example, the equation of Bernoulli as the simple equation. The interested reader is directed to [483] for more details on this.

6. The SEsM and the Method of Hirota and Inverse Scattering Transform Method

The SEsM can be connected to other famous methods for obtaining exact solutions of nonlinear partial differential equations such as the method of Hirota and the Inverse Scattering Transform (IST) method. Below, we discuss the connection between the SEsM and the method of Hirota and between the IST and the SEsM for the cases of the Korteweg-de Vries equation and the nonlinear Scrödinger equation.

6.1. The SEsM and the Method of Hirota

The Hirota method [136,137,138,139,140,141] is a simple and popular method for obtaining multisoliton solutions of nonlinear PDEs, which are integrable. Usually, in Step 1 of this method, a transformation of the solved NPDE is made. Frequently, the form of this transformation is
u ( x , t ) = 2 2 x 2 f ( x , t ) = 2 f 2 f x 2 f x 2 f 2 .
In Step 2 of the method, one searches for a solution in the form
f = α + ϵ f 1 + ϵ 2 f 2 + ϵ 3 f 3 + .
where α and ϵ are parameters. Solution (72) is substituted in the solved nonlinear PDE. In Step 3 of the method, one solves the obtained equations for the orders ϵ , ϵ 2 , ϵ 3 , …. In Step 4 of the method, the solution is constructed based on the solutions of the equations for the different powers of ϵ .
Hirota introduced the bilinear operators, which can be written as [497]
D x m ( a . b ) = j = 0 m ( 1 ) ( m j ) m ! j ! ( m j ) ! j a x j m j b x m j , D x m D t n ( a . b ) = j = 0 m i = 0 n ( 1 ) ( m + n j i ) m ! j ! ( m j ) ! n ! i ! ( n i ) ! i + j a t i x j m + n i j b t n i x m j .
These operators allow bilinear operators to have useful properties such as finding a solution for the sequence of equations for the powers of ϵ in an effective way.
The Hirota method is connected to the SEsM in the following way [460]:
Proposition 3.
The method of Hirota is connected to a specific case of the SEsM where the transformation in Step 1 of the SEsM is the same as the transformation used in the Hirota method, the representation of function f using f 1 , f 2 , from Step 2 of the SEsM is (72), the differential equations for f 1 , f 2 , from Step 3 of the SEsM are the chain of equations obtained for the orders ϵ, ϵ 2 , … within the scope of the Hirota method, and the simple equations that are used in the construction of the solutions for f 1 , f 2 , are differential equations for exponential functions.
One example of the demonstration of the connection between the SEsM and the method of Hirota is in obtaining the three-soliton solution of the Korteweg-de Vries equation:
u t + 6 u u x + 3 u x 3 = 0 .
Step 1 of the SEsM is to transform the nonlinearity of this equation using (71). We obtain
f 2 f x t + f x f t + f 4 f x 4 4 f x 3 f x 3 + 3 2 f x 2 2 = 0 ,
which can be written as ( D x D t + D x 4 ) ( f · f ) = 0 . Step 2 of the SEsM is to fix the composite function for the solution of the equation. This is made by (72). The substitution of (72) in (75) leads to the differential equations for f 1 , f 2 , …:
2 x t + 3 x 3 f 1 = 0 ,
2 x t + 3 x 3 f 2 = D x ( D t + D x 3 ) ( f 1 · f 1 ) ,
2 x t + 3 x 3 f 3 = D x ( D t + D x 3 ) ( f 1 · f 2 + f 2 · f 1 ) ,
etc. The next step of the SEsM is to choose the simple equations for f 1 , f 2 , …. These are equations for exponential functions. One starts with d g 1 d η 1 = g 1 ; g 1 = exp ( η 1 ) , where η 1 = λ 1 x + ω 1 t + σ 1 and λ 1 , ω 1 and σ 1 are parameters and chooses f 1 as
f 1 = exp ( η 1 ) .
The substitution of (79) in (76) leads to the dispersion relation ω 1 + λ 1 3 = 0 . f 1 satisfies (77) and f 2 can be taken to be 0. f 3 , f 4 , can also be set to 0. One obtains
f = 1 + ϵ f 1 .
The parameter ϵ is absorbed by σ 1 and the solution for f corresponds to the single soliton solution for u by the use of (71).
For obtaining the two-soliton solution of the Korteweg-de Vries equation, one uses the solution of (76) constructed by the solutions of two simple equations d g 1 d η 1 = g 1 ; g 1 = exp ( η 1 ) ; d g 2 d η 2 = g 2 ; g 2 = exp ( η 2 ) . Here, η i = λ i x + ω i t + σ i and λ i , ω i and σ i , i = 1 , 2 are parameters. Then,
f 1 = exp ( η 1 ) + exp ( η 2 ) .
Two algebraic equations occur from the substitution of (81) in (76): ω i + λ i 3 = 0 , i = 1 , 2 . The substitution of (81) in (77) leads to an equation for f 2 , for which the solution is f 2 = a 12 exp ( η 1 + η 2 ) , with a 12 = ( λ 1 λ 2 ) 2 ( λ 1 + λ 2 ) 2 . f 3 , f 4 , etc., can be set to zero and one arrives at the two-soliton solution of the KdV equation
f = 1 + ϵ [ exp ( η 1 ) + exp ( η 2 ) ] + ϵ 2 a 12 exp ( η 1 + η 2 ) .
Here, ϵ and ϵ 2 are absorbed by σ 1 , 2 and a 1 , 2 .
For obtaining the three-soliton solution of the Korteweg-de Vries equation, one needs a solution for f 1 constructed using the solutions of three simple equations for the exponential functions d g i d η i = g i ; g i = exp ( η i ) ; i = 1 , 2 , 3 . Here, η i = λ i x + ω i t + σ i and λ i , ω i and σ i , i = 1 , 2 , 3 are parameters. f 1 is the sum of the solutions of the three simple equations
f 1 = exp ( η 1 ) + exp ( η 2 ) + exp ( η 3 ) .
The substitution of (83) in (76) leads to three dispersion relations ω i + λ i 3 = 0 , i = 1 , 2 , 3 . For f 2 , after a substitution of (83) in (77), one obtains f 2 = a 12 exp ( η 1 + η 2 ) + a 13 exp ( η 1 + η 3 ) + a 23 exp ( η 2 + η 3 ) , with a i j = ( λ i λ j ) 2 ( λ i + λ j ) 2 , i , j = 1 , 2 , 3 , i < j . The substitution of f 1 and f 2 in (78) leads to the following relationship for f 3 : f 3 = b 123 exp ( η 1 + η 2 + η 3 ) , b 123 = a 12 a 13 a 23 . This leads to f 4 = 0 , f 5 = f 6 = = 0 . The solution for f is
f = 1 + ϵ [ exp ( η 1 ) + exp ( η 2 ) + exp ( η 3 ) ] + ϵ 2 [ a 12 exp ( η 1 + η 2 ) + a 13 exp ( η 1 + η 3 ) + a 23 exp ( η 2 + η 3 ] + ϵ 3 b 123 exp ( η 1 + η 2 + η 3 ) ,
and this leads to the three-soliton solution of the Korteweg-de Vries equation.

6.2. The SEsM and the Inverse Scattering Transform Method: The Case of the Korteweg-de Vries Equation

Due to the excellent work of Rosales [498,499], we can easily connect the SEsM to various equations that are solvable by the Inverse Scattering Transform method. We demonstrate this for the Korteweg-de Vries equation and the nonlinear Schrödinger equation. We consider the KdV equation
u t + 6 u u x + 3 u x 3 = 0 .
We skip Step 1 of the SEsM (the transformation of the nonlinearity). In Step 2 of the SEsM, we present u ( x , t ) as a composite function of other functions u 1 , u 2 ,
u ( x , t ) = n = 1 ϵ n u n ( x , t ) .
Initially, one treats ϵ as a small parameter. The substitution of (86) in (85) leads to equations for u 1 , u 2 , as follows:
u n t + 3 u n x 3 = 6 x j = 1 n 1 u j u n j , n = 2 , 3 ,
and the equation for u 1 is
u 1 t + 3 u 1 x 3 = 0 .
The solutions of the obtained equations for u 1 , u 2 , … are connected through a Fourier series to the solution of the simple equation for exponential functions. For u 1 we use
u 1 = C d λ ( k ) ( k ) exp [ i k x + k 3 t ] ,
where d λ ( k ) is an appropriate measure in the complex plane C and the term ( k ) is introduced for convenience. The form of u 1 influences the relationships for the equations for u 2 , u 3 , …. One obtains
u 2 t + 3 u 2 x 3 = 6 i C 2 d λ ( k 1 ) d λ ( k 2 ) k 1 k 2 ( k 1 + k 2 ) exp { i [ ( k 1 + k 2 ) x + ( k 1 3 + k 2 3 ) t ] } ,
u 3 t + 3 u 3 x 3 = 6 x ( u 1 u 2 + u 2 u 3 ) = = 6 i C 3 d λ ( k 1 ) d λ ( k 2 ) d λ ( k 3 ) ( k 1 + k 2 ) ( k 1 + k 2 + k 3 ) exp { i [ ( k 1 + k 2 + k 3 ) x + ( k 1 3 + k 2 3 + k 3 3 ) t ] } .
The solutions for u 2 , u 3 , are assumed to be of a form similar to that of u 1 :
u n ( x , t ) = C n [ d λ ( k ) ] n Φ n ( k 1 , , k n ) exp { i Ω n } .
The substitution of (92) in the equations for u 2 , u 3 , leads to the fixation of Φ n and Ω n . The fixation leads to
Φ 2 ( k 1 , k 2 ) = 1 ; Ω 2 = ( k 1 + k 2 ) x + ( k 1 3 + k 2 3 ) t , Φ 3 ( K 1 , k 2 , k 3 ) = k 1 + k 2 + k 3 ( k 1 + k 2 ) ( k 2 + k 3 ) ; Ω 3 = ( k 1 + k 2 + k 3 ) x + ( k 1 2 + k 2 2 + k 3 2 ) t , Φ 4 = k 1 + k 2 + k 3 + k 4 ( k 1 + k 2 ) ( k 2 + k 3 ) ( k 3 + k 4 ) ; Ω 4 = ( k 1 + k 2 + k 3 + k 4 ) x + ( k 1 2 + k 2 2 + k 3 2 + k 4 2 ) t .
This leads to the following relationship for u n
u n = ( 1 ) n + 1 i x C n [ d λ ( k ) ] n exp ( i Ω n ) j = 1 n 1 ( k j + k j + 1 ) ,
and the solution of the Korteweg-de Vries equation is
u = i x n = 1 ( ϵ ) n C n [ d λ ( k ) ] n exp ( i Ω n ) j = 1 n 1 ( k j + k j + 1 ) .
One rewrites (95) as
u ( x , t ) = x n = 1 n ( ϵ ) n C n d [ λ ( k ) ] n p ^ ( k 1 ) P ^ ( k 1 , k 2 ) P ^ ( k 2 , k 3 ) , P ^ ( k n 1 , k n ) p ^ ( k n ) ,
where p ^ k = exp i 2 [ k x ω ( k ) t ] ; P ^ ( k , q ) = i p ^ ( k ) p ^ ( q ) k + q ; ω ( k ) = k 3 . Two possibilities exist for the measure d λ ( k ) , discrete or continuous. For the case of discrete d λ ( k ) , the integral can be replaced by a sum: C d λ ( k ) f ( k ) = m a m 2 f ( i k m ) . Thus,
u = x n = 1 ( ϵ ) n m j p m 1 P m 1 , m 2 P m 2 , m 3 P m n 1 , m n p m n .
Here, P m q = p m p q k m + k q ; p m = a m p ^ ( i k m ) = a m exp 1 2 ( k m x + k m 3 t ) , and (97) can be written as
u = ϵ x p T ( I + ϵ P ) 1 p ,
where p is the column vector of all p m , P is the square matrix of all P m q , and p T is the transpose of p. The matrix P is real symmetric and positive definite for k m > 0 and real for a m . Thus, u is nonsingular for ϵ > 0 and ϵ can be absorbed into the coefficients a m . Because of this, Solution (98) is not limited to only small values of ϵ .
Solution (98) can be rewritten in a more common form. Observing that P x = 1 2 p p T , we arrive at
u = 2 x Tr [ ( I + P ) 1 P x ] = 2 2 x 2 Tr [ ln ( I + P ) ] = 2 2 x 2 ln det ( I + P ) .
This is the relationship for the multisoliton solution of the Korteweg-de Vries equation.
For the case of continuous measure d λ ( k ) , we obtain the Gelfand–Levitan–Marchenko equation from the methodology of the Inverse Scattering Transform method as follows. P is considered more a general operator than a square matrix: ( P f ) ( k ) = C d λ ( l ) P ^ ( k , l ) f ( l ) . Thus, the solution of the KdV equation is
u = ϵ x C d λ ( k ) p ^ ( k ) [ ( I + ϵ P ) 1 p ( k ) .
Solution (100) can be written as
u = 2 x K ( x , x ) .
Here, K ( x , y ) = ϵ 2 p T ( x ) [ I + ϵ P ( x ) ] 1 p ( y ) , and after a calculation, one obtains
K ( x , y ) = ϵ 2 p T ( x ) p ( y ) ϵ 2 x d z K ( x , z ) p T ( z ) p ( y ) .
This is the GLM equation with
B ( x + y ) = ϵ 2 p T ( x ) p ( y ) = ϵ 2 d λ ( k ) exp i k 2 ( x + y ) + k 3 t .

6.3. The SEsM and the Inverse Scattering Transform Method: The Case of the Nonlinear Schrödinger Equation

We consider the nonlinear Schrödinger equation
i φ t + 2 φ x 2 + 2 σ φ 2 φ = 0 .
Step 1 of the SEsM (the transformation of the nonlinearity) is skipped. In Step 2 of the SEsM, we present φ ( x , t ) as a composite function similar to (86). Performing the calculations for the solutions of the equations for the components of φ using differential equations for exponential functions as simple equations (Steps 3 and 4 of the SEsM) we arrive at
φ = n = 0 ; m = 2 n + 1 d λ 1 d μ 2 d λ 3 d μ 4 d λ m σ n ϵ m C m exp ( i Ω m ) j = 1 2 n ( k j + k j + 1 ) .
In (105), Ω n 1 n [ k j x + ( 1 ) j k j 2 t ] . d λ ( k ) and d μ ( k ) are measures such that d ν * ( k * ) = d λ ( k ) ( * denotes complex conjugation). d λ j = d λ ( k j ) ; d μ j = d μ ( k j ) . There is a quantity analogous to P from the case of the Korteweg-de Vries equation above. The definition of this quantity for the case of the nonlinear Schrödinger equation is
P ^ ( k , l ) = i p ^ * ( k ) p ^ ( l ) k * + l .
In Ω n , one has two dispersion relationships: ω 1 ( k ) = k 2 and ω 2 ( k ) = k 2 . In p ^ , we take the first of them. Denoting the conjugate operator of P as P ¯ , one can rewrite (105) as
φ = ϵ C d λ ( k ) p ^ ( k ) ( I + σ ϵ 2 P ¯ P ) 1 ( k ) = ϵ p T ( I + σ ϵ 2 P ¯ P ) 1 p .
For φ 2 , one obtains
φ 2 = ϵ 2 x p T ( I + σ ϵ 2 P ¯ P ) 1 P ¯ p ¯ .
This is reduced to the relationship of Zakharov and Shabat [500] for the multiple-envelope solutions of the nonlinear Schrödinger equation for the specific case of (108) of discrete measure (discrete λ ( k ) ) and when p is a vector and P is a square matrix. In this case,
φ 2 = 1 σ 2 x 2 ln det ( I + σ ϵ 2 P ¯ P ) .

7. Special Application: The SEsM and the SIR Model of Epidemics

Below we discuss the exact solutions of an equation. This equation is obtained based on the SIR model in epidemiology. The obtained solutions can describe epidemic waves caused by different diseases (COVID-19-inclusive).
The equation is obtained using the classic idea of Kermak and McKendrick [501] for the transformation of the SIR model with constant coefficients to a single nonlinear differential equation. We consider an epidemic in a population. The population is divided into three groups: susceptible individuals—S; infected individuals—I; recovered individuals—R. The model equations for the time changes in the numbers of individuals in the above three groups are
d S d t = τ N S I , d I d t = τ N S I ρ I , d R d t = ρ I .
In (refsir1), τ is the transmission rate and ρ is the recovery rate. These rates are assumed to be constants.
N = S + I + R ,
is the total population, which is assumed to be constant. The model (110) can be reduced to a single equation for R as follows. From the last equation of (110) we have
I = 1 ρ d R d t .
The substitution of (112) in the first equation of (110) leads to the relationship
S = S ( 0 ) exp τ ρ N [ R R ( 0 ) ] .
Here, S ( 0 ) and R ( 0 ) are the numbers of susceptible and recovered individuals at the time t = 0 . The substitution of (111) and (113) in the last equation of (110) leads to the differential equation for R
d R d t = ρ N R S ( 0 ) exp τ ρ N ( R R ( 0 ) ) .
Below, we assume R ( 0 ) = 0 (no recovered individuals at t = 0 ). Let us consider the ratio τ R ρ N . The maximum of the ratio R / N is 1 and if τ < < ρ , then τ R ρ N < < 1 . Let us consider a kind of epidemic, where τ < < ρ (transmission rate is much lower than the recovery rate). Then, exp τ ρ N R can be represented as a Taylor series
exp τ ρ N R = j = 0 M τ ρ N R j .
M has infinite value in the full Taylor series but we can truncate it at M = 2 , M = 3 ,..., if τ ρ N R is small enough. From (114) we obtain
d R d t = ρ N R S ( 0 ) j = 0 M τ ρ N R j , M = 2 , 3 ,
We apply the most simple version of the SEsM to (116). The solution is searched as a composite function of a single function U and the form of the composite function is a polynomial
R = k = 0 K μ k U k ,
where μ l and K are parameters. We also fix the form of the simple equation for U
d U d t = l = 0 L ν l U l ,
where ν l and L are parameters.
In order to proceed with the application of the SEsM, we have to obtain a balance equation that connects the parameters M, K, and L. The form of this balance equation is
L 1 = K ( M 1 ) .
Then, we proceed as follows. We fix the value of M. This fixes the form of (116). Then, we fix the value of K. This fixes the form of the composite function for the solution of the equation. Then, from (119) we obtain the value of L. The fixation of the parameters M, K, and L and the substitution of (117) and (118) in (116) lead us to the system of nonlinear algebraic equations as required by the algorithm of the SEsM. The nontrivial solutions of this system lead to solutions of R. The substitution of these solutions in (112) leads to the corresponding solutions for I. Finally, the solutions for R and I are substituted in (111) and we obtain the solution for S.
Below, we consider the simplest case M = 2 . This is the classical Kermak–McKendrick case. From (116), the equation we have to solve is
d R d t = α R 2 + β R + γ ,
where γ = ρ [ N S ( 0 ) ] ; β = τ S ( 0 ) N ρ ; α = τ 2 S ( 0 ) ρ N 2 , and (120) is an equation of the Riccati kind. From (119),
L = K + 1 .
Let us consider the case K = 1 , L = 2 . In this case,
R = μ 0 + μ 1 U ,
and the simple equation is
d U d t = ν 0 + ν 1 U + ν 2 U 2 .
The substitution of (122) and (123) in (120) leads to a system of nonlinear algebraic equations. One nontrivial solution of this system is μ 0 = 0 , μ 1 = 1 ; ν 0 = γ ; ν 1 = β , ν 2 = α . The simple equation, Equation (123), coincides with (120). As we consider the case τ < < ρ then, β 2 4 α γ ρ [ ρ 2 τ S ( 0 ) N ] > 0 . This is because S ( 0 ) < < N . In addition, we have β 2 4 α γ = τ 2 S ( 0 ) 2 N 2 1 + 4 N S ( 0 ) S ( 0 ) > τ 2 S ( 0 ) 2 N 2 1 2 τ R ρ N 2 = ( 2 α R + β ) 2 . A particular solution of (120) is
R ( t ) = β 2 α θ 2 α tanh θ ( t + C ) 2 ,
where θ 2 = β 2 4 α γ > 0 and C is a constant of integration. The substitution of (124) in (112) leads to the well-known bell-shaped curve for the evolution of the number of infected individuals I
I = θ 2 2 α ρ sech 2 θ ( t + C ) 2 .
Note that θ 2 2 α ρ is positive as α is negative.
We know a particular solution (124) of (120). Then, we can write the general solution of (120) as R = β 2 α θ 2 α tanh θ ( t + C ) 2 + D v , where D is a constant and v ( t ) is the solution of the linear differential equation
d v d t θ tanh θ ( t + C ) 2 v = α D .
The solution of (126) is
v = cosh 2 θ ( t + C 2 E 2 α D θ tanh θ ( t + C 2 ,
where E is a constant of integration. Then, the general solution of Equation (120) is
R ( t ) = β 2 α θ 2 α tanh θ ( t + C ) 2 + D cosh 2 θ ( t + C 2 E 2 α D θ tanh θ ( t + C 2 .
This solution introduces corrections to the well-known curve for the time evolution of recovered persons within the scope of the SIR model. The substitution of (128) in (112) leads to the following relationship for the time evolution of the number of infected persons:
I = θ 2 2 α ρ sech 2 θ ( t + C ) 2 D ρ θ E sinh θ ( t + C ) α D cosh [ θ ( t + C ) ] cosh 4 θ ( t + C ) 2 E 2 α D θ tanh θ ( t + C ) 2 2 .
This solution leads to corrections to the well-known bell-shaped curve for the time evolution of infected individuals in the SIR model. For example, let t = C . Then, from (125),
I = θ 2 2 α ρ ,
and from (129),
I = θ 2 2 α ρ α D 2 ρ E .
The presence of the general solution (128) of (116) makes the further consideration of the subcases of the case M = 2 unnecessary. What we can obtain by considering these cases are the specific cases of the general solution (128).

8. Other Examples of the Application of the SEsM

The SEsM and its most simple version, the MMSE, have numerous applications. Below, we mention several results. In most cases, we search for traveling wave solutions of the kind u ( x , t ) = u ( ξ ) = u ( x v t ) , and the composite function from Step 2 of the SEsM is given by the power series of the solution V of the simple equation
u ( ξ ) = μ = ν η θ μ [ V ( ξ ) ] μ .
In (132), ν > 0 , μ > 0 , θ μ are parameters. In most cases, we use as the simple equation
d V d ξ = α = 0 β γ α V ( ξ ) α ,
where γ α is a parameter. Equation (133) has two important specific cases: the equations of Bernoulli and Riccati. Their solutions are expressed by elementary functions. For the Bernoulli equation γ 0 = 0 ; γ 1 = a ; γ 2 = = γ k 1 = 0 ; γ k = b ; γ k + 1 = = γ β = 0 . The equation is
d V d ξ = a V ( ξ ) + b [ V ( ξ ) ] k ,
where k is an integer and k > 1 . The solutions that are used are
V ( ξ ) = a exp [ a ( k 1 ) ( ξ + ξ 0 ) ] 1 b exp [ a ( k 1 ) ( ξ + ξ 0 ) ] k 1 , b < 0 , a > 0 ; V ( ξ ) = a exp [ a ( k 1 ) ( ξ + ξ 0 ) ] 1 + b exp [ a ( k 1 ) ( ξ + ξ 0 ) ] k 1 , b > 0 , a < 0 .
ξ 0 is a constant of integration.
For the Riccati equation, γ 0 = c ; γ 1 = b ; γ 2 = a , γ 3 = = γ β = 0 . The equation is
d V d ξ = a [ V ( ξ ) ] 2 + b V ( ξ ) + c .
The solution that is used is
V ( ξ ) = b 2 a θ 2 a tanh θ ( ξ + ξ 0 ) 2 .
Here, θ 2 = b 2 4 a c > 0 and ξ 0 is a constant of integration.
As a specific case of the use of the equation of Riccati, we also consider the extended tanh-function equation
d V d ξ = c ¯ 2 a ¯ 2 V 2 .
Equation (138) is obtained from (136) when b = 0 , a = a ¯ 2 , and c = c ¯ 2 . The solution of (138) that is used is
V ( ξ ) = c ¯ a ¯ tanh [ a ¯ c ¯ ( ξ + ξ 0 ) ] ,
where a ¯ 2 V ( ξ ) 2 < c ¯ 2 and ξ 0 is a constant of integration.
The first illustration of the methodology of the SEsM is for obtaining traveling wave solutions of the class of equations
p = 1 N 1 α p p Q t p + q = 1 N 2 β q q Q x q + m = 1 M μ m Q m = 0 .
In (140), α p , β q , and μ m are parameters. Specific cases of (140) are used as model equations in the theory of the migration of populations. If N 1 = 1 , N 2 = 2 , α 1 = 1 , β 1 = 0 , β 2 = 1 , (140) is reduced to an equation of the kind Q t 2 Q x 2 = f ( Q ) , with f ( Q ) = m = 1 M μ m Q m . Such kinds of equations are widely used in the modeling of genetic diffusion.
One searches for traveling wave solutions Q ( x , t ) = Q ( ξ ) = Q ( x v t ) . Taking into account that N = max ( N 1 , N 2 ) , one can extend (140) to
n = 1 N α n n Q t n + β n n Q x n + m = 1 M μ m Q m = 0 ,
where the coefficients of the additional terms are equal to 0. The traveling wave solution and ν n = β n + α n ( v ) n reduce (141) to
n = 1 N ν n d n Q d ξ n + m = 1 M μ m Q m = 0 .
As the simple equation, we use the equation of Bernoulli (134). The substitution of (132) and (134) in (142) leads to a polynomial P = κ r ϕ r + κ r 1 ϕ r 1 + + κ 0 = 0 . In P, r is an integer. The coefficients κ depend on the parameters of the solved equation and the parameters of the solution. The solution is obtained when all coefficients are set to 0. The result is a system of nonlinear algebraic relationships
κ l = 0 , l = r , r 1 , , 0
Next, we move to Step 3 of the SEsM, ensuring that each equation from (143) has at least two terms by balancing the highest powers arising in (142). This leads to the following balance equation (for the case when the Bernoulli equation is used as the simple equation): L M = L + N ( k 1 ) . Thus, the choice of the highest power L of the solution polynomial leads to the fixation of the parameter k for the Bernoulli equation
k = 1 + L ( M 1 ) N .
If the Riccati equation is used as the simple equation, then the balance equation is L ( M 1 ) = N .
Several examples of exact solutions follow for cases where the Bernoulli equation is used as the simple equation. We first set N = 2 , M = 2 , and (142) becomes
ν 1 d 2 Q d ξ 2 + ν 2 d Q d ξ + μ 2 Q 2 + μ 1 Q = 0 .
A specific case of this equation is the equation studied by Ablowitz and Zeppetela [502] d 2 Q d ξ 2 + v d Q d ξ + Q ( 1 Q ) = 0 . Ablowitz and Zeppetela obtained the solution Q ( ξ ) = 1 [ 1 r exp ( ξ / 6 ) ] 2 , r < 0 , v = 5 / 6 . Another solution of the last equation is obtained using the tanh-method. This solution is [503] Q ( ξ ) = Q ( x v t ) = 1 4 1 tanh 1 2 6 ξ 2 , v = 5 6 .
We search for a solution of (145) of the kind Q ( ξ ) = i = 0 L a i ϕ i ; d ϕ d ξ = b ϕ ( 1 + L / 2 ) + a ϕ . L = 2 is assumed and for this case, in Step 4 of the SEsM one obtains the system of nonlinear algebraic relationships between the parameters of the solutions and the parameters of (142)
6 ν 1 a 2 b 2 + μ 2 a 2 2 = 0 , 2 ν 2 a 2 b + 10 ν 1 a 2 b a + 2 μ 2 a 2 a 1 + 2 ν 1 a 1 b 2 = 0 , μ 1 a 2 + μ 2 a 1 2 + 2 ν 2 a 2 a + ν 2 a 1 b + 4 ν 1 a 2 a 2 + 2 μ 2 a 2 a 0 + 3 ν 1 a 1 a b = 0 , μ 1 a 1 + 2 μ 2 a 1 a 0 + ν 1 a 1 a 2 + ν 2 a 1 a = 0 , μ 1 a 0 + μ 2 a 0 2 = 0 .
A solution of (146) is a 0 = 6 25 ν 2 2 ν 1 μ 2 , a 1 = 0 , a 2 = 6 ν 1 b 2 μ 2 , μ 1 = 6 25 ν 2 2 ν 1 . The corresponding solution of (145) is
Q ( ξ ) = 6 25 ν 2 2 ν 1 μ 2 6 ν 1 b 2 μ 2 a exp [ a ( ξ + ξ 0 ) ] 1 b exp [ a ( ξ + ξ 0 ) ] 2 .
for the case when b < 0 and a > 0 . Another solution of (146) is a 0 = a 1 2 4 a 2 , ν 2 = 5 μ 1 6 a , μ 2 = 4 μ 1 a 2 a 1 2 , b = 2 a a 2 a 1 . The corresponding solution of (142) is
Q ( ξ ) = a 2 a exp ( a ( ξ + ξ 0 ) ) 1 2 a a 2 a 1 exp ( ξ + ξ 0 ) 2 + a 1 a exp ( a ( ξ + ξ 0 ) ) 1 2 a a 2 a 1 exp ( ξ + ξ 0 ) + a 1 2 4 a 2 .
Next, we set L = 4 ( L = 3 is impossible as k has to be an integer). In Step 3 of the SEsM, we obtain the system of nonlinear algebraic relationships between the parameters
24 ν 1 a 4 b 2 + μ 2 a 4 2 = 0 , 2 μ 2 a 4 a 3 + 15 ν 1 a 3 b 2 = 0 , 4 ν 2 a 4 b + 40 ν 1 a 4 b a + 8 ν 1 a 2 b 2 + 2 μ 2 a 4 a 2 + μ 2 a 3 2 = 0 , 24 ν 1 a 3 b a + 2 μ 2 a 4 a 1 + 3 ν 1 a 1 b 2 + 3 ν 2 a 3 b + 2 μ 2 a 3 a 2 = 0 , 2 μ 2 a 4 a 0 + 16 ν 1 a 4 a 2 + 4 ν 2 a 4 a + μ 1 a 4 + 2 μ 2 a 3 a 1 + μ 2 a 2 2 + 2 ν 2 a 2 b + 12 ν 1 a 2 b a = 0 , μ 1 a 3 + 9 ν 1 a 3 a 2 + 3 ν 2 a 3 a + ν 2 a 1 b + 2 μ 2 a 2 a 1 + 2 μ 2 a 3 a 0 + 4 ν 1 a 1 b a = 0 , 4 ν 1 a 2 a 2 + μ 1 a 2 + 2 μ 2 a 2 a 0 + 2 ν 2 a 2 a + μ 2 a 1 2 = 0 , ν 2 a 1 a + ν 1 a 1 a 2 + 2 μ 2 a 1 a 0 + μ 1 a 1 = 0 , μ 1 a 0 + μ 2 a 0 2 = 0 .
One solution of this system is a 0 = 6 25 ν 2 2 ν 1 μ 2 , a 1 = a 2 = a 3 = 0 , a 4 = 24 ν 1 b 2 μ 2 , a = 1 10 ν 2 ν 1 , μ 1 = 6 25 ν 2 2 ν 1 . The corresponding solution of (145) ( b > 0 , a < 0 ) is
Q ( ξ ) = 6 25 ν 2 2 ν 1 μ 2 1 b exp [ ν 2 ( ξ + ξ 0 ) 5 ν 1 ] 1 + b exp [ ν 2 ( ξ + ξ 0 ) 5 ν 1 ] 2 .
Another solution of (149) is a 4 = b a 2 2 a , a 3 = 0 , a 1 = 0 , a 0 = a a 2 2 b , μ 1 = 12 5 a ν 2 , μ 2 = 24 5 b ν 2 a 2 , ν 1 = ν 2 10 a . The corresponding solution of (145) is
Q ( ξ ) = b a 2 2 a a exp [ 2 a ( ξ + ξ 0 ) ] 1 + b exp [ 2 a ( ξ + ξ 0 ) ] 2 + a 2 a exp [ 2 a ( ξ + ξ 0 ) ] 1 + b exp [ 2 a ( ξ + ξ 0 ) ] + a a 2 2 b .
Here, b > 0 , a < 0 .
Next, one can set N = 2 , M = 3 . The equation that is to be solved is
ν 1 d 2 Q d ξ 2 + ν 2 d Q d ξ + μ 3 Q 3 + μ 2 Q 2 + μ 1 Q = 0 .
This equation is connected to the Huxley equation for gene propagation [504]. Other equations having polynomial nonlinearity of the third order are the Burgers–Huxley equation and equations that model aspects of neuronal activity [505].
Several solutions of (152) are below. For the case L = 1 , the SEsM leads to a solution of the kind
Q ( ξ ) = i = 0 L a i ϕ ( ξ ) i , d Q d ξ = b ϕ 1 + L + a ϕ .
For the case b > 0 , a < 0 , the solution is
Q ( ξ ) = b 2 μ 3 μ 1 μ 3 ( 2 μ 2 ν 1 ν 2 2 μ 3 ν 1 3 ν 1 2 μ 3 ν 1 ) exp [ ( 2 μ 2 ν 1 ν 2 2 μ 3 ν 1 3 ν 1 2 μ 3 ν 1 ) ( ξ + ξ 0 ) ] 1 + b exp [ ( 2 μ 2 ν 1 ν 2 2 μ 3 ν 1 3 ν 1 2 μ 3 ν 1 ) ( ξ + ξ 0 ) ] .
For L = 2 , the solution of (152) for the case a > 0 , b < 0 is
Q ( ξ ) = 2 b 2 ν 1 μ 3 μ 3 2 ν 1 μ 1 μ 3 ν 2 ν 1 μ 3 6 ν 1 ν 1 μ 3 exp [ 2 2 ν 1 μ 1 μ 3 ν 2 ν 1 μ 3 6 ν 1 ν 1 μ 3 ( ξ + ξ 0 ) ] 1 b exp [ 2 2 ν 1 μ 1 μ 3 ν 2 ν 1 μ 3 6 ν 1 ν 1 μ 3 ( ξ + ξ 0 ) ] μ 1 μ 3 μ 3 .
One can continue by assuming N = 2 , M = 4 . The equation that has to be solved is
ν 1 d 2 Q d ξ 2 + ν 2 d Q d ξ + μ 4 Q 4 + μ 3 Q 3 + μ 2 Q 2 + μ 1 Q = 0 .
Equations of this kind arise in the theory of the migration of populations [476,477]. The application of the SEsM leads to the conclusion that this equation has solutions of the form Q ( ξ ) = i = 0 L a i [ ϕ ( ξ ) ] i , d ϕ d ξ = b ϕ ( 1 + 3 L / 2 ) + a ϕ . For L = 2 , one possible solution for the case a < 0 , b > 0 is
Q ( ξ ) = 10 ν 1 b 2 3 μ 3 ν 1 μ 3 8 10 ν 1 μ 4 exp 3 μ 3 ν 1 μ 3 8 10 ν 1 μ 4 ( ξ + ξ 0 ) 1 + b exp 3 μ 3 ν 1 μ 3 8 10 ν 1 μ 4 ( ξ + ξ 0 ) 2 / 3 μ 3 4 μ 4 .
for the case b > 0 and a < 0 .
Next, we can set N = 2 , M = 5 . The equation that has to be solved is
ν 1 d 2 Q d ξ 2 + ν 2 d Q d ξ + μ 5 Q 5 + μ 4 Q 4 + μ 3 Q 3 + μ 2 Q 2 + μ 1 Q = 0 .
The SEsM leads to the conclusion that this equation has solutions of the kind Q ( ξ ) = i = 0 L a i [ ϕ ( ξ ) ] i , d ϕ d ξ = b ϕ ( 1 + 2 L ) + a ϕ . For the case L = 1 and for a > 0 , b < 0 , the solution is [478]
Q ( ξ ) = a 1 Q 1 Q 2 + a 0 ,
where
Q 1 = 4 ν 1 b a 0 2 + ν 2 a 1 2 + 4 ν 1 a 0 b ( ν 1 b a 0 3 + 2 a 0 ν 2 a 1 2 ) + ν 2 a 1 ( ν 2 a 1 3 + 4 ν 1 b a 0 3 ) 2 ν 1 a 1 2 × exp 2 4 ν 1 b a 0 2 + ν 2 a 1 2 + 4 ν 1 a 0 b ( ν 1 b a 0 3 + 2 a 0 ν 2 a 1 2 ) + ν 2 a 1 ( ν 2 a 1 3 + 4 ν 1 b a 0 3 ) 2 ν 1 a 1 2 ( ξ + ξ 0 ) .
Q 2 = 1
b exp 2 4 ν 1 b a 0 2 + ν 2 a 1 2 + 4 ν 1 a 0 b ( ν 1 b a 0 3 + 2 a 0 ν 2 a 1 2 ) + ν 2 a 1 ( ν 2 a 1 3 + 4 ν 1 b a 0 3 ) 2 ν 1 a 1 2 ( ξ + ξ 0 ) .
The solutions above are for the use of the equation of Bernoulli as the simple equation in the SEsM. One can also use other equations as simple equations. One example is the equation of Riccati. Numerous solutions of the discussed equations obtained based on the application of the SEsM based on the Riccati equation as the simple equation can be seen in [478].
The SEsM has been applied to obtain exact solutions of more complicated classes of equations. One such class is connected to modeling PDEs from ecology and population dynamics [480]
p = 1 N 1 A p ( Q ) p Q t p + r = 2 N 2 B r ( Q ) Q t r + s = 1 N 3 C s ( Q ) s Q x s + u = 2 N 4 D u ( Q ) Q x u + F ( Q ) = 0 .
A p ( Q ) , B r ( Q ) , C s ( Q ) , D u ( Q ) , and F ( Q ) are polynomials of Q. Specific cases of (160) have been obtained as model equations in the theory of the migration of populations [476,477]. We search for traveling wave solutions Q ( x , t ) = Q ( ξ ) = Q ( x v t ) . This transforms (160) into
p = 0 n 1 A ¯ p ( Q ) d p Q d ξ p + r = 2 n 2 B ¯ r ( Q ) d Q d ξ r + F ( Q ) = 0 .
In (161), n 1 = max ( N 1 , N 3 ) , n 2 = max ( N 2 , N 4 ) and A ¯ p ( Q ) = A p ( Q ) ( v ) p + C p ( Q ) , B ¯ r ( Q ) = B r ( Q ) ( v ) r + D r ( Q ) . A ¯ p , B ¯ r , and F are assumed as follows: A ¯ p ( Q ) = i = 0 M 1 p q i Q i , B ¯ r ( Q ) = i = 0 M 2 r q * Q i , F ( Q ) = i = 0 M q ¯ Q i . Here, q i , q i * , and q ¯ i are coefficients. We discuss obtaining the balance equations for several specific cases below. The corresponding exact solutions can be seen in [480].
We discuss the case of the application of the SEsM for obtaining the exact solution of (161) when the simple equation is the equation of Bernoulli and the solution is searched in the form
Q ( ξ ) = i = 0 L a i [ ϕ ( ξ ) ] i .
ϕ ( ξ ) is a solution of the equation of Bernoulli, which is used here as the simple equation. Thus, we fix the form of the composite function for the solution and we fix the simple equation. Next, we have to go to Step 3 of the SEsM to obtain the balance equations. In order to obtain the balance, we have to estimate the maximum powers of the term of the solved equation for the case when the equation of Bernoulli is used as the simple equation. We name the maximum powers BL 1 p , BL 2 r , and BL 3 . The estimations are as follows:
Term A ¯ p ( Q ) d p Q d ξ p BL 1 p = L M 1 p + L + p ( k 1 ) , p = 0 , 1 , , n 1
L M 1 p is the maximum power of ϕ ( ξ ) arising from the polynomial A ¯ p ( Q ) .
Term B ¯ r ( Q ) d Q d ξ r BL 2 r = L M 2 r + r ( L + k 1 ) , r = 2 , , n 1
L M 2 r is the maximum power of ϕ ( ξ ) arising from the polynomial B ¯ r ( Q ) .
Term F ( Q ) BL 3 = L M
We denote as max( BL 1 p ) the largest value among all BL 1 p and with max( BL 2 r ), the largest value among all BL 2 r . The theoretically possible balance equations are
max ( BL 1 p ) = max ( BL 2 r ) ; max ( BL 1 p ) = L M ; max ( BL 2 r ) = L M .
Often, one of these three terms is smaller than the other two.
Several examples of the determination of the balance equations are below. For the equation
Q t = D 2 Q x 2 + d 2 Ψ d Q 2 Q x 2 + d Ψ d Q 2 Q x 2 ,
where D is the constant coefficient of diffusion, Ψ ( Q ) is a polynomial of Q as follows Ψ ( Q ) = m 1 = 0 M 1 a m 1 Q m 1 . The largest powers of ϕ ( ξ ) are
BL 11 = L + k 1 ; BL 12 = L M 1 + 2 ( k 1 ) ; BL 22 = L M 1 + 2 ( k 1 ) .
Equation (164) has to be a nonlinear PDE. Then, M 1 > 2 . This leads to BL 12 > BL 11 and BL 22 > BL 11 . The balance equation becomes
BL 12 = BL 22 .
The next example is the balance equation for the reaction–diffusion equation with a diffusion coefficient that depends on population density
Q t = d D d Q Q x 2 + D ( Q ) 2 Q x 2 + F ( Q ) = 0 ,
where we assume D ( Q ) = m 1 = 0 M 1 a m 1 Q m 1 ; F ( Q ) = m = 0 M b m Q m . For the case of the use of the Bernoulli equation as the simple equation, the largest powers of ϕ ( ξ ) arising from the four terms of (166) are BL 11 = L + k 1 ; BL 12 = L M 1 + L + 2 ( k 1 ) , BL 22 = L ( M 1 1 ) + 2 ( L + k 1 ) ; BL 3 = L M . BL 11 is smaller than BL 12 = BL 22 . Two possibilities remain. If M is such that BL 3 is smaller than BL 12 , then the balance equation is BL 12 = BL 22 , and no requirements are imposed on L and M 1 . The second possibility is to balance BL 12 and BL 3 . Thus, BL 12 = BL 3 L ( M 1 + 1 ) + 2 ( k 1 ) = L M .
The balance equations for the other differential equations of the studied class, as well as the corresponding exact solutions, can be seen in [480].
A more complicated example is connected to the extended Camassa–Holm equation [484]. The Camassa–Holm equation is a model equation for the x component u of the fluid velocity at a certain depth z 0 below the fluid surface. The Camassa–Holm equation was introduced as a completely integrable bi-Hamiltonian dispersive shallow water equation [506]. The Camassa–Holm equation is [484,507]
U T + 2 κ U X + 3 ϵ U U X ϵ δ 2 3 U X 2 T = ϵ 2 δ 2 2 U X 2 U X 2 + U 3 U X 3 .
ϵ and δ are the amplitude parameter and the shallowness parameter and U is a quantity connected to the height of the wave η over the water surface.
η 5 3 U + 1 4 5 3 ϵ U 2 1 5 ϵ δ 2 2 U X 2 .
Below, we present an exact solution of the generalized Camassa–Holm equation [484]
U T + p 1 U X p 2 3 U X 2 T p 3 3 X 3 f ( U ) + p 4 X 1 2 g ( U ) + p 5 2 d 2 f d U 2 U X 2 = 0 .
p 1 ; p 2 ; p 3 ; p 4 ; p 5 are parameters. Specific cases of (169) are as below. For f ( U ) = 1 2 U 2 ; g ( U ) = C + 3 U 2 (C is a constant) and p 3 = p 4 p 5 , one obtains
U T + p 1 U X p 2 3 U X 2 T + 3 p 4 U U X 2 p 3 U X 2 U X 2 p 3 U 3 U X 3 = 0 .
Equation (169) contains the Camassa–Holm equation as a particular case when p 2 = 2 κ ; p 4 = ϵ ; p 2 = ϵ δ 2 ; p 3 = ϵ 2 δ 2 . In general, in Equation (169), f ( U ) and g ( U ) can be arbitrary polynomials of U.
We skip Step 1 of the SEsM (no transformation of the nonlinearity). We search for traveling wave solutions and introduce the coordinate ξ = X v T , where v is the velocity of the traveling wave. In Step 2 of the SEsM, we fix the composite function for the solution to be a power series of the solution of the simple equation. In Step 3 of the SEsM, we have to determine the balance equations for the corresponding simple equations. For the case when the simple equation is the equation of Bernoulli, the balance equation is
ν 1 = 2 ( k 1 ) J I .
where k > 1 and J > I > 1 . For I > 1 , another balance equation is possible. It is
I ν 1 + 3 ( k 1 ) = I ν 1 + 3 ( k 1 ) ,
when ( J I ) ν 1 < 2 ( k 1 ) . This balance is possible because of the existence of two terms in (169), which lead to the maximum power I ν 1 + 3 ( k 1 ) of the terms of the resulting polynomial, and these powers are maximum powers for the entire polynomial when ( J I ) ν 1 < 2 ( k 1 ) .
A third balance equation is possible in the case of the use of the Bernoulli equation as the simple equation. This balance is when I = 1 . The corresponding balance equation is
ν 1 = 2 ( k 1 ) J 1 .
Another possibility for the simple equation is the equation of Riccati (136). The solution of Equation (169) is, again, of the kind (132). The possible balance equations in this case are
ν 1 = I ; J < ν 1 + 2 ; J < I + 2 ,
ν 1 + 2 = J ; I < ν 1 ; I + 2 < J ,
J = I + 2 ; ν 1 < I ; ν 1 + 2 < J .
The possible balance equations for the case when the extended tanh-equation is used as the simplest equation are the same as for the case of the Riccati equation.
We return to the case of the use of the SEsM with the equation of Bernoulli as the simple equation. The values of parameters in the balance equation, Equation (171), are set to J = 3 , I = 2 , k = 2 . Then, from (171), ν 1 = 2 . The equation we solve is
U T + p 1 U X p 2 3 U X 2 T p 3 3 X 3 ( q 0 + q 1 U + q 2 U 2 ) + p 4 X 1 2 ( r 0 + r 1 U + r 2 U 2 + r 3 U 3 ) + p 5 q 2 U X 2 = 0 .
The solution of Equation (170) is searched in the form
U ( ξ ) = θ 0 + θ 1 V ( ξ ) + θ 2 V ( ξ ) 2 .
V ( ξ ) is the solution of the equation of Bernoulli for k = 2 . Substituting (178) in (170), we obtain a system of seven nonlinear relationships between the parameters of the solved equation and the parameters of the solution (178). Step 4 of the SEsM is to solve this system of algebraic equations. One nontrivial solution is
θ 2 = 8 q 2 a c ( 5 p 3 p 4 p 5 ) p 4 r 3 ; θ 1 = 8 q 2 c 2 ( 5 p 3 p 4 p 5 ) p 4 r 3 ; θ 0 = Θ 0 6 q 2 p 4 r 3 ( 4 p 3 p 4 p 5 ) ; r 1 = R 1 12 p 2 p 4 r 3 ( 4 p 3 p 4 p 5 ) ; Θ 0 = 36 q 2 2 a 2 p 4 p 5 p 3 + 4 q 2 2 a 2 p 4 2 p 5 2 3 p 2 v p 4 r 3 + 3 p 3 q 1 p 4 r 3 + 80 q 2 2 a 2 p 3 2 10 q 2 p 3 r 2 p 4 + 2 q 2 p 4 2 p 5 r 2 ; R 1 = ( 1520 q 2 4 a 4 p 4 2 p 5 2 p 3 2 256 q 2 4 a 4 p 4 3 p 5 3 p 3 4 q 2 2 p 4 4 p 5 2 r 2 2 + 32 q 2 2 p 3 r 2 2 p 4 3 p 5 60 q 2 2 p 3 2 r 2 2 p 4 2 + 3840 q 2 4 a 4 p 3 4 12 p 3 2 q 1 p 4 2 r 3 q 2 r 2 + 12 p 2 v p 4 2 r 3 q 2 p 3 r 2 18 p 2 v p 4 2 r 3 2 p 3 q 1 + 9 p 3 2 q 1 2 p 4 2 r 3 2 + 16 q 2 4 a 4 p 4 4 p 5 4 + 9 p 2 2 v 2 p 4 2 r 3 2 3968 q 2 4 a 4 p 4 p 5 p 3 3 24 v q 2 2 r 3 p 4 3 p 5 2 + 192 v q 2 2 r 3 p 4 2 p 3 p 5 384 v q 2 2 r 3 p 4 p 3 2 192 p 1 q 2 2 r 3 p 4 2 p 3 p 5 + 384 p 1 q 2 2 r 3 p 4 p 3 2 + 24 p 1 q 2 2 r 3 p 4 3 p 5 2 ) .
This leads to the following solution of Equation (177):
U ( ξ ) = Θ 0 6 q 2 p 4 r 3 ( 4 p 3 p 4 p 5 ) + 8 q 2 c 2 ( 5 p 3 p 4 p 5 ) p 4 r 3 a exp [ a ( ξ + ξ 0 ) ] 1 c exp [ a ( ξ + ξ 0 ) ] + 8 q 2 a c ( 5 p 3 p 4 p 5 ) p 4 r 3 a exp [ a ( ξ + ξ 0 ) ] 1 c exp [ a ( ξ + ξ 0 ) ] 2 ,
for c < 0 and a > 0 . For c > 0 and a < 0 , the solution is
U ( ξ ) = Θ 0 6 q 2 p 4 r 3 ( 4 p 3 p 4 p 5 ) 8 q 2 c 2 ( 5 p 3 p 4 p 5 ) p 4 r 3 a exp [ a ( ξ + ξ 0 ) ] 1 + c exp [ a ( ξ + ξ 0 ) ] + 8 q 2 a c ( 5 p 3 p 4 p 5 ) p 4 r 3 a exp [ a ( ξ + ξ 0 ) ] 1 + c exp [ a ( ξ + ξ 0 ) ] 2 .
Another possibility is to use the equation of Riccati (136) as the simple equation. We use (174) for the balance equation and set the parameters to J = 1 ; ν 1 = 2 ; I = 2 . The equation that is to be solved is
U T + p 1 U X p 2 3 U X 2 T p 3 3 X 3 ( q 0 + q 1 U + q 2 U 2 ) + p 4 X 1 2 ( r 0 + r 1 U ) + p 5 q 2 U X 2 = 0 .
In Step 4 of the SEsM, we obtain a nonlinear algebraic system of eight relationships between the parameters of the solved equation and the parameters of the solution. A nontrivial solution of this system is
θ 2 = 15 a 2 ( p 4 r 1 + 2 v 2 p 1 ) q 2 p 4 p 4 ( 8 c 2 d a + 16 a 2 d 2 + c 4 ) , θ 1 = 15 a c ( p 4 r 1 + 2 v 2 p 1 ) q 2 p 4 p 4 ( 8 c 2 d a + 16 a 2 d 2 + c 4 ) , θ 0 = Θ 0 q 2 p 4 p 5 ( 8 c 2 d a + 16 a 2 d 2 + c 4 ) , Θ 0 = 1 4 ( 40 a d p 4 r 1 80 a d v + 80 a d p 1 + 5 c 2 p 4 r 1 10 c 2 v + 10 c 2 p 1 + 16 p 4 p 5 q 1 c 2 d a 32 p 4 p 5 q 1 a 2 d 2 2 p 4 p 5 q 1 c 4 80 p 2 v c 2 d a + 160 p 2 v a 2 d 2 + 10 p 2 v c 4 ) .
The corresponding solution of Equation (182) based on the solution (137) of the Riccati equation is
U ( ξ ) = Θ 0 q 2 p 4 p 5 ( 8 c 2 d a + 16 a 2 d 2 + c 4 ) + 15 a c ( p 4 r 1 + 2 v 2 p 1 ) q 2 p 4 p 4 ( 8 c 2 d a + 16 a 2 d 2 + c 4 ) c 2 a + θ 2 a tanh θ ( ξ + ξ 0 ) 2 15 a 2 ( p 4 r 1 + 2 v 2 p 1 ) q 2 p 4 p 4 ( 8 c 2 d a + 16 a 2 d 2 + c 4 ) c 2 a + θ 2 a tanh θ ( ξ + ξ 0 ) 2 2 .
In (184), θ 2 = c 2 4 a d > 0 .
Finally, we can use the extended tanh-equation, Equation (138), as the simple equation. From the corresponding balance equation, we choose the parameters I = J = 2 and ν 1 = 2 . The equation to be solved is
U T + p 1 U X p 2 3 U X 2 T p 3 3 X 3 ( q 0 + q 1 U + q 2 U 2 ) + p 4 X 1 2 ( r 0 + r 1 U + r 2 U 2 ) + p 5 q 2 U X 2 = 0 .
In Step 4 of the SEsM, one obtains a system of nonlinear algebraic equations that contains eight relationships. A nontrivial solution of this system is
θ 2 = 60 a ¯ 4 ( 2 v p 5 q 2 + 2 p 1 p 5 q 2 + r 1 p 4 p 5 q 2 + 5 r 2 p 2 v r 2 p 4 p 5 q 1 ) p 4 ( 64 p 5 2 q 2 2 d ¯ 4 a ¯ 2 25 r 2 2 ) , θ 1 = 0 , θ 0 = Θ 0 p 4 ( 25 r 2 2 + 64 p 5 2 q 2 2 d ¯ 4 a ¯ 4 ) , Θ 0 = 1 2 ( 160 p 5 q 2 a ¯ 2 d ¯ 2 v 160 p 5 q 2 a ¯ 2 d ¯ 2 p 1 80 p 5 q 2 a ¯ 2 d ¯ 2 r 1 p 4 400 a ¯ 2 d ¯ 2 r 2 p 2 v + 80 p 5 a ¯ 2 d ¯ 2 r 2 p 4 q 1 + 320 p 2 v p 5 q 2 d ¯ 4 a ¯ 4 50 r 2 v + 50 r 2 p 1 + 25 r 2 r 1 p 4 64 p 4 p 5 2 q 1 q 2 d ¯ 4 a ¯ 4 ) p 3 = 1 5 p 4 p 5 .
The corresponding solution of the solved equation, Equation (185), is
U ( ξ ) = Θ 0 p 4 ( 25 r 2 2 + 64 p 5 2 q 2 2 d ¯ 4 a ¯ 4 ) + 60 a ¯ 2 d ¯ 2 ( 2 v p 5 q 2 + 2 p 1 p 5 q 2 + r 1 p 4 p 5 q 2 + 5 r 2 p 2 v r 2 p 4 p 5 q 1 ) p 4 ( 64 p 5 2 q 2 2 d ¯ 4 a ¯ 2 25 r 2 2 ) tanh 2 [ a ¯ d ¯ ( ξ + ξ 0 ) ] .
The next example we consider is connected to the generalized Swift–Hohenberg equation. The 1+1-dimensional version is discussed:
u t = r u 1 + 2 2 u + N ( u ) .
Here, r is a parameter and N ( u ) is a polynomial nonlinearity N ( u ) s = 0 s * c s * u s , where c s * , s = 0 , 1 , 2 , are parameters. Traveling wave solutions are searched. We set c 0 = c 0 * ; c 1 = c 1 * + r 1 ; …, c s = c s * , and the solved equation becomes
u t + 2 2 u x 2 + u 4 x 4 = s = 0 s * c s u s .
The balance equation that is obtained in Step 3 of the SEsM is
η = 4 ( β 1 ) s * 1 .
We consider the specific case s * = 5 . The corresponding Swift–Hohenberg equation becomes
u t + 2 2 u x 2 + 4 u x 4 = c 0 + c 1 u + c 2 u 2 + c 3 u 3 + c 4 u 4 + c 5 u 5 .
The balance equation is
η = β 1 .
For β = 2 , η = 1 . β = 2 means that we use the equation of Riccati (136) or the extended tanh-Equation (138) as the simple equation. For the case of the use of the extended tanh-equation as the simple equation, in Step 4 of the SEsM, one obtains the system of nonlinear algebraic relationships.
24 a 4 c 5 θ 1 4 = 0 , c 4 + 5 c 5 θ 0 = 0 , 40 a 3 c + 4 a 2 c 3 θ 1 2 10 c 5 θ 0 2 θ 1 2 4 c 4 θ 0 θ 1 2 = 0 , v a c 2 θ 1 3 c 3 θ 0 θ 1 10 c 5 θ 0 3 θ 1 6 c 4 θ 0 2 θ 1 = 0 , 2 c 2 θ 0 + 4 c a + 16 a 2 c 2 5 c 5 θ 0 4 3 c 3 θ 0 2 4 c 4 θ 0 3 c 1 = 0 , c 0 c 1 θ 0 v θ 1 c c 2 θ 0 2 c 4 θ 0 4 c 5 θ 0 5 c 3 θ 0 3 = 0 .
A nontrivial solution of this system is
c 5 = 24 a 4 θ 1 4 , c 4 = 120 a 4 θ 0 θ 1 4 , c 3 = 4 a 2 10 θ 1 2 a c + θ 1 2 + 60 a 2 θ 0 2 θ 1 4 , c 2 = a v θ 1 3 + 120 a 2 θ 0 θ 1 2 c + 12 a θ 0 θ 1 2 + 240 a 3 θ 0 3 θ 1 4 , c 1 = 2 a θ 0 v θ 1 3 + 60 a 2 θ 0 2 θ 1 2 c + 6 a θ 0 2 θ 1 2 + 60 a 3 θ 0 4 + 2 θ 1 4 c + 8 θ 1 4 a c 2 θ 1 4 , c 0 = 24 a 4 θ 0 5 + a θ 0 2 v θ 1 3 + 40 a 3 θ 0 3 θ 1 2 c + 4 a 2 θ 0 3 θ 1 2 + 4 a θ 0 θ 1 4 c + 16 a 2 θ 0 θ 1 4 c 2 + v θ 1 5 c θ 1 4 .
The corresponding solution of (191) is
u ( ξ ) = θ 0 + θ 1 c ¯ a ¯ tanh [ a ¯ c ¯ ( ξ + ξ 0 ) ] ,
where a = a ¯ 2 and c = c ¯ 2 .
For the case β = 3 , one obtains from (192) η = 2 . The simple equation is
d V d ξ = γ 0 + γ 1 V + γ 2 V 2 + γ 3 V 3 ,
which for γ 0 = γ 2 = 0 and γ 1 = a , as well as γ 3 = b , is reduced to the equation of Bernoulli (134). A nontrivial solution of the corresponding system of nonlinear algebraic relationships obtained in Step 4 of the SEsM is
θ 1 = 0 , c 5 = 384 b 4 θ 2 4 , c 4 = 960 b 3 ( 2 b θ 0 + θ 2 a ) θ 2 4 , c 3 = 16 b 2 ( 240 b 2 θ 0 2 240 b θ 0 θ 2 a + 50 θ 2 2 a 2 + θ 2 2 ) θ 2 4 , c 2 = 2 b θ 2 4 ( v θ 2 3 + 120 θ 2 3 a 3 + 12 θ 2 3 a 1920 b 3 θ 0 3 + 2880 b 2 θ 0 2 θ 2 a 1200 b θ 0 θ 2 2 a 2 24 b θ 0 θ 2 2 ) , c 1 = 2 θ 2 4 ( 960 b 4 θ 0 4 1920 b 3 θ 0 3 θ 2 a + 1200 b 2 θ 0 2 θ 2 2 a 2 + 24 b 2 θ 0 2 θ 2 2 v a θ 2 4 + 8 a 4 θ 2 4 + 2 b θ 0 v θ 2 3 240 b θ 0 θ 2 3 a 3 24 b θ 0 θ 2 3 a + 4 a 2 θ 2 4 ) , c 0 = 2 θ 0 θ 2 4 ( 192 b 4 θ 0 4 480 b 3 θ 0 3 θ 2 a + 400 b 2 θ 0 2 θ 2 2 a 2 + 8 b 2 θ 0 2 θ 2 2 v a θ 2 4 + 8 a 4 θ 2 4 + b θ 0 v θ 2 3 120 b θ 0 θ 2 3 a 3 12 b θ 0 θ 2 3 a + 4 a 2 θ 2 4 ) .
The solution of the Swift–Hohenberg Equation (191) is
u ( ξ ) = θ 0 + θ 2 a exp [ 2 a ( ξ + ξ 0 ) ] 1 b exp [ 2 a ( ξ + ξ 0 ) ] ,
for the case b < 0 , a > 0 and
u ( ξ ) = θ 0 θ 2 a exp [ 2 a ( ξ + ξ 0 ) ] 1 + b exp [ 2 a ( ξ + ξ 0 ) ] ,
for the case b > 0 , a < 0 . θ 0 and θ 2 are free parameters.
The next example is connected to the generalized Rayleigh equation
2 u t 2 2 u x 2 + F u t + G ( u ) = 0 .
F and G are the functions of the corresponding terms. For u = u ( t ) , G ( u ) = u , one obtains the Rayleigh equation. For F ( d u d t ) = ϵ d u d t d u d t 2 , one obtains the specific case of the Rayleigh equation
d 2 y d t 2 ϵ d u d t d u d t 2 + u = 0 .
Equation (201) occurs in the theory of sound.
We consider the case where F and G are polynomials
F u t = j = 0 m a j u t j ; G ( u ) = s = 0 s * c s u s .
Traveling wave solutions u ( x , t ) = u ( ξ ) = u ( x v t ) are searched for the solved equation
2 u t 2 + j = 0 m a j u t j 2 u x 2 = s = 0 s * c s u s .
The balance equation that occurs in Step 3 of the SEsM is
η = m ( β 1 ) s * m .
For the case m = 2 , we solve the equation
2 u t 2 + j = 0 2 a j u t j 2 u x 2 = s = 0 s * c s u s ,
and the balance equation is
η = 2 ( β 1 ) s * 2 .
Here, we set β = 2 . This fixes the simple equation to be the equation of Riccati or the extended tanh-equation. In addition, we set s * = 3 and the solved equation becomes
2 u t 2 + j = 0 2 a j u t j 2 u x 2 = c 0 + c 1 u + c 2 u 2 + c 3 u 3 .
For the case of the extended tanh-equation as the simple equation, one obtains a system of nonlinear relationships in Step 4 of the SEsM. A nontrivial solution of this system is
θ 2 = 4 a 2 v 2 a 2 c 3 , θ 1 = 0 , θ 0 = 1 6 a 2 c 3 v 2 ( 3 v 2 c 3 3 c 3 2 c 2 a 2 v 2 + 2 ( 18 v 2 c 3 2 12 c 1 a 2 2 v 4 c 3 + 4 c 2 2 a 2 2 v 4 9 v 4 c 3 2 9 c 3 2 ) 1 / 2 ) , a 1 = 0 , c = 1 8 v 4 a a 2 2 ( 18 v 2 c 3 2 12 c 1 a 2 2 v 4 c 3 + 4 c 2 2 a 2 2 v 4 9 v 4 c 3 2 9 c 3 2 ) 1 / 2 , c 0 = 1 108 c 3 2 a 2 3 v 6 [ ( 18 v 2 c 3 2 12 c 1 a 2 2 v 4 c 3 + 4 c 2 2 a 2 2 v 4 9 v 4 c 3 2 9 c 3 2 ) 1 / 2 ( 9 v 4 c 3 2 18 v 2 c 3 2 + 9 c 3 2 + 12 c 1 a 2 2 v 4 c 3 4 c 2 2 a 2 2 v 4 ) 8 c 2 3 a 2 3 v 6 108 a 0 c 3 2 a 2 3 v 6 108 v 6 c 3 c 2 2 a 2 2 + 324 v 6 c 3 2 c 1 a 2 2 + 108 c 3 c 2 2 a 2 2 v 4 324 c 1 a 2 2 v 4 c 3 2 + 270 v 6 c 3 3 810 v 4 c 3 3 + 810 v 2 c 3 3 270 c 3 3 + 36 c 1 a 2 3 v 6 c 3 c 2 ] .
The solution of (207) is
u ( ξ ) = 1 6 a 2 c 3 v 2 ( 3 v 2 c 3 3 c 3 2 c 2 a 2 v 2 + 2 ( 18 v 2 c 3 2 12 c 1 a 2 2 v 4 c 3 + 4 c 2 2 a 2 2 v 4 9 v 4 c 3 2 9 c 3 2 ) 1 / 2 ) + θ 2 { [ 1 8 v 4 a 2 a 2 2 ( 18 v 2 c 3 2 12 c 1 a 2 2 v 4 c 3 + 4 c 2 2 a 2 2 v 4 9 v 4 c 3 2 9 c 3 2 ) 1 / 2 ] 1 / 2 tanh [ ( 1 8 v 4 a 2 2 ( 18 v 2 c 3 2 12 c 1 a 2 2 v 4 c 3 + 4 c 2 2 a 2 2 v 4 9 v 4 c 3 2 9 c 3 2 ) 1 / 2 ) 1 / 2 ( ξ + ξ 0 ) ] } .
Next, we set β = 3 . The simple equation, Equation (133), is again reduced to the equation of Bernoulli (134) with k = 3 . Moreover, we set s * = 3 . Then, from the balance equation, η = 4 . The solved equation is (207), and in Step 4 of the SEsM, one arrives at the system of nonlinear algebraic relationships. A nontrivial solution of this system is
θ 4 = 16 a 2 v 2 b 2 c 3 , θ 3 = 0 , θ 2 = 16 a 2 v 2 b a c 3 , θ 1 = 0 , θ 0 = 1 6 a 2 v 2 c 3 ( 8 a 2 2 v 4 a 2 2 c 2 a 2 v 2 + 3 v 2 c 3 3 c 3 ) , a 1 = 0 , c 1 = 1 12 c 3 a 2 2 v 4 ( 4 c 2 2 a 2 2 v 4 + 64 a 2 4 v 8 a 4 + 9 v 4 c 3 2 18 v 2 c 3 2 + 9 c 3 2 ) , c 0 = 1 108 c 3 2 a 2 3 v 6 ( 27 c 3 3 + 27 v 6 c 3 3 81 v 4 c 3 3 + 81 v 2 c 3 3 + 512 a 2 6 v 12 a 6 + 4 c 2 3 a 2 3 v 6 108 c 0 c 3 2 a 2 3 v 6 192 c 2 a 2 5 v 10 a 4 27 c 2 a 2 v 2 c 3 2 + 54 c 2 a 2 v 4 c 3 2 27 c 2 a 2 v 6 c 3 2 ) .
The corresponding solution of the solved equation is
u ( ξ ) = 1 6 a 2 v 2 c 3 ( 8 a 2 2 v 4 a 2 2 c 2 a 2 v 2 + 3 v 2 c 3 3 c 3 ) + 16 a 2 v 2 b a 2 c 3 exp [ 2 a ( ξ + ξ 0 ) ] 1 b exp [ 2 a ( ξ + ξ 0 ) ] + 16 a 2 v 2 b 2 c 3 a exp [ 2 a ( ξ + ξ 0 ) ] 1 b exp [ 2 a ( ξ + ξ 0 ) ] 2 ,
for the case b < 0 , a > 0 and
u ( ξ ) = 1 6 a 2 v 2 c 3 ( 8 a 2 2 v 4 a 2 2 c 2 a 2 v 2 + 3 v 2 c 3 3 c 3 ) 16 a 2 v 2 b a 2 c 3 exp [ 2 a ( ξ + ξ 0 ) ] 1 + b exp [ 2 a ( ξ + ξ 0 ) ] + 16 a 2 v 2 b 2 c 3 a exp [ 2 a ( ξ + ξ 0 ) ] 1 + b exp [ 2 a ( ξ + ξ 0 ) ] 2 ,
for the case b > 0 , a < 0 .
The last example illustrates the application of the SEsM to the case of the use of two simple equations. The set of equations that is discussed contains equations of the nonlinear Schrödinger kind [488]:
i q t + a q x x + q k = A k 1 B b k q 2 k = 0 .
In order to obtain a solution of Equation (213) one needs the two simplest equations. We skip Step 1 of the SEsM (no transformation of the nonlinearity, as the nonlinearity is a polynomial one). In Step 2 of the SEsM, we can choose a composite function for the solution. The choice is q ( x , t ) = g ( ξ ) h ( x , t ) . Here, g ( ξ ) is a real function ( ξ = α x + β t ). h ( x , t ) is a complex function. The two simple equations are for g ( ξ ) and h ( x , t ) . Below, h * ( x , t ) denotes the complex conjugate function of h ( x , t ) . The substitution of the form of q ( x , t ) in Equation (213) leads to (we denote g ξ as g )
i β g h + i g h t + α 2 a g h + 2 α a g h x + a g h x x + ( g h ) k = A k 1 B g 2 k ( h h * ) k = 0 .
Here, g ξ g denotes g ξ . The simple equation for h ( x , t ) is h ζ = i h , ζ = κ x + ω t + θ . The solution is h ( ζ ) = exp ( i ζ ) = exp [ i ( κ x + ω t + θ ) ] . This choice of h leads to the equation for g ( ξ )
α 2 a g + ( 2 α κ a + β ) i g ( ω + κ 2 a ) g + g k = A k 1 B b k g 2 k = 0 .
g ( ξ ) must be a real function and because of this, β = 2 α κ a . We substitute this in (215) and multiply the result by g . What is obtained is integrated with respect to ξ . The result is
α 2 a g 2 ( ω + κ 2 a ) g 2 c + k = A k 1 B b k k + 1 g 2 k + 2 = 0 .
c is a constant of integration. We introduce a new variable u = g σ . Here, σ is a parameter. The equation for u ( ξ ) becomes
u 2 = σ 2 ( ω + κ 2 a ) α 2 a u 2 + σ 2 c α 2 a u 2 ( σ 1 ) / σ k = A k 1 B σ 2 b k α 2 a ( k + 1 ) u 2 ( k + σ ) / σ .
Equation (217) has various specific cases. Three of them are:
  • σ = 1 , k = 0 , , B . Then, (217) becomes the following equation:
    u 2 = k = 0 B c k u 2 ( k + 1 ) , c 0 = ω + κ 2 a b 0 α 2 a ; c k = b k α 2 a ( k + 1 ) , k = 1 , , B .
  • σ = 2 , k = 2 , , B . Then, (217) becomes
    u 2 = k = 2 B c k u k + 2 , c 2 = b 2 α 2 a ; c 1 = 4 c α 2 a ; c 0 = 4 ω + 4 κ 2 a 4 b 0 α 2 a ; c k = 4 b k α 2 a ( k + 1 ) , k = 1 , , B .
  • c = 0 , σ : arbitrary positive integer that is different from 1, B = D σ , k = n σ , n = 1 , , D . Then, (217) becomes
    u 2 = n = 1 D c n u 2 ( n + 1 ) , c 1 = σ 2 b σ α 2 a ( 1 σ ) , c 0 = σ 2 ( ω + κ 2 a b 0 ) α 2 a , c n = σ 2 b n σ α 2 a ( σ n + 1 ) , n = 1 , , D .
The solutions of (218)–(220) can be expressed by the special function V a ( ξ ; q , l , m ) from [486]. Remember that this special function is the solution of d q u d ξ q l = j = 0 m a j u j . Here, q, l, and m are positive integers and a = ( a 0 , , a m ) . The solutions of (213) are:
  • σ = 1 , k = 0 , , B .
    q ( x , t ) = V a 1 ( ξ ; 1 , 2 , 2 B + 2 ) exp [ i ( k x + ω t + θ ) ] .
    Here, a 1 = ( a 0 , , a 2 B + 2 ) and a 0 = 0 , a 2 i = c i 1 , a 2 i 1 = 0 , i = 1 , , B + 1 .
  • σ = 2 , k = 2 , , B .
    q ( x , t ) = [ V a 2 ( ξ ; 1 , 2 , B + 2 ) ] 1 / 2 exp [ i ( k x + ω t + θ ) ] .
    Here, a 2 = ( a 0 , , a B + 2 ) and a i = c i 2 , i = 0 , , B + 2 .
  • c = 0 , σ : an arbitrary positive integer that is different from 1, B = D σ , k = n σ , n = 1 , , D .
    q ( x , t ) = [ V a 3 ( ξ ; 1 , 2 , 2 D + 2 ) ] 1 / σ exp [ i ( k x + ω t + θ ) ] .
    Here, a 3 = ( a 0 , , a 2 D + 2 ) and a 2 i = c i 1 , a 2 i + 1 = 0 , i = 0 , D + 1 .
The above solutions can be expressed by the specific case of the V-functions for the selected values of the parameters. One example is for the solutions given by Equation (222). These solutions are:
  • Case B = 0 . Equation (213) is the solution, which is
    q ( x , t ) = { cosh [ ± ( b 2 / a ) 1 / 2 ( x 2 κ a t ) ] } 1 / 2 exp { i [ κ x + 1 4 ( b 2 4 κ 2 a + 4 b 0 ) t + θ ] } .
  • Case B = 1 . Equation (213) is
    i q t + a q x x + b 2 q q 4 + b 0 q + b 1 q q 2 = 0 .
    The solution of Equation (225) is
    q ( x , t ) = { [ ± ( b 1 / ( 2 a ) ) 1 / 2 ( x 2 κ a t ) ; 0 , 2 b 2 / b 1 ] } 1 / 2 exp { i [ κ x + ( b 0 κ 2 a ) t + θ ] } .
    Here, is the elliptic function of Weierstrass.
    Case B = 2 . Equation (213) is
    i q t + a q x x + b 2 q q 4 + b 0 q + b 1 q q 2 + b 2 q q 4 = 0 ,
    The solution is expressed by the Jacobi elliptic functions [488].
The list of examples of the application of the SEsM could continue. For more applications, see, for example, [460,463,467,468,469,470,484,485,487,488,508,509,510].

9. The SEsM and Other Methods for Obtaining Specific Exact Solutions of Nonlinear Nonintegrable Differential Equations

An interesting direction for future research is the relationship between the SEsM and other methods for obtaining exact solutions of nonlinear differential equations. We have shown above that by using the appropriate composite function of Step 2 of the SEsM and simple equations for exponential functions, the SEsM can be connected to the method of Hirota and the Inverse Scattering Transform method. These two methods are often used for obtaining exact solutions of integrable nonlinear differential equations. The situation with respect to the methods for obtaining particular solutions of nonintegrable nonlinear differential equations is as follows. It was shown that the SEsM contains as specific cases [511] the Jacobi Elliptic Function Expansion method [512], F-Expansion method [513,514,515], Modified Simple Equation method [516], Trial Function method [517,518], General Projective Riccati Equations method [519], and First Integral method [520]. Other methods that demonstrate specific cases of the SEsM are [495] the Homogeneous Balance method [521,522,523,524] and the Auxiliary Equation method [525]. The list continues [466] with the G’/G method [526,527], Exp-function method [528], tanh method [529], and the method of the Fourier series. We note that the SEsM can also deal with differential equations containing fractional derivatives [465]. The list can be continued and we present here the following conjecture.
Conjecture 1.
Any method for obtaining exact analytical solutions of nonintegrable nonlinear differential equations is a specific case of the SEsM.
The motivation for this conjecture is as follows. The SEsM is based on known solutions of simple equations and the solution of the more complicated solved equation is constructed as a composite function of these solutions. Any method for obtaining exact analytical solutions of nonlinear differential equations has to construct these solutions using known functions. These known functions are the solutions of the simple equations. The solution of the solved equation is a function of these known functions. This is the composite function used in the SEsM. Then, if one wants to find a method for obtaining the exact solution of a nonlinear differential equation that is not a specific case of the SEsM, one has to do one of the following things:
  • Use the forms of the known functions that are not solutions of any simple equations. This will be extremely difficult as one has to use functions that are not solutions of differential equations.
  • Construct the searched solution as a function of the known functions, which is not a composite function. This is also an extremely difficult task.

10. Concluding Remarks

In this review article, we discuss an effective method for obtaining exact analytical solutions of nonlinear differential equations, the Simple Equations Method (SEsM). Through a relatively large overview of the literature, we show that the research area connected to the search for exact solutions of nonlinear differential equations is of much interest because these equations are used in the mathematical models of numerous complex systems, and the mathematical aspect of the search of exact solutions is also quite interesting. The SEsM has a simple algorithm based on the appropriate transformation of the nonlinearity of the solved equation (if necessary), and then the solution is constructed as a composite function of solutions of simpler differential equations. The SEsM has four steps. In Step 1, one may use a transformation to convert the nonlinearity of the solved equation to a more treatable kind of nonlinearity. A treatable kind of nonlinearity is polynomial nonlinearity. Step 2 of the SEsM is connected with the use of a composite function in order to build the solution of the solved equation. We show that the composite function favors the occurrence of a certain class of simple equations that define a special function, the V-function. This function contains as specific cases numerous elementary and special functions and leads to the occurrence of systems of polynomials in the process of the application of the SEsM to the solved equation. The SEsM transforms the solved equation into a sum of functions multiplied by some coefficients, which are relationships between the parameters of the solved equation and the parameters of the solution. The setting of these coefficients to 0 leads to a system of nonlinear algebraic equations. In order to obtain a nontrivial solution of the solved equation, this system of nonlinear algebraic relationships must be balanced. This leads to additional relationships called balance equations. The balance equations are important as they may lead to a fixation of the form of the simple equation and the form of the composite function for the solution of the solved equation. Thus, the first three steps of the SEsM convert the solved equation to a system of nonlinear algebraic relationships. The solution of this system is the task of Step 4 of the SEsM. Each obtained nontrivial solution of the system of nonlinear algebraic relationships leads to a nontrivial solution of the solved equation.
We note that in general, the SEsM is constructed to be able to solve systems of nonlinear differential equations. Above, we presented the version of the method for the solution of one nonlinear differential equation. We demonstrated that this version of the SEsM is connected to well-known methods for obtaining multisoliton exact analytical solutions of integrable differential equations, the Inverse Scattering Transform method, and the method of Hirota. We showed how the SEsM led to the equation of Gelfand–Levitan–Marchenko for the case of the Korteweg-de Vries equation and to the relationships of Zakharov and Shabat for the case of the nonlinear Schrödinger equation. We discussed numerous transformations of nonpolynomial nonlinearities to polynomial nonlinearities in Step 1 of the SEsM. The need to use composite functions arose naturally in the SEsM as the solution had to be constructed as a composite function of the solution of the simple equations. This led to the occurrence of the Faa di Bruno formula for the derivatives of the composite functions in the process of the application of the SEsM.
The use of composite functions in the SEsM favors simple equations with polynomial nonlinearity. Thus, simple equations, such as the equations of Bernoulli and Riccati and the equations for the elliptic functions of Weierstrass and Jacobi, are often used in the SEsM. The connection between the choice of the simple equation and the form of the composite function for the solution of the solved equation was discussed with special attention in the text. Numerous examples of the application of the SEsM were presented. These applications illustrated the methodology and presented exact analytical solutions of many equations obtained over the years with the use of the SEsM in our research. Special attention was given to the application of the SEsM to differential equations arising in connection with the SIR model of epidemic spreading. The obtained solutions can be used for modeling COVID-19 waves. Other examples were connected to exact analytical solutions of the equations of Fisher and Huxley, as well as the solutions of the generalized Camassa–Holm, Swift–Hohenberg, and Rayleigh equations. Special attention was also devoted to obtaining exact solutions of nonlinear differential equations such as the nonlinear equation of Schrödinger. The reason is that we showed the application of the SEsM for the case of two simple equations in the process of obtaining exact analytical solutions of these equations.
As we have shown above, the methodology of the SEsM is effective. It is connected to well-known methods for obtaining exact multisoliton solutions of nonlinear differential equations and contains as specific cases many methods for obtaining specific exact analytical solutions of nonintegrable equations. Future research on the SEsM will be in two directions: (i) an extension of the area of the application of the methodology by treating additional classes of nonlinearities and equations; and (ii) the application of the methodology to systematic methods for obtaining exact solutions of nonlinear differential equations.

Funding

This research was partially supported by the project BG05 M2OP001-1.001-0008 “National Center for Mechatronics and Clean Technologies”, funded by the Operating Program “Science and Education for Intelligent Growth” of Republic of Bulgaria.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Latora, V.; Nicosia, V.; Russo, G. Complex Networks. Principles, Methods, and Applications; Cambridge University Press: Cambridge, UK, 2017; ISBN 978-1-107-10318-4. [Google Scholar]
  2. Chian, A.C.-L. Complex Systems Approach to Economic Dynamics; Springer: Berlin, Germany, 2007; ISBN 978-3-540-39752-6. [Google Scholar]
  3. Vitanov, N.K. Science Dynamics and Research Production. Indicators, Indexes, Statistical Laws and Mathematical Models; Springer: Cham, Switzerland, 2016; ISBN 978-3-319-41629-8. [Google Scholar]
  4. Treiber, M.; Kesting, A. Traffic Flow Dynamics: Data, Models, and Simulation; Springer: Berlin, Germany, 2013; ISBN 978-3-642-32460-4. [Google Scholar]
  5. May, R.M.; Levin, S.A.; Sugihara, G. Complex Systems: Ecology for Bankers. Nature 2008, 451, 893–895. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Ivanova, K.; Ausloos, M. Application of the Detrended Fluctuation Analysis (DFA) Method for Describing Cloud Breaking. Phys. A 1999, 274, 349–354. [Google Scholar] [CrossRef]
  7. Kutner, R.; Ausloos, M.; Grech, D.; Di Matteo, T.; Schinckus, C.; Stanley, H.E. Manifesto for a Post-Pandemic Modeling. Phys. A 2019, 516, 240–253. [Google Scholar] [CrossRef] [Green Version]
  8. Simon, J.H. The Economic Consequences of Immigration; The University of Michigan Press: Ann Arbor, MI, USA, 1999; ISBN 978-0472086160. [Google Scholar]
  9. Drazin, P.G. Nonlinear Systems; Cambridge University Press: Cambridge, UK, 1992; ISBN 0-521-40489-4. [Google Scholar]
  10. Dimitrova, Z.I. Numerical Investigation of Nonlinear Waves Connected to Blood Flow in an Elastic Tube with Variable Radius. J. Theor. Appl. Mech. 2015, 45, 79–92. [Google Scholar] [CrossRef] [Green Version]
  11. Kawasaki, K.; Ohta, T. Kink Dynamics in One-Dimensional Nonlinear Systems. Phys. A 1982, 116, 573–593. [Google Scholar] [CrossRef]
  12. Dimitrova, Z. On Traveling Waves in Lattices: The Case of Riccati Lattices. J. Theor. Appl. Mech. 2012, 42, 3–22. [Google Scholar] [CrossRef]
  13. Ganji, D.D.; Sabzehmeidani, Y.; Sedighiamiri, A. Nonlinear Systems in Heat Transfer; Elsevier: Amsterdam, The Netherlands, 2018; ISBN 978-0-12-812024-8. [Google Scholar]
  14. Kantz, H.; Schreiber, T. Nonlinear Time Series Analysis; Cambridge University Press: Cambridge, UK, 2004; ISBN 978-0511755798. [Google Scholar]
  15. Verhulst, F. Nonlinear Differential Equations and Dynamical Systems; Springer: Berlin, Germany, 2006; ISBN 978-3-540-60934-6. [Google Scholar]
  16. Mills, T. Applied Time Series Analysis; Academic Press: London, UK, 2019; ISBN 978-012-813117-6. [Google Scholar]
  17. Struble, R. Nonlinear Differential Equations; Dover: New York, NY, USA, 2018; ISBN 978-0486817545. [Google Scholar]
  18. Vitanov, N.K.; Dimitrova, Z.I.; Ausloos, M. Verhulst-Lotka-Volterra Model of Ideological Struggle. Phys. A 2010, 389, 4970–4980. [Google Scholar] [CrossRef] [Green Version]
  19. Grossberg, S. Nonlinear Neural Networks: Principles, Mechanisms, and Architectures. Neural Netw. 1981, 1, 17–61. [Google Scholar] [CrossRef] [Green Version]
  20. Brezis, H.; Browder, F. Partial Differential Equations in the 20th Century. Adv. Math. 1998, 135, 76–144. [Google Scholar] [CrossRef] [Green Version]
  21. Kruskal, M. Nonlinear Wave Equations. In Dynamical Systems, Theory and Applications; Moser, J., Ed.; Springer: Berlin, Germany, 1975; pp. 310–354. [Google Scholar] [CrossRef]
  22. Witting, J. On the Highest and Other Solitary Waves. SIAM J. Appl. Math. 1975, 28, 700–719. [Google Scholar] [CrossRef]
  23. Ablowitz, M.J. Nonlinear Evolution Equations-Continuous and Discrete. SIAM Rev. 1977, 19, 663–684. [Google Scholar] [CrossRef]
  24. Taha, T.R.; Ablowitz, M.J. Analytical and Numerical Aspects of Certain Nonlinear Evolution Equations. I. Analytical. J. Comput. Phys. 1984, 55, 192–202. [Google Scholar] [CrossRef]
  25. Camassa, R.; Hyman, J.M.; Luce, B.P. Nonlinear Waves and Solitons in Physical Systems. Phys. D Nonlinear Phenom. 1998, 123, 1–20. [Google Scholar] [CrossRef]
  26. Zabusky, N.J. Fermi–Pasta–Ulam, Solitons and the Fabric of Nonlinear and Computational Science: History, Synergetics, and Visiometrics. Chaos Interdiscip. J. Nonlinear Sci. 2005, 15, 015102. [Google Scholar] [CrossRef]
  27. Zabusky, N.J.; Kruskal, M.D. Interaction of “Solitons” in a Collisionless Plasma and the Recurrence of Initial States. Phys. Rev. Lett. 1965, 15, 240–243. [Google Scholar] [CrossRef] [Green Version]
  28. Kadomtsev, B.B.; Karpman, V.I. Nonlinear Waves. Sov. Phys. Uspekhi 1971, 14, 40–60. [Google Scholar] [CrossRef]
  29. Hasegawa, A.; Tappert, F. Transmission of Stationary Nonlinear Optical Pulses in Dispersive Dielectric Fibers. I. Anomalous Dispersion. Appl. Phys. Lett. 1973, 23, 142–144. [Google Scholar] [CrossRef]
  30. Boyd, J.P. Equatorial Solitary Waves. Part I: Rossby Solitons. J. Phys. Oceanogr. 1980, 10, 1699–1717. [Google Scholar] [CrossRef]
  31. Boyd, J.P. Equatorial Solitary Waves. Part 2: Envelope Solitons. J. Phys. Oceanogr. 1983, 13, 428–449. [Google Scholar] [CrossRef]
  32. Polturak, E.; deVegvar, P.G.N.; Zeise, E.K.; Lee, D.M. Solitonlike Propagation of Zero Sound in Superfluid 3He. Phys. Rev. Lett. 1981, 46, 1588–1591. [Google Scholar] [CrossRef]
  33. Lonngren, K.E. Soliton Experiments in Plasmas. Plasma Phys. 1983, 25, 943–982. [Google Scholar] [CrossRef]
  34. Burger, S.; Bongs, K.; Dettmer, S.; Ertmer, W.; Sengstock, K.; Sanpera, A.; Shlyapnikov, G.V.; Lewenstein, M. Dark Solitons in Bose–Einstein Condensates. Phys. Rev. Lett. 1999, 83, 5198–5201. [Google Scholar] [CrossRef] [Green Version]
  35. Whitham, G.B. Linear and Nonlinear Waves; Wiley: New York, NY, USA, 1999; ISBN 0-471-35942-4. [Google Scholar]
  36. Tanaka, M. Physics of Nonlinear Waves; Morgan & Claypool: London, UK, 2020; ISBN 9781681737140. [Google Scholar]
  37. Debnath, L. (Ed.) Nonlinear Waves; Cambridge University Press: Cambridge, UK, 1983; ISBN 0-521-25468-X. [Google Scholar]
  38. Infeld, E.; Rowlands, G. Nonlinear Waves, Solitons and Chaos; Cambridge University Press: Cambridge, UK, 2000; ISBN 0-521-63212-9. [Google Scholar]
  39. Popivanov, P.; Slavova, A. Nonlinear Waves: An Introduction; World Scientific: Singapore, 2010; ISBN 9789813107953. [Google Scholar]
  40. Popivanov, P.; Slavova, A. Nonlinear Waves: A Geometric Approach; World Scientific: Singapore, 2018; ISBN 9789813271623. [Google Scholar]
  41. Jang, J.K.; Erkintalo, M.; Murdoch, S.G.; Coen, S. Ultraweak Long-range Interactions of Solitons Observed Over Astronomical Distances. Nat. Photonics 2013, 7, 657–663. [Google Scholar] [CrossRef] [Green Version]
  42. Wadati, M. Introduction to Solitons. Pramana 2001, 57, 841–847. [Google Scholar] [CrossRef]
  43. Craig, W.; Guyenne, P.; Hammack, J.; Henderson, D.; Sulem, C. Solitary Water Wave Interactions. Phys. Fluids 2006, 18, 057106. [Google Scholar] [CrossRef] [Green Version]
  44. Kulikovskiii, A.; Sveshnikova, E. Nonlinear Waves in Elastic Media; CRC Press: Boca Raton, FL, USA, 1995; ISBN 0-8493-8643-8. [Google Scholar]
  45. Maugin, G. Nonlinear Waves in Elastic Crystals; Oxford University Press: Oxford, UK, 1999; ISBN 0-19-853484-1. [Google Scholar]
  46. Mielke, A.; Kirchgässner, K. (Eds.) Structure and Dynamics of Nonlinear Waves in Fluids; World Scientific: Singapore, 1995; ISBN 981-02-2124-X. [Google Scholar]
  47. Osborne, A.R. Nonlinear Topics in Ocean Physics; North-Holland: Amsterdam, The Netherlands, 1991; ISBN 9780444597823. [Google Scholar]
  48. Kundu, A. (Ed.) Tsunami and Nonlinear Waves; Springer: Berlin, Germany, 2007; ISBN 9783540712565. [Google Scholar]
  49. Kluwick, A. (Ed.) Nonlinear Waves in Real Fluids; Springer: Wien, Austria, 2014; ISBN 9783709126080. [Google Scholar]
  50. Ma, Q. Advances in Numerical Simulation of Nonlinear Water Waves; World Scientific: Singapore, 2010; ISBN 9789812836502. [Google Scholar]
  51. Guo, B.; Tian, L.; Yan, Z.; Ling, L.; Wang, Y.-F. Rogue Waves: Mathematical Theory and Applications in Physics; De Gryuiter: Berlin, Germany, 2017; ISBN 9783110470574. [Google Scholar]
  52. Grimshaw, R. (Ed.) Nonlinear Waves in Fluids: Recent Advances and Modern Applications; Springer: Berlin, Germany, 2007; ISBN 9783211380253. [Google Scholar]
  53. Kim, C.-H. Nonlinear Waves and Offshore Structures; World Scientific: Singapore, 2008; ISBN 9789813102484. [Google Scholar]
  54. Mishin, E.; Streltsov, A. Nonlinear Wave and Plasma Structures in the Auroral and Subauroral Geospace; Elsevier: Amsterdam, The Netherlands, 2021; ISBN 9780128209318. [Google Scholar]
  55. Jeffery, A.; Engelbrecht, J. (Eds.) Nonlinear Waves in Solids; Springer: Wien, Austria, 1994; ISBN 978-3-211-82558-7. [Google Scholar]
  56. Nazarov, V.; Radostin, A. Nonlinear Acoustic Waves in Micro-Inhomogeneous Solids; Wiley: Chchester, UK, 2005; ISBN 9781118456088. [Google Scholar]
  57. Cottam, M.G. (Ed.) Linear and Nonlinear Spin Waves in Magnetic Films and Superlattices; World Scientific: Singapore, 1994; ISBN 9789810210069. [Google Scholar]
  58. Maimistov, A.I.; Basharov, A.M. Nonlinear Optical Waves; Kluwer: Dordrecht, The Netherlands, 1999; ISBN 9780792357520. [Google Scholar]
  59. Leble, S.B. Nonlinear Waves in Waveguides; Springer: Berlin, Germany, 2013; ISBN 9783642754203. [Google Scholar]
  60. Epstein, I.R.; Pojman, J.A. An Introduction to Nonlinear Chemical Dynamics; Oxford University Press: Oxford, UK, 1998; ISBN 0-19-509670-3. [Google Scholar]
  61. Bhatnagar, P.L. Nonlinear Waves in One-Dimensional Dispersive Systems; Clarendon Press: Oxford, UK, 1980; ISSN 0964-9174. [Google Scholar]
  62. Belashov, V.Y.; Vladimirov, S.V. Solitary Waves in Dispersive Complex Media; Springer: Berlin, Germany, 2005; ISBN 978-3-540-23376-3. [Google Scholar]
  63. Ablowitz, M.J. Nonlinear Dispersive Waves: Asymptotic Analysis and Solitons; Cambridge University Press: Cambridge, UK, 2011; ISBN 9781107012547. [Google Scholar]
  64. Benjamin, T.B.; Bona, J.L.; Mahony, J.J. Model Equations for Long Waves in Nonlinear Dispersive Systems. Philos. Trans. R. Soc. Lond. Ser. A Math. Phys. Sci. 1972, 272, 47–78. [Google Scholar] [CrossRef]
  65. Karpman, V.I. Non-Linear Waves in Dispersive Media; Pergamon Press: Oxford, UK, 1975; ISBN 0-08-017720-4. [Google Scholar]
  66. Korpel, A.; Banerjee, P.P. A Heuristic Guide to Nonlinear Dispersive Wave Equations and Soliton-Type Solutions. Proc. IEEE 1984, 72, 1109–1130. [Google Scholar] [CrossRef]
  67. Fillipov, A.T. The Versatile Soliton; Springer: New York, NY, USA, 2010; ISBN 9780817649746. [Google Scholar]
  68. Bullough, R.K.; Caudrey, P.J. (Eds.) Solitons; Springer: Berlin, Germany, 2011; ISBN 9783642814501. [Google Scholar]
  69. Longren, K.E.; Scott, A.C. Solitons in Action; Academic Press: New York, NY, USA, 1978; ISBN 9780124555808. [Google Scholar]
  70. Lamb, G.L. Elements of Soliton Theory; Wiley: New York, NY, USA, 1980; ISBN 9780471045595. [Google Scholar]
  71. Drazin, P.G. Solitons; Cambridge University Press: Cambridge, UK, 1983; ISSN 0076-0552. [Google Scholar]
  72. Ablowitz, M.J.; Fuchssteiner, B.; Kruskal, M. Topics in Soliton Theory and Exactly Solvable Nonlinear Equations; World Scientific: Singapore, 1987; ISBN 9789813237957. [Google Scholar]
  73. Rogers, C.; Schief, W.K. Bäcklund and Darboux Transformations: Geometry and Modern Applications in Soliton Theory; Cambridge University Press: Cambridge, UK, 2002; ISBN 9780521012881. [Google Scholar]
  74. Iliev, I.D.; Khristov, E.; Kirchev, K.P. Spectral Methods in Soliton Equations; Kongman: Harlow, UK, 1994; ISBN 0-582-23963-X. [Google Scholar]
  75. Manton, N.; Sutcliffe, P. Topological Solitons; Cambridge University Press: Cambridge, UK, 2004; ISBN 9781139454698. [Google Scholar]
  76. Cercignani, C. Solitons-Theory and Application. Nuovo C. Riv. 1977, 7, 429–469. [Google Scholar] [CrossRef]
  77. Akhmediev, N.; Ankiewicz, A. Dissipative Solitons; Springer: Berlin, Germany, 2005; ISBN 9783540233732. [Google Scholar]
  78. Belinski, V.A.; Verdaguer, E. Gravitational Solitons; Cambridge University Press: Cambridge, UK, 2001; ISBN 9780521805865. [Google Scholar]
  79. Weigel, H. Chiral Soliton Models for Baryons; Springer: Berlin, Germany, 2007; ISBN 9783540754367. [Google Scholar]
  80. Ricketts, D.S.; Ham, D. Electrical Solitons: Theory, Design, and Applications; CRC Press: Boca Raton, FL, USA, 2011; ISBN 9781351833691. [Google Scholar]
  81. Zakharov, V.E.; Wabnitz, S. (Eds.) Optical Solitons: Theoretical Challenges and Industrial Perspectives; Springer: Berlin, Germany, 2013; ISBN 9783662038079. [Google Scholar]
  82. Drazin, P.G.; Johnson, R.S. Solitons: An Introduction; Cambridge University Press: Cambridge, UK, 1989; ISBN 9780521336550. [Google Scholar]
  83. Novikov, S.; Manakov, S.V.; Pitaevskii, L.P.; Zakharov, V.E. Theory of Solitons: The Inverse Scattering Method; Springer: Berlin, Germany, 1984; ISBN 9780306109775. [Google Scholar]
  84. Dickey, L.A. Soliton Equations and Hamiltonian Systems; World Scientific: Singapore, 2003; ISBN 9789812794512. [Google Scholar]
  85. Malomed, B.A. Soliton Management in Periodic Systems; Springer: New York, NY, USA, 2006; ISBN 9780387293349. [Google Scholar]
  86. Davydov, A.S. Solitons in Molecular Systems; Springer: Dordrecht, The Netherlands, 2013; ISBN 9789401730259. [Google Scholar]
  87. Olver, P.J.; Sattiger, D.H. Solitons in Physics, Mathematics, and Nonlinear Optics; Springer: New York, NY, USA, 2012; ISBN 9781461390336. [Google Scholar]
  88. Dauxois, T.; Peyrard, M. Physics of Solitons; Cambridge University Press: Cambridge, UK, 2006; ISBN 9780521854214. [Google Scholar]
  89. Scott, A.C.; Chu, F.Y.F.; McLaughlin, D.W. The Soliton: A New Concept in Applied Science. Proc. IEEE 1973, 61, 1443–1483. [Google Scholar] [CrossRef]
  90. Karpman, V.I. Soliton Evolution in the Presence of Perturbation. Phys. Scr. 1979, 20, 462–478. [Google Scholar] [CrossRef]
  91. Whitam, G.B. Lectures on Wave Propagation; Springer: Berlin, Germany, 1979; ISBN 3-540-08945-4. [Google Scholar]
  92. Jimbo, M.; Miwa, T. Solitons and Infinite Dimensional Lie Algebras. Publ. Res. Inst. Math. Sci. Kioto Univ. 1983, 19, 943–1001. [Google Scholar] [CrossRef] [Green Version]
  93. Gibbon, J.D. A Survey of the Origins and Physical Importance of Soliton Equations. Philos. Trans. R. Soc. Lond. Ser. A Math. Phys. Sci. 1985, 315, 335–365. [Google Scholar] [CrossRef]
  94. Newell, A.C. Solitons in Mathematics and Physics; Society for Industrial and Applied Mathematics: Philadelphia, PA, USA, 1985; ISBN 978-0-898711-96-7. [Google Scholar]
  95. Kivshar, Y.S.; Malomed, B.A. Dynamics of Solitons in Nearly Integrable Systems. Rev. Mod. Phys. 1989, 61, 763–915. [Google Scholar] [CrossRef] [Green Version]
  96. Remoissenet, M. Waves Called Solitons: Concepts and Experiments; Springer: Berlin, Germany, 2013; ISBN 978-3-642-05819-2. [Google Scholar]
  97. Yu, L. Solitons & Polarons in Conducting Polymers; World Scientific: Singapore, 1988; ISBN 9789971500542. [Google Scholar]
  98. Maxon, S. Cylindrical Solitons in a Warm, Multi-Ion Plasma. Phys. Fluids 1976, 19, 266–271. [Google Scholar] [CrossRef]
  99. Yajima, T.; Wadati, M. Solitons in Electron Beam Plasma. J. Phys. Soc. Jpn. 1990, 59, 3237–3248. [Google Scholar] [CrossRef]
  100. Huberman, B.A. Superfluid Solitons in Helium Films. Phys. Rev. Lett. 1978, 41, 1389–1393. [Google Scholar] [CrossRef]
  101. Kodama, Y.; Maruta, A.; Hasegawa, A. Long Distance Communications with Solitons. Quantum Opt. B 1994, 6, 463–516. [Google Scholar] [CrossRef]
  102. Gardner, C.S.; Greene, J.M.; Kruskal, M.D.; Miura, R.M. Method for Solving the Korteweg-deVries Equation. Phys. Rev. Lett. 1967, 19, 1095–1097. [Google Scholar] [CrossRef]
  103. Ablowitz, M.J.; Kaup, D.J.; Newell, A.C.; Segur, H. The Inverse Scattering Transform-Fourier Analysis for Nonlinear Problems. Stud. Appl. Math. 1974, 53, 249–315. [Google Scholar] [CrossRef]
  104. Ablowitz, M.J. Lectures on the Inverse Scattering Transform. Stud. Appl. Math. 1978, 58, 17–94. [Google Scholar] [CrossRef]
  105. Wadati, M.; Kamijo, T. On the Extension of Inverse Scattering Method. Prog. Theor. Phys. 1974, 52, 397–414. [Google Scholar] [CrossRef] [Green Version]
  106. Lax, P.D. Integrals of Nonlinear Equations of Evolution and Solitary Waves. Commun. Pure Appl. Math. 1968, 21, 467–490. [Google Scholar] [CrossRef]
  107. Zakharov, V.E.E.; Shabat, A.B. A Scheme for Integrating the Nonlinear Equations of Mathematical Physics by the Method of the Inverse Scattering Problem. I. Funktsional’nyi Anal. I Ego Prilozheniya 1974, 8, 43–53. [Google Scholar] [CrossRef]
  108. Zakharov, V.E.E.; Shabat, A.B. Integration of Nonlinear Equations of Mathematical Physics by the Method of Inverse Scattering. II. Funktsional’nyi Anal. I Ego Prilozheniya 1979, 13, 13–22. [Google Scholar] [CrossRef]
  109. Fokas, A.S. A Unified Transform Method for Solving Linear and Certain Nonlinear PDEs. Proc. R. Soc. Lond. Ser. A Math. Phys. Eng. Sci. 1997, 453, 1411–1443. [Google Scholar] [CrossRef]
  110. Fokas, A.S.; Ablowitz, M.J. Method of Solution for a Class of Multidimensional Nonlinear Evolution Equations. Phys. Rev. Lett. 1983, 51, 7–10. [Google Scholar] [CrossRef]
  111. Manakov, S.V. The Method of the Inverse Scattering Problem, and Two-dimensional Evolution Equations. Uspekhi Mat. Nauk 1976, 31, 245–246. [Google Scholar]
  112. Manakov, S.V. The Inverse Scattering Transform for the Time-dependent Schrödinger Equation and Kadomtsev-Petviashvili Equation. Physica D 1981, 3, 420–427. [Google Scholar] [CrossRef]
  113. Dodd, R.K.; Bullough, R.K. The Generalised Marchenko Equation and the Canonical Structure of the AKNS-ZS Inverse Method. Phys. Scr. 1979, 20, 514–530. [Google Scholar] [CrossRef]
  114. Ablowitz, M.J.; Kruskal, M.; Segur, H. A Note on Miura’s Transformation. J. Math. Phys. 1979, 20, 999–1003. [Google Scholar] [CrossRef]
  115. Kaup, D.J.; Newell, A.C. Evolution Equations, Singular Dispersion Relations, and Moving Eigenvalues. Adv. Math. 1979, 31, 67–100. [Google Scholar] [CrossRef] [Green Version]
  116. Newell, A.C.; Redekopp, L.G. Breakdown of Zakharov-Shabat Theory and Soliton Creation. Phys. Rev. Lett. 1977, 38, 377–380. [Google Scholar] [CrossRef]
  117. Beals, R.; Coifman, R.R. Scattering and Inverse Scattering for First Order Systems. Commun. Pure Appl. Math. 1984, 37, 39–90. [Google Scholar] [CrossRef]
  118. Kaup, D.J. A perturbation Expansion for the Zakharov–Shabat Inverse Scattering Transform. SIAM J. Appl. Math. 1976, 31, 121–133. [Google Scholar] [CrossRef]
  119. Kaup, D.J. On the inverse scattering problem for cubic eigenvalue problems of the class ψxxx + 6x + 6 = λψ. Stud. Appl. Math. 1980, 62, 189–216. [Google Scholar] [CrossRef]
  120. Beals, R.; Coifman, R.R. The D-Bar Approach to Inverse Scattering and Nonlinear Evolutions. Phys. D Nonlinear Phenom. 1986, 18, 242–249. [Google Scholar] [CrossRef]
  121. Ablowitz, M.J. Applications of Slowly Varying Nonlinear Disperive Wave Theories. Stud. Appl. Math. 1971, 50, 329–344. [Google Scholar] [CrossRef]
  122. Goriely, A. Integrability and Nonintegrability of Dynamical Systems; World Scientific: Singapore, 2001; ISBN 981-02-3533-X. [Google Scholar]
  123. Calogero, F. A Method to Generate Solvable Nonlinear Evolution Equations. Lett. Al Nuovo Cimento 1975, 14, 443–447. [Google Scholar] [CrossRef]
  124. Newell, A.C. The General Structure of Integrable Evolution Equations. Proc. R. Soc. Lond. A Math. Phys. Sci. 1979, 365, 283–311. [Google Scholar] [CrossRef]
  125. Gupta, M.R. Exact Inverse Scattering Solution of a Non-Linear Evolution Equation in a Non-Uniform Medium. Phys. Lett. A 1979, 72, 420–422. [Google Scholar] [CrossRef]
  126. Shimizu, T.; Wadati, M. A New Integrable Nonlinear Evolution Equation. Prog. Theor. Phys. 1980, 63, 808–820. [Google Scholar] [CrossRef] [Green Version]
  127. Fokas, A.S. On the Integrability of Linear and Nonlinear Partial Differential Equations. J. Math. Phys. 2000, 41, 4188–4237. [Google Scholar] [CrossRef]
  128. Degasperis, A.; Lombardo, S. Multicomponent Integrable Wave Equations: I. Darboux- Dressing Transformation. J. Phys. A Math. Theor. 2007, 40, 961–977. [Google Scholar] [CrossRef] [Green Version]
  129. Kupershmidt, B.A. Integrable and Superintegrable Systems; World Scientific: Singapore, 1990; ISBN 9789810203160. [Google Scholar]
  130. Yang, J. Nonlinear Waves in Integrable and Nonintegrable Systems; SIAM: Philadelphia, PA, USA, 2010; ISBN 9780898719680. [Google Scholar]
  131. Fokas, A.S. Integrable Nonlinear Evolution Partial Differential Equations in 4+2 and 3+1 Dimensions. Phys. Rev. Lett. 2006, 96, 190201. [Google Scholar] [CrossRef] [PubMed]
  132. Zakharov, V.E.E.; Manakov, S.V. Construction of Higher-dimensional Nonlinear Integrable Systems and of Their Solutions. Funct. Anal. Its Appl. 1985, 19, 89–101. [Google Scholar] [CrossRef]
  133. Calogero, F.; Degasperis, A. Extension of the Spectral Transform Method for Solving Nonlinear Evolution Equations. Lett. Al Nuovo Cimento 1978, 22, 131–137. [Google Scholar] [CrossRef]
  134. Kichenassamy, S. Nonlinear Wave Equations; CRC Press: Boca Raton, FL, USA, 2021; ISBN 9781003210276. [Google Scholar]
  135. Deift, P.; Lund, F.; Trubowitz, E. Nonlinear Wave Equations and Constrained Harmonic Motion. Proc. Natl. Acad. Sci. USA 1980, 77, 716–719. [Google Scholar] [CrossRef] [Green Version]
  136. Hirota, R. Exact Solution of the Korteweg-de Vries Equation for Multiple Collisions of Solitons. Phys. Rev. Lett. 1971, 27, 1192–1194. [Google Scholar] [CrossRef]
  137. Hirota, R.; Satsuma, J. Soliton Solutions of a Coupled Korteweg-de Vries Equation. Phys. Lett. A 1981, 85, 407–408. [Google Scholar] [CrossRef]
  138. Gibbon, J.D.; Radmore, P.; Tabor, M.; Wood, D. The Painlev’e Property and Hirota’s Method. Stud. Appl. Math. 1985, 72, 39–63. [Google Scholar] [CrossRef]
  139. Matsukidaira, J.; Satsuma, J.; Strampp, W. Soliton Equations Expressed by Trilinear Forms and Their Solutions. Phys. Lett. A 1990, 147, 467–471. [Google Scholar] [CrossRef]
  140. Satsuma, J.; Kajiwara, K.; Matsukidaira, J.; Hietarinta, J. Solutions of the Broer-Kaup System Through Its Trilinear Form. J. Phys. Soc. Jpn. 1992, 61, 3096–3102. [Google Scholar] [CrossRef]
  141. Hirota, R. The Direct Method in Soliton Theory; Cambridge University Press: Cambridge, UK, 2004; ISBN 0-521-83660-3. [Google Scholar]
  142. Miura, R.M.; Gardner, C.S.; Kruskal, M.D. Korteweg-de Vries Equation and Generalizations. II. Existence of Conservation Laws and Constants of Motion. J. Math. Phys. 1968, 9, 1204–1209. [Google Scholar] [CrossRef]
  143. Gardner, C.S.; Greene, J.M.; Kruskal, M.D.; Miura, R.M. Korteweg-de Vries Equation and Generalizations. VI. Methods for Exact Solution. Commun. Pure Appl. Math. 1974, 27, 97–133. [Google Scholar] [CrossRef]
  144. Bona, J.L.; Smith, R. The Initial-Value Problem for the Korteweg-de Vries Equation. Philos. Trans. R. Soc. Lond. Ser. A Math. Phys. Sci. 1975, 278, 555–601. [Google Scholar] [CrossRef] [Green Version]
  145. Novikov, S.P. The Periodic Problem for the Korteweg–de Vries Equation. Funktsional’nyi Anal. I Ego Prilozheniya 1974, 8, 54–66. [Google Scholar] [CrossRef]
  146. Zabusky, N.J.; Galvin, C.J. Shallow-water Waves, the Korteweg-de Vries Equation and Solitons. J. Fluid Mech. 1971, 47, 811–824. [Google Scholar] [CrossRef]
  147. Wadati, M. The Exact Solution of the Modified Korteweg-de Vries Equation. J. Phys. Soc. Jpn. 1972, 32, 1681. [Google Scholar] [CrossRef]
  148. Hirota, R. Nonlinear Partial Difference Equations. I. A Difference Analogue of the Korteweg-de Vries Equation. J. Phys. Soc. Jpn. 1977, 43, 1424–1433. [Google Scholar] [CrossRef]
  149. Wahlquist, H.D.; Estabrook, F.B. Bäcklund Transformation for Solutions of the Korteweg-de Vries Equation. Phys. Rev. Lett. 1973, 31, 1386–1388. [Google Scholar] [CrossRef]
  150. Wadati, M. The Modified Korteweg-de Vries Equation. J. Phys. Soc. Jpn. 1973, 34, 1289–1296. [Google Scholar] [CrossRef]
  151. Sawada, K.; Kotera, T. A Method for Finding N-soliton Solutions of the KdV Equation and KdV-like Equation. Prog. Theor. Phys. 1974, 51, 1355–1367. [Google Scholar] [CrossRef] [Green Version]
  152. Lax, P.D. Periodic Solutions of the KdV Equation. Commun. Pure Appl. Math. 1975, 28, 141–188. [Google Scholar] [CrossRef]
  153. Miura, R.M. The Korteweg–deVries Equation: A Survey of Results. SIAM Rev. 1976, 18, 412–459. [Google Scholar] [CrossRef]
  154. Caudrey, P.J.; Dodd, R.K.; Gibbon, J.D. A New Hierarchy of Korteweg–de Vries Equations. Proc. R. Soc. Lond. A Math. Phys. Sci. 1976, 351, 407–422. [Google Scholar] [CrossRef]
  155. Satsuma, J. N-soliton Solution of the Two-dimensional Korteweg-deVries Equation. J. Phys. Soc. Jpn. 1976, 40, 286–290. [Google Scholar] [CrossRef]
  156. Ablowitz, M.J.; Segur, H. Asymptotic Solutions of the Korteweg-deVries Equation. Stud. Appl. Math. 1977, 57, 13–44. [Google Scholar] [CrossRef]
  157. Hammack, J.L.; Segur, H. The Korteweg-de Vries Equation and Water Waves. Part 3. Oscillatory Waves. J. Fluid Mech. 1978, 84, 337–358. [Google Scholar] [CrossRef] [Green Version]
  158. Grimshaw, R. Slowly Varying Solitary Waves. I. Korteweg-de Vries Equation. Proc. R. Soc. Lond. A Math. Phys. Sci. 1979, 368, 359–375. [Google Scholar] [CrossRef]
  159. Ito, M. An Extension of Nonlinear Evolution Equations of the K-dV (mK-dV) Type to Higher Orders. J. Phys. Soc. Jpn. 1980, 49, 771–778. [Google Scholar] [CrossRef]
  160. Miles, J.W. The Korteweg-de Vries Equation: A Historical Essay. J. Fluid Mech. 1981, 106, 131–147. [Google Scholar] [CrossRef]
  161. Lax, P.D.; David Levermore, C. The Small Dispersion Limit of the Korteweg-de Vries Equation. I. Commun. Pure Appl. Math. 1983, 36, 253–290. [Google Scholar] [CrossRef]
  162. Wadati, M. Stochastic Korteweg-de Vries Equation. J. Phys. Soc. Jpn. 1983, 52, 2642–2648. [Google Scholar] [CrossRef]
  163. Bona, J.L.; Souganidis, P.E.; Strauss, W.A. Stability and Instability of Solitary Waves of Korteweg-de Vries Type. Proc. R. Soc. Lond. A Math. Phys. Sci. 1987, 411, 395–412. [Google Scholar] [CrossRef]
  164. Boyd, J.P. Theta Functions, Gaussian Series, and Spatially Periodic Solutions of the Korteweg–de Vries Equation. J. Math. Phys. 1982, 23, 375–387. [Google Scholar] [CrossRef] [Green Version]
  165. Nijhoff, F.; Capel, H. The Discrete Korteweg-de Vries Equation. Acta Appl. Math. 1995, 39, 133–158. [Google Scholar] [CrossRef] [Green Version]
  166. Crighton, D.G. Applications of KdV. Acta Appl. Math. 1995, 39, 39–67. [Google Scholar] [CrossRef]
  167. Kenig, C.E.; Ponce, G.; Vega, L. A Bilinear Estimate with Applications to the KdV Equation. J. Am. Math. Soc. 1996, 9, 573–603. [Google Scholar] [CrossRef]
  168. Goncharenko, V.M. Multisoliton Solutions of the Matrix KdV Equation. Theor. Math. Phys. 2001, 126, 81–91. [Google Scholar] [CrossRef]
  169. Linares, F.; Ponce, G. Introduction to Nonlinear Dispersive Equations; Springer: New York, NY, USA, 2014; ISBN 978-0-38784-898-3. [Google Scholar]
  170. Dubrovin, B.A.E.; Matveev, V.B.; Novikov, S.P. Non-linear Equations of Korteweg-de Vries Type, Finite-zone Linear Operators, and Abelian Varieties. Russ. Math. Surv. 1976, 31, 59–146. [Google Scholar] [CrossRef]
  171. Ablowitz, M.J.; Kaup, D.J.; Newell, A.C.; Segur, H. Nonlinear Evolution Equations of Physical Significance. Phys. Rev. Lett. 1973, 31, 125–127. [Google Scholar] [CrossRef]
  172. Ting, A.C.; Chen, H.H.; Lee, Y.C. Exact Solutions of a Nonlinear Boundary Value Problem: The Vortices of the Two-Dimensional sinh-Poisson Equation. Phys. D Nonlinear Phenom. 1987, 26, 37–66. [Google Scholar] [CrossRef]
  173. Rubinstein, J. Sine-Gordon Equation. J. Math. Phys. 1970, 11, 258–266. [Google Scholar] [CrossRef]
  174. Hirota, R. Exact Solution of the Sine-Gordon Equation for Multiple Collisions of Solitons. J. Phys. Soc. Jpn. 1972, 33, 1459–1463. [Google Scholar] [CrossRef]
  175. Scott, A.C.; Chu, F.Y.; Reible, S.A. Magnetic-flux Propagation on a Josephson Transmission Line. J. Appl. Phys. 1976, 47, 3272–3286. [Google Scholar] [CrossRef]
  176. Dodd, R.K.; Bullough, R.K. Bäcklund Transformations for the Sine–Gordon Equations. Proc. R. Soc. Lond. A Math. Phys. Sci. 1976, 351, 499–523. [Google Scholar] [CrossRef]
  177. Orfanidis, S.J. Sine-Gordon Equation and Nonlinear σ Model on a Lattice. Phys. Rev. D 1978, 18, 3828–3832. [Google Scholar] [CrossRef]
  178. Pöppe, C. Construction of Solutions of The Sine-Gordon Equation by Means of Fredholm Determinants. Phys. D Nonlinear Phenom. 1983, 9, 103–139. [Google Scholar] [CrossRef]
  179. Weiss, J. The Sine-Gordon Equations: Complete and Partial Integrability. J. Math. Phys. 1984, 25, 2226–2235. [Google Scholar] [CrossRef]
  180. Malomed, B.A. Emission From, Quasi-Classical Quantization, and Stochastic Decay of Sine-Gordon Solitons in External Fields. Phys. D Nonlinear Phenom. 1987, 27, 113–157. [Google Scholar] [CrossRef]
  181. Ablowitz, M.J.; Kaup, D.J.; Newell, A.C.; Segur, H. Method for Solving the Sine- Gordon Equation. Phys. Rev. Lett. 1973, 30, 1262–1264. [Google Scholar] [CrossRef]
  182. Ablowitz, M.J.; Prinari, B.; Trubatch, A.D. Discrete and Continuous Nonlinear Schrödinger Systems; Cambridge University Press: Cambridge, UK, 2004; ISBN 0-521-52437-2. [Google Scholar]
  183. Zabusky, N.J. Solitons and Bound States of the Time-independent Schrödinger Equation. Phys. Rev. 1968, 168, 124–128. [Google Scholar] [CrossRef]
  184. Zakharov, V.E.; Manakov, S.V. On the Complete Integrability of a Nonlinear Schrödinger Equation. Theor. Math. Phys. 1974, 19, 551–559. [Google Scholar] [CrossRef]
  185. Segur, H. Asymptotic Solutions and Conservation Laws for the Nonlinear Schrödinger Equation. II. J. Math. Phys. 1976, 17, 714–716. [Google Scholar] [CrossRef]
  186. Nogami, Y.; Warke, C.S. Soliton Solutions of Multicomponent Nonlinear Schrödinger Equation. Phys. Lett. A 1976, 59, 251–253. [Google Scholar] [CrossRef]
  187. Jaulent, M.; Miodek, I. Nonlinear Evolution Equations Associated with ‘Enegry-Dependent Schrödinger potentials’. Lett. Math. Phys. 1976, 1, 243–250. [Google Scholar] [CrossRef]
  188. Pereira, N.R. Soliton in the Damped Nonlinear Schrödinger Equation. Phys. Fluids 1977, 20, 1735–1743. [Google Scholar] [CrossRef] [Green Version]
  189. Fukuda, I.; Tsutsumi, M. On Coupled Klein-Gordon-Schrödinger Equations, II. J. Math. Anal. Appl. 1978, 66, 358–378. [Google Scholar] [CrossRef] [Green Version]
  190. Kawata, T.; Inoue, H. Exact Solutions of the Derivative Nonlinear Schrödinger Equation Under the Nonvanishing Conditions. J. Phys. Soc. Jpn. 1978, 44, 1968–1976. [Google Scholar] [CrossRef]
  191. Kaup, D.J.; Newell, A.C. An Exact Solution for a Derivative Nonlinear Schrödinger Equation. J. Math. Phys. 1978, 19, 798–801. [Google Scholar] [CrossRef]
  192. Ma, Y.C. The Perturbed Plane-Wave Solutions of the Cubic Schrödinger Equation. Stud. Appl. Math. 1979, 60, 43–58. [Google Scholar] [CrossRef]
  193. Boiti, M.; Pempinelli, F. Nonlinear Schrödinger Equation, Bäcklund Transformations and Painlev’e Transcendents. Il Nuovo C. B 1980, 59, 40–58. [Google Scholar] [CrossRef]
  194. Balakrishnan, R. Dynamics of a Generalised Classical Heisenberg Chain. Phys. Lett. A 1982, 92, 243–246. [Google Scholar] [CrossRef]
  195. Cohen, A.; Kappeler, T. Scattering and Inverse Scattering for Steplike Potentials in the Schrödinger Equation. Indiana Univ. Math. J. 1985, 34, 127–180. [Google Scholar] [CrossRef]
  196. Hasse, R.W. A General Method for the Solution of Nonlinear Soliton and Kink Schrödinger Equations. Z. Phys. B Condens. Matter 1980, 37, 83–87. [Google Scholar] [CrossRef]
  197. Ma, Y.C.; Ablowitz, M.J. The Periodic Cubic Schrödinger Equation. Stud. Appl. Math. 1981, 65, 113–158. [Google Scholar] [CrossRef]
  198. Clarkson, P.A.; Cosgrove, C.M. Painlevé Analysis of the Non-Linear Schrödinger Family of Equations. J. Phys. A Math. Gen. 1987, 20, 2003–2024. [Google Scholar] [CrossRef]
  199. Olmedilla, E. Multiple Pole Solutions of The Non-Linear Schrödinger Equation. Phys. D Nonlinear Phenom. 1987, 25, 330–346. [Google Scholar] [CrossRef]
  200. Tracy, E.R.; Chen, H.H. Nonlinear Self-modulation: An Exactly Solvable Model. Phys. Rev. A 1988, 37, 815–839. [Google Scholar] [CrossRef]
  201. Davies, B. Higher Conservation Laws for the Quantum Non-linear Schrödinger Equation. Phys. A Stat. Mech. Its Appl. 1990, 167, 433–456. [Google Scholar] [CrossRef]
  202. Potasek, M.J.; Tabor, M. Exact Solutions for an Extended Nonlinear Schrödinger Equation. Phys. Lett. A 1991, 154, 449–452. [Google Scholar] [CrossRef]
  203. Clarkson, P.A. Dimensional Reductions and Exact Solutions of a Generalized Nonlinear Schrödinger Equation. Nonlinearity 1992, 5, 453–472. [Google Scholar] [CrossRef]
  204. Tasgal, R.S.; Potasek, M.J. Soliton Solutions to Coupled Higher-Order Nonlinear Schrödinger Equations. J. Math. Phys. 1992, 33, 1208–1215. [Google Scholar] [CrossRef]
  205. Boffetta, G.; Osborne, A.R. Computation of the Direct Scattering Transform for the Nonlinear Schrödinger Equation. J. Comput. Phys. 1992, 102, 252–264. [Google Scholar] [CrossRef]
  206. Manas, M. Darboux Transformations for the Nonlinear Schrödinger Equations. J. Phys. A Math. Gen. 1996, 29, 7721–7737. [Google Scholar] [CrossRef]
  207. Nakkeeran, K.; Porsezian, K.; Sundaram, P.S.; Mahalingam, A. Optical Solitons in N-Coupled Higher Order Nonlinear Schrödinger Equations. Phys. Rev. Lett. 1998, 80, 1425–1428. [Google Scholar] [CrossRef]
  208. Calogero, F.; Degasperis, A.; Xiaoda, J. Nonlinear Schrödinger-type Equations from Multiscale Reduction of PDEs. I. Systematic Derivation. J. Math. Phys. 2000, 41, 6399–6443. [Google Scholar] [CrossRef]
  209. Serkin, V.N.; Belyaeva, T.Y.L. High-Energy Optical Schrödinger Solitons. J. Exp. Theor. Phys. Lett. 2001, 74, 573–577. [Google Scholar] [CrossRef]
  210. Tenorio, C.H.; Vargas, E.V.; Serkin, V.N.; Granados, M.A.; Belyaeva, T.L.; Moreno, R.P.; Lara, L.M. Dynamics of Solitons in the Model of Nonlinear Schrödinger Equation with an External Harmonic Potential: II. Dark Solitons. Quantum Electron. 2005, 35, 929–937. [Google Scholar] [CrossRef]
  211. Aktosun, T.; Demontis, F.; Van Der Mee, C. Exact Solutions to the Focusing Nonlinear Schrödinger Equation. Inverse Probl. 2007, 23, 2171–2195. [Google Scholar] [CrossRef]
  212. Sulem, C.; Sulem, P.L. The Nonlinear Schrödinger Equation: Self-Focusing and Wave Collapse; Springer: Berlin, Germany, 2007; ISBN 0-387-98611-1. [Google Scholar]
  213. Bottman, N.; Deconinck, B.; Nivala, M. Elliptic Solutions of the Defocusing NLS Equation are Stable. J. Phys. A Math. Theor. 2011, 44, 285201. [Google Scholar] [CrossRef]
  214. Suslov, S. On Integrability of Nonautonomous Nonlinear Schrödinger Equations. Proc. Am. Math. Soc. 2012, 140, 3067–3082. [Google Scholar] [CrossRef]
  215. Zhai, B.G.; Zhang, W.G.; Wang, X.L.; Zhang, H.Q. Multi-Rogue Waves and Rational Solutions of The Coupled Nonlinear Schrödinger Equations. Nonlinear Anal. Real World Appl. 2013, 14, 14–27. [Google Scholar] [CrossRef]
  216. Ling, L.; Zhao, L.C. Simple Determinant Representation for Rogue Waves of the Nonlinear Schrödinger Equation. Phys. Rev. E 2013, 88, 043201. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  217. Chowdury, A.; Kedziora, D.J.; Ankiewicz, A.; Akhmediev, N. Soliton Solutions of an Integrable Nonlinear Schrödinger Equation With Quintic Terms. Phys. Rev. E 2014, 90, 032922. [Google Scholar] [CrossRef]
  218. Kedziora, D.J.; Ankiewicz, A.; Chowdury, A.; Akhmediev, N. Integrable Equations of the Infinite Nonlinear Schrödinger Equation Hierarchy with Time Variable Coefficients. Chaos Interdiscip. J. Nonlinear Sci. 2015, 25, 103114. [Google Scholar] [CrossRef]
  219. Ling, L.; Zhao, L.C.; Guo, B. Darboux Transformation and Multi-Dark Soliton for N-Component Nonlinear Schrödinger Equations. Nonlinearity 2015, 28, 3243–3271. [Google Scholar] [CrossRef]
  220. Fokas, A.S. Integrable Multidimensional Versions of the Nonlocal Nonlinear Schrödinger Equation. Nonlinearity 2016, 29, 319. [Google Scholar] [CrossRef]
  221. Biondini, G.; Mantzavinos, D. Long-Time Asymptotics for the Focusing Nonlinear Schrödinger Equation with Nonzero Boundary Conditions at Infinity and Asymptotic Stage of Modulational Instability. Commun. Pure Appl. Math. 2017, 70, 2300–2365. [Google Scholar] [CrossRef] [Green Version]
  222. Feng, B.F.; Luo, X.D.; Ablowitz, M.J.; Musslimani, Z.H. General Soliton Solution to a Nonlocal Nonlinear Schrödinger Equation with Zero and Nonzero Boundary Conditions. Nonlinearity 2018, 31, 5385–5409. [Google Scholar] [CrossRef]
  223. Ablowitz, M.J.; Feng, B.F.; Luo, X.D.; Musslimani, Z.H. Inverse Scattering Transform For The Nonlocal Reverse Space–Time Nonlinear Schrödinger Equation. Theor. Math. Phys. 2018, 196, 1241–1267. [Google Scholar] [CrossRef]
  224. Bilman, D.; Miller, P.D. A Robust Inverse Scattering Transform for the Focusing Nonlinear Schrödinger Equation. Commun. Pure Appl. Math. 2019, 72, 1722–1805. [Google Scholar] [CrossRef] [Green Version]
  225. Fedele, R.; Schamel, H. Solitary Waves in the Madelung’s Fluid: Connection Between the Nonlinear Schrödinger Equation and the Korteweg-de Vries Equation. Eur. Phys. J. B-Condens. Matter Complex Syst. 2002, 27, 313–320. [Google Scholar] [CrossRef]
  226. Fedele, R.; Schamel, H.; Shukla, P.K. Solitons in the Madelung’s Fluid. Phys. Scr. 2002, T98, 18–23. [Google Scholar] [CrossRef]
  227. Fedele, R. Envelope Solitons versus Solitons. Phys. Scr. 2002, 65, 502–508. [Google Scholar] [CrossRef]
  228. Serkin, V.N.; Hasegawa, A. Exactly Integrable Nonlinear Schrodinger Equation Models with Varying Dispersion, Nonlinearity and Gain: Application for Soliton Dispersion. IEEE J. Sel. Top. Quantum Electron. 2002, 8, 418–431. [Google Scholar] [CrossRef]
  229. Hyman, J.M.; Nicolaenko, B. The Kuramoto-Sivashinsky Equation: A Bridge Between PDE’s and Dynamical Systems. Phys. D Nonlinear Phenom. 1986, 18, 113–126. [Google Scholar] [CrossRef]
  230. Quispel, G.R.W.; Nijhoff, F.W.; Capel, H.W. Linearization of the Boussinesq Equation and the Modified Boussinesq Equation. Phys. Lett. A 1982, 91, 143–145. [Google Scholar] [CrossRef]
  231. Clarkson, P.A. Nonclassical Symmetry Reductions of the Boussinesq Equation. Chaos Solitons Fractals 1995, 5, 2261–2301. [Google Scholar] [CrossRef]
  232. Clarkson, P.A.; Dowie, E. Rational Solutions of the Boussinesq Equation and Applications to Rogue Waves. Trans. Math. Its Appl. 2017, 1, tnx003. [Google Scholar] [CrossRef]
  233. Johnson, R.S. A Two-Dimensional Boussinesq Equation for Water Waves and Some of its Solutions. J. Fluid Mech. 1996, 323, 65–78. [Google Scholar] [CrossRef]
  234. Tsuzuki, T. Nonlinear Waves in the Pitaevskii-Gross Equation. J. Low Temp. Phys. 1971, 4, 441–457. [Google Scholar] [CrossRef]
  235. Clarkson, P.A. Painlevé Analysis and the Complete Integrability of a Generalized Variable- Coefficient Kadomtsev-Petviashvili Equation. IMA J. Appl. Math. 1990, 44, 27–53. [Google Scholar] [CrossRef]
  236. Zhao, H.Q.; Ma, W.X. Mixed Lump–kink Solutions to the KP Equation. Comput. Math. Appl. 2017, 74, 1399–1405. [Google Scholar] [CrossRef]
  237. Ma, Y.C. The Complete Solution of the Long-wave–Short-wave Resonance Equations. Stud. Appl. Math. 1978, 59, 201–221. [Google Scholar] [CrossRef]
  238. Fokas, A.S.; Sung, L.Y. On the Solvability of the N-wave, Davey-Stewartson and Kadomtsev-Petviashvili equations. Inverse Probl. 1992, 8, 673–708. [Google Scholar] [CrossRef]
  239. Rao, J.; Zhang, Y.; Fokas, A.S.; He, J. Rogue Waves of the Nonlocal Davey–Stewartson I Equation. Nonlinearity 2018, 31, 4090–4104. [Google Scholar] [CrossRef]
  240. Bock, T.L.; Kruskal, M.D. A Two-Parameter Miura Transformation of the Benjamin-Ono Equation. Phys. Lett. A 1979, 74, 173–176. [Google Scholar] [CrossRef]
  241. Kaup, D.J.; Matsuno, Y. The Inverse Scattering Transform for the Benjamin–Ono Equation. Stud. Appl. Math. 1998, 101, 73–98. [Google Scholar] [CrossRef]
  242. Fokas, A.S.; Fuchssteiner, B. The Hierarchy of the Benjamin-Ono Equation. Phys. Lett. A 1981, 86, 341–345. [Google Scholar] [CrossRef]
  243. Parker, A. On the Camassa-Holm Equation and a Direct Method of Solution I. Bilinear Form and Solitary Waves. Proc. R. Soc. Lond. Ser. A Math. Phys. Eng. Sci. 2004, 460, 2929–2957. [Google Scholar] [CrossRef]
  244. Parker, A. On the Camassa–Holm Equation and a Direct Method of Solution. III. N-soliton Solutions. Proc. R. Soc. A Math. Phys. Eng. Sci. 2005, 461, 3893–3911. [Google Scholar] [CrossRef]
  245. Hone, A.N.W. The Associated Camassa-Holm Equation and the KdV Equation. J. Phys. A Math. Gen. 1999, 32, L307–L314. [Google Scholar] [CrossRef]
  246. Matsuno, Y. Cusp and Loop Soliton Solutions of Short-Wave Models for the Camassa–Holm and Degasperis–Procesi Equations. Phys. Lett. A 2006, 359, 451–457. [Google Scholar] [CrossRef]
  247. Miura, R.M. Conservation Laws for the Fully Nonlinear Long Wave Equations. Stud. Appl. Math. 1974, 53, 45–56. [Google Scholar] [CrossRef]
  248. Dye, J.M.; Parker, A. An Inverse Scattering Scheme for the Regularized Long-Wave Equation. J. Math. Phys. 2000, 41, 2889–2904. [Google Scholar] [CrossRef]
  249. Verheest, F.; Hereman, W. Conservations Laws and Solitary Wave Solutions for Generalized Schamel Equations. Phys. Scr. 1994, 50, 611–614. [Google Scholar] [CrossRef]
  250. Fordy, A.P.; Gibbons, J. Integrable Nonlinear Klein-Gordon Equations and Toda Lattices. Commun. Math. Phys. 1980, 77, 21–30. [Google Scholar] [CrossRef]
  251. Ablowitz, M.J.; Haberman, R. Nonlinear Evolution Equations-Two and Three Dimensions. Phys. Rev. Lett. 1975, 35, 1185–1188. [Google Scholar] [CrossRef]
  252. Matsuno, Y. Multisoliton Solutions of the Degasperis–Procesi Equation and Their Peakon Limit. Inverse Probl. 2005, 21, 1553–1570. [Google Scholar] [CrossRef]
  253. Constantin, A.; Ivanov, R.I.; Lenells, J. Inverse Scattering Transform for the Degasperis–Procesi Equation. Nonlinearity 2010, 23, 2559–2575. [Google Scholar] [CrossRef] [Green Version]
  254. Satsuma, J.; Ablowitz, M.J.; Kodama, Y. On An Internal Wave Equation Describing a Stratified Fluid With Finite Depth. Phys. Lett. A 1979, 73, 283–286. [Google Scholar] [CrossRef]
  255. Kako, F.; Mugibayashi, N. Complete Integrability of General Nonlinear Differential-Difference Equations Solvable by the Inverse Method. II. Prog. Theor. Phys. 1979, 61, 776–790. [Google Scholar] [CrossRef] [Green Version]
  256. Kaup, D. A Higher-order Water-Wave equation and the Method for Solving It. Prog. Theor. Phys. 1975, 54, 396–408. [Google Scholar] [CrossRef]
  257. Qiao, Z. New Integrable Hierarchy, Its Parametric Solutions, Cuspons, One-Peak Solitons, and M/W-Shape Peak Solitons. J. Math. Phys. 2007, 48, 082701. [Google Scholar] [CrossRef] [Green Version]
  258. Lamb, G.L., Jr. Analytical Descriptions of Ultrashort Optical Pulse Propagation in a Resonant Medium. Rev. Mod. Phys. 1971, 43, 99–124. [Google Scholar] [CrossRef]
  259. Haus, H.A.; Wong, W.S. Solitons in Optical Communications. Rev. Mod. Phys. 1996, 68, 423–444. [Google Scholar] [CrossRef]
  260. Moloney, J.V.; Newell, A.C. Nonlinear Optics. Phys. D Nonlinear Phenom. 1990, 44, 1–37. [Google Scholar] [CrossRef]
  261. Kumar, A. Soliton Dynamics in a Monomode Optical Fibre. Phys. Rep. 1990, 187, 63–108. [Google Scholar] [CrossRef]
  262. Wabnitz, S.; Kodama, Y.; Aceves, A.B. Control of Optical Soliton Interactions. Opt. FiberTechnol. 1995, 1, 187–217. [Google Scholar] [CrossRef]
  263. Hasegawa, A. An Historical Review of Application of Optical Solitons for High Speed Communications. Chaos Interdiscip. J. Nonlinear Sci. 2000, 10, 475–485. [Google Scholar] [CrossRef]
  264. Leblond, H. Half-Cycle Optical Soliton in Quadratic Nonlinear Media. Phys. Rev. A 2008, 78, 013807. [Google Scholar] [CrossRef] [Green Version]
  265. Zakharov, V.E.; Kuznetsov, E.A. Optical Solitons and Quasisolitons. J. Exp. Theor. Phys. 1998, 86, 1035–1046. [Google Scholar] [CrossRef]
  266. Hasegawa, A. Optical solitons in fibers. In Optical Solitons in Fibers; Hasegawa, A., Ed.; Springer: Berlin, Germany, 1989; pp. 1–74. [Google Scholar] [CrossRef]
  267. Kivshar, Y.S.; Agrawal, G. Optical Solitons: From Fibers to Photonic Crystals; Academic Press: Amsterdam, The Netherlands, 2003; ISBN 97800805380. [Google Scholar]
  268. Bullough, R.K.; Jack, P.M.; Kitchenside, P.W.; Saunders, R. Solitons in Laser Physics. Phys. Scr. 1979, 20, 364–381. [Google Scholar] [CrossRef]
  269. Haus, H.A. Optical Fiber Solitons, Their Properties and Uses. Proc. IEEE 1993, 81, 970–983. [Google Scholar] [CrossRef]
  270. Potasek, M.J. Novel Femtosecond Solitons in Optical Fibers, Photonic Switching, and Computing. J. Appl. Phys. 1989, 65, 941–953. [Google Scholar] [CrossRef]
  271. Ferreira, M.F.S. Solitons in Optical Fiber Systems; Wiley: Hoboken, NJ, USA, 2022; ISBN 9781119506676. [Google Scholar]
  272. Mollenauer, L.F.; Gordon, J.P. Solitons in Optical Fibers: Fundamentals and Applications; Elsevier: Amsterdam, The Netherlands, 2006; ISBN 9780080465067. [Google Scholar]
  273. Makhankov, V.G. Soliton Phenomenology; Kluwer: Dordrecht, The Netherlands, 1990; ISBN 9789401074940. [Google Scholar]
  274. Kodama, Y. Solitons in Two-Dimensional Shallow Water; SIAM: Philadelphia, PA, USA, 2018; ISBN 9781611975529. [Google Scholar]
  275. Song, Y.; Shi, X.; Wu, C.; Tang, D.; Zhang, H. Recent Progress of Study on Optical Solitons in Fiber Lasers. Appl. Phys. Rev. 2019, 6, 021313. [Google Scholar] [CrossRef]
  276. Lederer, F.; Stegeman, G.I.; Christodoulides, D.N.; Assanto, G.; Segev, M.; Silberberg, Y. Discrete Solitons in Optics. Phys. Rep. 2008, 463, 1–126. [Google Scholar] [CrossRef]
  277. Moloney, J.; Newell, A. Nonlinear Optics; CRC Press: Boca Raton, FL, USA, 2018; ISBN 978-0-8133-4118-7. [Google Scholar]
  278. Bishop, A.R. Solitons in Condensed Matter Physics. Phys. Scr. 1979, 20, 409–423. [Google Scholar] [CrossRef]
  279. Toda, M. Studies of a Non-linear Lattice. Phys. Rep. 1975, 18, 1–123. [Google Scholar] [CrossRef]
  280. Flaschka, H. The Toda Lattice. II. Existence of Integrals. Phys. Rev. B 1974, 9, 1924–1925. [Google Scholar] [CrossRef]
  281. Flaschka, H. On the Toda Lattice. II: Inverse-Scattering Solution. Prog. Theor. Phys. 1974, 51, 703–716. [Google Scholar] [CrossRef] [Green Version]
  282. Hirota, R.; Suzuki, K. Theoretical and Experimental Studies of Lattice Solitons in Nonlinear Lumped Networks. Proc. IEEE 1973, 61, 1483–1491. [Google Scholar] [CrossRef]
  283. Toda, M.; Wadati, M. A Soliton and Two Solitons in an Exponential Lattice and Related Equations. J. Phys. Soc. Jpn. 1973, 34, 18–25. [Google Scholar] [CrossRef]
  284. Wadati, M. Transformation Theories for Nonlinear Discrete Systems. Prog. Theor. Phys. Suppl. 1976, 59, 36–63. [Google Scholar] [CrossRef] [Green Version]
  285. Wilson, G. The Modified Lax and Two-Dimensional Toda Lattice Equations Associated with Simple Lie algebras. Ergod. Theory Dyn. Syst. 1981, 1, 361–380. [Google Scholar] [CrossRef] [Green Version]
  286. Kaup, D.J. The Forced Toda Lattice: An Example of an Almost Integrable System. J. Math. Phys. 1984, 25, 277–281. [Google Scholar] [CrossRef]
  287. Ueno, K.; Takasaki, K. Toda Lattice Hierarchy. Adv. Stud. Pure Math. 1984, 4, 1–95. [Google Scholar] [CrossRef]
  288. Ablowitz, M.J.; Ladik, J.F. A Nonlinear Difference Scheme and Inverse Scattering. Stud. Appl. Math. 1976, 55, 213–229. [Google Scholar] [CrossRef]
  289. Kac, M.; van Moerbeke, P. On an Explicitly Soluble System of Nonlinear Differential Equations Related to Certain Toda Lattices. Adv. Math. 1975, 16, 160–169. [Google Scholar] [CrossRef] [Green Version]
  290. Chen, H.H.; Liu, C.S. Bäcklund Transformation Solutions of the Toda Lattice Equation. J. Math. Phys. 1975, 16, 1428–1430. [Google Scholar] [CrossRef]
  291. Varma, C.M. Dynamics of Anharmonic Lattices: Solitons and the Central-Peak Problem in One Dimension. Phys. Rev. B 1976, 14, 244–254. [Google Scholar] [CrossRef]
  292. Hirota, R. Discrete Analogue of Generalized Toda Equation. J. Phys. Soc. Jpn. 1981, 50, 3785–3791. [Google Scholar] [CrossRef]
  293. Levi, D.; Pilloni, L.; Santini, P.M. Integrable Three-dimensional Lattices. J. Phys. A Math. Gen. 1981, 14, 1567–1575. [Google Scholar] [CrossRef]
  294. Date, E.; Jinbo, M.; Miwa, T. Method for Generating Discrete Soliton Equations. I. J. Phys. Soc. Jpn. 1982, 51, 4116–4124. [Google Scholar] [CrossRef]
  295. Toda, M. Nonlinear Lattice and Soliton Theory. IEEE Trans. Circuits Syst. 1983, 30, 542–554. [Google Scholar] [CrossRef]
  296. Wiersma, G.L.; Capel, H.W. Lattice Equations, Hierarchies and Hamiltonian Structures. Phys. A Stat. Mech. Its Appl. 1987, 142, 199–244. [Google Scholar] [CrossRef]
  297. Papageorgiou, V.G.; Nijhoff, F.W.; Capel, H.W. Integrable Mappings and Nonlinear Integrable Lattice Equations. Phys. Lett. A 1990, 147, 106–114. [Google Scholar] [CrossRef]
  298. Veselov, A.P. Growth and Integrability in the Dynamics of Mappings. Commun. Math. Phys. 1992, 145, 181–193. [Google Scholar] [CrossRef]
  299. Braun, O.M.; Kivshar, Y.S. Nonlinear Dynamics of the Frenkel–Kontorova Model. Phys. Rep. 1998, 306, 1–108. [Google Scholar] [CrossRef]
  300. Takhtadzhan, L.A.; Faddeev, L.D. The Quantum Method of the Inverse Problem and The Heisenberg XYZ Model. Russ. Math. Surv. 1979, 34, 11–68. [Google Scholar] [CrossRef]
  301. Balakrishnan, R. On the Inhomogeneous Heisenberg Chain. J. Phys. C Solid State Phys. 1982, 15, L1305–L1308. [Google Scholar] [CrossRef]
  302. Latha, M.M.; Vasanthi, C.C. An Integrable Model of (2 + 1)-Dimensional Heisenberg Ferromagnetic Spin Chain and Soliton Excitations. Phys. Scr. 2014, 89, 065204. [Google Scholar] [CrossRef]
  303. Porsezian, K.; Daniel, M.; Lakshmanan, M. On the Integrability Aspects of the One- Dimensional Classical Continuum Isotropic Biquadratic Heisenberg Spin Chain. J. Math. Phys. 1992, 33, 1807–1816. [Google Scholar] [CrossRef]
  304. Fogedby, H.C. Solitons and Magnons in the Classical Heisenberg Chain. J. Phys. Math. Gen. 1980, 13, 1467–1499. [Google Scholar] [CrossRef]
  305. Slavin, A.N.; Rojdestvenski, I.V. “Bright” and “Dark” Spin Wave Envelope Solitons in Magnetic Films. IEEE Trans. Magn. 1994, 30, 37–45. [Google Scholar] [CrossRef]
  306. Lam, L.; Prost, J. (Eds.) Solitons in Liquid Crystals; Springer: New York, NY, USA, 1992; ISBN 978-1-4612-6946-5. [Google Scholar]
  307. Maugin, G.A. Solitons in Elastic Solids (1938–2010). Mech. Res. Commun. 2011, 38, 341–349. [Google Scholar] [CrossRef]
  308. Erofeev, V.I.; Klyueva, N.V. Solitons and Nonlinear Periodic Strain Waves in Rods, Plates, and Shells (A Review). Acoust. Phys. 2002, 48, 643–655. [Google Scholar] [CrossRef]
  309. Slavin, A.N.; Büttner, O.; Bauer, M.; Demokritov, S.O.; Hillebr, S.B.; Kostylev, M.P.; Kalinikos, B.A.; Shevchenko, T.; Rapoport, Y. Collision Properties of Quasi-One-Dimensional Spin Wave Solitons and Two-Dimensional Spin Wave Bullets. Chaos Interdiscip. J. Nonlinear Sci. 2003, 13, 693–701. [Google Scholar] [CrossRef]
  310. Ichikawa, Y.H. Topics on Solitons in Plasmas. Phys. Scr. 1979, 20, 296–305. [Google Scholar] [CrossRef]
  311. Ikezi, H.; Taylor, R.J.; Baker, D.R. Formation and Interaction of Ion-Acoustic Solitions. Phys. Rev. Lett. 1970, 25, 11–14. [Google Scholar] [CrossRef]
  312. Tran, M.Q. Ion Acoustic Solitons in a Plasma: A Review of Their Experimental Properties and Related Theories. Phys. Scr. 1979, 20, 317–327. [Google Scholar] [CrossRef]
  313. Sakanaka, P.H. Formation and Interaction of Ion-Acoustic Solitary Waves in a Collisionless Warm Plasma. Phys. Fluids 1972, 15, 304–310. [Google Scholar] [CrossRef]
  314. Ablowitz, M.J.; Segur, H. On The Evolution of Packets of Water Waves. J. Fluid Mech. 1979, 92, 691–715. [Google Scholar] [CrossRef]
  315. Nakamura, Y. Experiments on Ion-acoustic Solitons in Plasmas. IEEE Trans. Plasma Sci. 1982, 10, 180–195. [Google Scholar] [CrossRef]
  316. Yajima, N.; Oikawa, M. Formation and Interaction of Sonic-Langmuir Solitons: Inverse Scattering Method. Prog. Theor. Phys. 1976, 56, 1719–1739. [Google Scholar] [CrossRef] [Green Version]
  317. Yuen, H.C.; Lake, B.M. Nonlinear Dynamics of Deep-water Gravity Waves. Adv. Appl. Mech. 1982, 22, 67–229. [Google Scholar] [CrossRef]
  318. Segur, H. The Korteweg-de Vries Equation and Water Waves. Solutions of The Equation. Part 1. J. Fluid Mech. 1973, 59, 721–736. [Google Scholar] [CrossRef]
  319. Hirota, R.; Satsuma, J. N-Soliton Solutions of Model Equations for Shallow Water Waves. J. Phys. Soc. Jpn. 1976, 40, 611–612. [Google Scholar] [CrossRef]
  320. Zeytounian, R.K. Nonlinear Long Waves on Water and Solitons. Sov. Phys.-Uspekhi 1995, 38, 1333–1382. [Google Scholar] [CrossRef]
  321. Tzirtzilakis, E.; Xenos, M.; Marinakis, V.; Bountis, T.C. Interactions and Stability of Solitary Waves in Shallow Water. Chaos Solitons Fractals 2002, 14, 87–95. [Google Scholar] [CrossRef]
  322. Geyer, A. Solitary Traveling Water Waves of Moderate Amplitude. J. Nonlinear Math. Phys. 2012, 19 (Suppl. S1), 1240010. [Google Scholar] [CrossRef] [Green Version]
  323. Constantin, A. Nonlinear Water Waves with Applications to Wave-Current Interactions and Tsunamis; SIAM: Philadelphia, PA, USA, 2011; ISBN 978-1-61344-434-4. [Google Scholar]
  324. Madsen, P.A.; Fuhrman, D.R.; Schäffer, H.A. On the Solitary Wave Paradigm for Tsunamis. J. Geophys. Res. Ocean. 2008, 113, C12012. [Google Scholar] [CrossRef]
  325. Osborne, A.R.; Burch, T.L. Internal Solitons in the Andaman Sea. Science 1980, 208, 451–460. [Google Scholar] [CrossRef] [PubMed]
  326. Helfrich, K.R.; Melville, W.K. Long Nonlinear Internal Waves. Annu. Rev. Fluid Mech. 2006, 38, 395–425. [Google Scholar] [CrossRef]
  327. Miles, J.W. On Internal Solitary Waves. Tellus 1979, 31, 456–462. [Google Scholar] [CrossRef] [Green Version]
  328. Christie, D.R.; Muirhead, K.J.; Hales, A.L. Intrusive Density Flows in the Lower Troposphere: A Source of Atmospheric Solitons. J. Geophys. Res. Ocean. 1979, 84, 4959–4970. [Google Scholar] [CrossRef]
  329. Vlachogiannis, M.; Bontozoglou, V. Observations of Solitary Wave Dynamics of Film Flows. J. Fluid Mech. 2001, 435, 191–215. [Google Scholar] [CrossRef]
  330. Ludu, A.; Draayer, J.P. Nonlinear Modes of Liquid Drops as Solitary Waves. Phys. Rev. Lett. 1998, 80, 2125–2128. [Google Scholar] [CrossRef]
  331. Ablowitz, M.J.; Baldwin, D.E. Nonlinear Shallow Ocean-Wave Soliton Interactions on Flat Beaches. Phys. Rev. E 2012, 86, 036305. [Google Scholar] [CrossRef] [Green Version]
  332. Ablowitz, M.J.; Segur, H. Long Internal Waves in Fluids of Great Depth. Stud. Appl. Math. 1980, 62, 249–262. [Google Scholar] [CrossRef]
  333. Gibbon, J.D.; James, I.N.; Moroz, I.M. An Example of Soliton Behaviour in a Rotating Baroclinic Fluid. Proc. R. Soc. Lond. A Math. Phys. Sci. 1979, 367, 219–237. [Google Scholar] [CrossRef]
  334. Osborne, A.R.; Segre, E.; Boffetta, G.; Cavaleri, L. Soliton Basis States in Shallow- Water Ocean Surface Waves. Phys. Rev. Lett. 1991, 67, 592–595. [Google Scholar] [CrossRef] [PubMed]
  335. Walker, L.R. Interfacial Solitary Waves in a Two-fluid Medium. Phys. Fluids 1973, 16, 1796–1804. [Google Scholar] [CrossRef]
  336. Bona, J.L.; Rajopadhye, S.V.; Schonbek, M.E. Models for Propagation of Bores. I. Two- Dimensional Theory. Differ. Integral Equ. 1994, 7, 699–734. [Google Scholar]
  337. Rajopadhye, S.V. Some Models for the Propagation of Bores. J. Differ. Equ. 2005, 217, 179–203. [Google Scholar] [CrossRef] [Green Version]
  338. Shrira, V.I.; Geogjaev, V.V. What Makes the Peregrine Soliton so Special as a Prototype of Freak Waves? J. Eng. Math. 2010, 67, 11–22. [Google Scholar] [CrossRef]
  339. Redekopp, L.G. On The Theory of Solitary Rossby Waves. J. Fluid Mech. 1977, 82, 725–745. [Google Scholar] [CrossRef]
  340. Slunyaev, A.; Kharif, C.; Pelinovsky, E.; Talipova, T. Nonlinear Wave Focusing on Water of Finite Depth. Phys. D Nonlinear Phenom. 2002, 173, 77–96. [Google Scholar] [CrossRef]
  341. Davydov, A.S. Solitons in Quasi-One-Dimensional Molecular Structures. Sov. Phys. Uspekhi 1982, 25, 603–613. [Google Scholar] [CrossRef]
  342. Hyman, J.M.; McLaughlin, D.W.; Scott, A.C. On Davydov’s Alpha-Helix Solitons. Phys. D Nonlinear Phenom. 1981, 3, 23–44. [Google Scholar] [CrossRef]
  343. Scott, A.C. Dynamics of Davydov Solitons. Phys. Rev. A 1982, 26, 578–595. [Google Scholar] [CrossRef]
  344. Yomosa, S. Solitary Waves in Large Blood Vessels. J. Phys. Soc. Jpn. 1987, 56, 506–520. [Google Scholar] [CrossRef]
  345. Purwins, H.G.; Bödeker, H.U.; Amiranashvili, S. Dissipative Solitons. Adv. Phys. 2010, 59, 485–701. [Google Scholar] [CrossRef]
  346. Yomosa, S. Solitary Excitations in Muscle Proteins. Phys. Rev. A 1985, 32, 1752–1758. [Google Scholar] [CrossRef] [PubMed]
  347. Orfanidis, S.J. Soliton Solutions of the Massive Thirring Model and the Inverse Scattering Transform. Phys. Rev. D 1976, 14, 472–478. [Google Scholar] [CrossRef]
  348. Kuznetsov, E.A.; Mikhailov, A.V. On the Complete Integrability of the Two-dimensional Classical Thirring Model. Theor. Math. Phys. 1977, 30, 193–200. [Google Scholar] [CrossRef]
  349. Matsutani, S.; Tsuru, H. Physical Relation Between Quantum Mechanics and Solitons on a Thin Elastic Rod. Phys. Rev. A 1992, 46, 1144–1147. [Google Scholar] [CrossRef]
  350. Kaup, D.J.; Newell, A.C. Solitons as Particles, Oscillators, and in Slowly Changing Media: A Singular Perturbation Theory. Proc. R. Soc. Lond. A Math. Phys. Sci. 1978, 361, 413–446. [Google Scholar] [CrossRef]
  351. Bogoyavlenskii, O.I. Breaking Solitons in 2+ 1-Dimensional Integrable Equations. Russ. Math. Surv. 1990, 45, 1–86. [Google Scholar] [CrossRef]
  352. Zakharov, V.E.; Kuznetsov, E.A. Solitons and Collapses: Two Evolution Scenarios of Nonlinear Wave Systems. Sov. Phys.-Uspekhi 2012, 55, 535–556. [Google Scholar] [CrossRef] [Green Version]
  353. Balakrishnan, R. Soliton Propagation in Nonuniform Media. Phys. Rev. A 1985, 32, 1144–1149. [Google Scholar] [CrossRef] [PubMed]
  354. Sanchez, A.; Vazquez, L. Nonlinear Wave Propagation in Disordered Media. Int. J. Mod. Phys. B 1991, 5, 2825–2882. [Google Scholar] [CrossRef] [Green Version]
  355. Kuznetsov, E.A.; Dias, F. Bifurcations of Solitons and Their Stability. Phys. Rep. 2011, 507, 43–105. [Google Scholar] [CrossRef]
  356. Bass, F.G.; Kivshar, Y.S.; Konotop, V.V.; Sinitsyn, Y.A. Dynamics of Solitons Under Random Perturbations. Phys. Rep. 1988, 157, 63–181. [Google Scholar] [CrossRef]
  357. Keener, J.P.; McLaughlin, D.W. Solitons Under Perturbations. Phys. Rev. A 1977, 16, 777–790. [Google Scholar] [CrossRef]
  358. Kodama, Y.; Ablowitz, M.J. Perturbations of Solitons and Solitary Waves. Stud. Appl. Math. 1981, 64, 225–245. [Google Scholar] [CrossRef]
  359. Malomed, B.A. Vortex Solitons: Old Results and New Perspectives. Phys. D Nonlinear Phenom. 2019, 399, 108–137. [Google Scholar] [CrossRef] [Green Version]
  360. Matveev, V.B. Positons: Slowly Decreasing Analogues of Solitons. Theor. Math. Phys. 2002, 131, 483–497. [Google Scholar] [CrossRef]
  361. Appert, K.; Vaclavik, J. Dynamics of Coupled Solitons. Phys. Fluids 1977, 20, 1845–1849. [Google Scholar] [CrossRef]
  362. Boiti, M.; Martina, L.; Pempinelli, F. Multidimensional Localized Solitons. Chaos Solitons Fractals 1995, 5, 2377–2417. [Google Scholar] [CrossRef] [Green Version]
  363. Hammack, J.; Henderson, D.; Guyenne, P.; Yi, M. Solitary-Wave Collisions. In Advances in Engineering Mechanics—Reflections and Outlooks: In Honor of Theodore YT Wu; Chwang, A.T., Teng, M.H., Valentine, D.T., Eds.; World Scientific: Singapore, 2005; pp. 173–194. [Google Scholar] [CrossRef] [Green Version]
  364. Lund, F.; Regge, T. Unified Approach to Strings and Vortices with Soliton Solutions. Phys. Rev. D 1976, 14, 1524–1535. [Google Scholar] [CrossRef]
  365. Benney, D.J. A General Theory for Interactions Between Short and Long Waves. Stud. Appl. Math. 1977, 56, 81–94. [Google Scholar] [CrossRef]
  366. Chu, F.Y.; Scott, A.C. Inverse Scattering Transform for Wave-Wave Scattering. Phys. Rev. A 1975, 12, 2060–2064. [Google Scholar] [CrossRef]
  367. Pohlmeyer, K. Integrable Hamiltonian Systems and Interactions through Quadratic Constraints. Commun. Math. Phys. 1976, 46, 207–221. [Google Scholar] [CrossRef]
  368. Olver, P.J.; Rosenau, P. Tri-Hamiltonian Duality Between Solitons and Solitary-Wave Solutions Having Compact Support. Phys. Rev. E 1996, 53, 1900–1906. [Google Scholar] [CrossRef]
  369. Gerdjikov, V.S.; Vilasi, G.; Yanovski, A.B. Integrable Hamiltonian Hierarchies: Spectral and Geometric Methods; Springer: Berlin, Germany, 2008; ISBN 978-3-540-77054-1. [Google Scholar]
  370. Reshetikhin, N.Y.; Faddeev, L.D. Hamiltonian Structures for Integrable Models of Field Theory. In Fifty Years of Mathematical Physics: Selected Works of Ludwig Faddeev; Ge, M., Niemi, A.J., Eds.; World Scientific: Singapore, 2016; pp. 323–338. ISBN 978-981-4340-95-3. [Google Scholar]
  371. Degasperis, A.; Procesi, M. Asymptotic Integrability. Symmetry Perturbation Theory 1999, 1, 23–37. [Google Scholar]
  372. Babelon, O.; Bernard, D.; Talon, M. Introduction to Classical Integrable Systems; Cambridge University Press: Cambridge, UK, 2003; ISBN 978-0-511-05571-3. [Google Scholar]
  373. Ablowitz, M.; Muslimani, Z.H. Integrable Nonlocal Nonlinear Equations. Stud. Appl. Math. 2016, 139, 7–59. [Google Scholar] [CrossRef] [Green Version]
  374. Wadati, M.; Olmedilla, E.; Akutsu, Y. Lax Pair for the One-Dimensional Hubbard Model. J. Phys. Soc. Jpn. 1987, 56, 1340–1347. [Google Scholar] [CrossRef]
  375. Ponce, G. Lax Pairs and Higher Order Models for Water Waves. J. Differ. Equ. 1993, 102, 360–381. [Google Scholar] [CrossRef]
  376. Ramani, A.; Grammaticos, B.; Satsuma, J. Integrability of Multidimensional Discrete Systems. Phys. Lett. A 1992, 169, 323–328. [Google Scholar] [CrossRef]
  377. Terng, C.L.; Uhlenbeck, K. Geometry of Solitons. Not. Am. Math. Soc. 2000, 47, 17–25. [Google Scholar]
  378. Terng, C.L. Soliton Equations and Differential Geometry. J. Differ. Geom. 1997, 45, 407–445. [Google Scholar] [CrossRef]
  379. Cavalcante, J.A.; Tenenblat, K. Conservation Laws for Nonlinear Evolution Equations. J. Math. Phys. 1988, 29, 1044–1049. [Google Scholar] [CrossRef]
  380. Kamchatnov, A.M. On Whitham Theory for Perturbed Integrable Equations. Phys. D Nonlinear Phenom. 2004, 188, 247–261. [Google Scholar] [CrossRef] [Green Version]
  381. Wadati, M.; Sanuki, H.; Konno, K. Relationships Among Inverse Method, Bäcklund Transformation and an Infinite Number of Conservation Laws. Prog. Theor. Phys. 1975, 53, 419–436. [Google Scholar] [CrossRef] [Green Version]
  382. Benjamin, T.B. The Stability of Solitary Waves. Proc. R. Soc. Lond. A Math. Phys. Sci. 1972, 328, 153–183. [Google Scholar] [CrossRef]
  383. Pego, R.L.; Weinstein, M.I. On Asymptotic Stability of Solitary Waves. Phys. Lett. A 1992, 162, 263–268. [Google Scholar] [CrossRef] [Green Version]
  384. Tao, T. Why Are Solitons Stable? Bull. Am. Math. Soc. 2009, 46, 1–33. [Google Scholar] [CrossRef] [Green Version]
  385. Wahlquist, H.D.; Estabrook, F.B. Prolongation Structures of Nonlinear Evolution Equations. J. Math. Phys. 1975, 16, 1–7. [Google Scholar] [CrossRef]
  386. Corones, J. Solitons and Simple Pseudopotentials. J. Math. Phys. 1976, 17, 756–759. [Google Scholar] [CrossRef]
  387. Olshanetsky, M.A.; Perelomov, A.M. Classical Integrable Finite-Dimensional Systems Related to Lie Algebras. Phys. Rep. 1981, 71, 313–400. [Google Scholar] [CrossRef]
  388. Krichever, I.M. Methods of Algebraic Geometry in the Theory of Non-linear Equations. Russ. Math. Surv. 1077, 32, 185–213. [Google Scholar] [CrossRef]
  389. Fokas, A.S. Symmetries and Integrability. Stud. Appl. Math. 1987, 77, 253–299. [Google Scholar] [CrossRef]
  390. Wadati, M.; Konno, K.; Ichikawa, Y.H. New Integrable Nonlinear Evolution Equations. J. Phys. Soc. Jpn. 1979, 47, 1698–1700. [Google Scholar] [CrossRef]
  391. Palais, R. The Symmetries of Solitons. Bull. Am. Math. Soc. 1997, 34, 339–403. [Google Scholar] [CrossRef] [Green Version]
  392. Clarkson, P.A.; Kruskal, M.D. New Similarity Reductions of the Boussinesq Equation. J. Math. Phys. 1989, 30, 2201–2213. [Google Scholar] [CrossRef]
  393. Ma, W.X.; Strampp, W. An Explicit Symmetry Constraint for the Lax Pairs and the Adjoint Lax Pairs of AKNS Systems. Phys. Lett. A 1994, 185, 277–286. [Google Scholar] [CrossRef]
  394. Correa, F.; Dunne, G.V.; Plyushchay, M.S. The Bogoliubov–de Gennes System, the AKNS Hierarchy, and Nonlinear Quantum Mechanical Supersymmetry. Ann. Phys. 2009, 324, 2522–2547. [Google Scholar] [CrossRef] [Green Version]
  395. Neugebauer, G.; Meinel, R. General N-Soliton Solution of the Class on Arbitrary Background. Phys. Lett. A 1984, 100, 467–470. [Google Scholar] [CrossRef]
  396. Flaschka, H.; Newell, A.C.; Ratiu, T. Kac-Moody Lie Algebras and Soliton Equations: II. Lax Equations Associated with A1 (1). Phys. D Nonlinear Phenom. 1983, 9, 300–323. [Google Scholar] [CrossRef]
  397. Ramani, A.; Grammaticos, B.; Bountis, T. The Painlevé Property and Singularity Analysis of Integrable and Non-integrable Systems. Phys. Rep. 1989, 180, 159–245. [Google Scholar] [CrossRef]
  398. Hone, A.N. Painlevé Tests, Singularity Structure and Integrability. In Integrability; Mikhailov, A.V., Ed.; Springer: Berlin, Germany, 2009; pp. 245–277. ISBN 978-3-540-88111-7. [Google Scholar]
  399. Conte, R.; Musette, M. The Painlevé Handbook; Springer: Cham, Switzerland, 2008; ISBN 978-3-030-53339-7. [Google Scholar]
  400. Newell, A.C.; Tabor, M.; Zeng, Y.B. A Unified Approach to Painlevé Expansions. Phys. D Nonlinear Phenom. 1987, 29, 1–68. [Google Scholar] [CrossRef]
  401. Kruskal, M.D.; Clarkson, P.A. The Painlevé-Kowalevski and Poly-Painlevé Tests for Integrability. Stud. Appl. Math. 1992, 86, 87–165. [Google Scholar] [CrossRef]
  402. Ablowitz, M.J.; Ramani, A.; Segur, H. A Connection Between Nonlinear Evolution Equations and Ordinary Differential Equations of P-Type. I. J. Math. Phys. 1980, 21, 715–721. [Google Scholar] [CrossRef]
  403. Ablowitz, M.J.; Ramani, A.; Segur, H. A Connection Between Nonlinear Evolution Equations and Ordinary Differential Equations of P-type. II. J. Math. Phys. 1980, 21, 1006–1015. [Google Scholar] [CrossRef]
  404. Karpman, V.I.; Maslov, E.M. Perturbation Theory for Solitons. Sov. Phys.-JETP 1977, 46, 281–291. [Google Scholar] [CrossRef]
  405. Ablowitz, M.J.; Satsuma, J. Solitons and Rational Solutions of Nonlinear Evolution Equations. J. Math. Phys. 1978, 19, 2180–2186. [Google Scholar] [CrossRef]
  406. Nakamura, A.A. Direct Method of Calculating Periodic Wave Solutions to Nonlinear Evolution Equations. I. Exact Two-periodic Wave Solution. J. Phys. Soc. Jpn. 1979, 47, 1701–1705. [Google Scholar] [CrossRef]
  407. Konopelchenko, B.; Strampp, W. The AKNS Hierarchy as Symmetry Constraint of the KP Hierarchy. Inverse Probl. 1991, 7, L17–L24. [Google Scholar] [CrossRef]
  408. Nijhoff, F.W.; Quispel, G.R.W.; Capel, H.W. Direct Linearization of Nonlinear Difference-Difference Equations. Phys. Lett. A 1983, 97, 125–128. [Google Scholar] [CrossRef]
  409. Wadati, M.; Sogo, K. Gauge Transformations in Soliton Theory. J. Phys. Soc. Jpn. 1983, 52, 394–398. [Google Scholar] [CrossRef]
  410. Chen, H.H. General Derivation of Bäcklund Transformations from Inverse Scattering Problems. Phys. Rev. Lett. 1974, 33, 925–928. [Google Scholar] [CrossRef]
  411. Lambert, F.; Springael, J. Soliton Equations and Simple Combinatorics. Acta Appl. Math. 2008, 102, 147–178. [Google Scholar] [CrossRef]
  412. Calogero, F.; Degasperis, A. Nonlinear Evolution Equations Solvable by the Inverse Spectral Transform.—I. Il Nuovo Cimento B 1976, 32, 201–242. [Google Scholar] [CrossRef]
  413. Gerdjikov, V.S. Generalised Fourier Transforms for The Soliton Equations. Gauge-Covariant formulation. Inverse Probl. 1986, 2, 51–74. [Google Scholar] [CrossRef]
  414. Pickering, A. The Singular Manifold Method Revisited. J. Math. Phys. 1996, 37, 1894–1927. [Google Scholar] [CrossRef]
  415. Sakhnovich, A. Generalized Backlund-Darboux Transformation: Spectral Properties and Nonlinear Equations. J. Math. Anal. Appl. 2001, 262, 274–306. [Google Scholar] [CrossRef] [Green Version]
  416. Satsuma, J. A Wronskian Representation of N-soliton Solutions of Nonlinear Evolution Equations. J. Phys. Soc. Jpn. 1979, 46, 359–360. [Google Scholar] [CrossRef]
  417. Segur, H.; Ablowitz, M.J. Asymptotic Solutions of Nonlinear Evolution Equations and a Painlevé Transcedent. Phys. D Nonlinear Phenom. 1981, 3, 165–184. [Google Scholar] [CrossRef]
  418. Makhankov, V.G. Dynamics of Classical Solitons (in Non-integrable Systems). Phys. Rep. 1978, 35, 1–128. [Google Scholar] [CrossRef]
  419. Fornberg, B.; Whitham, G.B. A Numerical and Theoretical Study of Certain Nonlinear Wave Phenomena. Philos. Trans. R. Soc. Lond. Ser. A Math. Phys. Sci. 1978, 289, 373–404. [Google Scholar] [CrossRef]
  420. Channell, P.J.; Scovel, C. Symplectic Integration of Hamiltonian systems. Nonlinearity 1990, 3, 231–259. [Google Scholar] [CrossRef]
  421. Vliegenthart, A.C. On Finite-difference Methods for the Korteweg-de Vries Equation. J. Eng. Math. 1971, 5, 137–155. [Google Scholar] [CrossRef]
  422. Argyris, J.; Haase, M. An Engineer’s Guide to Soliton Phenomena: Application of the Finite Element Method. Comput. Methods Appl. Mech. Eng. 1987, 61, 71–122. [Google Scholar] [CrossRef]
  423. Feng, B.F.; Mitsui, T. A Finite Difference Method for the Korteweg-de Vries and the Kadomtsev-Petviashvili Equations. J. Comput. Appl. Math. 1998, 90, 95–116. [Google Scholar] [CrossRef] [Green Version]
  424. Argyris, J.; Haase, M.; Heinrich, J.C. Finite Element Approximation to Two-Dimensional Sine-Gordon Solitons. Comput. Methods Appl. Mech. Eng. 1991, 86, 1–26. [Google Scholar] [CrossRef]
  425. Alexander, M.E.; Morris, J.L. Galerkin Methods Applied to Some Model Equations for Non- linear Dispersive Waves. J. Comput. Phys. 1979, 30, 428–451. [Google Scholar] [CrossRef]
  426. Griffiths, D.F.; Mitchell, A.R.; Morris, J.L. A Numerical Study of the Nonlinear Schrödinger equation. Comput. Methods Appl. Mech. Eng. 1984, 45, 177–215. [Google Scholar] [CrossRef]
  427. Bhatta, D.D.; Bhatti, M.I. Numerical Solution of KdV Equation Using Modified Bernstein Polynomials. Appl. Math. Comput. 2006, 174, 1255–1268. [Google Scholar] [CrossRef]
  428. Bratsos, A.G.; Tsitouras, C.; Natsis, D.G. Linearized Numerical Schemes for the Boussinesq Equation. Appl. Numer. Anal. Comput. Math. 2005, 2, 34–53. [Google Scholar] [CrossRef]
  429. Liu, H.; Yan, J. A Local Discontinuous Galerkin Method for the Korteweg–de Vries Equation With Boundary Effect. J. Comput. Phys. 2006, 215, 197–218. [Google Scholar] [CrossRef]
  430. Dehghan, M.; Shokri, A. A Numerical Method for KdV Equation Using Collocation and Radial Basis Functions. Nonlinear Dyn. 2007, 50, 111–120. [Google Scholar] [CrossRef]
  431. Liu, P.L.F.; Cheng, Y. A Numerical Study of the Evolution of a Solitary Wave over a Shelf. Phys. Fluids 2001, 13, 1660–1667. [Google Scholar] [CrossRef]
  432. Bridgman, T.; Hereman, W.; Quispel, G.R.W.; van der Kamp, P.H. Symbolic Computation of Lax Pairs of Partial Difference Equations Using Consistency Around the Cube. Found. Comput. Math. 2013, 13, 517–544. [Google Scholar] [CrossRef] [Green Version]
  433. Gordoa, P.R.; Pickering, A. Nonisospectral Scattering Problems: A Key to Integrable Hierarchies. J. Math. Phys. 1999, 40, 5749–5786. [Google Scholar] [CrossRef]
  434. Herbst, B.M.; Morris, J.L.; Mitchell, A.R. Numerical Experience With the Nonlinear Schrödinger Equation. J. Comput. Phys. 1985, 60, 282–305. [Google Scholar] [CrossRef]
  435. Meiss, J.D.; Pereira, N.R. Internal Wave Solitons. Phys. Fluids 1978, 21, 700–702. [Google Scholar] [CrossRef]
  436. Ponce Dawson, S.; Ferro Fontan, C. Soliton Decay of Nonlinear Alfven waves: Numerical Studies. Phys. Fluids 1988, 31, 83–89. [Google Scholar] [CrossRef]
  437. Hopf, E. The Partial Differential Equation: ut + uux = ϵuxx. Commun. Pure Appl. Math. 1950, 3, 201–230. [Google Scholar] [CrossRef]
  438. Cole, J.D. On a Quasi-Linear Parabolic Equation Occurring in Aerodynamics. Q. Appl. Math. 1951, 9, 225–236. [Google Scholar] [CrossRef]
  439. Ablowitz, M.J.; Clarkson, P.A. Solitons, Nonlinear Evolution Equations and Inverse Scattering; Cambridge University Press: Cambridge, UK, 1991; ISBN 978-0511623998. [Google Scholar]
  440. Tabor, M. Chaos and Integrability in Dynamical Systems; Wiley: New York, NY, USA, 1989; ISBN 978-0471827283. [Google Scholar]
  441. Carrielo, F.; Tabor, M. Similarity Reductions from Extended Painlevé Expansions for Nonintegrable Evolution Equations. Phys. D 1991, 53, 59–70. [Google Scholar] [CrossRef]
  442. Carrielo, F.; Tabor, M. Painlev’e Expansions for Nonintegrable Evolution Equations. Phys. D 1989, 39, 77–94. [Google Scholar] [CrossRef]
  443. Weiss, J.; Tabor, M. Carnevalle, G. The Painlevé Property for Partial Differential Equations. J. Math. Phys. 1983, 24, 522–526. [Google Scholar] [CrossRef]
  444. Kudryashov, N.A. On Types of Nonlinear Nonintegrable Equations with Exact Solutions. Phys. Lett. A 1991, 155, 269–275. [Google Scholar] [CrossRef]
  445. Kudryashov, N.A. Simplest Equation Method to Look for Exact Solutions of Nonlinear Differential Equations. Chaos Solitons Fractals 2005, 24, 1217–1231. [Google Scholar] [CrossRef] [Green Version]
  446. Kudryashov, N.A.; Loguinova, N.B. Extended Simplest Equation Method for Nonlinear Differential Equations. Appl. Math. Comput. 2008, 205, 361–365. [Google Scholar] [CrossRef]
  447. Kudryashov, N.A. Partial Differential Equations with Solutions Having Movable First-Order Singularities. Phys. Lett. A 1992, 169, 237–242. [Google Scholar] [CrossRef]
  448. Kudryashov, N.A. Exact Solitary Waves of the Fisher Equation. Phys. Lett. A 2005, 342, 99–106. [Google Scholar] [CrossRef]
  449. Kudryashov, N.A. One Method for Finding Exact Solutions of Nonlinear Differential Equations. Commun. Nonlinear Sci. Numer. Simul. 2012, 17, 2248–2253. [Google Scholar] [CrossRef] [Green Version]
  450. Kudryashov, N.A. Exact Soliton Solutions of the Generalized Evolution Equation of Wave Dynamics. J. Appl. Math. Mech. 1988, 52, 361–365. [Google Scholar] [CrossRef]
  451. Kudryashov, N.A. Exact Solutions of Nonlinear Wave Equations Arising in Mechanics. J. Appl. Math. Mech. 1990, 54, 372–375. [Google Scholar] [CrossRef]
  452. Kudryashov, N.A. Exact Solutions and Integrability of the Duffing–Van der Pol Equation. Regul. Chaotic Dyn. 2018, 23, 471–479. [Google Scholar] [CrossRef]
  453. Kudryashov, N.A. Exact Solutions of the Equation for Surface waves in a Convecting Fluid. Appl. Math. Comput. 2019, 344, 97–106. [Google Scholar] [CrossRef]
  454. Kudryashov, N.A. A Generalized Model for Description of Propagation Pulses in Optical Fiber. Optik 2019, 189, 42–52. [Google Scholar] [CrossRef]
  455. Kudryashov, N.A. First Integrals and Solutions of the Traveling Wave Reduction for the Triki–Biswas Equation. Optik 2019, 185, 275–281. [Google Scholar] [CrossRef]
  456. Kudryashov, N.A. Highly Dispersive Optical Solitons of the Generalized Nonlinear Eighth-Order Schrödinger Equation. Optik 2020, 206, 164335. [Google Scholar] [CrossRef]
  457. Kudryashov, N.A. The Generalized Duffing Oscillator. Commun. Nonlinear Sci. Numer. Simul. 2021, 93, 105526. [Google Scholar] [CrossRef]
  458. Urbain, F.; Kudryashov, N.A.; Tala-Tebue, E.; Hubert, M.B.; Doka, S.Y.; Crepin, K.T. Exact Solutions of the KdV Equation with Dual-Power Law Nonlinearity. Comput. Math. Math. Phys. 2021, 61, 431–435. [Google Scholar] [CrossRef]
  459. Kudryashov, N.A. Solitary waves of the generalized Sasa-Satsuma equation with arbitrary refractive index. Optik 2021, 232, 166540. [Google Scholar] [CrossRef]
  460. Vitanov, N.K.; Dimitrova, Z.I.; Vitanov, K.N. Simple Equations Method (SEsM): Algorithm, Connection with Hirota Method, Inverse Scattering Transform Method, and Several Other Methods. Entropy 2021, 23, 10. [Google Scholar] [CrossRef] [PubMed]
  461. Vitanov, N.K. Recent Developments of the Methodology of the Modified Method of Simplest Equation with Application. Pliska Stud. Math. Bulg. 2019, 30, 29–42. [Google Scholar]
  462. Vitanov, N.K. Modified Method of Simplest Equation for Obtaining Exact Solutions of Nonlinear Partial Differential Equations: History, recent development and studied classes of equations. J. Theor. Appl. Mech. 2019, 49, 107–122. [Google Scholar] [CrossRef]
  463. Vitanov, N.K. The Simple Equations Method (SEsM) For Obtaining Exact Solutions of Nonlinear PDEs: Opportunities Connected to the Exponential Functions. AIP Conf. Proc. 2019, 2159, 030038. [Google Scholar] [CrossRef]
  464. Vitanov, N.K. Simple Equations Method (SEsM): Review and New Results. AIP Conf. Ser. 2022, 2459, 020003. [Google Scholar] [CrossRef]
  465. Vitanov, N.K. Simple Equations Method (SEsM) and Nonlinear PDEs with Fractional Derivatives. AIP Conf. Ser. 2022, 2459, 030040. [Google Scholar] [CrossRef]
  466. Vitanov, N.K.; Dimitrova, Z.I. Simple Equations Method (SEsM) and Other Direct Methods for Obtaining Exact Solutions of Nonlinear PDEs. AIP Conf. Proc. 2019, 2159, 030039. [Google Scholar] [CrossRef]
  467. Vitanov, N.K. Schrödinger Equation and Nonlinear Waves. In Understanding the Schrödinger Equation; Simpao, V., Little, H., Eds.; Nova Science Publishers: New York, NY, USA, 2020; pp. 37–92. ISBN 978-1-53617-662-9. [Google Scholar]
  468. Dimitrova, Z.I.; Vitanov, N.K. Travelling Waves Connected to Blood Flow and Motion of Arterial Walls. In Water in Biomechanical and Related Systems; Gadomski, A., Ed.; Springer: Cham, Switzerland, 2021; pp. 243–263. ISBN 978-3-030-67226-3. [Google Scholar]
  469. Vitanov, N.K.; Dimitrova, Z.I. Simple Equations Method and Non-linear Differential Equations with Non-polynomial Non-linearity. Entropy 2021, 23, 1624. [Google Scholar] [CrossRef]
  470. Vitanov, N.K.; Dimitrova, Z.I.; Vitanov, K.N. On the Use of Composite Functions in the Simple Equations Method to Obtain Exact Solutions of Nonlinear Differential Equations. Computation 2021, 9, 104. [Google Scholar] [CrossRef]
  471. Vitanov, N.K. Simple Equations Method (SEsM) for Obtaining Exact Solutions of Nonlinear Differential Equations. In Advances in Mathematics Research; Baswell, A.R., Ed.; Nova Science Publishers: New York, NY, USA, 2022; pp. 105–138. ISBN 978-1-68507-892-8. [Google Scholar]
  472. Martinov, N.; Vitanov, N. On the Correspondence Between the Self-consistent 2D Poisson-Boltzmann Structures and the Sine-Gordon Waves. J. Phys. A Math. Gen. 1992, 25, L51–L56. [Google Scholar] [CrossRef]
  473. Martinov, N. Vitanov. On Some Solutions of the Two-Dimensional Sine-Gordon Equation. J. Phys. A Math. Gen. 1992, 25, L419–L426. [Google Scholar] [CrossRef]
  474. Vitanov, N.K. On Travelling Waves and Double-Periodic Structures in Two-Dimensional Sine–Gordon Systems. J. Phys. A Math. Gen. 1996, 29, 5195–5207. [Google Scholar] [CrossRef]
  475. Vitanov, N.K. Breather and Soliton Wave Families for the Sine-Gordon Equation. Proc. Roy. Soc. Lond. A 1998, 454, 2409–2423. [Google Scholar] [CrossRef]
  476. Vitanov, N.K.; Jordanov, I.P.; Dimitrova, Z.I. On Nonlinear Dynamics of Interacting Populations: Coupled Kink Waves in a System of Two Populations. Commun. Nonlinear Sci. Numer. Simul. 2009, 14, 2379–2388. [Google Scholar] [CrossRef]
  477. Vitanov, N.K.; Jordanov, I.P.; Dimitrova, Z.I. On Nonlinear Population Waves. Appl. Math. Comput. 2009, 215, 2950–2964. [Google Scholar] [CrossRef]
  478. Vitanov, N.K. Application of Simplest Equations of Bernoulli and Riccati Kind for Obtaining Exact Traveling-Wave Solutions for a Class of PDEs with Polynomial Nonlinearity. Commun. Nonlinear Sci. Numer. Simul. 2010, 15, 2050–2060. [Google Scholar] [CrossRef]
  479. Vitanov, N.K.; Dimitrova, Z.I. Application of The Method of Simplest Equation for Obtaining Exact Traveling-Wave Solutions for Two Classes of Model PDEs from Ecology and Population Dynamics. Commun. Nonlinear Sci. Numer. Simul. 2010, 15, 2836–2845. [Google Scholar] [CrossRef]
  480. Vitanov, N.K.; Dimitrova, Z.I.; Kantz, H. Modified Method of Simplest Equation and its Application to Nonlinear PDEs. Appl. Math. Comput. 2010, 216, 2587–2595. [Google Scholar] [CrossRef]
  481. Vitanov, N.K. Modified Method of Simplest Equation: Powerful Tool for Obtaining Exact and Approximate Traveling-Wave Solutions of Nonlinear PDEs. Commun. Nonlinear Sci. Numer. Simul. 2011, 16, 1176–1185. [Google Scholar] [CrossRef]
  482. Vitanov, N.K.; Dimitrova, Z.I.; Vitanov, K.N. On the Class of Nonlinear PDEs That Can be Treated by the Modified Method of Simplest Equation. Application to Generalized Degasperis–Processi Equation and B-Equation. Commun. Nonlinear Sci. Numer. Simul. 2011, 16, 3033–3044. [Google Scholar] [CrossRef]
  483. Vitanov, N.K. On Modified Method of Simplest Equation for Obtaining Exact and Approximate Solutions of Nonlinear PDEs: The Role of the Simplest Equation. Commun. Nonlinear Sci. Numer. Simul. 2011, 16, 4215–4231. [Google Scholar] [CrossRef]
  484. Vitanov, N.K.; Dimitrova, Z.I.; Kantz, H. Application of the Method of Simplest Equation for Obtaining Exact Traveling-Wave Solutions for the Extended Korteweg–de Vries Equation and Generalized Camassa–Holm Equation. Appl. Math. Comput. 2013, 219, 7480–7492. [Google Scholar] [CrossRef]
  485. Vitanov, N.K.; Dimitrova, Z.I. Solitary Wave Solutions for Nonlinear Partial Differential Equations that Contain Monomials of Odd and Even Grades with Respect to Participating Derivatives. Appl. Math. Comput. 2014, 247, 213–217. [Google Scholar] [CrossRef]
  486. Vitanov, N.K.; Dimitrova, Z.I.; Vitanov, K.N. Modified Method of Simplest Equation for Obtaining Exact Analytical Solutions of Nonlinear Partial Differential Equations: Further Development of the Methodology with Applications. Appl. Math. Comput. 2015, 269, 363–378. [Google Scholar] [CrossRef] [Green Version]
  487. Vitanov, N.K.; Dimitrova, Z.I.; Ivanova, T.I. On Solitary Wave Solutions of a Class of Nonlinear Partial Differential Equations Based on the Function 1/cosh(αx+βt)n. Appl. Math. Comput. 2017, 315, 372–380. [Google Scholar] [CrossRef]
  488. Vitanov, N.K.; Dimitrova, Z.I. Modified Method of Simplest Equation Applied to the Nonlinear Schrödinger Equation. J. Theor. Appl. Mech. Sofia 2018, 48, 59–68. [Google Scholar] [CrossRef] [Green Version]
  489. Vitanov, N.K. Simple Equations Method (SEsM) and Its Connection with the Inverse Scattering Transform Method. AIP Conf. Proceedibgs 2021, 2321, 030035. [Google Scholar] [CrossRef]
  490. Nikolova, E.V.; Jordanov, I.P.; Dimitrova, Z.I.; Vitanov, N.K. Evolution of non-linear Waves in a Blood-Filled Artery with an Aneurysm. AIP Conf. Proc. 2017, 1895, 07002. [Google Scholar] [CrossRef] [Green Version]
  491. Jordanov, I.P.; Vitanov, N.K. On the Exact Traveling Wave Solutions of a Hyperbolic Reaction- Diffusion Equation. Stud. Comput. Intell. 2019, 793, 199–210. [Google Scholar] [CrossRef]
  492. Nikolova, E.V.; Chilikova-Lubomirova, M.; Vitanov, N.K. Exact Solutions of a Fifth-Order Korteweg–de Vries–type Equation Modeling non-linear Long Waves in Several Natural Phenomena. AIP Conf. 2021, 2321, 030026. [Google Scholar] [CrossRef]
  493. Vitanov, N.K.; Vitanov, K.N. Population Dynamics in Presence of State Dependent Fluctuations. Comput. Math. Appl. 2013, 68, 962–971. [Google Scholar] [CrossRef]
  494. Vitanov, N.K.; Dimitrova, Z.I. Simple Equations Method (SEsM) and Its Particular Cases: Hirota Method. AIP Conf. Proc. 2021, 2321, 030036. [Google Scholar] [CrossRef]
  495. Dimitrova, Z.I.; Vitanov, K.N. Homogeneous Balance Method and Auxiliary Equation Method as Particular Cases of Simple Equations Method (SEsM). AIP Conf. Proc. 2021, 2321, 030004. [Google Scholar] [CrossRef]
  496. Constantine, G.M.; Savits, T.H. A Multivariate Faa di Bruno Formula with Applications. Trans. Am. Math. Soc. 1996, 348, 503–520. [Google Scholar] [CrossRef]
  497. Hereman Zhuang, W.W. Symbolic Computation of Solitons via Hirota’s Bilinear Method; Department of Mathematical and Computer Sciences Colorado School of Mines: Golden, CO, USA, 1994; preprint. [Google Scholar]
  498. Rosales, R.R. Exact Solution of Some Nonlinear Evolution Equations. Stud. Appl. Math. 1978, 59, 117–151. [Google Scholar] [CrossRef]
  499. Rosales, R.R. Exact Solution of Some Nonlinear Evolution Equations. Ph.D. Thesis, California Institute of Technology, Pasadena, CA, USA, 1977. [Google Scholar]
  500. Zakharov, V.; Shabat, A. Exact Theory of Two-dimensional Self-focusing and One-dimensional Self-modulation of Waves in Nonlinear Media. J. Exp. Theor. Phys. 1971, 61, 118–134. [Google Scholar]
  501. Kermack, W.O.; McKendrick, A.G. A Contribution to the Mathematical Theory of Epidemics. Proc. R. Soc. Lond. Ser. A 1927, 115, 700–721. [Google Scholar] [CrossRef] [Green Version]
  502. Ablowitz, M.J.; Zeppetella, A. Explicit Solutions of Fisher’s Equation for a Special Wave Speed. Bull. Math. Biol. 1979, 41, 835–840. [Google Scholar] [CrossRef]
  503. Malfliet, W. Solitary Wave Solutions of Nonlinear Wave Equations. Am. J. Phys. 1992, 60, 650–654. [Google Scholar] [CrossRef]
  504. Broadbridge, P.; Bradshaw, B.H.; Fulford, G.R.; Aldis, G.K. Huxley and Fisher Equations for Gene Propagation: An Exact Solution. ANZIAM J. 2002, 44, 11–20. [Google Scholar] [CrossRef] [Green Version]
  505. Scott, A.C. Nonlinear Science. Emergence and Dynamics of Coherent Structures; Oxford University Press: Oxford, UK, 1999; ISBN 9780198528524. [Google Scholar]
  506. Camassa, R.; Holm, D.D. An Integrable Shallow Water Equation with Peaked Solitons. Phys. Rev. Lett. 1993, 71, 1661–1664. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  507. Johnson, R.S. The Classical Problem of Water Waves: A Reservoir of Integrable and Nearly Integrable Equations. J. Nonlinear Math. Phys. 2003, 10 (Suppl. S1), 72–92. [Google Scholar] [CrossRef] [Green Version]
  508. Nikolova, E.V. Exact Travelling-Wave Solutions of the Extended Fifth-Order Korteweg–deVries Equation via Simple Equations Method (SEsM): The Case of Two Simple Equations. Entropy 2022, 24, 1288. [Google Scholar] [CrossRef]
  509. Vitanov, N.; Dimitrova, Z. On Waves and Distributions in Population Dynamics. Biomath 2012, 1, 1209253. [Google Scholar] [CrossRef] [Green Version]
  510. Dimitrova, Z.I. Several Examples of Application of the Simple Equations Method (SEsM) for Obtaining Exact Solutions of Nonlinear PDEs. AIP Conf. Proc. 2022, 2459, 030005. [Google Scholar] [CrossRef]
  511. Dimitrova, Z.I. On Several Specific Cases of the Simple Equations Method (SEsM): Jacobi Elliptic Function Expansion method, F-Expansion Method, Modified Simple Equation Method, Trial Function Method, General Projective Riccati Equations Method, and First Intergal Method. AIP Conf. Proc. 2022, 2459, 030006. [Google Scholar] [CrossRef]
  512. Liu, S.; Fu, Z.; Liu, S.; Zhao, Q. Jacobi elliptic function expansion method and periodic wave solutions of nonlinear wave equations. Phys. Lett. A 2001, 289, 69–74. [Google Scholar] [CrossRef]
  513. Wang, M.; Zhou, Y. The Periodic Wave Solutions for the Klein–Gordon–Schrödinger Equations. Phys. Lett. A 2003, 318, 84–92. [Google Scholar] [CrossRef]
  514. Wang, M.; Li, X. Applications of F-expansion to Periodic Wave Solutions for a new Hamiltonian Amplitude Equation. Chaos Solitons Fractals 2005, 24, 1257–1268. [Google Scholar] [CrossRef]
  515. Ren, Y.-J.; Zhang, K.-Q. A Generalized F-expansion Method to find Abundant Families of Jacobi Elliptic Function solutions of the (2 + 1)-Dimensional Nizhnik–Novikov–Veselov Equation. Chaos Solitons Fractals 2006, 27, 959–979. [Google Scholar] [CrossRef]
  516. Jawad, A.J.M.; Petkovic, M.D.; Biswas, A. Modified Simple Equation Method for Nonlinear Evolution Equations. Appl. Math. Comput. 2010, 217, 869–877. [Google Scholar] [CrossRef]
  517. Liu, S.; Fu, Z.; Liu, S.-D.; Zhao, Q. A Simple Fast Method in Finding Particular Solutions of Some Nonlinear PDE. Appl. Math. Mech. 2001, 22, 326–331. [Google Scholar] [CrossRef]
  518. Xie, Y.; Tang, J. A Unified Trial Function Method in Finding the Explicit and Exact Solutions to Three NPDEs. Phys. Scr. 2006, 74, 197–200. [Google Scholar] [CrossRef]
  519. Li, B.; Chen, Y. General Projective Riccati Equation Method and Exact Solutions for Generalized KdV-type and KdV–Burgers-type Equations with Nonlinear Terms of Any Order. Chaos Solitons Fractals 2004, 19, 977–984. [Google Scholar] [CrossRef]
  520. Feng, Z.S. The First Integer Method to Study the Burgers-Korteweg-de Vries Equation. J. Phys. A 2002, 35, 343–349. [Google Scholar] [CrossRef]
  521. Wang, M.-L. Solitary Wave Solutions for Variant Boussinesq equations. Phys. Lett. A 1995, 199, 169–172. [Google Scholar] [CrossRef]
  522. Wang, M.-L. Exact Solutions for a Compound KdV-Burgers Equation. Phys. Lett. A 1996, 213, 279–287. [Google Scholar] [CrossRef]
  523. Wang, M.-L.; Zhou, Y.; Li, Z. Application of a Homogeneous Balance Method to Exact Solutions of Nonlinear Equations in Mathematical Physics. Phys. Lett. A 1996, 216, 67–75. [Google Scholar] [CrossRef]
  524. Fan, E.; Zhang, H. A note on the homogeneous balance method. Phys. Lett. A 1998, 246, 403–406. [Google Scholar] [CrossRef]
  525. Sirendaoreji; Jiong, S. Auxiliary Equation Method for Solving Nonlinear Partial Differential Equations. Phys. Lett. A 2003, 309, 387–396. [Google Scholar] [CrossRef]
  526. Dimitrova, Z.I. Relation Between G’/G-expansion Method and the Modified Method of Simplest Equation. Comptes Rendus L’Academie Bulg. Des Sci. 2012, 65, 1513–1520. [Google Scholar]
  527. Wang, M.L.; Li, X.Z.; Zhang, J.L. The (G’/G)-Expansion Method and Traveling Wave Solutions of Nonlinear Evolution Equations in Mathematical Physics. Phys. Lett. A 2008, 372, 417–423. [Google Scholar] [CrossRef]
  528. He, J.H.; Wu, X.H. Exp-Function Method for Nonlinear Wave Equations. Chaos Solitons Fractals 2006, 30, 700–708. [Google Scholar] [CrossRef]
  529. Malfliet, W.; Hereman, W. The tahn method I: Exact solutions of nonlinear evolution and wave equations. Phys. Scr. 1996, 54, 563–568. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Vitanov, N.K. Simple Equations Method (SEsM): An Effective Algorithm for Obtaining Exact Solutions of Nonlinear Differential Equations. Entropy 2022, 24, 1653. https://doi.org/10.3390/e24111653

AMA Style

Vitanov NK. Simple Equations Method (SEsM): An Effective Algorithm for Obtaining Exact Solutions of Nonlinear Differential Equations. Entropy. 2022; 24(11):1653. https://doi.org/10.3390/e24111653

Chicago/Turabian Style

Vitanov, Nikolay K. 2022. "Simple Equations Method (SEsM): An Effective Algorithm for Obtaining Exact Solutions of Nonlinear Differential Equations" Entropy 24, no. 11: 1653. https://doi.org/10.3390/e24111653

APA Style

Vitanov, N. K. (2022). Simple Equations Method (SEsM): An Effective Algorithm for Obtaining Exact Solutions of Nonlinear Differential Equations. Entropy, 24(11), 1653. https://doi.org/10.3390/e24111653

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

Article Metrics

Back to TopTop