Generalized ODE reduction algorithm for bounded degree transformation
Abstract
The integrability problem of rational first-order ODEs y=M(x,y)N(x,y), where M,N ∈ 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 systemduarte2021efficientcheze2020symbolic and becomes impractical for ODE of large size. In huang2025algorithm, we have proposed an algorithm to find the inverse of a local rational transformation y A(x,y)B(x,y) that transforms a rational ODE to a simpler and more tractable structure y=Σi=0nfi(x)yi, whose integrability under linear transformation \x F(x),y P(x)y+Q(x)\ can be detected by Maple efficiently CHEBTERRAB2000204cheb2000first. 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.
Turn this paper into a full lesson
ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.