# Homogeneous-Like Generalized Cubic Systems

Authors: G. R. Nicklason
Journal: International Journal of Differential Equations
 (2016)

---

## Abstract

We consider properties and center conditions for plane polynomial systems of the formsx˙=-y-p1(x,y)-p2(x,y), y˙=x+q1(x,y)+q2(x,y) where p1, q1 and p2, q2 are polynomials of degrees n and 2n-1, respectively, for integers n≥2. We restrict our attention to those systems for which yp2(x,y)+xq2(x,y)=0. In this case the system can be transformed to a trigonometric Abel equation which is similar in form to the one obtained for homogeneous systems (p2=q2=0). From this we show that any center condition of a homogeneous system for a given n can be transformed to a center condition of the corresponding generalized cubic system and we use a similar idea to obtain center conditions for several other related systems. As in the case of the homogeneous system, these systems can also be transformed to Abel equations having rational coefficients and we briefly discuss an application of this to a particular Abel equation.


---

## Body

## 1. Introduction

In this work we consider differential polynomial systems in the plane having the form of a linear center perturbed by homogeneous polynomials of degreesn and 2n-1 where n≥2 is an integer. We refer to these asgeneralized cubic systems since they contain the cubic system (n=2) as a particular case. Specifically, we assume(1)dxdt=-y-p1x,y-p2x,y,dydt=x+q1x,y+q2x,y,where p1, q1 and p2, q2 are homogeneous polynomials of degrees n and 2n-1, respectively. We will also have occasion to consider the reduced (homogeneous) problem(2)dxdt=-y-p1x,y,dydt=x+q1x,y,in which the perturbation consists of a single polynomial. Corresponding to (1), (2) are the first-order differential equations: (3)dydx=-x+q1x,y+q2x,yy+p1x,y+p2x,y,(4)dydx=-x+q1x,yy+p1x,y. In his original work [1] Poincaré developed a method for determining if the origin is a center by seeking an analytic solution to the equation y′=Q/P where P, Q are polynomials satisfying P(0,0),Q(0,0)=0. For (4) it takes the form(5)Ux,y=12x2+y2+∑k=2∞Uαkx,y,where Uαk(x,y) is a homogeneous polynomial of degree αk=k(n-1)+2. This solution is required to satisfy the condition(6)dUdt=∂U∂xdxdt+∂U∂ydydt=∑k=1∞Vαkx2+y2αk≡∑l=2∞V~2lx2+y2l.Here, Vαk, V~2l are calledLyapunov coefficients and they are homogeneous polynomials in the coefficients of the system.Most of the known center conditions for systems of type (1) are for cubic systems and for these there is an extensive literature. We note particularly the works of Lloyd et al. [2–9], Żołądek [10], Alwash [11], Cherkas and Romanovski [12], and the many references therein. A great number of these results were obtained by an exhaustive analysis of the Lyapunov coefficients. In contrast the main results in this work were found by studying the differential equations for these systems and determining what we can learn from them. Our study will consider systems which satisfy the condition(7)yp2x,y+xq2x,y=0.The only complete set of center conditions [3] known for systems of this type is for the cubic system. We will show that all generalized cubic systems which satisfy (7) can be transformed to Abel differential equations of various types and from these we will deduce a probable set (at this time we do not consider the completeness of the center conditions because it would involve an extensive analysis of the Lyapunov coefficients) of complete center conditions for the quintic system and rederive those for the cubic system. Complete sets of center conditions are known for these two cases for systems which also satisfy yp1(x,y)+xq1(x,y)=0. These are discussed in [11, 13] where it is shown that the only center conditions possible for these systems are symmetric centers. In Section 3 we use trigonometric forms for the Lyapunov coefficients to obtain a simple rederivation of these conditions.For the purposes of this work an integrating factorμ of (3) is a function such that(8)∂∂yμx,yx+q1x,y+q2x,y-∂∂xμx,yy+p1x,y+p2x,y=0.We make use of the fact by Reeb [14] that if μ is analytic and nonzero on a neighborhood of the critical point (0,0) then the corresponding system is a center of (1) and (2). The only integrating factors that we mention in this paper are of this type.The Cartesian forms provided by (3) and (4) are difficult to work with so we consider various transformations of these equations. In particular it is well known that (4) can be transformed to an Abel equation of the first kind in which the coefficients are trigonometric polynomials in cos⁡θ, sin⁡θ. We show that the particular subcase (7) of (3) that we consider is also transformable to an Abel equation having similar form. By this equation we can easily relate the homogeneous form with its corresponding generalized cubic form and this leads directly to many of the general results that we will establish for these systems. We use it to show that any center condition for a homogeneous system can be transformed to a center condition of the corresponding generalized cubic system.In Section3 we present most of the main results and in Section 4 we apply these to various forms of generalized cubic systems to obtain several new center conditions for these systems. Sections 5 and 6 are devoted primarily to a brief discussion of simple closed invariant curves which can occur in these systems and by considering particular examples determine whether or not they are limit cycles. In the final section we present an example of a solvable Abel differential equation which cannot currently be solved by the Computer Algebra System Maple.


## 2. Derivation of Related Equations and Aspects of the Abel Differential Equation

Here we obtain the equations related to (3) and (4) which we will use in this work. Since this development is intended for both odd and even values of n, we need to pay close attention to the forms of the trigonometric polynomials which arise and to forms which various substitutions can take.Using a polar coordinate transformationx=rcos⁡θ, y=rsin⁡θ in (3) and (4), respectively, we obtain (9)drdθ=ξ1θrn+ξ2θr2n-11+η1θrn-1+η2θr2n-2,(10)drdθ=ξ1θrn1+η1θrn-1,where(11)ξiθ=sin⁡θqicos⁡θ,sin⁡θ-cos⁡θpicos⁡θ,sin⁡θ,ηiθ=sin⁡θpicos⁡θ,sin⁡θ+cos⁡θqicos⁡θ,sin⁡θfor i=1,2. Here ξ1, η1 and ξ2, η2 are homogeneous trigonometric polynomials of degrees n+1 and 2n, respectively, in cos⁡θ, sin⁡θ. If center conditions are defined in terms of these trigonometric polynomials, as is frequently the case, these equations can be used to define the corresponding coefficient functions p1,…,q2. Inverting expressions (11) for the case i=1 we obtain(12)p1x,y=-rncos⁡θξ1θ+rnsin⁡θη1θ=rn-1-xξ1θ+yη1θ.The degrees of ξ1, η1 are such that this expression might not be a polynomial, although only the highest degree terms n+1 will contribute to this possibility. Writing ξ1=ξ1H+ξ1~, η1=η1H+η1~ where ξ1H, η1H are the degree n+1 terms and ξ1~, η1~ contain all other terms, we see that this highest degree term can be written as(13)p1Hx,y=rn-1-xξ1Hθ+yη1Hθ=rn-1-xacos⁡n+1θ+bsin⁡n+1θ+rn-1yccos⁡n+1θ+dsin⁡n+1θ,where a, b, c, d are constants. For n=2 this is(14)p1Hx,y=-ax4+3b-cx3y-3a+dx2y2-b-3cxy3+dy4x2+y2.The remainder of this with respect to x is 4(b-c)xy3-4(a+d)y4 which will be 0 if d=-a, b=c. The same conditions ensure that q1H is a polynomial and also holds for the case n=3. We can now show in an inductive fashion that these conditions are sufficient for all n≥2. Suppose it is true for all n≤N where N≥3 is an integer. Then for n=N+1 we have(15)p1Hx,y=rN-xacos⁡N+2θ+bsin⁡N+2θ+rNyccos⁡N+2θ+dsin⁡N+2θ.Since the trigonometric polynomials are degree N+2, the denominator is of order rN+2 when expressed in Cartesian form. We can express this in terms of ξ1H for n=N by expanding the trigonometric functions and setting c=b, d=-a. We obtain(16)p1Hx,y=-rNxcos⁡θ+ysin⁡θacos⁡n+1θ+bsin⁡n+1θ,which is a homogeneous polynomial of degree N+1 when evaluated at cos⁡θ=x/r, sin⁡θ=y/r. The same consideration also hold for q1H and by obvious extension to the nonlinearities of degree 2n-1. Since (n+1)ξ1H+η1H′=(n+1)a+dcos⁡n+1θ+b-csin⁡n+1θ these conditions can be obtained by the vanishing of this expression. In general, the highest order terms must satisfy ((n-1)i+2)ξiH+ηiH′=0 for i=1,2 This basic restriction will apply to several of our results and it is sometimes satisfied by requiring that the highest order terms vanish. In the following we will simply indicate the highest allowable degrees of the trigonometric polynomials based on this condition. We will also drop the subscript 1 for the homogeneous system and simply refer to p, q, ξ, η.Equations (9) and (10) (with η1=η) can be further transformed in a number of ways. We continue by using the transformation given in [15] (17)ρθ=rθn-11+ηθrθn-1.With this (9) and (10), respectively, become (18)dρdθ=ψ4θρ4+ψ3θρ3+ψ2θρ21-η1θρ+η2θρ2,(19)dρdθ=-n-1ξθηθρ3+n-1ξθ-η′θρ2,where in (18) (20)ψ2θ=n-1ξ1θ-η1′θ,ψ3θ=n-1ξ2θ-2ξ1θη1θ+η1θη1′θ,ψ4θ=n-1ξ1θη12θ-ξ2θη1θ-η2θη1′θ. Condition (7) we are considering is equivalent to η2=0. In this case we can show that ψ4ρ2+ψ3ρ+ψ2=(1-η1ρ)((n-1)(ξ2-ξ1η1)ρ+ψ2) so (18) reduces to the trigonometric Abel equation(21)dρdθ=n-1ξ2θ-ξ1θη1θρ3+n-1ξ1θ-η1′θρ2.This equation is very interesting in the fact that it shares many of the same analytical properties as (19) but describes a much larger class of systems (1). It is the equation upon which most of the results obtained in this paper are based. Since η2=0, the maximum degree of ξ2 in (21) is 2n-2 whereas both ξ1, η1 can have maximal degree n+1. If we further set η1=0 which is the case discussed in [11, 13] for n=2,3 we obtain (22)dρdθ=n-1ξ2θρ3+n-1ξ1θρ2which is essentially the same equation as the corresponding form of (9).Since several of the classes of solutions which we obtain are based on (21), we briefly review certain properties of Abel equations. In this work we will consider both types of Abel equations. An Abel equation of the first kind has the form (23)dydx=f3xy3+f2xy2+f1xy+f0xand an Abel equation of the second kind has the form(24)dydx=f3xy3+f2xy2+f1xy+f0xg1xy+g0x,where the coefficient functions are assumed to be suitably differentiable functions of x. Form (24) can always be transformed to an Abel equation of the first kind by the variable change (25)yx=1g1xux-g0xg1x. For a general Abel equation of the first kind it is possible to define recursively an infinite sequence of relative invariants by [16] (26)s3x=f0xf32x+227f23x+13f3xf2′x-f2xf3′x-f1xf2xf3x,(27)s2k+1x=f3xs2k-1′x+2k-113f22x-f3′x-f1xf3xs2k-1xfor k≥2. From these, a sequence of absolute invariants can be formed. If the first invariant I1=s53/s35 is constant, the Abel equation can be transformed to a separable equation. This is the only general class of Abel equations which is integrable by quadrature. We note that if f0=f2=0 then s3=0 and the Abel equation is a Bernoulli equation. In cases more general than that just indicated, the vanishing of s3 means the Abel equation is transformable to a Bernoulli equation.A sufficient condition that an Abel equation of form (23) with f0=0 has a constant first invariant I1 is that the coefficient functions satisfy a relation of the form (28)f2xf3xf2x′+f1xf3x=Cf22x,where C is a constant. This gives (29)I1=7291-3C39C-22. In [17] we showed that a homogeneous system could be transformed to an Abel equation of the first kind having rational coefficients. This is also true for systems (4) which satisfy (7). Setting y=vx, interchanging the roles of x and v, and letting u=xn-1, we obtain an Abel equation of the second kind(30)dudv=-n-1p21,vu3+p11,vu2+vuRvu+v2+1,where R(v)=q1(1,v)+vp1(1,v). An Abel equation of the first kind is obtained by using (25). Renaming the variables this gives form (23) where (31)f3x=-n-1x2+13p21,xR2x+n-1x2+12p11,xRx-n-1xx2+1,f2x=3n-1x2+12p21,xR2x-2n-1x2+1p11,xRx+x2+1R′xRx+n-3x,f1x=-3n-1x2+1p21,xR2x+n-1p11,xRx-R′xRx,f0x=n-1p21,xR2x.In the following we will have cause to refer to arbitrary trigonometric polynomials. The basic form for these is(32)α0+∑k=1Nα2kcos⁡2kθ+β2ksin⁡2kθif n is odd and (33)∑k=1Nα2k-1cos⁡2k-1θ+β2k-1sin⁡2k-1θif n is even. The value of N depends upon n and will vary according to the circumstance.


## 3. Basic Results for Generalized Cubic Systems

In this section we develop some of the basic results which are standard for generalized cubic systems. We begin with the main result which shows that any center condition for a homogeneous system of degreen can be transformed into a center condition of the generalized cubic system having the same value of n. In this way we can truly think of the homogeneous systems as being nontrivial particular cases (p2,q2≠0) of the corresponding generalized cubic systems. We present several applications of the theorem and some of the following propositions in the next section.Theorem 1.
Letξ, η given by (11) define a center of (2) and (4) for some integer n≥2. Then this condition can be transformed to a center condition for the same n of the generalized cubic system defined by ξ1, η1, ξ2 with η2=0.Proof.
We will show that there exist generalized cubic systems which satisfy the differential equation (21) which is the same as (19) for the homogeneous system. The two equations will be the same if (n-1)ξ1-η1′=(n-1)ξ-η′ and ξ2-ξ1η1=-ξη. Set ξ1=ξ+ξ~ and η1=η+η~ where ξ~, η~ are trigonometric polynomials of degree not greater than n-1 defined by either (32) or (33) according to n being odd or even. The required identities will be satisfied if we take (34)ξ1θ=ξθ+ξ~θ=ξθ+1n-1η~′θ,η1θ=ηθ+η~θ,ξ2θ=ξθη~θ+1n-1η~′θηθ+1n-1η~′θη~θ.In this η~ can be arbitrarily chosen.A couple of points are worth mentioning regarding these results. Sinceξ, η are assumed to define a homogeneous polynomial system, the manner in which ξ1, η1 are defined ensures that they too define polynomials. Also, the maximum degree of ξ2 is 2n-2 but the terms ξη~, η~′η could be of degree 2n if the degree of ξ, η is n+1. However, in accordance with the definition of ξ~, η~ we can show that these highest degree terms will cancel leaving an expression for ξ2 having degree 2n-2.Another significant consequence of Theorem1 concerns the integrability of the systems involved. It follows directly that if the original homogeneous system is integrable, then so is the resulting generalized cubic system. This is because (19) and (21) are the same and (19) must be solvable if the homogeneous system is integrable.We now present several more results which help characterize the nature of certain generalized cubic systems. The first two are valid for general systems and the remainder for the caseη2=0.Proposition 2.
Letψ2, ψ3, ψ4 and η1, η2 be odd and even trigonometric polynomials, respectively, defined by (11) and (20) (or there exists a translation θ→θ+θ0 for which it is true). Then the solution of (18) is an even function of θ and the origin is a center for system (1) and (3).Proof.
It is straightforward to show thatρ(θ) and ρ(-θ) satisfy the same differential equation and the evenness of the solution gives ρ(-π)=ρ(π).The conditions of Proposition2 are clearly satisfied if ξ1, ξ2 are odd and η1, η2 are even. In this case the solution r(θ) of (9) is also an even function. These conditions definesymmetric ortime-reversible centers because the x-axis is a line of symmetry for the phase portrait.The next result shows that a center condition of a generalized cubic system for a givenn is also a center condition for any generalized cubic system with n1=n+2k where k≥1 is an integer.Proposition 3.
Letn1≥2 be an integer and n2>n1 an integer having the same parity as n1. Suppose that ξ1,n1, η1,n1, ξ2,n1, η2,n1 defined by (11) define a center for (1) and (3) for n1. Then ξ1,n2=(n1-1)/(n2-1)ξ1,n1, ξ2,n2=(n1-1)/(n2-1)ξ2,n1, η1,n2=η1,n1, η2,n2=η2,n1 define a center for n2.Proof.
Sincen1 gives a center the solution of (9) is 2π-periodic. Setting R=rn-1 in (9), we have for n=n2 form of it(35)n2-1ξ1,n2θ=n2-1n1-1n2-1ξ1,n1θ=n1-1ξ1,n1θwith a similar result for ξ2,n2. Thus the transformed equation is unchanged by the substitution and n2 system is also a center.The following is similar to Theorem1 in structure but instead relates to generalized cubic systems having a particular form. It is established in the same manner as Theorem 1.Proposition 4.
Letn≥2 be an integer and let ξ1, η1, ξ2 with η2=0 define a center condition of (1) and (21). Then the system defined by ξ1∗, η1∗, ξ2∗ with η2∗=0, where (36)ξ1∗θ=ξ1θ+ξ~θ=ξ1θ+1n-1η~′θ,η1∗θ=η1θ+η~θ,ξ2∗θ=ξ2θ+ξ1θη~θ+1n-1η~′θη1θ+1n-1η~′θη~θ,satisfies the same equation (21) as the original system and defines a center for (1). In this η~(θ) is a trigonometric polynomial of degree not greater than n-1 defined by either (32) or (33) according to n being odd or even.We will refer to any pair of systems which satisfy the conditions of Proposition4 asconjugate systems. Conjugate systems also exist for the homogeneous case but it is much more difficult to find them because we do not have the freedom allowed in generalized cubic systems by the simple definition of ξ2∗. We gave such a pair of homogeneous systems in [18], but unfortunately one of the systems was presented incorrectly. They should be (37)px,y=a+29n-1b2axn-1y,qx,y=axn+bxn-1y-n-1a-2n-29n-1b2axn-2y2-13n-2xn-3y3,px,y=b3n-1xn-a-29n-1b2axn-1y-2n-3b3n-1xn-2y2,qx,y=2n-3b3n-1xn-1y-n-2a-29n-1b2axn-2y2-n-22b3n-1xn-3y3for n≥2 and for arbitrary parameters a≠0 and b. These systems can be mapped to center conditions for the corresponding generalized cubic systems which will lead to solvable Abel equations (30) and (31), although we do not pursue this at this time. An example of such a calculation is given in the final section of the paper. It is also possible to find other conjugate homogeneous systems and, of those we know, they produce centers of generalized cubic systems characterized by part (2) of Proposition 8 for the particular case of quintic systems.The case whereη1=0 as well has been extensively studied and complete sets of center conditions [11, 13] are known for n=2,3 systems. In the following result we consider the general case of these systems and relate them to certain homogeneous systems.Proposition 5.
Letn≥4 be an integer and suppose that a generalized cubic system satisfies η1=η2=0. Then any center condition of the homogeneous system of degree n-2 is a center condition of this system.Proof.
In this case the system satisfies (22) where ξ1, ξ2 are arbitrary and have degrees n-1 and 2(n-1), respectively. The result is established by noting that we can simply select ξ1=((n-3)ξ-η′)/(n-1) and ξ2=-(n-3)ξη/(n-1), where ξ, η define a center of the homogeneous system.
n = 2,3 cases are not covered by this result; however, we can easily show that they can have only symmetric centers. For n=3 the most general form is given by (32) with ξ1(θ)=α2cos⁡2θ+β2sin⁡2θ+α0 and ξ2(θ)=γ4cos⁡4θ+δ4sin⁡4θ+γ2cos⁡2θ+δ2sin⁡2θ+γ0. Since there always exists a transformation θ→θ+θ0 such that α2cos⁡2(θ+θ0)+β2sin⁡2(θ+θ0) is odd, we can take α2=0. In [19] we used trigonometric integrals to calculate the first five Lyapunov coefficients of a certain homogeneous system. Adapting them to this system (ξ-η′/(n-1)→ξ1,-ξη→ξ2) we see that the first four Lyapunov coefficients will be zero if (38)∫02πξ1θdθ=∫02πξ2θdθ=∫02πξ¯1θξ2θdθ=∫02πξ¯1θ2ξ2θdθ=0,where ξ¯1′=ξ1. Evaluating these integrals and assuming ξ1≠0 (otherwise it is just a homogeneous system) we find sequentially that α0=γ0=γ2=γ4=0. This leaves ξ1(θ)=β2sin⁡2θ and ξ2(θ)=δ4sin⁡4θ+δ2sin⁡2θ and in view of Proposition 2 these are symmetric centers. The proof for n=2 case can be carried out in exactly the same fashion. We will say more about n≥4 cases in the next section.


## 4. Center Conditions for Generalized Cubic Systems Such Thatη2=0

It is generally accepted that two of the independent center conditions for homogeneous systems are the Hamiltonian and symmetric systems. On the basis of Theorem1 we will determine how each of these conditions transform to center conditions of generalized cubic systems which satisfy η2=0. The Hamiltonian condition for homogeneous systems is most easily given in terms of polar representation as ξ=-η′/(n+1). Applying the results of Theorem 1 and rescaling ρ transforms (21) to (39)dρdη=n2-14n2ηρ3+ρ2which satisfies (28) with a value C=(n2-1)/(4n2). Thus we always obtain a constant invariant equation with I1=729(n2+3)3/(4n2(n2-9)2) except in the case n=3 when C=2/9. In this last case we obtain a Bernoulli rather than constant invariant form.For homogeneous systems the basic condition for symmetric centers is given byξ, η being, respectively, odd and even. In this case the coefficients of ρ3, ρ2 in (19) are both odd and the solutions of both (10) and (19) are even functions. The condition r(-θ)=r(θ) shows that these solutions are symmetric with respect to x-axis. Application of Theorem 1 to this system gives a generalized cubic system in which the coefficient functions ξ1, η1, ξ2 are defined in terms of an arbitrary trigonometric polynomial of type (32) or (33) having maximum degree n-1. With ξ, η having the stated parity, even for n=2 or n=3 cases there is no reason to expect any of the coefficient functions of the generalized cubic system to be either even or odd. By construction (19) and (21) are still the same equation and have an even solution ρ(θ) which is the condition which guarantees that they are centers. On the other hand there is no reason to expect (9) to have an even solution (or a translated r(θ+θ0) even solution), so these are not generally symmetric centers. A similar situation also exists for homogeneous systems and we gave examples of these in [20].We now apply Theorem1 to the two homogeneous systems for which the center conditions are fully known, that is, n=2,3 cases. The general quadratic system can be written as (the Kapteyn form [21]) (40)px,y=a0x2+a1xy+a2y2,qx,y=b0x2+b1xy-b0y2and it is known that there are 4 independent center conditions for this system.Theorem 6.
The homogeneous quadratic system written in Kapteyn form has a center at the origin if and only if one of the following conditions is satisfied:(1)
b1-2a0=a1+2b0=0, the Hamiltonian system;(2)
a1=b0=0, the reversible or symmetric case;(3)
a1+2b0=b1+3a0+5a2=b02+2a22+a0a2=0, the critical case;(4)
a0+a2=0, the Lotka-Volterra case.The case of center conditions for the cubic system satisfyingη2=0 was considered in both [3, 22]. In each case 4 separate conditions are found and we can easily show that they correspond exactly to the systems obtained by applying Theorem 1 and Proposition 4 to the known conditions given in Theorem 6. In this case there is a one-to-one correspondence between the center conditions for the two systems and it gives rise to the question if this is true for all n. It is not directly true as there are several more parameters in the cubic system than in the homogeneous system. However, if we eliminate some of these parameters from the conditions for the cubic system we can recover the above conditions. It may be possible that some type of restricted converse of this type is true, but we also doubt that it will ever be proved. On the other hand we think it is a real possibility that we can use results from the generalized cubic system to find as yet unknown center conditions for the corresponding homogeneous problem. We also believe that it is quite likely that application of Theorem 1 and Proposition 4 to the homogeneous degree 3 system will produce a complete set of conditions for the quintic (n=3) system.Lunkevich, Sibiriskii, and Malkin [23, 24] have shown that the center conditions for the homogeneous system of degree 3 have three independent components.Theorem 7.
The general homogeneous cubic system can be written as(41)px,y=a0x3+a1x2y+a2xy2+a3y3,qx,y=b0x3+b1x2y+a1xy2+b3y3.This system has a center at the origin if and only if one of the following conditions is satisfied:(1)
3a1-b1=a2-3b3=0, the Hamiltonian system;(2)
a2-b1=a3-b0=a0-b3=0, the reversible or symmetric case;(3)
a2-b1+3(a0-b3)=a3-b0=5a0+2a2-b3=a1+3a3=0 and (3a2+b3)(a2+7b3)+100a32=0, the critical case.As with the case of the quadratic system, each of these conditions transforms directly to a center condition forn=3 quintic system. We are not aware of these having been previously given, so we formally present them below in Proposition 8. As we have seen the Hamiltonian case transforms to a Bernoulli equation with s3=0 as defined by (26) and the symmetric case transforms to a system in which the coefficients of (21) are odd. Since these ideas were developed with respect to a polar representation we will continue to present them in this manner. In order to do so we need to express each of the homogeneous systems in polar form. We can make the conditions of the third case solvable without radicals by setting 3a2+b3=-Ka3 and a2+7b3=100a3/K. From this we find (42)ξθ=a3sin⁡4θ-K2-100a320Kcos⁡4θ-K2+100a310Kcos⁡2θ,ηθ=a3cos⁡4θ+K2-100a320Ksin⁡4θ-K2+100a320Ksin⁡2θ.If a3=0 we can replace these forms by one of (43)ξθ=-a2cos⁡4θ+2a2cos⁡2θ,ηθ=a2sin⁡4θ+a2sin⁡2θor(44)ξθ=-b3cos⁡4θ-2b3cos⁡2θ,ηθ=b3sin⁡4θ-b3sin⁡2θ.Proposition 8.
If any of the following conditions are satisfied, then the corresponding system is a center ofn=3 (quintic) generalized cubic system defined by ξ1, η1, ξ2 with η2=0. In these ξ1=ξ+ξ~=ξ+η~′/2, η1=η+η~, where η~ is an arbitrary trigonometric polynomial having maximum degree 2 as defined by (32) and ξ, η define a center of the homogeneous system. The systems given by (1), (2)(a), and (3) are based on Theorem 1.(1)
Transformation of Hamiltonian condition is as follows: letη be an arbitrary trigonometric polynomial of degree 4 given by (32). Define ξ=-η′/4 and set (45)ξ2θ=-14η′θη~θ+12η~′θηθ+12η~′θη~θ.(2)
Transformation of symmetric condition is as follows:
(a)
Basic form: letξ, η be respectively odd and even arbitrary trigonometric polynomials of degree 4 defined by (32). Define ξ2 by (46)ξ2θ=ξθη~θ+12η~′θηθ+12η~′θη~θ.Then the coefficient functions-2ξη=2(ξ2-ξ1η1) and 2ξ-η′=2ξ1-η1′ of (19) and (21) are odd functions and the solutions of these equations are even functions such that ρ(-π)=ρ(π). These forms do not generally define symmetric centers unless further conditions are imposed on ξ1, η1.(b)
Alternate form: letξ, η be respectively odd and even arbitrary trigonometric polynomials of degree 4 defined by (32). Define ξ2 to be any trigonometric polynomial given by (32) of degree 4 such that its odd part (ξ2)o≠0 and its even part (ξ2)e is equal to the even part of ξ1η1. That is, (ξ2)e=(ξ1η1)e. Then the coefficient functions -2ξη, 2(ξ2-ξ1η1) and 2ξ-η′=2ξ1-η1′ of (19) and (21) are once again odd functions and the subsequent analysis is the same as in the previous part.(3)
Transformation of critical case is as follows: letξ, η be given by any one of (42)–(44) and define ξ2 as in part (a) of (2).We give Cartesian forms of the systems defined by (1) and (2)(a). For the Bernoulli case we rederived the result in a Cartesian format in order to obtain a simpler form in which some of the constants have been suitably redefined. The transformation of the Hamiltonian case forn=3 gives the following quintic system: (47)p1x,y=a0x3+a1x2y+3b3xy2+a3y3,q1x,y=b0x3+3a0x2y+b2xy2+b3y3,p2x,y=c0x5+c1x4y+c2x3y2+c3x2y3+c4xy4and q2(x,y)=-yp2(x,y)/x, where (48)c1=13a1-b2b0+132a1-3b0+b2c0a0,c2=a1-b2a0-3c0+3b3c0a0,c3=19a12+a1b2-2b22-13a1-3a3+2b2c0a0,c4=13a1-b2b3-b3c0a0.The symmetric form is developed from the forms ξ(θ)=c4sin⁡4θ+b2sin⁡2θ and η(θ)=c4cos⁡4θ+c2cos⁡2θ+c0 with η~ given by (32) with N=1. It can be expressed as(49)p1x,y=β2x3+a1x2y-3β2xy2+a3y3,q1x,y=a3+2α2+2c2x3-3β2x2y+a1-6α2+4b2-2c2xy2+β2y3,p2x,y=β2a3+2α2+2c2x5+c1x4y-2β23α2-2b2+c2x3y2+c3x2y3-β2a3xy4and q2(x,y)=-yp2(x,y)/x, where (50)a1=α0+3α2-2b2+c0+c2-3c4,a3=α0-α2+c0-c2+c4,c1=2α0α2-2b2+2c4+2α22-2b2-c0-c2+c4-2β22,c3=2α0α2-2b2-2c4-2α22+2b2+c0-c2-c4+2β22.For K≠0 the homogeneous system upon which case (3) of Proposition 8 is based has an integrating factor given by (51)μx,y=Aa33T16-30a32KT13T2+150K2a3Ax2+y2+40Kxy-20000K4a32T14-20Ka3T1T2+100K24,where T1=10x+Ky, T2=Kx+10y, A=K2+100. The corresponding integrating factor for the quintic system is a massive expression having the same basic form but with an additional factor which clearly arises from η~. It can be given as μ(x,y)=(α0+α2)x2+2β2xy+(α0-α2)y2+12P(x,y)/Q4(x,y) where P, Q are polynomials having degrees 6, 4, respectively, such that P(0,0),Q(0,0)≠0.The conditions provided by Theorems6 and 7 can also be applied directly to systems which satisfy η1=0 as well. By Theorem 1 and Proposition 5, Theorem 6 provides a set of 4 center conditions for the septic (n=4) system and Theorem 7 gives a set of 3 conditions for the nonic (n=5) system. Once again we think these conditions are probably complete, but we do not attempt to establish that herein. For systems satisfying (22) a type of converse can be found. That is, for general values of n and using arbitrary (i.e., not necessarily center producing) forms for the coefficient functions of (22), we can transform it to an equation of the form (21) in which ξ1, η1, ξ2 have the proper degrees n-1, n-1, and 2n-6, respectively. However, we are not able to satisfy the remaining conditions on ξ1, η1 such that they would define a polynomial system, so this form of (21) would instead produce a rational system. These additional conditions would depend upon the actual forms of the coefficient functions of (22) and in this regard it is clear that they must exist for any of the center cases which are a consequence of Theorem 1.


## 5. Constant Invariant Solutions and Limit Cycles

The origin of the modified system(52)dxdt=λx-y-p1x,y-p2x,y,dydt=x+λy+q1x,y+q2x,ywill be a focus if λ≠0 or a center or focus if λ=0. In [25] Giné and Llibre considered the problem for the case when η1=η2=0 for general values of n. With the help of (28) they constructed systems with either a center or a focus and from this determined the existence of certain limit cycles. In this and the following section we carry out the same type of analysis for our η2=0 systems and also extend it to certain more general cases for which η2≠0 as well. It is straightforward to show that (52) can be transformed to (53)drdθ=λr+ξ1θrn+ξ2θr2n-11+η1θrn-1+η2θr2n-2and if η2=0, the counterpart of (21) is(54)dρdθ=n-1ξ2θ-ξ1θη1θ+λη12θρ3+n-1ξ1θ-2λn-1η1θ-η1′θρ2+λn-1ρ.The form of (54) having η1=0 was used in [25] and it was shown that this equation has constant invariant solutions. We show that the general form of (54) has similar solutions which for particular values of λ and C subsume those in [25]. All the systems that we obtain in this section areDarboux integrable since they lead to constant invariant Abel equations.Using the coefficient functions of (54) we can write (28) as (55)dudθ+λn-1u=Cψθ,where u=f3/f2 and ψ=(n-1)ξ1-2λ(n-1)η1-η1′. Solving this equation, substituting for u, and isolating ξ2 give(56)ξ2θ=Ce-λn-1θn-1∫eλn-1θψθdθ+Aψθ-λη12θ+ξ1θη1θ.The integration constant A is nonzero only in the case for which λ=0 and n is odd. The result does produce a trigonometric polynomial and it does give constant invariant forms for (56); however, the only way to satisfy the maximum degree requirement of ξ2 being 2n-2 is to restrict the maximum degrees of ξ1, η1 to n-1.In [18] we showed that a constant invariant Abel equation (23) always has particular solutions of the form(57)yx=1f3xK′s3x1/3-f2x3,where s3 is given by (26) and K′ is a constant. In the case where f0=0 and (28) is satisfied we can write f3′f2-f3f2′+f1f2f3=Cf23 and s3 becomes (58)s3x=227-C3f23x.With this the particular solution is just (59)yx=Kf2xf3xfor some constant K. For (54) K is defined by the equation K2+K+C=0 which only has real solutions if C≤1/4. We can use the coefficient functions of (54) to find the particular solutions, but it is easier to use the equation obtained from it by using the usual rationalizing substitution z=tan⁡θ. There are no useful results when n is even, so we will restrict our attention to odd values by considering the case n=3. Also, the expressions obtained become very large so we will further restrict our attention to those systems in which C=-2 and λ=1. There is little loss of generality by making these specific choices, except for the values C=0,1/4 which must be dealt with separately. We will briefly consider these cases at the end of the section. Also, we do not consider λ=0 because these systems are centers. For the indicated choices there are two separate solutions arising from (59) given by (60)ρz=-4z2+1KC2a2+8a3+B-b3z2+8a2+2B+10b3z+2a2+7b0-3b1+b2+b3,where B=b0-b1-b2 and K/C=-1/2,1. Substituting for z will then give the appropriate forms for (54).The Cartesian form of the system which is obtained from (56) can be given as (61)p1x,y=a2-b1+b3x3+a3+b0-b2x2y+a2xy2+a3y3,q1x,y=b0x3+b1x2y+b2xy2+b3y3,p2x,y=c0x5+c1x4y+c2x3y2+c3x2y3+c4xy4and q2(x,y)=-yp2(x,y)/x, where (62)c0=143a2-2b1+3b32a2-3b1+b2+b3+8b02+1433a2-30b1+4b2+29b3b0,c1=1218a2-19b1-3b2+37b3a2+b12+52b1b2-212b1b3-12b22-b2b3+192b32-32b02+8a3b0+1235a2-3b1-14b2+40b3b0+2a2-3b1+b2+b3a3,c2=10a2-3b1-8b2+19b3a2+12b12+32b1b2+b22-192b2b3+192b32+19a3b0-122a2+b1+2b2+19b3b0+18a2-10b1-b2+19b3a3,c3=-12b1+3b2+17b3a2+12b1b2+12b1b3+12b22+b2b3-192b32+9a32-12b02-a3b0-12a2-b1b0+20a2-3b1-8b2+19b3a3,c4=-14a2+3b32a2+b0-b1-b2-b3+9a32+a3b0-b1+b2+8b3a3.Reversing the transformations which led to (54) and its equivalent form in terms of z gives us two invariant curves of the original system. From (60) and (61) we obtain a pair of conics C1, C2 having general form Akx2+Bkxy+Cky2+Dk=0 for k=1,2. The coefficients of these are linear expressions given in terms of the six arbitrary parameters a2, a3, b0, b1, b2, b3 of the system. We have (63)2a2+5b0-3b1+b2+b3x2+23a2+b0-b1-b2+4b3xy+2a2+6a3+b0-b1-b2-b3y2-2=0,2a2+11b0-3b1+b2+b3x2+26a2+b0-b1-b2+7b3xy+2a2+12a3+b0-b1-b2-b3y2+4=0.Since the determinant of the coefficient matrix of the linear system defined by the coefficients of these relations is nonzero, we can uniquely solve for any set of values {A1,B1,C1,A2,B2,C2} in terms of the system parameters. That is, any combination of hyperbolae, pairs of lines, circles, ellipses, or null conics is possible in this system. The situation is similar for any other choice of λ, C except for those values explicitly mentioned, so nothing is lost by considering the specific values λ=1, C=-2. Also, in determining the form of (63) we find that the transformation from (60) shows that x2+y2=0 is another invariant of these systems.A necessary condition that a curve be a limit cycle is that it is a closed curve encircling the critical point. For the invariant curves defined by (63) this means that we are dealing with ellipses or circles which have their centers at the origin. However, not all such curves are limit cycles because they are dependent upon the particular configuration. If one invariant curve (circle or ellipse) lies entirely inside another, then they are both limit cycles provided that no critical points of the system lie on either of the curves. However, if they intersect the flow pattern is markedly different with the points of intersection being saddle points (see (64) below and the results in the next section). In the following we give examples of both of these situations. The system (64)dxdt=x-y-49144x3+19x2y-37144xy2+19y3+13576x5+1144x4y+132x3y2+1144x2y3+5576xy4,dydt=x+y-536x3-53144x2y-536xy2-41144y3+13576x4y+1144x3y2+132x2y3+1144xy4+5576y5has a solution given by U(x,y)=C, where(65)Ux,y=3ln⁡x2+y2+6arctan⁡xy-ln⁡2x2+xy+y2-24-2ln⁡x2+y2-4.Since the circle x2+y2=4 lies entirely inside the ellipse 2x2+xy+y2=24 and the system has no other real-valued critical points, both are limit cycles. The origin is an unstable focus (λ=1) so the circle is stable and the ellipse is unstable as t→∞. The system (66)dxdt=x-y-516x3-316x2y-54xy2+34y3+164x5+564x4y+1764x3y2+564x2y3+14xy4,dydt=x+y-18x3-516x2y-1716xy2-54y3+164x4y+564x3y2+1764x2y3+564xy+14y5has the circle x2+y2=4 and the ellipse x2+16y2=16 as invariant curves. It has nine real-valued critical points; one at the origin and four others at the intersections of the invariant curves. These last four are saddle points which are given by (x∗,y∗)=(±4/5,±2/5) in which we take all combinations of the signs. Translating each of these points to the origin we find that the eigenvalues of the linear portion of the resulting systems are given by the two sets of values r=±1,∓2. The remaining four critical points are centers which are located in the four regions outside one of the curves and inside the other. They are given by (x′,y′)=(-50α3+14α,2α) where α is a root of 125x4-45x2+2=0. Translating each to the origin, we find that the eigenvalues of the linear parts satisfy one of the two relations r2+(123±1341)/50=0. Hence they are pure imaginary and a rotation and rescaling of the time will allow us to write the systems in standard form as X˙=-Y+⋯, Y˙=X+⋯. The critical points are centers because the systems are integrable.In presenting system (66) our original intention was to give a system in which the invariant curves intersect at some general set of points, but we found the complexity of the resulting behaviour somewhat surprising. For example, let D be the region containing the critical point at the origin. The boundary ∂D of D consists of four arcs, two from the ellipse and two from the circle, which meet at the four critical points (x∗,y∗). Clearly, these are points where a unique tangent fails to exist on ∂D. Any trajectory originating in D approaches ∂D in a fashion similar to that of a general limit cycle but is displaced from this boundary as it approaches one of the critical (saddle) points. Thus ∂D acts somewhat like a stable limit cycle for trajectories in D except that in this case it is not smooth. The behaviour in the four bounded regions outside D as well as in the unbounded region exterior to both curves is also quite interesting. We believe that it would be of interest to analyze this system (or ones like it) more fully than we can do here.The final case we present is the simple system(67)p1x,y=13R12R22x2+y23R12+R22x-2R12+R22y,q1x,y=-13R12R22x2+y22R12+R22x+3R12+R22y,p2x,y=-1R12R22xx2+y22,q2x,y=1R12R22yx2+y22.It has only a single critical point at the origin if R2≠R1 and has two circular limit cycles x2+y2=R12, x2+y2=R22. For R12=1 and R22=3 a somewhat different system than this, based on the values λ=-3, C=3/16 and the condition η1=0, was given in [23]. See also the results in the next section.The preceding discussion does not apply directly to the two casesC=0,1/4. If C=0 we can see from (56) that the coefficient of ρ3 in the trigonometric Abel equation (54) vanishes and the equation reduces to a Bernoulli equation. The roots of K2+K+C=0 are K=0,-1 and the main consequence of this is that there is only one conic C1 rather than the two which occur in the general case. For C=1/4 the roots are equal and the two conics coalesce, so we again have a situation in which there is only a single conic.


## 6. Concentric Circular Invariant Curves

In [26] Llibre and Rodríguez construct a vector field having an arbitrary distribution of circular limit cycles and show that this system has a Darboux first integral. The solution of (67) is very close to the general form given in [26] but with some minor differences in the coefficients. Here we consider somewhat more general conditions for which a generalized cubic system for which n being odd can have two isolated, concentric circular invariant curves by removing the condition that η2 should be zero. This will produce systems which have (at least) two isolated closed trajectories, but in many cases these are not limit cycles. For the forms that we consider these systems are transformable to Riccati equations which frequently can be solved in terms of special functions. Clearly, the easiest way to search for such solutions is to use a polar representation. Let ξ1, ξ2 be given by (32) and write ξ2=a+ξ2~, ξ1=b+ξ1~ where a, b are constants and ξ2~, ξ1~ have no constant parts. If r=R0≠0 is a constant solution of (53) then we must have λR0+R0nξ1+R02n-1ξ2=R0(λ+R0n-1(b+ξ1~)+R02(n-1)(a+ξ2~))=0. Since the trigonometric functions (as in (32)) are independent this splits into two equations: (68)aR02n-1+bR0n-1+λ=0,ξ1~θ+R0n-1ξ2~θ=0.If R1, R2 are distinct roots of the quadratic this gives a=λ/(R1R2)2 and b=-λ(R12+R22)/(R1R2)2. Then from the second equation of (68) we have (69)ξ2~θ=-1R1n-1ξ1~θ=-1R2n-1ξ1~θwhich can be satisfied only if ξ1~=ξ2~=0. Hence, in order to have two distinct solutions of this type we assume in the following that ξ2=λ/(R1R2)2, ξ1=-λ(R12+R22)/(R1R2)2. This can be seen to be the case for system (67). At this point we can choose η1, η2 arbitrarily subject to the necessary degree constraints. If we take η1=η2=0 and R1=1, R2=3, λ=-3 we obtain the system given in [25] mentioned in the last section.In what follows we will takeη1(θ)=c2cos⁡2θ+d2sin⁡2θ+c0 and η2(θ)=α2cos⁡2θ+β2sin⁡2θ+α0. These forms can be assumed for any odd value of n, but we will once again restrict our attention to the particular case n=3. With r2=R, (53) becomes(70)dRdθ=2λR12R22RR-R12R-R221+c2cos⁡2θ+d2sin⁡2θ+c0R+α2cos⁡2θ+β2sin⁡2θ+α0R2.Eliminating the trigonometric functions by setting z=tan⁡θ and then interchanging the roles of z and R lead to the Riccati equation (71)dzdR=R12R222λα0-α2R2+c0-c2R+1RR-R12R-R22z2+R12R22λβ2R+d2R-R12R-R22z+R12R222λα0+α2R2+c0+c2R+1RR-R12R-R22.The usual substitution y=-u′/(P(x)u) for y′=P(x)y2+Q(x)y+R(x) will then further convert this to a linear, second-order equation. We do not give the result of this conversion except to mention that it has certain similarities to the generalHeun differential equation [27] (72)y′′+γz+δz-1+ϵz-ay′+αβz-qzz-1z-ay=0and its confluent forms. These equations are very general and include the hypergeometric cases 2F1, 1F1, 0F1 as well as other classes of functions. The solutions of (71) are frequently expressible in terms ofHeun functions with R=0, R12, R22 reflecting the singularities which appear in (72). In addition to the solutions given in [25], it is possible for (71) to have other elementary solutions. For example, the pair of systems given by (73)η1θ=±R12-R22R12R22cos⁡2θ±λR12-R22R12R22sin⁡2θ-R12+R22R12R22,η2θ=∓R14-R242R12R22cos⁡2θ+R14+R242R12R22reduce the linearized form of (71) to the Euler equation z′′+z′/R-z/(4λ2R2)=0.Each of the systems defined by (71) has the property that the circles r=R1,R2 are invariant curves, although most of the solutions are non-Liouvillian. By considering certain numerical examples it can be seen that there are systems for which both, one, or neither of these circles is a limit cycle. As we saw in the cases of systems (64)–(67) the exact nature of the curves is dependent upon the existence and location of other critical points of the system.We have not considered the remaining possibility in (68) when the quadratic has equal roots. Instead we will look at the case of the center-focus form (λ=0) when the system has a single invariant circle r=R1. Similar calculations show that a relation, ξ2=-ξ1/R1n-1, of the same type as the second equation of (68) must be satisfied. Due to parity considerations, such an equality is only possible for odd values of n and shows that circular invariant curves centered at the origin cannot exist if n is even. Whether or not the circle is a limit cycle is once again dependent upon the structure of the system. The system (74)dxdt=-y+45x3-225x2y+85xy2+45y3-45xx2+y2x2-4xy-4y2,dydt=x+125x3-85x2y-145xy2-45y3-25x2+y23x3+8x2y-5xy2-2y3has x2+y2=1 as an invariant curve and in this case it is not a limit cycle. The origin is asymptotically stable (as t→∞) and apart from this the system has eight other real-valued critical points. Four of these lie on the circle with two of these being unstable nodes and two saddle points. The saddle points are given by (75)x,y=±3813±31426,±33826∓1413≈±0.906,±0.423and the nodes by(76)x,y=±3813∓31426,±33826±1413≈±0.042,±0.999.All trajectories inside the circle appear to emanate (as t→-∞) from one of the nodes or the other, but it does not seem to be possible to specify from which one a particular trajectory originates. As a trajectory approaches a neighbourhood of one of the saddle points, it is deflected away in a direction towards the origin. The remaining four critical points lie just outside the circle and in conjunction with those on the circle give rise to a fairly complex flow pattern in that region.


## 7. Application to Abel Differential Equations

Many physical systems can be converted to Abel differential equations and because of this it is always of interest when solvable Abel equations are encountered. In an earlier work [17] we showed that a homogeneous system can be transformed to such an equation (set p2=0, p1=p, q1=q in (30)) having rational coefficients and herein we have shown that η2=0 generalized cubic systems can also be similarly transformed. One consequence of this is that any integrable system of this type can be transformed to a solvable Abel equation. In the papers by Cheb-Terrab and Roche [16, 28] the authors clearly demonstrate the need for a convenient method of classifying Abel equations. Two Abel equations are said to belong to the sameequivalence class if one can be transformed into the other using a transformation of a specific type (see [28]). If this is possible then if one equation of a particular class is solvable so are all the other members of that class. Many of the ideas in these papers have been incorporated in the computer algebra system Maple. All symbolic computations in this paper were carried out in the most recent version Maple 2016 which has an excellent suite of routines for solving such equations; however, we encountered several equations which could not be solved. These include (21) and its corresponding rational form given by setting z=tan⁡θ as well as (30) and its first kind form defined by (31). In the following we will present an example of a system which produces Abel equations which were not solved by the software but which can be transformed to a solvable equation using a nonstandard transformation.One simple cubic system defined from case 4 of Theorem6 by the functions (77)ξθ=142a1-b2sin⁡3θ-14a2+2b1cos⁡3θ-142a1+b2sin⁡θ+14a2-2b1cos⁡θ,ηθ=14a2+2b1sin⁡3θ+142a1-b2cos⁡3θ+14a2-2b1sin⁡θ+142a1+b2cos⁡θand η~(θ)=α1cos⁡θ+β1sin⁡θ with a1=5/2,  a2=b1=-b2=3,  α1=5,  β1=9/2 produced a system having four invariant lines. The phase plane equation (3) is (78)dydx=-4x+30x2+48xy-30y2+75x2y+95xy2-75y34y-6x2+28xy+6y2-75x3-95x2y+75xy2.This is easily solved by Maple using symmetries and the resulting rational (Darboux) solution can be expressed as(79)Ux,y=5y+26621x-2y+23913x+18y+27710x+9y+2182=C.This cubic system produces an Abel equation (30) given by (80)dydx=-515x2-19x-15y3+23x2+14x-3y2+4xy23x3-x2+21x+15y+4x2+4and the corresponding homogeneous system produces the equation(81)dydx=-6x2+x-1y2-2xyx-12x+53x+1y-2x2-2,neither of which was solved by Maple. Since the transformation y=-2u/((9x+10)u-2) transforms (80) into (81) these equations belong to the same equivalence class. Moreover, the equation arising from (21) has the same first absolute invariant I1 (when suitably modified) as (80) and (81), so it too belongs to the same equivalence class. Thus failure of the algorithm to solve one means that the others would also not be solved. However, each of these equations is solvable since solutions can be obtained by suitably transforming (79).If we transform (81) by setting (82)x=-3v-5u+139v-2u+13,y=313v-239u+13we obtain the Abel equation (83)dvdu=-2333v2-29u-91v30uv-11u2+52u.Since (82) is not of a type which will preserve the equivalence class of the original equation (or even guarantee that the new equation is of Abel type), (83) belongs to a different class than (81). In the terminology of [28] it is of type AIA (Abel, Inverse-Abel) which means that if we interchange the roles of u and v we obtain a new equation having the same form, but of a different equivalence class. In the case of this example, both of these new equations are of type constant invariant and hence are solvable.The primary reason for presenting this particular example is because of its simplicity and clarity. In this study we have encountered systems whose first integrals are much more complicated than the one just given but will still generate solvable Abel equations. At this time we have not studied the nature of these new Abel equations, although it is highly unlikely that they are as simple as (83) with regard to solvability.



---

*Source: 7640340-2016-09-05.xml*