Three New Classes of Solvable N-Body Problems of Goldfish Type with Many Arbitrary Coupling Constants

Abstract: Three new classes of N-body problems of goldfish type are identified, with N an arbitrary positive integer (N ě 2). These models are characterized by nonlinear Newtonian (“accelerations equal forces”) equations of motion describing N equal point-particles moving in the complex z-plane. These highly nonlinear equations feature many arbitrary coupling constants, yet they can be solved by algebraic operations. Some of these N-body problems are isochronous, their generic solutions being all completely periodic with an overall period T independent of the initial data (but quite a few of these solutions are actually periodic with smaller periods T{p with p a positive integer); other models are isochronous for an open region of initial data, while the motions for other initial data are not periodic, featuring instead scattering phenomena with some of the particles incoming from, or escaping to, infinity in the remote past or future.


Introduction
Recently, a new technique to identify many-body problems solvable by algebraic operations has been introduced [1,2], and several examples of such models have been discussed [1][2][3][4][5][6].In the present paper, three new classes of such models are introduced and discussed (the impatient reader may immediately glance at these findings reported in the next section).We term these models "of goldfish type" because their Newtonian ("accelerations equal forces") equations of motion read as follows: : z n " p2r `1q iω 9 z n `r pr `1q ω 2 z n ´N ÿ "1, ‰n " p 9 z n `irωz n q p 9 z `irωz q z n ´z with the nonlinear functions f n ˆ z, ¨ z ˙appropriately defined (see below).(The original goldfish model is the special case of these equations of motion with r " 0 and f n ˆ z, ¨ z ˙" 0; after its first identification as a solvable model [7], and its tentative recognition as a "goldfish" [8], this N-body problem and some of its extensions have been investigated in several publications (see, for instance, [9][10][11][12][13][14][15][16][17][18][19][20]).
Notation 1.1.Above, and hereafter, N is an arbitrary positive integer (N ě 2); indices such as n, m, run over the positive integers from 1 to N (unless otherwise indicated: see, for instance, the restriction on the values of in (1)); the N complex coordinates z n " z n ptq identify the positions of N points moving in the complex z-plane as functions of the time t; below, we also introduce other complex variables ζ n " ζ n pτq that depend on the complex variable τ related to the real variable t ("time") as follows: τ " τ ptq " exp piωtq ´1 iω .(2a) Above, and hereafter, i is the imaginary unit (i 2 " ´1), ω is an arbitrary (nonvanishing) real constant, so that τ ptq vanishes at the initial time t " 0 and is periodic in time with period T 0 , τ p0q " 0 , τ pt `T0 q " τ ptq , T 0 " 2π |ω| .
Above, and hereafter, we adopt the standard notation according to which superimposed dots denote differentiations with respect to the time variable t, so that, for instance, while primes appended to a function indicate differentiations with respect to the argument of that function (see, for instance, (3b) and (3c), where the relevant variable is of course τ).The auxiliary complex coordinates ζ n " ζ n pτq are related to the N complex coordinates z n " z n ptq-which are the main protagonists of our N-body problems, as they identify the positions of N equal unit-mass point-particles moving in the complex z-plane according to the Newtonian ("accelerations equal forces") equations of motion of type (1)-by the following relations: Above, and hereafter, the number r is required to be real and rational, where, of course, the numerator q and the denominator p are coprime integers, and, for definiteness, we hereafter assume that p is positive, p ě 1.
Finally, note that above we often omitted indicating explicitly the argument of functions, and we will do so throughout whenever this can be done without causing confusion.
In this paper, a key role is played by the following monic τ-dependent polynomial of degree N in the (complex) variable ζ, Above, and hereafter, the symbol ζ denotes the unordered set of the N complex numbers ζ n (and likewise z is the unordered set of the N complex numbers z n ), while γ (see (4a) and below) is the N-vector with components γ m (and likewise w is the N-vector with components w m ).The symbol ř N n 1 ,n 2 ,...,n m ˚denotes the sum from 1 to N over the m indices n 1 , n 2 , ..., n m with the restriction that these indices be all different among themselves, so that for N " 2, for N " 3, and so on.Let us also display the corresponding formulas for the τ-derivative of γ m " γ m pτq, where the symbol ř N n 1 ,n 2 ,...,n m ˚n˚d enotes the sum from 1 to N over the m indices n 1 , n 2 , ..., n m with the restriction that these indices be all different among themselves and all different from the index n, so that for N " 2, for N " 3, and so on.Analogous formulas hold, of course, for the relations among the N coefficients w m " w m ptq and the N zeros z n " z n ptq of the polynomial (4b), and their time-derivatives; for instance, Let us also report the relations-implied by (3a)-among the N coefficients w m ptq of the polynomial (4b) and the N coefficients γ m pτq of the polynomial (4a) with τ related to t via (2) and r defined by (3d), implying (i) Of course, the assignment of the N coefficients of a polynomial defines uniquely the corresponding unordered set of its N zeros, but generally it only allows to compute explicitly these N zeros for N ď 4. (ii) Moreover-and quite relevantly in our context (see below)-if a polynomial features a dependence on an additional variable (as, for instance, the dependence of the polynomial P N pz; z ptq ; w ptqq on the real variable t ("time") (see (4b)), then the unordered character of the set of its N zeros z n ptq is generally only relevant at one value of time, say at the "initial" time t " 0, since, at other values of time, the ordering gets generally determined by the natural requirement that the functions z n ptq evolve continuously over time.This prescription then fixes, for all time, the ordering of the zeros z n ptq-i.e., the assignment of the value n of its index to each zero z n ptq-as long as the coefficients w m ptq evolve themselves continuously and unambiguously over time and moreover no "collision" of two or more zeros occurs over the time evolution, i.e., for all time z n ptq ‰ z ptq if n ‰ (since clearly such collisions imply a loss of identity of the coinciding zeros).However, this identification requires an analysis of the time evolution of the N zeros z n ptq not only in the complex z-plane, but in fact over the N-sheeted Riemann surface associated to the N roots of the polynomial P N pz; z ptq ; w ptqq , and/or over the evolution of each coefficient w m ptq if its time evolution takes itself place on a Riemann surface (as it indeed happens in the cases discussed below).
A key formula for the identification and investigation of solvable Newtonian N-body problems is the following relation [1] among the τ-evolutions of the N zeros and the N coefficients of the monic polynomial p N pζ; τq (see 4a)) of degree N in its argument ζ and depending on the extra variable τ: In the present paper-which is a sequel to [6], so that we dispense below from some of the remarks reported there about the significance of the solvable many-body models treated herein-we focus on the N-body models that obtain via (9) for the coordinates ζ n pτq-and especially for the corresponding coordinates z n ptq (see ( 3))-when the coefficients γ m pτq evolve according to the following system of decoupled nonlinear Ordinary Differential Equations (ODEs): where the N "coupling constants" g m are arbitrary (possibly complex), while the two (sets of) parameters a m and b m are instead hereafter required to be real rational numbers and to be related to each other (and to the parameter r, (see (3d)) by the single relation: which-provided neither r nor 1 `mr vanish, i.e., r ‰ 0 and r ‰ ´1{m for m in its range from 1 to N, as we now assume (the case r " 0 shall be discussed below separately)-is necessary and sufficient to guarantee that the corresponding system of ODEs satisfied by the coefficients w m ptq be autonomous, reading as follows (as the diligent reader will easily verify via (10a) and ( 8)) : In fact, we will focus below only on certain specific assignments of the parameters r, a m , b m which allow the explicit solution of the ODEs (10a) in terms of elementary functions.The corresponding solvable N-body models satisfied by the coordinates z n " z n ptq are displayed-and their properties discussed-in the following Section 2, with the corresponding proofs provided in Section 3, while the special case with a m " 2 and b m " ´1 (and with g m an arbitrary rational number) is treated in [6].

Results
The Newtonian equations of motion of the first class of N-body models treated in this paper read as follows: with r an arbitrary rational number (see (3d); of course r ‰ ´1{m for m " 1, 2, ..., N if g m ‰ 0), ω an arbitrary nonvanishing real number, the N coupling constants g m arbitrary complex numbers (not all vanishing) and (see Notation 1.2) - We then assert that this N-body model is solvable by algebraic operations, its remarkable properties being detailed by the following Proposition 2.1.The N complex coordinates z n ptq providing the solution at time t of the initial-values problem of the Newtonian equations of motion (11) are the N zeros of the monic polynomial (4b), the N coefficients w m ptq of which are given, in terms of the initial data z n p0q , 9 z n p0q , by the following formulas: - These functions w m ptq are, of course, defined by continuity in t from their initial values w m p0q (i.e., this prescription identifies the determination of all the rational roots appearing in the above formulas), and note that, if g m " 0, the formula (12b) must be replaced by Let us also report the properties of these solutions, (12), when they are generated by generic initial data, z n p0q , 9 z n p0q, excluding the nongeneric initial data identified in Proposition 2.2 for which the system of evolution equations (11) runs into a singularity at a finite time.
These solutions are all nonsingular for all time, remaining in a finite region-the size of which depends on the initial data-of the complex z-plane and featuring no particle collisions, i.e., for all time t, z n ptq ‰ z ptq if n ‰ ; and they are all completely periodic with a period T which is an integer multiple of the basic period T 0 (see (2b)), z n pt `Tq " z n ptq , T " KT 0 , (13) with the positive integer K restricted as follows: p ď K ď p N!, where, of course, p is the denominator of the rational number r (see (3d)).Note the arbitrariness of the real (nonvanishing) rational number r, and especially of the N complex coupling constants g m .In addition, note that Proposition 2.1 implies that all generic solutions of the N-body model characterized by the Newtonian equations of motion (11a)-excluding the nongeneric solutions which are singular (see below Proposition 2.2)-are completely periodic with the same period T MAX " pp N!q T 0 .However, there are lots of solutions that are completely periodic with periods which are integer submultiples of T MAX .The detailed identification of these solutions and their periods is a nontrivial matter, as shown, for instance, by the discussion of this phenomenology in the paper [11]-that treats the "periodic goldfish model" (for this terminology, see [8]), which is in fact characterized by the same equations of motions (11a), but with all coupling constants vanishing, g m " 0-and by the detailed investigation of the structure of the Riemann surfaces associated with other analogous many-body models [21][22][23][24][25].
Proposition 2.2.The solutions of the Newtonian equations of motion (11) may feature singularities via two phenomena, both of which correspond to nongeneric initial data.
The first phenomenon is characterized by initial data satisfying-for at least one value m of the index m in its range from 1 to N such that g m ‰ 0 and ´m r is not a positive integer-the equality where τp1q m is defined in terms of the initial data z n p0q, 9 z n p0q as above (see (12c)).
This singularity occurs at the time t " ťp1q defined as follows: ťp1q " min " ťp1q m ı (with the minimum taken over the values of the index m satisfying the condition ( 14)), with where θ p1q m is defined mod r2πs as follows: Note that the fact that θ p1q m is real is implied by (14), and that certainly 0 ă ťp1q m ă T 0 , hence 0 ă ťp1q ă T 0 .
The second phenomenon causing the equations of motion (11a) to hit a singularity is the occurrence of a collision of two (or possibly more) particles at some time t " t c ; so that, for some indices n, with n ‰ (of course both in their interval from 1 to N) there holds the equality z n pt c q " z pt c q, causing the term z n ptq ´z ptq appearing in the denominators in the right-hand side of (11a) to vanish.This phenomenon corresponds to the fact that two of the zeros of the polynomial (4b) with (12) coincide, a fact that clearly only happens for nongeneric initial data (in the complex z-plane), although the condition on the initial data that would cause this phenomenon to happen-and the time at which this phenomenon would happen-can be computed explicitly only for small values of N (and even then the result is not very enlightening).
Let us complete our discussion of the first class of solvable N-body problems by displaying the equations of motion (11) in the (simplest) N " 2 case: with r an arbitrary rational number (r ‰ 0, r ‰ ´1{2, r ‰ ´1), and g 1 , g 2 two arbitrary complex numbers (not both vanishing).Note that, for r " ´2{3, both exponents on the right-hand side of these ODEs are integers, since then 1{ p1 `rq " 3, 1{ p1 `2rq " ´3, and this is as well the case for r " ´3{4, since then 1{ p1 `rq " 4, 1{ p1 `2rq " ´2.
Let us end our treatment of the first class of solvable N-body problems characterized by the Newtonian equations of motion (11) by pointing out that, for r " ´2{3 and N " 3, the exponents on the right-hand side of (11a) are all integers, indeed for r " ´2{3, the three exponents p2 `mrq { p1 `mrq take, for m " 1, 2, 3, the three values 4, ´2, 0.
The Newtonian equations of motion of the second class of solvable N-body models read as follows: with the various quantities defined as above (in particular 9 w m defined in terms of 9 z n and z n by (11c) with Notation 1.2).Note that this might be considered the special case of the first class of models (see above) with r " 0, a m " 2 and b m " 0, which was previously excluded because it requires a special treatment (see Section 3).
The solvable character of this N-body system is demonstrated by the following Proposition 2.3.The N complex coordinates z n ptq providing the solution at time t of the initial-values problem of the Newtonian equations of motion ( 17) are the N zeros of the monic polynomial (4b), the N coefficients w m ptq of which are given, in terms of the initial data z n p0q , 9 z n p0q, by the following formula: of course with τ ptq " rexp piωtq ´1s { piωq (see (2a)).
Here and below, the function ln is defined by continuity in t from its vanishing value at t " 0 where τ ptq " 0, and, of course, w m p0q, respectively, 9 w m p0q are defined in terms of the initial data z n p0q and 9 z n p0q by (12f), respectively, (12g), and if g m " 0 then (18) becomes w m ptq " w m p0q `9 w m p0q τ ptq, while if 9 w m p0q " 0 it yields w m ptq " w m p0q.
Clearly, these coefficients w m ptq are completely periodic in the time t with period T 0 , see (2a), w m pt `T0 q " w m ptq , iff the initial data satisfy the inequality iff instead the initial data satisfy the opposite inequality, they are periodic in t except for a constant shift (independent of the initial data!) over each period T 0 , so that w m pt `kT 0 q " w m ptq ´2πik g m , k " ˘1, ˘2, ˘3, ..., hence, in this second case, they diverge as t Ñ ˘8.In addition-in the intermediate, nongeneric case in which the initial data imply the equality -the coefficient w m ptq diverges at the finite times ťp2q m " pi ωq ´1 ln ´1 `iω τp2q m ¯mod pT 0 q , (21b) which are, of course, real thanks to (21a).Correspondingly, the particles coordinates z n ptq-being the N zeros of the polynomial (4b) with the coefficients w m ptq (see (18))-are periodic with period K T 0 -K being a positive integer in the range from 1 to N! (see below Remark 3.1)-iff the initial data satisfy the inequality (19b), while iff instead the initial data satisfy the opposite inequality (20a) at least one of the particle coordinates z n ptq comes from or escapes to infinity in the remote past and future: see, for instance, the relevant discussion in Appendix G ("Asymptotic behavior of the zeros of a polynomial whose coefficients diverge exponentially") of the book [9].In addition, of course, if the (nongeneric) initial data imply validity of the equality (21a), the equations of motion run into a singularity at t " ťp2q m (see (21b)).Other nongeneric initial data causing the equations of motion (17) to run into a singularity at a finite time are those leading to particle collisions.Note that generally the nongeneric initial data causing singularities are also those that separate the regions of initial data associated to different behaviors of the model, including the emergence of the higher periodicities associated to values of K larger than unity as well as the periodic and nonperiodic cases.
We do not display explicitly the equations of motion of this second model in the simplest N " 2 case because they can be immediately obtained by setting r " 0 in those of the first model (see ( 16)).
We conclude our report on the properties of the second class of solvable N-body problems with the following remark (proven in Section 3): Remark 2.1.The system of Newtonian equations of motion ( 17) is Hamiltonian, albeit with a time-dependent Hamiltonian.
The third class of solvable N-body problems is characterized by the Newtonian equations of motion where ω is an arbitrary nonvanishing real number, r is an arbitrary rational number (of course if g m ‰ 0 then r ‰ 0 and r ‰ 1{m for the integer m in its range from 1 to N), the N coupling constants g m are arbitrary complex numbers (not all vanishing), and the quantities w m " w m ptq , respectively, 9 w m " 9 w m ptq are expressed in terms of the coordinates z n " z n ptq and 9 z n " 9 z n ptq via (11b), respectively (11c).The solvable character of this N-body problem is demonstrated by the following Proposition 2.4.The N complex coordinates z n ptq providing the solution at time t of the initial-values problem of the Newtonian equations of motion (22) are the N zeros of the monic polynomial (4b), the N coefficients w m ptq of which are given, in terms of the initial data z n p0q , 9 z n p0q, by the following formulas: where w m p0q , respectively, 9 w m p0q , are defined in terms of the initial data z n p0q and 9 z n p0q by (12f), respectively, (12g), and of course the determinations of the rational roots are implied by the requirement that these formulas be valid at the initial time t " 0 and thereafter by continuity in t, and again if g m " 0, the expression (23b) must be replaced by (12h).
It is plain that, if neither mr nor 1{ pmrq are integers, for generic initial data the function γ m pτq features, in the complex τ-plane, two rational branch points at τ " τp3q m , respectively, at τ " τp3q m , and that the periodicity of γ m ptq " γ m pτ ptqq as function of t (see (23b)), is determined by the location in the complex τ-plane of these two branch points with respect to the circle C centered at τ " i{ω and with radius 1{ |ω| on which rotates the point τ ptq " rexp pi ω tq ´1s { pi ωq as function of the time t.If both branch points are located outside the circle C-and the condition on the initial data determining this outcome is clearly validity of both inequalities with τp3q m , respectively, τp3q m defined in terms of the initial data by (23d), respectively, (23e)-then clearly γ m ptq " γ m pτ ptqq as a function of t is periodic with period T 0 , see (2a), hence w m ptq is periodic in t with period T pm " pm T 0 , with pm " p{m if this number is integer; otherwise, pm " p (where of course p is the denominator of r: see (23a) and (3d)).In addition, as a consequence-if the inequalities (24a) are valid for all values of m in the range from 1 to N-then the polynomial (4b) is periodic in t with period where the MinimumCommonMultiple must be evaluated for the values of m " m such that g m ‰ 0; hence (see below Remark 3.1), the N coordinates z n ptq are periodic with period Q T, where Q is a positive integer in the range from 1 to N!.
It is also plain that the solutions z n ptq are periodic in t with a period which is a positive integer multiple Q of T 0 even if the initial data imply instead that some of the inequalities (24a) are reversed, but, in these cases, the determination of the outcome-in particular, of the value of Q-requires, to begin with, a standard analysis of the structure of the Riemann surface associated to the function of the complex variable τ, and consequently of the periodicity of ϕ m ptq " ϕ m pτ ptqq as a function of the real variable t ("time") when τ ptq " rexp piωtq ´1s { piωq so that τ ptq travels on the circle C with The next integration can be performed in terms of elementary functions only for special assignments (satisfying the restriction (10b)) of the parameters a m and b m , to which we restrict attention in the present paper.
Our first assignment is b m " 0, (29a) implying via (10b) and with A m , respectively, B p1q m defined by (28b) (or, equivalently, (12d)), respectively, (12e).Here, we assume, of course, that r ‰ 0 (implying a m ‰ 2) and moreover that r ‰ ´1{m with m " 1, 2, ..., N (implying that a m is a finite rational number for all values of m in its range from 1 to N).
This ODE (29c) can now be easily integrated, yielding (12b) with (12c).These developments clearly prove the first part of Proposition 2.1.
To prove the second part of Proposition 2.1 we ascertain, to begin with, the periodicity properties as functions of the time variable t of the coefficients γ m ptq " γ m pτ ptqq psee (12b)).The starting point is the observation that τ ptq-see (2) or (12a)-is a periodic function of t with period T 0 , rotating in the complex τ-plane on the circle C centered at the point i{ω and having radius 1{ |ω|.Hence, any holomorphic function of τ is as well periodic in t with period T 0 ; this clearly is (for all values of m in its range from 1 to N) the case of the functions γ m pτq p see (12b)), if r is a negative integer.If instead r is a positive integer, the functions γ m pτq are meromorphic in τ, featuring a polar singularity at τ " τp1q m p see (12b)).In this case, γ m ptq " γ m pτ ptqq is again generally periodic in t with period T 0 , γ m pt `T0 q " γ m ptq, but for the nongeneric assignments of the initial data such that τp1q m falls on the circle C-note that τp1q m does depend on the initial data (see (12c) and (12e)), and that the condition for this to happen is validity of the equality then the function γ m ptq " γ m pτ ptqq diverges at the times ťp1q m " pi ωq ´1 ln ´1 `iω τp1q m ¯mod pT 0 q .(30b) Finally, if r is rational but not integer, i.e., p ą 1 (see (3d)), and m is not an integer multiple of p, γ m pτq features a rational branch point at τ " τp1q m (see (12b)), then the evolution of γ m ptq " γ m pτ ptqq as function of the time t depends on the location of the branch point τp1q m in the complex τ-plane, whether it falls outside, inside or just on the circle C. The latter case requires that the initial data satisfy the condition (30a), implying again that they are not generic.While clearly the condition that the branch point be located outside the circle C is validity of the inequality and for the corresponding nongeneric initial data the function γ m ptq " γ m pτ ptqq is again periodic with period T 0 .If instead the branch point falls inside the circle C, and the condition on the initial data for this to happen is validity of the opposite inequality the periodicity of γ m ptq " γ m pτ ptqq gets modified: the period is then Tp " pT 0 (see (12b) and (3d)), unless p{m " pm is an integer, in which case the period is Ťp m " pm T 0 .
Next, let us discuss the t-periodicity of w ptq , related to γ ptq " γ pτ ptqq by (8a) or (12a).We then note that the prefactor exp pirmωtq in (8a) or (12a) is also periodic in t (see (3d)), with the same period Ťp m " pm T 0 , where pm " p unless p{m is an integer, in which case pm " p{m.We may therefore conclude that, for all generic initial data, the functions w m ptq are periodic with period Ťp m , w m `t `Ť pm ˘" w m ptq . (33) Next, let us discuss the periodicity of the zeros z n ptq of the polynomial (4b) with coefficients w m ptq periodic as indicated just above.In this context, the following Remark 3.1 is relevant.Remark 3.1.If a time-dependent polynomial P N pz; tq , of degree N in z, is time-periodic with period T, P N `z; t `T ˘" P N pz; tq, the unordered set z ptq of its N zeros z n ptq is of course periodic with the same period T, z `t `T ˘" z ptq (since after a period the polynomial is unchanged); however, due to the possibility that these zeros, as it were, "exchange their places" over their time evolution, the period of each individual zero z n ptq, considered as a continuous function of time, may be a positive integer multiple of T; although of course that multiple cannot exceed the number N! of permutations of the N elements of the unordered set z ptq (for a detailed discussion of this phenomenology in analogous many-body contexts see [11,[21][22][23][24][25]).
We can therefore conclude that the N zeros z n ptq of the monic polynomial (4b), the N coefficients w m ptq of which are periodic as described above, are also periodic, for generic initial data, with periods KT 0 , with the positive integer K restricted as follows: since the very definition of pm (see above) implies that p is the Minimum Common Multiple of the N parameters pm .This ends the proofs of the findings reported in Propositions 2.1 and 2.2.The proofs of the other two Propositions reported in Section 2 follow, but below we omit the proofs of some aspects of the results reported in Section 2, which are analogous to those detailed above in the formulations and proofs of Propositions 2.1 and 2.2.
Hence, the only aspect relevant to prove Proposition 2.3 that we do report is the derivation of the formula (18).The starting point is the ODE which is the special case of (10a) with r " 0, a m " 2 and b m " 0. The integration of this ODE is a trivial task, yielding γ m pτq " γ m p0q ´ln and the fact that this implies ( 18) is an immediate consequence of (8) with r " 0.
As for the proof of Remark 2.1, it is based on the observation that the decoupled nonlinear system of N ODEs : w m " iω 9 w m `g p 9 w m q 2 (36a) is Hamiltonian, since it is implied by the standard Hamiltonian equations of motion

Outlook
In this Section 4 we outline tersely possible extensions of the results obtained in this paper.A tentative list of such possible developments is provided in the last Section of [6]: the first of those reported there is essentially implemented in the present paper, while the others are still open.Also open is the possibility to investigate a combination of the three classes of models treated above, and of the class of models treated in [6], by assuming that, for different values of the parameter m, some coefficients γ m evolve according to the three different solvable subcases of the ODEs (28) discussed above in Section 3 or to the case discussed in [6].In addition, of course, also open is the investigation of the hierarchies of solvable models associated with each of those mentioned herein via the notion of generations of monic polynomials such that the coefficients of the polynomials of the next generation coincide with the zeros of a polynomial of the current generation [2].
Other N-body problems-more general than those treated in [6] and in the present paper-can be investigated by analogous techniques to those employed above, while allowing each γ m ptq to satisfy a solvable second-order ODE still belonging to the class (10) but forsaking the restriction on the parameters a m and b m necessary and sufficient to allow the solution of these ODEs in terms of elementary functions (for instance, a natural extension might include special functions of elliptic and hyperelliptic type).

Notation 1 . 2 .
which features the N zeros ζ n pτq and the N coefficients γ m pτq; and by its counterpart, P N pz; z ptq ; w ptqq " the N zeros z n ptq and the N coefficients w m ptq (see Notation 1.1).The N coefficients γ m " γ m pτq of the polynomial (4a) are of course expressed in terms of its N zeros ζ n " ζ n pτq via the formulas Two comments on the problem to determine the N zeros ζ n , respectively, z n of a monic polynomial of degree N in ζ, respectively, z from its N coefficients γ m , respectively, w m (see (4a), respectively, ( (23)e, finally, via (38a) with (39) and (8), we arrive at the formulas(23), thereby proving Proposition 2.4.