Generalized ODE reduction algorithm for bounded degree transformation
By: Shaoxuan Huang
Potential Business Impact:
Simplifies hard math problems for computers.
The integrability problem of rational first-order ODEs $y^{\prime}=\frac{M(x,y)}{N(x,y)}$, where $M,N \in \mathbb{R}[x,y]$ is a long-term research focus in the area of dynamical systems, physics, etc. Although the computer algebra system such as Mathematica, Maple has developed standard algorithms to tackle its first integral expressed by Liouvillian or special function, this problem is quite difficult and the general method requires specifying a tight degree bound for the Darboux polynomial. Computing the bounded degree first integral, in general, is very expensive for a computer algebra system\cite{duarte2021efficient}\cite{cheze2020symbolic} and becomes impractical for ODE of large size. In \cite{huang2025algorithm}, we have proposed an algorithm to find the inverse of a local rational transformation $y \to \frac{A(x,y)}{B(x,y)}$ that transforms a rational ODE to a simpler and more tractable structure $y^{\prime}=\sum_{i=0}^nf_i(x)y^i$, whose integrability under linear transformation $\left\{x \to F(x),y \to P(x)y+Q(x)\right\}$ can be detected by Maple efficiently \cite{CHEBTERRAB2000204}\cite{cheb2000first}. In that paper we have also mentioned when $M(x,y),N(x,y)$ of the reducible structure are not coprime, canceling the common factors in $y$ will alter the structure which makes that algorithm fail. In this paper, we consider this issue. We conclude that with the exact tight degree bound for the polynomial $A(x,y)$ given, we have an efficient algorithm to compute such transformation and the reduced ODE for "quite a lot of" cases where $M,N$ are not coprime. We have also implemented this algorithm in Maple and the code is available in researchgate.
Similar Papers
Generalized ODE reduction algorithm with bounded degree transformation
Symbolic Computation
Simplifies complex math equations for easier study.
Faster multivariate integration in D-modules
Symbolic Computation
Solves hard math problems for counting special graphs.
Support bound for differential elimination in polynomial dynamical systems
Symbolic Computation
Finds hidden math rules in changing systems.