REVIEW 3 major objections 5 minor 7 references
A system of 2 nonlinearly coupled ODEs which is explicitly solvable and possibly isochronous provided its coefficients are suitably restricted
T0 review · 3 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read A system of two coupled ODEs with 23 coefficients is explicitly solvable when 12 coefficient constraints hold, and becomes isochronous under two further conditions.
desk verdict A genuine but incremental extension of a known linearization trick; the forward construction is solid, the sufficiency of the 12 constraints needs external verification. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing device is the polynomial substitution (6) together with the linear system (2). The substitution is deliberately chosen so that differentiating $y$ along (1) yields expressions that collapse to linear functions of $y$ when the coefficients obey (9). The inverse formulas (7)--(8) invert the map through a square root, and the discriminant $S(t)$ controls the singularity line. The 12 constraints (17)--(18), and their versions in the special cases, are obtained by eliminating the 15 auxiliary parameters from the 23 formulas (9), and they certify that a given system belongs to the explicitly solvable class.
What would settle it
A reader can settle the claim by picking a concrete instance of the class, for instance the isochronous Example 2 with initial data $x_1(0)=1$, $x_2(0)=0.2$; the explicit formula gives a real trajectory only up to the time $t_S$ when $x_1(t)=0$, so a numerical integration of the original ODEs that continued smoothly past that time would show the formula misses a branch. Conversely, substituting the formulas (7) and (3) into (1) by computer algebra for generic parameters would expose any algebraic error in the coefficient formulas (9).
Extended reading notes
Core claim
The central discovery is that the class of systems (1) with $f_k,g_k,h_\ell$ given by the explicit formulas (9) in terms of the 15 parameters $\eta_{nm},\xi_n,\alpha_\ell,\beta_n,\gamma,\lambda_n,\mu$ is exactly solvable: the polynomial map (6), $y_1=\alpha_0x_1^2+\alpha_1x_1x_2+\alpha_2x_2^2+\beta_1x_1+\beta_2x_2+\gamma$, $y_2=\lambda_1x_1+\lambda_2x_2+\mu$, transforms the nonlinear equations (1) into the linear system (2), whose solution is the explicit exponential formula (3). The inverse map (7)--(8) then gives $x_1(t),x_2(t)$ explicitly, up to a square root whose sign is fixed by the initial data and whose vanishing marks the only obstruction to continued uniqueness. Inverting the coefficient formulas produces 12 constraints on the 23 coefficients, leaving 11 coefficients freely assignable; when in addition $\eta_{22}=-\eta_{11}$ and $(\eta_{11})^2+\eta_{12}\eta_{21}=-\omega^2<0$, all nonsingular solutions are periodic with common period $T=2\pi/\omega$.
Load-bearing premise
The construction assumes the map (6) can be inverted on the whole time interval where the formula is used: in the generic case $C_2\neq 0$, or $C_1(t)\neq 0$ when $C_2=0$, and the chosen square-root branch must not cross $S(t)=0$; the paper leaves open exactly which initial data meet these conditions.
Editorial extensions
If this is right
- Any system whose coefficients satisfy the 12 constraints has explicit solution formulas, so its solutions can be written without numerical integration, and the formulas contain all dependence on initial data.
- When the two conditions (4) hold, every solution that stays nonsingular for all time is periodic with the same period $2\pi/\omega$; the paper's examples show that some trajectories instead terminate at the singularity line before one period.
- The polynomial subcase $P_1\equiv 1$ yields explicitly solvable systems of two coupled cubic ODEs, including non-homogeneous cubics not covered by the previous homogeneous treatment.
- The class contains systems with stable, unstable, and neutral equilibria, as well as systems with no real equilibria; the equilibrium type is inherited from the associated linear system because the map (6) sends the linear equilibrium to the two preimages (33).
- The transformation reduces the number of independent coefficients: a linear change of variables puts the system in a form with 19 coefficients satisfying 8 constraints, with one equation becoming purely polynomial.
Reading between the lines
- This suggests the same linearization mechanism is likely to work for higher-degree numerators or more variables: any polynomial map from $x$ to $y$ that linearizes the flow imposes matching conditions on the coefficients, and the 12-constraint structure is one instance of that matching.
- The constraints are shown to be sufficient, not necessary; a natural question is whether a wider set of solvable systems exists whose coefficients do not satisfy (17)--(18), which would amount to finding a different linearizing map.
- The explicit formulas provide a cheap laboratory for singularity studies: one can choose the linear system to be periodic and the parameters so that the image trajectory hits $P_1=0$, yielding computable values of $t_S$ and a test case for numerical methods for rational vector fields with movable singularities.
- Because the period $T$ is set by $\omega$ alone, the class can be used to design two-dimensional ODE models with a prescribed common oscillation frequency, which may be useful when synchronized periodic outputs are observed.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies the planar system (1) of two first-order ODEs with cubic numerators and a common linear denominator, containing 23 coefficients. The authors introduce a quadratic-linear ansatz (6) relating the variables x_n to variables y_n that evolve according to the linear system (2). Substitution yields 23 explicit formulas (9) expressing the coefficients of (1) in terms of 15 parameters, and hence explicit solutions via (3) and (7)–(8). The paper then claims to invert these formulas and obtains 12 explicit constraints—eqs. (17)–(18), (19), and (25) in the three subcases of Section 2.4—that are asserted to be sufficient for the system to belong to the explicitly solvable class. It also identifies conditions (4) that imply isochrony, and presents five worked examples illustrating periodic, singular, and hyperbolic behaviors. Throughout, the manuscript is candid about possible loss of determinacy when the square-root S(t) in (7c) vanishes, but the formal claims in the abstract and Section 4 are stated without such qualifications.
Significance. If the elimination algebra is correct, the result is a useful and interesting contribution to the literature on explicitly solvable and isochronous planar polynomial/rational ODE systems, in the spirit of earlier work by Calogero and collaborators. The constructive ansatz is transparent, the explicit solution formulas are valuable, and the examples in Section 3 provide concrete illustrations of the different qualitative regimes. The paper does not ship machine-checked proofs or a verification artifact, which is a real limitation because the central sufficiency claim rests on lengthy asserted elimination steps. The branch-choice and finite-time-singularity caveats are acknowledged in Remark 2.2-1, but they are not integrated into a precise theorem, so the phrase 'arbitrary initial data' in the abstract overstates what is actually proved.
major comments (3)
- [§2.4.1, eqs. (13)–(18) and eq. (7)] The inversion in Section 2.4.1 uses the transformation (7), which is only valid when C2 ≠ 0, but the case C2 = 0 with h1 and h2 both nonzero is not treated. This is not merely hypothetical: for α0 = 1, α1 = 3, α2 = 2, λ1 = λ2 = 1, one obtains h1 = h2 = −1 while C2 = 0. In this regime the transformation must be replaced by (8), yet no inversion formulas are supplied and the text does not state that such coefficient sets are excluded. Since the paper claims that the 12 constraints are sufficient whenever they hold (with nonvanishing denominators), the proof is incomplete for an open subset of the image of (9). Please add this case explicitly or restrict the sufficiency statement to C2 ≠ 0, and check whether the constraints themselves force C2 ≠ 0.
- [Abstract, Remark 1-1, and Remark 2.2-1] The abstract and Remark 1-1 promise that the initial-value problem can be explicitly solved 'with arbitrary initial data', but the construction requires a choice of branch in (7a), continuation past values where S(t) in (7c) vanishes, and nonvanishing of the denominators appearing in the inversion formulas. Remark 2.2-1 concedes that solutions may cease to be determined at a finite time tS when S(tS) = 0. The paper never delimits the set of initial data for which the explicit formula defines a genuine solution on a maximal interval, nor does it state the branch-selection rule as part of a formal theorem. This ambiguity affects the central claim and should be resolved by stating a precise theorem with the maximal interval, the branch rule, and the admissible initial data.
- [§2.4, eqs. (17)–(18), (19), (25)] The derivation of the 12 constraints is described as 'trivial if tedious algebra' performed with Mathematica, but no derivation, intermediate algebraic system, or machine-readable verification is included. Because the sufficiency claim is the main result and the formulas are long, I could not independently verify from the text that eqs. (17)–(18) are exactly the elimination conditions, and a single algebraic error there would invalidate the theorem. Please provide a reproducible computation—for example, a short Mathematica or SymPy script that checks both the substitution of (9) into (1) and the elimination leading to (17)–(18)—or list the reduced intermediate equations, and explicitly state all denominator conditions such as g6 ≠ 0, g8 ≠ 0, g9 ≠ 0, h1 ≠ 0, h2 ≠ 0, g8h1 − 3g9h2 ≠ 0, and 9g6g9h1 − (g8)^2 h2 ≠ 0.
minor comments (5)
- [§3.2, eq. (49a)] The displayed expression for y1(0) contains an apparent typo: the term '−3 [x2(0)]^2' appears twice, likely intended to be '−3 x2(0)'.
- [§3.3, eq. (55b) and §3.4, eq. (63a)] In eq. (55b) the second equilibrium point is labeled with X(+)_2 instead of X(−)_2, and in eq. (63a) both equilibria are labeled with the superscript (+); the second should be (−).
- [§2.5, eq. (27)] The symbol λ is used both as a parameter in eqs. (6) and (9) and as the norm λ = sqrt(λ1^2 + λ2^2) in eq. (27); please use a different symbol for the norm to avoid confusion.
- [Throughout] There are numerous typographical errors, including 'isochonous', 'equilibra', 'sy stem', and 'ODES'; a careful copyedit is recommended.
- [§3.1–3.5] Several displayed solutions, such as eqs. (48a), (72a), and (74a), are extremely long and contain unbalanced parentheses in places (e.g., eq. (48a)); moving the full expressions to an appendix or providing a supplementary computation file would improve readability and checkability.
Circularity Check
No significant circularity: the solvable class is constructed from a solvable linear system and an explicit invertible change of variables; the 12 constraints are derived compatibility conditions, not fitted or assumed conclusions.
full rationale
The derivation chain starts from the explicitly solvable linear system (2) and the explicit transformation (6)-(7). The paper claims that if y solves (2) and x,y are related by (6), then x solves (1) provided the coefficients are given by the 23 formulas (9); this is a direct substitution statement and is the construction itself, not a fitted conclusion. The 12 constraints in Section 2.4 are obtained by eliminating the 15 parameters from the 23 formulas (9), i.e., they are compatibility conditions for membership in the image of a parametrization; they are not assumed inputs. Nothing in the paper defines the target solvability in terms of the constraints, nor does it use the conclusion as an input. The only self-references are to background work ([5], [7]); Remark 3.1-1 explicitly says the earlier homogeneous case is not included, so that citation is not load-bearing. The unresolved issues flagged by the skeptical reading—that the Mathematica elimination is not displayed and the C2=0 with h1,h2 nonzero regime is not inverted—are verification/completeness concerns, not circularity. Remark 2.2-1 also flags branch and singularity limitations, which concern the domain of validity of the solutions, not circularity. Therefore the paper is self-contained in its construction and no circular step is present. Score 0.
Assumptions & free parameters
free parameters (7)
- alpha_0, alpha_1, alpha_2 (quadratic coefficients in y1 ansatz, eq. (6a))
- beta_1, beta_2 (linear coefficients in y1 ansatz, eq. (6a))
- gamma (constant in y1 ansatz, eq. (6a))
- lambda_1, lambda_2 (linear coefficients in y2 ansatz, eq. (6b))
- mu (constant in y2 ansatz, eq. (6b))
- eta_11, eta_12, eta_21, eta_22 (matrix of linear system (2))
- xi_1, xi_2 (inhomogeneity of linear system (2))
assumptions (4)
- standard math Explicit solution of the 2x2 linear system (2) via eigenvalues, eq. (3).
- domain assumption Generic invertibility of the change of variables: denominators such as C2, h1, h2, alpha_1^2 - 4 alpha_0 alpha_2, g6, g8, g9 and the combined denominators in eqs. (13)-(18) are nonzero.
- domain assumption For real initial data and real parameters, the square root S(t) in (7c) remains real on the interval where the solution is defined, and the sign chosen at t=0 is preserved by continuity.
- standard math The lengthy polynomial identities in Section 2.4, verified with Mathematica, are correct as stated.
Cite this review
Pith. "Pith review of A system of 2 nonlinearly coupled ODEs which is explicitly solvable and possibly isochronous provided its coefficients are suitably restricted." pith.science (2026). https://pith.science/paper/KCM6YGXO
@misc{pith2026250524370,
author = {Pith},
title = {Pith review of: A system of 2 nonlinearly coupled ODEs which is explicitly solvable and possibly isochronous provided its coefficients are suitably restricted},
year = {2026},
howpublished = {\url{https://pith.science/paper/KCM6YGXO}},
note = {Machine review of arXiv:2505.24370}
}
abstract
In this paper we discuss some remarkable properties of the autonomous system of 2 first-order Ordinary Differential Equations (ODEs), which equates the derivatives $\dot{x}_n(t)$ ($n = 1, 2$) of the 2 dependent variables $x_n(t)$ to the ratios of polynomials (with constant coefficients) in the 2 variables $x_n (t)$: each of the 2 (a priori different) polynomials $P_3^{(n)}(x_1, x_2)$ in the 2 numerators is of degree 3; the 2 denominators are instead given by the same polynomial $P_1(x_1, x_2)$ of degree 1. Hence this system features 23 a priori arbitrary input numbers, namely the 23 coefficients defining these 3 polynomials. Our main finding is to show that if these 23 coefficients are given by 23 (explicitly provided) formulas in terms of 15 a priori arbitrary parameters, then the initial values problem (with arbitrary initial data $x_n (0)$) for this dynamical system can be explicitly solved. We also show that it is possible (with the help of Mathematica) to identify 12 explicit constraints on these 23 coefficients, which are sufficient to guarantee that this system belongs to the class of systems we are focusing on. Several such explicitly solvable systems of ODEs are treated (including the subcase with $P_1(x_1, x_2) = 1$, implying that the right-hand sides of the ODEs are just cubic polynomials: no denominators!). Examples of the solutions of several of these systems are reported and displayed, including cases in which the solutions are isochronous.
Reference graph
Works this paper leans on
-
[1]
R. Garnier, Sur des syst` emes diff´ erentielles du second ordre dont l’ int´ egral g´ en´ eral est uniforme, Annales scientifiques de l’ ´E. N. S., 3 eme s´ erie, 77, 2 (1960), 123-144
work page 1960
-
[2]
E. L. Ince, Ordinary Differential Equations , Dover, 1956
work page 1956
-
[3]
S. H. Strogatz, Nonlinear Dynamics and Chaos , Perseus Books, 1994
work page 1994
-
[4]
F. Dumortier, J. Llibre and J. C. Art´ es,Qualitative Theory of Planar Differential Systems, Springer, 2006
work page 2006
-
[5]
F. Calogero, Isochronous Systems , Oxford University Press, 2008 (hardback version), 2012 (updated paperback version)
work page 2008
-
[6]
A. D. Polyanin and V. F. Zaitsev, Handbook of Ordinary Differential Equations , Taylor & Francis, 2017
work page 2017
- [7]
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.