Open access peer-reviewed chapter

On the Generalized Simplest Equations: Toward the Solution of Nonlinear Differential Equations with Variable Coefficients

By Gunawan Nugroho, Purwadi Agus Darwito, Ruri Agung Wahyuono and Murry Raditya

Submitted: November 13th 2020Reviewed: December 22nd 2020Published: January 14th 2021

DOI: 10.5772/intechopen.95620

Downloaded: 189


The simplest equations with variable coefficients are considered in this research. The purpose of this study is to extend the procedure for solving the nonlinear differential equation with variable coefficients. In this case, the generalized Riccati equation is solved and becomes a basis to tackle the nonlinear differential equations with variable coefficients. The method shows that Jacobi and Weierstrass equations can be rearranged to become Riccati equation. It is also important to highlight that the solving procedure also involves the reduction of higher order polynomials with examples of Korteweg de Vries and elliptic-like equations. The generalization of the method is also explained for the case of first order polynomial differential equation.


  • the simplest equation
  • Riccati equation
  • nonlinear differential equations
  • reduction of polynomial

1. Introduction

Despite the advent of supercomputers in numerical methods, increasing activities are devoted to solving nonlinear differential equations by analytical method in recent years [1, 2, 3]. Analytical solutions have their own importance concerning the physical phenomena as they are often pave the way to the construction of right theory [4]. Many methods have been proposed concerning this important problem and generally, for the problems with constant coefficients. One of the useful methods is the method of simplest Equation [5] or for some authors, the auxiliary Equations [6]. The method is built by the utilization of the first integral of simplest nonlinear differential equations, such as Bernoulli and Riccati Equations [7]. The method had produced many new solutions of the considered nonlinear differential equations, generally with constant coefficients [8, 9].

For the more general cases, we have found that the method can be extended such as involving the solution of the nonlinear differential equations with variable coefficients. The nature of variable coefficients often arises in the equation describing the heterogenous media and composites [10] or in other cases are produced by the coordinate transformation of the partial differential Equations [11]. Those two categories are developed rapidly in recent years with the capacity of high-speed super computers which sufficient for computing nonlinear problem with complex geometries [12, 13], as sometimes desired by engineering design activities. The role of analytical solutions is as a benchmark to validate the computer algorithm with simpler geometries as it is usually performed [14].

In this chapter, the solution method of the simplest equations is different from the cases of constant coefficients except, on Bernoulli equation. Hence, we will start from Riccati equation instead of Bernoulli equation as the simplest equation to highlight the novelty of the procedure. The method is then followed by examples and conclusion.


2. The first integral of the simplest equations with variable coefficients

2.1 Riccati equation

Consider the Riccati equation with variable coefficients as follows,


Let A=β1β2and the above equation can be rearranged as,


and is separated as


The solutions for β1and β2are




The relation for A=β1β2is thus,


Without loss of generality, suppose that β1=eξγdξand thus the above relation is performed as,


Rearrange Eq. (5) and integrate once,




The solution for Ais then,


where the coefficients aiwill be determined later from the substitution into the considered nonlinear differential equations.

2.2 The Jacobi and Weierstrass equations

It is interesting to know that other simplest equations can also be rearranged into the Riccati-type equations. The famous examples are Jacobi [15] and Weierstrass Equations [16], which can solve a large class of nonlinear differential equations. Let us consider Jacobi type equation with variable coefficients,


and the Weierstrass equation as follows,


Here, the reader should not be confused by the coefficients which represent different functions with the same index. Take ϕ=1ν+aξand the Weierstrass equation becomes Jacobi equation which admits the similar method of solution.

Concerning the search for obtaining solution of (7) and (8), the balancing principle suggests the substitution of the first order series ϕ=b6+b7Aas in the following,


Performing the Riccati equation Aξ=a1b7A2+a2A+a3into (9) and we generate the following expression,


The coefficients of polynomial are then related with the coefficients in (9) in order to determine a1,a2,a3,b6,b7as functions of the known b1,b2,b3,b4and b5as follows,


Hence, the first equation gives,


and the second equation is then,


The next relation produces,


Eqs. (12) and (13) are thus substituted into the fourth relation of (12) to form the third order polynomial equation in term of b6as follows,


which the roots will determine the solution for b6as functions of b1,b2,b3,b4,b7, or b6=f1b7in simple unknown variable. The step now is to find the polynomial expression for b7from the last relation of (10) as,


Therefore, the last equation gives the expression for b7as polynomial equation of higher order, and the generated polynomial is,


In this case the higher order polynomial will be solved by reducing the order.


3. Reduction of higher order polynomial

Consider the sixth order polynomial equation as in the following,


First, multiply the above equation with the functionαand rearranged as,


where B=αb7. The polynomial equation is cut as in the following,


Note that, the coefficients biin this section is different from the previous section. Expanding for the new coefficients,


Hence, the relation for coefficients is,


The fifth coefficient relation is rearranged as,


and the roots are,


Also, the last relation is rewritten as,


Note that performing (19) into (21) will remove b5and α. Thus, it is necessary to take other relation, i.e. a6α6=α12+b5, which will produce the cubic equation as follows,


which has the roots as,


Substituting back into a6α6=α12+b5to get,


Therefore,αis also determined by (24) and so all the coefficients, b1,b2,b3,b4,b6,φ. The polynomial equation of sixth order is then re-expressed as,


as reduced into the quartic equation the roots can be obtained by radical solution.

The procedure described by (17–25) can be applied and iterated into (16) until the polynomial equation of b7is reduced into quartic equation. Hence, all the coefficients for Riccati equation the first order series, i.e. a1,a2,a3,b6,b7are determined and produce the solution as,


Thus, following the method explained by (2–6) and (17–25), we have arrived at the solution of Jacobi and Weierstrass equations with variable coefficients.


4. Solution examples

4.1 The elliptic-like equation

As an application, consider the elliptic-like equation with forcing function,


The balance principle suggests that the solution should be in the form,


Substituting into (27) will reproduce the following expression,


The next step is to differentiate the Riccati equation once,


Substituting into Eq. (29) and it will produce the polynomial equation as in the following,


or the next step is to relate the coefficients as,


In this case, the first equation gives,


For the second equation,


Thus, provide the expression for b5as,


The third and fourth equations produce,


Substituting (32) into (31),


Replace b5with (10),


which then solves a2regardless of b4. In this case, we take b4as the chosen fundamental variable and the resulted coefficients, b5,a1,a2,a3depend on b4and with the known coefficients b1,b2,b3. Therefore, the solution of (27) is generated as,


4.2 Korteweg de Vries equation

The next example is for the Korteweg de Vries type equation,


The balancing principle with application of Riccati equation will determined the ansatz,


Performing into (35) will produce,


Performing the Riccati equation into (37) will produce the following polynomial,


From this step on, there is a little hope to solve all the coefficients as they are equal to zero. As it has also to be reduced, it is important to note that the problem of reduction here is different from the case of Jacobi equation since all the coefficients are in principle solvable in algebraic form. In this case, it is not practical to reduce the fifth order polynomial as the even highest power, i.e., as a tenth order polynomial equation. The calculation will become too tedious as the detail expression is needed in the reduced polynomial equation. The next sub section will illustrate the reduction of an odd highest power polynomial equation.

4.3 Reduction of fifth order polynomial

Consider Eq. (38) as follows,


Multiply by the function, βand rearrange,


where, B=βA. Rearranged Eq. (39) as given by,


Expanding the all the coefficients as,


Relate the coefficients as in the following,


The fifth equation of (41) gives the roots as,


Moving to the last equation, the functions b5and βdisappear from the operation. In this case we will consider the test function, b5+β10=d6β5, and will perform as,


Substituting for b3, the expression for βis,


Substitute back to, b5+β10=d6β5as follows


which then solves b5, β, φand thus generates all the coefficients of bi. The polynomial is then rewritten as,


which is reduced as,


Eq. (46) dictates that the relations, d1=d2=d3=d4=0will satisfy for the solution. Hence, the coefficients are then,


The first equation gives,


The second equation is rewritten as,


The solution for b6is then,


The third equation will produce,


Take the expression for b5as,




The fourth relation of (47) will generate,


which then produce the solution of a3.

Note that a1and a2are the chosen fundamental variables and according to (48–51) and with the known coefficients b1,b2,b3,b4, they will define, b5,b6,b7,a3,β. Therefore, the solution of Korteweg de Vries equation is generated as,


Since only a few of the considered equation has a special polynomial to be solved by equating all the variable coefficients to zero, it is important to note that the reduction of polynomial order would be an important step. Solving all coefficients to zero often be an obstacle because the difficulty would be the same or even more than the original nonlinear ODEs. In this case, the reduction of polynomial manipulates and reduces the need for solving all coefficients.

However, it is possible not to search for the expression of variable coefficients, i.e., b5,b6,b7,a1,a2and a3. First the roots of Eq. (37) are determined first as ϕ, and then Eq. (1) is decomposed as,


with A=BD. The solution of (53) is then,


The definition for Bis determined by substituting the polynomial solution, ϕinto (54) as in the following,


Rearranging the above equation as,


Differentiating once,


Eq. (55) is a first order ODE in Band can be easily solved, which then prove that A=ϕwithout establishing the explicit expression for variable coefficients.


5. Generalized method

In this section, the method of solution to the Riccati equation is extended for the class of the first order polynomial differential equation as,


The above equation can be always re-expressed as,




which the coefficients will be reformulated as,


In this case, we will always obtain the new coefficients bi. Proceeding into the other equations andmultiply the equation by the function α, to get,


where B=αA. Then, all the new coefficients in biwill be determined. The step is now to solve the Riccati equation. Let B=β2β3, the equation can be rearranged as,


and is separated as,


The solutions for β2and β3are,


The relation for B=β2β3is thus,


Without loss of generality, suppose that β2=φeξγdξand the above relation is performed as,


Rearrange the above equation as,


Let e2ξb1αφ=αφ1b0and integrate the above equation to get,


The solution for Bis then reduced into the solution of the polynomial equation. Thus, let A=α1B=ϕ, whereϕis the expression from the solution of the resulting polynomial equation which is similar to (38). The expression for αcan be determined by the inverse method as in (53–55) for the first order polynomial differential Eq. (56).


6. Conclusion

In this chapter, we propose the method of the simplest or the auxiliary equation to solve the nonlinear differential equation with variable coefficients. The method is based on the solution of the generalized Riccati equation as the simplest equation. It is found that the other known simplest equations, i.e., Jacobi and Weierstrass equation, are also solved by the Riccati equation. The applications with the variable coefficients elliptic-like and Korteweg de Vries equations show that the problem of solving nonlinear differential equations with variable coefficients are simplified, especially by the reduction of the resulting polynomial equation in solving the Korteweg de Vries equation. The generalization of the method is also derived in detail.


Conflict of interest

Authors declare that there is no conflict of interest.

© 2021 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Gunawan Nugroho, Purwadi Agus Darwito, Ruri Agung Wahyuono and Murry Raditya (January 14th 2021). On the Generalized Simplest Equations: Toward the Solution of Nonlinear Differential Equations with Variable Coefficients, Recent Developments in the Solution of Nonlinear Differential Equations, Bruno Carpentieri, IntechOpen, DOI: 10.5772/intechopen.95620. Available from:

chapter statistics

189total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Thermodynamic Stability Conditions as an Eigenvalues Fundamental Problem

By Francisco Nogueira Lima

Related Book

First chapter

Text Mining for Industrial Machine Predictive Maintenance with Multiple Data Sources

By Giancarlo Nota and Alberto Postiglione

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us