A Systematic Analysis of the Properties of the Generalised Painlev\'e--Ince Equation

Andronikos Paliathanasis; P.G.L. Leach

REVIEW 3 major objections 5 minor 21 references

A Systematic Analysis of the Properties of the Generalised Painlev\'e--Ince Equation

T0 review · 3 major / 5 minor · reviewed 2026-08-14 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read The generalized Painlevé–Ince equation passes the Painlevé test under inversion exactly when α² = 9β, the same condition that gives maximal symmetry.

desk verdict A correct but routine symmetry analysis with a confident concluding inference that does not survive contact with the ARS algorithm. read the letter →

arxiv 1908.04563 v1 pith:CHI7BBJS submitted 2019-08-13 nlin.SI math.CA

Andronikos Paliathanasis , P.G.L. Leach This is my paper · ORCID

classification nlin.SImath.CA MSC 34A3434C1434M35

keywords generalizedPainlevé–InceequationPainlevétestsingularityanalysisLiepointsymmetriesintegrabilityLaurentseriescoordinatetransformationresonances

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper asks when the generalized Painlevé–Ince equation $\ddot{x}+\alpha x\dot{x}+\beta x^3=0$ is integrable. For generic values of $\alpha$ and $\beta$ the equation has only two Lie point symmetries and fails the Painlevé test. By inverting the dependent variable, $x=1/y$, the authors obtain a transformed equation whose singularity analysis succeeds exactly when $\alpha^2=9\beta$, the same condition that gives the maximally symmetric classical case. They conclude that the classical Painlevé–Ince equation is the only member of the family that passes the Painlevé test under this inversion, and that this route resolves the older awkwardness about its resonances.

What carries the argument

The load-bearing mechanism is the inversion $x(t)=1/y(t)$, which maps (2) to $y\ddot{y}-2\dot{y}^2+\alpha\dot{y}-\beta=0$. In this transformed equation the leading-order terms are only the first two, so the singularity analysis becomes clean: a simple pole with an arbitrary leading coefficient and resonances at $-1$ and $0$. The decisive step is the computation of the higher-order terms of the Laurent expansion; their coefficients vanish exactly at $\alpha^2=9\beta$, and this vanishing is what the paper uses to conclude that the original equation passes the Painlevé test under the transformation. A Right Painlevé series, meaning an expansion valid on a disc centered on the singularity, is the type of series that emerges here.

What would settle it

A concrete check: substitute a Right Painlevé series $y=a\tau^{-1}+b+c\tau+\cdots$ into the inverted equation (19) and solve for the higher coefficients symbolically for arbitrary $a,b,\alpha,\beta$. The paper's claim predicts that every higher coefficient is a nonzero multiple of $\alpha^2-9\beta$ (or a power of it). If any coefficient vanishes for some $\alpha^2\neq 9\beta$, or fails to vanish at $\alpha^2=9\beta$, the 'if and only if' claim is false.

Watch

Extended reading notes

Core claim

The paper's central claim is that equation (2) passes the Painlevé test under the coordinate transformation $x=1/y$ if and only if $\alpha^2=9\beta$. Under this inversion, (2) becomes $y\ddot{y}-2\dot{y}^2+\alpha\dot{y}-\beta=0$, for which only the first two terms are dominant: the singularity is a simple pole, the leading-order coefficient is arbitrary, and the resonances are $-1$ and $0$, matching a Right Painlevé series. The higher-order coefficients of that series vanish precisely when $\alpha^2=9\beta$, i.e., exactly for the Painlevé–Ince form with maximal symmetry. The authors therefore infer that this parameter relation is the exact integrability condition in the Painlevé sense.

Load-bearing premise

The argument assumes that the singularity analysis of the inverted equation (19) tells us whether the original equation (2) passes the Painlevé test, even though the two are related by a point transformation and the Painlevé property is not generally invariant under such transformations.

Editorial extensions

If this is right

For generic $\alpha,\beta$ the equation has exactly two Lie point symmetries and the direct singularity analysis of (2) fails, so no Laurent-series-based integrability is available in those cases.
The relation $\alpha^2=9\beta$ coincides with maximal symmetry (eight Lie point symmetries) and with the vanishing of the higher-order coefficients in the inverted expansion, unifying the symmetry and singularity criteria for this family.
The inversion $x=1/y$ provides a concrete template for applying singularity analysis to equations whose direct Laurent expansion is obstructed.
All members of the family reduce to an algebraic equation via the two symmetries, so reduction-based integrability holds generally, while Painlevé-type integrability is confined to the special relation.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

The transfer of the Painlevé-test result from (19) back to (2) is an inference the paper itself labels as such ('we can infer'); the Painlevé property is not generally invariant under point transformations, so a reader should treat the 'if and only if' as a conjecture about this specific transformation unless a direct proof of equivalence is supplied.
The same inversion trick could be tried on other two-symmetry polynomial ODEs with a single balance: whenever the inverted equation has an arbitrary leading coefficient and resonances $-1,0$, the vanishing of higher-order coefficients may single out integrable parameter values.
A direct test of the claim would be to compute the Laurent expansion of (2) itself for a parameter set with $\alpha^2\neq 9\beta$ and check whether any of the two possible leading-order branches admits a complete two-constant series; if one does, the 'only if' direction would need refinement.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Short version: the symmetry part of this paper is fine and the resonance calculation for (2) checks out, but the concluding iff claim is not supported, and the stress-test note is right: (19) passes the ARS test for every α and β.

What is actually there: for general α, β, equation (2) has two Lie point symmetries; at β = α²/9 the algebra jumps to eight and this is the classical Painlevé–Ince equation. That is known, and the paper cites the relevant literature. The reduction to an algebraic equation via invariants is a clean piece of standard symmetry analysis, and the leading-order coefficients and resonance formulas for (2) are correct. Fine, but old.

Where it goes wrong: Section 4. After x = 1/y, equation (19) has only two dominant terms, a simple pole with arbitrary leading coefficient, and resonances at -1 and 0. That means a two-parameter Laurent series exists for all α and β. The coefficient c3 = (β - α²/9)/(20a) is merely fixed by recursion; ARS does not require non-resonant coefficients to vanish. So the vanishing of the higher coefficients at α² = 9β has no bearing on whether the test is passed. Worse, a pole of y is a zero of x, which is an ordinary point of (2); the singularity analysis of (19) is not examining movable poles of the original equation. The inference in the Conclusion is therefore a non sequitur. The earlier caveat that the particular case α² = 9β is not being considered also sits oddly with using that same case as the target of the iff.

Other softness: equation (8) is misprinted, q in (11) is never defined (it must be β), and the abstract's closing phrase is hasty. These are minor but do not help.

Who this is for: someone collecting the symmetry properties of this ODE family might find a convenient summary, but there is no new result here, and the paper's distinctive claim is wrong. I would not send it for review. If a specialty venue insists on looking, one referee would catch the ARS point quickly; but as a desk decision, reject.

Referee Report

3 major / 5 minor

Summary. The manuscript studies the generalized Painlevé–Ince equation (2) via Lie point symmetry analysis and ARS singularity analysis. It reports that for generic α and β the equation possesses two Lie point symmetries, while for β = α²/9 it has eight symmetries, and it writes the reduction to a first-order algebraic equation. In the singularity analysis it finds a simple pole with two possible leading-order coefficients and computes the corresponding resonances. In Section 4 the authors introduce the transformation x = 1/y and, based on the vanishing of higher-order coefficients in the Laurent expansion of the transformed equation (19), conclude that (2) passes the Painlevé test under this transformation if and only if α² = 9β.

Significance. The claimed iff would be a useful criterion, since it would single out the classical Painlevé–Ince parameter from a two-parameter family. The paper also contains a correct and clearly presented computation of the Lie symmetries, including the maximal-symmetry case, and a correct leading-order balance for the original equation. The discussion of negative resonances connects to a body of literature in which the authors have participated. However, the central new result is not established: the argument in Section 4 does not follow from the ARS analysis of (19), and the transformation x=1/y does not map the singularity problem of (2) to that of (19).

major comments (3)

[Section 4, Eq. (19) and concluding sentence] The inference that (2) passes the Painlevé test under the coordinate transformation (18) if and only if α²=9β is not supported by the preceding analysis. For equation (19), the leading-order analysis gives a simple pole with an arbitrary coefficient, and the resonances are -1 and 0 for all values of α and β; there is no resonance at which a compatibility condition could single out α²=9β. The vanishing of the higher-order coefficient at α²=9β is a property of a coefficient at a non-resonant exponent, and in the ARS algorithm such coefficients are determined by the equation rather than required to vanish. Hence (19) passes the Painlevé test for every parameter pair, and the observed vanishing cannot be used to establish an iff for (2). Moreover, x=1/y maps poles of x (the movable singularities of (2)) to zeros of y; the Laurent expansion of (19) about a pole of y therefore describes the behaviour of x at an ordinary point, so this analysis does not probe the singularity structure of (2).
[Section 4] The sentence 'Naturally we are not considering the particular case in which α²=9β' is in direct contradiction with the concluding claim that the Painlevé test is passed under (18) if and only if α²=9β. If the special case is excluded from the transformation analysis, the conclusion cannot state anything about that case; if it is included, the sentence is contradicted. As written, the argument is either circular or vacuous for the only case it claims to detect.
[Section 3, Eqs. (16)-(17)] The resonance formulas (16) and (17) are misprinted and cannot be checked as written; they contain ambiguous expressions such as 'α sqrt(α²-8β) - α² + 8β' without proper parentheses. For the original equation (2) the second resonance is s = 4 - α a with a given in (14), a standard result that the paper does not state. Using this formula, integer resonances occur for parameter values other than β=α²/9; for example β=0 gives s=2 and β=-α² gives s=3 on one branch. Thus the singularity analysis of (2) itself does not single out β=α²/9, and the paper's classification of integrable cases is incomplete.

minor comments (5)

[Eq. (8)] The reduced equation is printed as 'vv′3 = 0', which is not a valid differential equation; it should be v v' + α u v + β u^3 = 0 in the variables (6)-(7).
[Eq. (11)] The quantity q is undefined; from the reduction it should be q = β.
[Abstract] The phrase 'integrable is terms of Lie symmetries' should read 'integrable in terms'.
[Section 4] The statement that 'the coefficients of the higher-order terms in the expansion vanish' should be made precise; for instance, the coefficient of τ² in the Laurent expansion of y is c3 = (β - α²/9)/(20a), which vanishes exactly when α²=9β.
[Acknowledgements and references] The acknowledgements contain 'Surananee University o Technology' (missing 'f') and the reference list has 'Mubarakzyanov GM !963' (should be 1963).

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the derivation is self-contained; the Section 4 iff is logically unsupported but is not a reduction of the conclusion to its inputs.

full rationale

The paper's main computations—the Lie symmetry classification in Section 2 and the Laurent expansions in Sections 3 and 4—are performed directly on equations (2) and (19) with no fitted parameters and no quantity defined in terms of the target result. The condition β = α²/9 is first obtained in equation (5) from the symmetry algebra before any singularity analysis, so the later observation that the higher-order coefficients in the expansion of (19) vanish at that same value is not a case of importing the conclusion into the premise. The self-citations [5,14,15] concern the interpretation of negative resonances; this material is background context and is not needed for the claimed iff, which rests on the expansion of (19). The actual defect in the paper is a logical one: from the fact that coefficients vanish at α²=9β the authors assert an iff, and they do not justify transferring the Painlevé property of (19) back to (2) under x = 1/y. That is a correctness/soundness objection, not circularity. Therefore no circular step can be exhibited under the standard of quoting a reduction of a prediction to its input.

Assumptions & free parameters 0 free parameters · 3 assumptions · 0 invented entities

The paper introduces no new fitted parameters or entities. It uses the equation's constants α and β as given, and the special relation β = α²/9 is a known integrability condition, not an ad hoc fit.

assumptions (3)

domain assumption The SYM Mathematica add-on correctly computes the Lie point symmetries of (2).
The paper uses the commercial CAS add-on SYM without independent verification; the two symmetries are plausible and can be checked by hand.
standard math The ARS algorithm (Ablowitz-Ramani-Segur) is a valid test for the Painlevé property.
The paper relies on the standard ARS procedure for singularity analysis as referenced in [1-3].
ad hoc to paper The change of variable x = 1/y preserves the property relevant to the Painlevé test.
The paper infers the original equation's Painlevé property from the transformed equation's resonances without proving invariance of the property under this transformation.

how reviews work

0 comments

Cite this review

Pith. "Pith review of A Systematic Analysis of the Properties of the Generalised Painlev\'e--Ince Equation." pith.science (2026). https://pith.science/paper/CHI7BBJS

@misc{pith2026190804563,
  author       = {Pith},
  title        = {Pith review of: A Systematic Analysis of the Properties of the Generalised Painlev\'e--Ince Equation},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/CHI7BBJS}},
  note         = {Machine review of arXiv:1908.04563}
}

abstract

We consider the generalized Painlev\'e--Ince equation, \begin{equation*} \ddot{x}+\alpha x\dot{x}+\beta x^{3}=0 \end{equation*} and we perform a detailed study in terms of symmetry analysis and of the singularity analysis. When the free parameters are related as $\beta =\alpha ^{2}/9~$the given differential equation is maximally symmetric and well-known that it pass the Painlev\'{e} test. For arbitrary parameters we find that there exists only two Lie point symmetries which can be used to reduce the differential equation into an algebraic equation. However, the generalized Painlev\'{e}--Ince equation fails at the Painlev\'{e} test, except if we apply the singularity analysis for the new second-order differential equation which follows from the change of variable $x=1/y.$ We conclude that the Painlev\'{e}--Ince equation is integrable is terms of Lie symmetries and of the Painlev\'{e} test.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

21 extracted references · 21 canonical work pages

[1]

Ablowitz M J, Ramani A & Segur H (1978) Nonlinear Evolution Equatio ns and Ordinary Diﬀerential Equations of Painlev´ e TypeLett Nuovo Cimento 23 333-337

work page 1978
[2]

Ablowitz M J, Ramani A & Segur H (1980) A connection between non linear evolution equations and ordinary diﬀerential equations of P type I J Math Phys 21 715-721

work page 1980
[3]

Ablowitz M J, Ramani A & Segur H (1980) A connection between non linear evolution equations and ordinary diﬀerential equations of P type II J Math Phys 21 1006-1015

work page 1980
[4]

On the sy stematic approach to the clas- siﬁcation of diﬀerential equations by group theoretical methods

Andriopoulos K, Dimas S, Leach PGL & Tsoubelis D (2009). On the sy stematic approach to the clas- siﬁcation of diﬀerential equations by group theoretical methods. Journal of Computational and Applied Mathematics, 230(1), 224-232

work page 2009
[5]

Andriopoulos K & Leach PGL (2006) An interpretation of the pres ence of both positive and negative nongeneric resonances in the singularity analysis Physics Letters A 359 199-203 4

work page 2006
[6]

Andriopoulos K, Leach PGL and Maharaj A (2011) On Diﬀerential S equences Applied Mathematics and Information Sciences 5(3) 484-499

work page 2011
[7]

SYM: A new symmetry-ﬁnd ing package for Mathematica

Dimas S & Tsoubelis D (2004, October). SYM: A new symmetry-ﬁnd ing package for Mathematica. In Proceedings of the 10th International Conference in Modern Gro up Analysis (Unmiversity of Cyprus, pp. 64-70)

work page 2004
[8]

A new Mathematica-based pro gram for solving overdetermined systems of PDEs

Dimas S & Tsoubelis D (2006, June). A new Mathematica-based pro gram for solving overdetermined systems of PDEs. In 8th International Mathematica Symposium (A vignon)

work page 2006

Show all 21 references

[9]

Partial Diﬀerential Equations, Algebraic Computing and N onlinear Systems

Dimas S “Partial Diﬀerential Equations, Algebraic Computing and N onlinear Systems.” PhD thesis, Uni- versity of Patras, Greece (2008)

2008
[10]

Euler M, Euler N and Leach PGL (2007) The Riccati and Ermakov- Pinney hierarchies Journal of Nonlinear Mathematical Physics 14 290-310

2007
[11]

Euler N and Leach PGL (2009) Aspects of proper diﬀerential se quences of ordinary diﬀerential equations Theoretical and Mathematical Physics 159 473-486 (0040-5779/09/1591-0473)

2009
[12]

Euler N and Leach PGL (2009) A novel Riccati Sequence Journal of Nonlinear Mathematical Physics 16 s01 157-164

2009
[13]

Euler M, Euler N and Leach PGL (2011) Properties of the Caloger o-Degasperis-Ibragimov-Shabat Diﬀer- ential Sequence Lobachevskii Journal of Mathematics 32 61-70

2011
[14]

Feix MR, G´ eronimi C, Cair´ o L, Leach PGL, Lemmer RL and Bouqu et S ´E (1997) On the singularity analysis of ordinary diﬀerential equations invariant under time tran slation and rescaling Journal of Physics A: Mathematical and General 30 7437-7461

1997
[15]

Lemmer RL & Leach PGL (1993) The Painlev´ e test, hidden symme tries and the equation y′′ + yy ′3 = 0 J Phys A: Math Gen 26 5017-5024

1993
[16]

Mahomed FM and Leach PGL (1985) The linear symmetries of a non linear diﬀerential equation Quæstiones Mathematicæ 8 241-274

1985
[17]

Morozov VV (1958) Classiﬁcation of six-dimensional nilpotent Lie algebras Izvestia Vysshikh Uchebn Za- vendeni ˘ ı Matematika5, 161-171

1958
[18]

Mubarakzyanov GM (1963) On solvable Lie algebras Izvestia Vysshikh Uchebn Zavendeni ˘ ı Matematika32 114-123

1963
[19]

Mubarakzyanov GM !963) Classiﬁcation of real structures of ﬁ ve-dimensional Lie algebras Izvestia Vysshikh Uchebn Zavendeni ˘ ı Matematika34 99-106
[20]

Mubarakzyanov GM (1963) Classiﬁcation of solvable six-dimensio nal Lie algebras with one nilpotent base element Izvestia Vysshikh Uchebn Zavendeni ˘ ı Matematika35 104-116

1963
[21]

Paliathanasis A & Leach PGL (2016) Nonlinear Ordinary Diﬀerentia l Equations: A discussion on Symme- tries and Singularities International Journal of Geometric Methods in Modern Physi cs 13 1630009 5

2016

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.

[1] [1]

Ablowitz M J, Ramani A & Segur H (1978) Nonlinear Evolution Equatio ns and Ordinary Diﬀerential Equations of Painlev´ e TypeLett Nuovo Cimento 23 333-337

work page 1978

[2] [2]

Ablowitz M J, Ramani A & Segur H (1980) A connection between non linear evolution equations and ordinary diﬀerential equations of P type I J Math Phys 21 715-721

work page 1980

[3] [3]

Ablowitz M J, Ramani A & Segur H (1980) A connection between non linear evolution equations and ordinary diﬀerential equations of P type II J Math Phys 21 1006-1015

work page 1980

[4] [4]

On the sy stematic approach to the clas- siﬁcation of diﬀerential equations by group theoretical methods

Andriopoulos K, Dimas S, Leach PGL & Tsoubelis D (2009). On the sy stematic approach to the clas- siﬁcation of diﬀerential equations by group theoretical methods. Journal of Computational and Applied Mathematics, 230(1), 224-232

work page 2009

[5] [5]

Andriopoulos K & Leach PGL (2006) An interpretation of the pres ence of both positive and negative nongeneric resonances in the singularity analysis Physics Letters A 359 199-203 4

work page 2006

[6] [6]

Andriopoulos K, Leach PGL and Maharaj A (2011) On Diﬀerential S equences Applied Mathematics and Information Sciences 5(3) 484-499

work page 2011

[7] [7]

SYM: A new symmetry-ﬁnd ing package for Mathematica

Dimas S & Tsoubelis D (2004, October). SYM: A new symmetry-ﬁnd ing package for Mathematica. In Proceedings of the 10th International Conference in Modern Gro up Analysis (Unmiversity of Cyprus, pp. 64-70)

work page 2004

[8] [8]

A new Mathematica-based pro gram for solving overdetermined systems of PDEs

Dimas S & Tsoubelis D (2006, June). A new Mathematica-based pro gram for solving overdetermined systems of PDEs. In 8th International Mathematica Symposium (A vignon)

work page 2006