---
title: Simple Equations Method (SEsM) Overview
url: https://www.emergentmind.com/topics/simple-equations-method-sesm
type: topic
---

# Simple Equations Method (SEsM) Overview

Simple Equations Method (SEsM) is a general methodology for constructing exact, and in principle approximate, solutions of nonlinear partial differential equations by representing the sought solution as a function of solutions of one or several simpler auxiliary differential equations, called simple equations. Its central operation is to replace direct treatment of a nonlinear PDE with a structured cascade: an optional transformation of the dependent variable, an ansatz in auxiliary functions, reduction through simple ODEs or PDEs, balance relations, and solution of a resulting nonlinear algebraic system. In the literature cited here, SEsM is presented both as a direct solution method and as a unifying framework that contains many named expansion, Riccati-based, elliptic-function, exponential, and Fourier constructions as particular cases [1908.07459], [1908.01075], [1909.00330], [2411.07333], [2504.16660].

## 1. Conceptual basis

SEsM starts from a nonlinear PDE, schematically
\[
DE(u,\ldots)=0,
\]
with \(u=u(x,\ldots,t)\), and seeks a representation of the form
\[
u(x,\ldots,t)=T\big(F(x,\ldots,t)\big),
\]
or, in multi-component settings,
\[
u_i(x,\ldots,t)=T_i\big[F_i(x,\ldots,t),G_i(x,\ldots,t),\dots\big].
\]
The purpose of the transformation \(T\) is to simplify the nonlinearity, to make it polynomial, or, in favorable cases, to remove it. The cited papers list Painlevé-type expansions, \(u=4\arctan F\) for sine-Gordon, and \(u=4\tanh^{-1}F\) for sinh-Gordon or Poisson–Boltzmann type equations among the admissible transformations [1908.07459], [2504.16660].

The transformed function \(F\) is then expressed in auxiliary functions \(f_1,\dots,f_N\), typically through a multivariate polynomial-type ansatz,
\[
\begin{aligned}
F &= a + \sum_{i_1=1}^N B_{i_1}f_{i_1}
+ \sum_{i_1=1}^N\sum_{i_2=1}^N Y_{i_1,i_2}f_{i_1}f_{i_2}
+ \cdots \\
&\quad + \sum_{i_1=1}^N\cdots\sum_{i_N=1}^N
\Theta_{i_1,\ldots,i_N}f_{i_1}\cdots f_{i_N},
\end{aligned}
\]
although rational forms, finite series, and more general composite expressions are also explicitly allowed. Each auxiliary function is linked to a simple equation, meaning an ODE or PDE that is more simple than the original nonlinear PDE and whose solutions and algebraic or differential properties are known or manageable [1908.07459].

Within this framework, “simple” is relative rather than absolute. The admissible class includes first-order linear equations, Riccati and Bernoulli equations, equations generating trigonometric, hyperbolic, Jacobi elliptic, and Weierstrass elliptic functions, as well as, in the most general formulation, PDEs, ODEs with variable coefficients, stochastic differential equations, and equations with fractional derivatives. This breadth is one reason SEsM is described as a methodology rather than a single fixed ansatz [1908.07459].

## 2. Algorithmic structure

The literature presents SEsM in two closely related organizational forms: a detailed seven-step scheme and a more compressed four-step algorithm. Both descriptions encode the same core logic: transform, represent, choose simple equations, substitute, balance, and solve [1908.07459], [2504.16660].

The first stage is the optional transformation of the dependent variable. When the original PDE already has a tractable polynomial nonlinearity, this stage may be skipped and one sets \(u=F\). Otherwise, a nontrivial \(T\) is used to alter the nonlinearity or bilinearize the problem. The transformed PDE is then written in terms of \(F\), or of several transformed functions \(F_i,G_i,\dots\) in the multi-component case [1908.07459], [2411.07333].

The second stage is the structural ansatz. One chooses how \(F\) depends on the auxiliary functions. In single-function reductions this often becomes
\[
u(\xi)=\sum_{n=0}^{N} H_n\,v(\xi)^n,
\]
while in genuinely multivariate settings one may use mixed products, rational expressions, or composite functions such as
\[
h(x,t)=f\big(g^{(1)}(x,t),g^{(2)}(x,t)\big).
\]
The 2024 development places particular emphasis on quadratic composite forms in two simple functions, including
\[
h=a+\beta_1 g^{(1)}+\beta_2 g^{(2)}+\gamma_1 g^{(1)2}+\gamma_2 g^{(2)2}+\gamma_3 g^{(1)}g^{(2)},
\]
and its specialization to a pure product \(h(x,t)=A\,g^{(1)}(x)\,g^{(2)}(\delta t)\) [2411.07333].

The third stage is the choice of simple equations and, when needed, their reduction. Auxiliary PDEs may be reduced to ODEs through traveling-wave variables such as
\[
\xi=\alpha x+\beta t,\qquad \eta=\rho y+\sigma t,
\]
or analogous coordinates. The reduced functions may then be expressed again through still simpler functions,
\[
a(\xi)=A[v(\xi)],\qquad b(\zeta)=B[w(\zeta)],
\]
often as finite series in powers of \(v\) and \(w\). In this sense SEsM is explicitly nested: the original PDE is rewritten through intermediate fields, which are themselves rewritten through solutions of simple equations [1908.07459], [2504.16660].

The final stage is substitution, balancing, and algebraic closure. After all replacements, the PDE is converted into a sum of linearly independent functions or monomials in the simple-function variables,
\[
DE(u,\ldots)=\sum_{\alpha} C_\alpha(\text{parameters})\,\Phi_\alpha(v,w,\ldots)=0.
\]
One then imposes balance conditions so that at least two different terms contribute to each dominant power, and sets all coefficients \(C_\alpha\) to zero. The resulting nonlinear algebraic system determines coefficients of the ansatz, parameters of the simple equations, and possibly wave numbers, frequencies, and PDE parameters. Exact parameter solutions yield exact PDE solutions; numerical parameter solutions yield numerical approximations in the same structural ansatz [1908.07459].

## 3. Mathematical structures and balance procedures

The choice of simple equation determines much of the analytic content of a SEsM construction. The cited papers repeatedly use Riccati-type equations,
\[
v'(\xi)=a_0+a_1v(\xi)+a_2v(\xi)^2,
\]
Bernoulli-type equations,
\[
v'(\xi)=\alpha v(\xi)+\beta v(\xi)^n,
\]
first-order linear exponential equations,
\[
f_\ell'(\xi)=\ell f_\ell(\xi),
\]
and second-order oscillatory equations,
\[
v_k''(\xi)+k^2v_k(\xi)=0,\qquad w_k''(\xi)+k^2w_k(\xi)=0.
\]
These choices respectively generate rational-hyperbolic, power-law, exponential, and trigonometric/Fourier structures [1908.07459].

A more general class used in earlier work has the schematic form
\[
\frac{d^k g}{d\xi^k}=\sum_{j=0}^{m} d_j g^j,
\]
and the squared-derivative class
\[
\left(\frac{da}{d\xi}\right)^2=\sum_{j=0}^{p} d_j a^j,
\]
which encompasses Weierstrass and Jacobi elliptic functions. This is how SEsM produces elliptic and periodic wave families for generalized Kawahara-type equations and related higher-order models [1908.01075].

A distinctive feature of the 2019 multisoliton paper is the explicit use of more than one balance equation. For the generalized Kawahara family
\[
u_t+\left(\sum_{k=0}^{l}\alpha_k u^k\right)u_x+\beta u_{xxx}+\gamma u^m u_{xxxxx}=0,
\]
one balance relation is reported as
\[
N(l-m)=4(p-1),
\]
and for the squared-derivative simplest equation formulation the balance changes to
\[
N=\frac{2(p-2)}{l-m}.
\]
These relations constrain truncation order and the degrees admissible in the simple equation [1908.01075].

The same paper singles out the fractional-power simple equation
\[
f_\xi=n\big[f^{(n-1)/n}-f^{(n+1)/n}\big],
\]
with solution \(f(\xi)=\tanh^n \xi\). Its role is to handle PDEs containing fractional powers. For the traveling-wave treatment of
\[
u_t=p(u^a)_x+q(u^b)_{xx},
\]
the balance procedure yields
\[
a=1+\frac{2}{n},\qquad b=1+\frac{1}{n},
\]
and the resulting exact solution has the form
\[
u(x,t)=\delta \tanh^n(\xi),
\]
with \(\xi\) determined by the corresponding algebraic parameter relations. This establishes a direct route from a fractional-power simplest equation to kink-type traveling waves with integer or fractional exponents [1908.01075].

## 4. SEsM as a unifying framework

The unifying claim of SEsM is explicit and cumulative across the cited papers. The 2019 overview identifies the Modified Method of Simplest Equation, the \(G'/G\)-method, the Exp-function method, the Tanh-method, and the Fourier-series method as particular cases. The 2025 paper adds the Jacobi Elliptic Function Expansion Method, the F-Expansion method, the Modified Simple Equation method, the Trial Function Method, the General Projective Riccati Equations Method, and the First Integral Method [1908.07459], [2504.16660].

| Method | SEsM specialization | Source |
|---|---|---|
| Modified Method of Simplest Equation | One simple equation; polynomial ansatz in its solution | [1908.07459] |
| \(G'/G\)-method and \((G'/G)_N\)-chain | One simple equation for \(w=G'/G\); polynomial ansatz in \(w\) | [1908.07459] |
| Exp-function method | Several first-order linear simple equations \(f_\ell'=\ell f_\ell\); rational ansatz in exponentials | [1908.07459] |
| Tanh-method | One simple equation \(v'=1-v^2\); polynomial or modified polynomial ansatz in \(v=\tanh\xi\) | [1908.07459] |
| Fourier-series method | Infinitely many linear second-order simple equations for trigonometric modes | [1908.07459] |
| Jacobi Elliptic Function Expansion Method | One Jacobi elliptic simple equation; traveling-wave polynomial expansion | [2504.16660] |
| F-Expansion method | Polynomial ansatz in a function satisfying an elliptic-type ODE | [2504.16660] |
| Trial Function Method | One implicit simple equation for the trial function; finite power-series structure | [2504.16660] |
| General Projective Riccati Equations Method | Projective Riccati simple equation with projective-type ansatz | [2504.16660] |
| First Integral Method | Restricted first-integral ansatz acting as an implicit simple-equation structure | [2504.16660] |

The identification is not merely terminological. In each case, the cited papers specify which SEsM step is skipped, how many simple equations remain, which ansatz is selected, and which auxiliary equation serves as the simple equation. For the classical \(G'/G\)-method, for example, one sets \(w=G'/G\), obtains
\[
w'=-w^2+\sum_{j=0}^{N} B_j w^j,
\]
and then uses
\[
u(\xi)=\sum_{m=0}^{M} a_m w(\xi)^m,
\]
which is exactly the SEsM one-simple-equation construction specialized to a polynomial Riccati-type equation for \(w\) [1908.07459].

The same logic applies to Jacobi elliptic expansion and F-expansion methods. Their distinction from SEsM lies not in a different reduction principle but in a narrower choice of simple equation and ansatz. This suggests that SEsM functions as an umbrella formalism for a large class of direct exact-solution techniques rather than as a competing isolated method.

## 5. Multisolitons, elliptic families, and composite-function extensions

One major application of SEsM is the reproduction of Hirota-type multisoliton structures. In the KdV case, the multisoliton papers use a logarithmic transformation of a tau-function-like quantity \(F\) and the finite ansatz
\[
F(x,t)=1+f_1(x,t)+f_2(x,t)+c\,f_1(x,t)f_2(x,t),
\]
with
\[
\frac{\partial f_i}{\partial x}=\alpha_i f_i,\qquad \frac{\partial f_i}{\partial t}=\beta_i f_i,
\]
so that each \(f_i\) is exponential. Substitution yields the algebraic relations
\[
\beta_i=-\alpha_i^3,\qquad
c=\frac{(\alpha_1-\alpha_2)^2}{(\alpha_1+\alpha_2)^2},
\]
and hence the classical two-soliton interaction coefficient emerges inside SEsM. The same line of argument is generalized to \(N\)-soliton tau functions built from many exponential simple equations [1908.01075], [1909.00330].

The nonintegrable side is represented by generalized Kawahara-type equations and related higher-order nonlinear dispersive models. There SEsM combines Painlevé-derived transformations such as
\[
u(x,t)=2[\ln F(x,t)]_{xx},
\]
single-function or polynomial ansätze, and Riccati or elliptic simple equations to obtain solitary, kink-type, and elliptic traveling-wave solutions. These examples are presented as evidence that the method retains the exact-solution capability of the Modified Method of Simplest Equation while extending it to more general constructions [1908.01075], [1909.00330].

The 2024 paper extends SEsM beyond standard traveling-wave reductions by focusing on derivatives of composite functions of two simple-equation solutions. For
\[
h(x_1,\dots,x_d)=f\big(g^{(1)}(x_1,\dots,x_d),\dots,g^{(m)}(x_1,\dots,x_d)\big),
\]
it invokes the multivariate Faa di Bruno formula, following Constantine and Savits, to organize higher derivatives. In the two-variable, two-function case this yields explicit chain-rule formulas for \(h_x\), \(h_t\), \(h_{xx}\), and \(h_{tt}\), which are then inserted into the nonlinear PDE
\[
(1+h^2)\,h_{xx}-2h(h_x)^2-h_{tt}=h(1-h^2).
\]
With simple equations of Jacobi elliptic type for \(g^{(1)}\) and \(g^{(2)}\), the paper derives explicit solution families including
\[
h(x,t)=A\,\mathrm{cn}(ax;k_1)\,\mathrm{cn}(dt;k_2),
\]
\[
h(x,t)=A\,\frac{\mathrm{sn}(ax;k_1)\,\mathrm{cn}(ax;k_1)}{\mathrm{dn}(\delta t;k_2)},
\]
\[
h(x,t)=A\,\mathrm{dn}(ax;k_1)\,\mathrm{sn}(\delta t;k_2),
\]
and
\[
h(x,t)=A\,\mathrm{dn}(ax;k_1)\,\frac{\mathrm{sn}(\delta t;k_2)}{\mathrm{cn}(\delta t;k_2)},
\]
together with trigonometric and hyperbolic limits such as \(A\,\mathrm{sech}(ax)\cos(dt)\) [2411.07333].

## 6. Scope, limitations, and interpretation

SEsM is primarily a constructive method for exact particular solutions of nonlinear PDEs, especially traveling waves, kinks, solitary waves, periodic waves, and multisolitons. At the same time, the Fourier-series embedding shows that it also encompasses approximate constructions for linear PDEs, because truncated Fourier expansions fit the SEsM pattern of representing the solution through many simple equations. The cited literature further remarks that the method of orthogonal functions is likewise a particular case [1908.07459].

A recurrent misconception is to equate SEsM with a single finite power series in one auxiliary function. The cited papers explicitly reject that restriction: the one-simple-equation polynomial expansion is only the MMSE or modified-simple-equation corner of a much broader formalism. Another misconception is to treat SEsM as intrinsically a traveling-wave method. Traveling-wave reduction is common and often convenient, but the composite-function treatment of \(h(x,t)=f(g^{(1)}(x,t),g^{(2)}(x,t))\) shows that SEsM also operates directly with multivariable composite structures without an explicit traveling coordinate [1908.07459], [2411.07333].

Its limitations are also stated, or directly implied, in the cited work. There is no general algorithm for selecting the optimal transformation \(T\), the most effective ansatz for \(F\), or the most productive simple equations. The method searches within a chosen ansatz class, so solutions outside that class are missed. The algebraic systems that arise after substitution can become large and strongly nonlinear, often requiring computer algebra and sometimes yielding only numerical parameter values rather than closed forms. No completeness theorem is claimed: SEsM is a constructive framework for obtaining some exact or approximate solutions, not all possible solutions of a given PDE [1908.07459], [2504.16660].

The subsequent trajectory of the method broadens rather than narrows its scope. The 2019 papers emphasized unification, multisolitons, Painlevé-type transformations, more than one simple equation, and fractional-power nonlinearities. The 2024 paper formalized composite derivatives through multivariate Faa di Bruno machinery. The 2025 paper expanded the catalog of named methods that can be embedded into SEsM. Taken together, these developments portray SEsM as a flexible methodological architecture whose defining feature is not a particular special function or ansatz, but the systematic reduction of nonlinear PDEs to algebraic consistency conditions through appropriately chosen simple equations [1908.01075], [2411.07333], [2504.16660].

Source: https://www.emergentmind.com/topics/simple-equations-method-sesm