arxiv: v1 [math.ds] 2 Apr 2014

PERIODIC SOLUTIONS OF THE EXTENDED DUFFING-VAN DER POL OSCILLATOR arxiv:144.68v1 [math.ds] 2 Apr 214 RODRIGO D. EUZÉBIO1,2 AND JAUME LLIBRE 1 Abstract. In this paper some aspects on the periodic solutions of the extended Duffing-Van der Pol oscillator are discussed. Doing different rescaling of the variables and parameters of the system associated to the extended Duffing-Van der Pol oscillator we show that it can bifurcate one or three periodic solutions from a 2 dimensional manifold filled by periodic solutions of the referred system. For each rescaling we exhibit concrete values for which these bounds are reached. Beyond that we characterize the stability of some periodic solutions. Our approach is analytical and the results are obtained using averaging theory and some algebraic techniques. 1. Introduction 1.1. Setting the problem. A large number of non autonomous chaotic phenomena in physics, engineering, mechanics, biology among others are described by second-order differential systems of the form (1) ẍ = g(x,ẋ,t)+γ(t), where g(x,ẋ,t) is a continuous function and γ(t) is some external force. For instance, in biology system (1) models the FHN neuron oscillator and in the engineering system (1) it is a model to the horizontal platform system. The specific topic addressed in this paper concerns with other particular case of equation (1), namely, an extension of the forced Van der Pol equation with external excitation. Van der Pol s system also plays an important role in many applications in areas like engineering, biology, physics and seismology (see [2] and references therein). The system associated with the forced Van der Pol equation with external excitation is characterized by system (1) with functions g and γ in the form (2) g(x,ẋ,t) = ρ (1 x 2 )ẋ d dx V(x), γ(t) = δ cos(ωt), where ρ is the damping parameter, V(x) is the potential function and δ and ω are the amplitude and angular frequency of the driving force γ(t), respectively. We assume that ρ is non negative and δ and ω are positive. 21 Mathematics Subject Classification. Primary 34C7, 34C15, 34C25, 34C29, 37C6. Key words and phrases. Extended Duffing-Van der Pol oscillator, periodic solution, non autonomous systems, averaging theory. 1

2 R.D. EUZÉBIO AND J. LLIBRE The potential V(x) can be approximated by a finite Taylor expansion in series V 2 (x) = 1 2 ω2 x 2, V 4 (x) = 1 2 ω2 x 2 + 1 4 α x 4, V 6 (x) = 1 2 ω2 x 2 + 1 4 α x 4 + 1 6 λ x 6, where ω and λ are non zero and α is a real number. Almost all papers on forced excited Van der Pol systems deal with the potential V 2. However some papers concerning the potential V 4 have shown a lot of interesting behaviours (see [7], [12], [25], [27] and [28]). This case is usually refereed as Duffing Van der Pol oscillator or φ 4 Van der Pol oscillator. Nevertheless more recently some papers were published taking into account the potential V 6 meanly addressing with the problem of chaos control (see [13], [17] and [22]). The dynamics considering the potential V 6 is more complex and rich than the corresponding cases considering the potentials V 2 and V 4 (see [24] and [26]). This is the case that we will regard. This case is quoted in the literature as extended Duffing-Van der Pol oscillator or φ 6 Van der Pol oscillator. In this paper we will give an analytical treatment to system (1) in order to study its periodic solutions considering functions (2) and the potential V 6. Indeed, calling y = ẋ we obtain (3) ẋ = y, ẏ = ω 2 x+ρ y α x 3 ρ x 2 y λ x 5 +δ cos(ωt). System (3) becomes simpler if we perform a rescaling s = ωt in the time t and another one y = ω Y in the spatial variable y. In fact, calling again the new time s by t and the variable Y by y, after the rescaling we have (4) ẋ = y, ẏ = x+ρy αx 3 ρx 2 y λx 5 +δcos t, wherethenewparametersρ, α, λandδ aretherespectiveoldonesdividedby ω 2. The study of the periodic solutions of the non linear non autonomous 2π periodic differential system (4) will be the objective of this paper. When α, λ and δ are zero equation (4) is referred as unforced Van der Pol equation and has a unique stable periodic solution for ρ positive. Furthermore, if ρ is large this periodic solution remains and it describes a periodic oscillatory behaviour called relaxation oscillation. On the other hand, by considering δ non zero, system (4) can present nonlinear attractors from simple periodic solutions up to multiperiodicity and chaos. In [1] we can find a didactical discussion about these objects and some results using classical techniques are obtained for a special case of system (4) using particular values of ρ, δ, ω and initial conditions. There exists an exhaustive list of papers in the literature studying the properties of system (4) when δ is zero. For δ positive many open questions remain mainly due to the difficulty of integrate the system. In particular, in [19] the authors study some aspects of robust practical synchronization

PERIODIC SOLUTIONS OF THE φ 6 -VAN DER POL OSCILLATOR 3 for system (1) and apply the results to a particular case of system (4). In the same direction, in [14] the author investigate synchronization processes between chaotic attractors considering α = λ = in system (4). In [2] and [11] it is studied many aspects of system (4) for the case where α and λ are zero and ρ is large. Furthermore in [8] the authors provide sufficient conditions to the existence of one periodic solution under some analytical hypotheses for a general forced Van der Pol system considering ρ equal to zero. Periodic solutions for system (4) were found in [17], where the authors investigate chaos control. Besides, in [5] and [6] were obtained results on the existence of periodic solutions for an autonomous special case of the extended Duffing-Van der Pol oscillator. In [16], [3] and [31] we also can find some results on systems similar to system (4). In this paper we are concerned with periodic solutions of system (4). We will present sufficient conditions in order that this system possess one or three periodic solution and we will provide conditions on the parameters ρ, α, λ and δ for which these bounds are realizable. In Addition we will prove that it is not possible to obtain different bounds of periodic solutions using the methodology presented in this paper. The phase space of our non autonomous differential system is (x, y, t) and itspoincarémapisdefinedinthe(x,y) space. Whenpossiblethestability of the periodic solutions will be studied. As usual a periodic solution is stable if the eigenvalues of the fixed point associated to its Poincaré map have negative real part, otherwise there periodic solution is unstable. Inside the unstable periodic solutions the are two types. The unstable saddle periodic solutions having an eigenvalue with negative real part and the other with positive real part, and the repellor with both eigenvalues having negative real part. We should note that the method used here for studying periodic solutions can be applied to any periodic non autonomous differential system as done in [9] and [18]. In these papers the authors applied the method used in this paper in order to guarantee the existence of periodic solutions in a periodic FitzHugh Nagumo system and in the Vallis system, respectively. The paper is organized as follows. In subsection 1.2 the main results are stated and compared with other results. Also, some important points on the results are clarified. In subsection 2.1 we prove the results. In subsection 2.2 some aspects on the results are pointed out. Lastly, section 3 is devoted to give a brief summary of the results that we use from averaging theory. 1.2. Statement of the main results. In this section we present our results. We will behave under an analytical approach and for this reason the results are valid for a large range of the parameters of system (4), different from the major part of the results dealing with numerical techniques. In

4 R.D. EUZÉBIO AND J. LLIBRE addition, it is important to note that we are concerned with harmonic solutions in the sense that they do not bifurcate from periodic solutions with multiple period. We have the following results. Theorem 1. Consider ε > sufficiently small, (ρ,δ,α,λ) = (εr,εd,ε n 2 a, ε n 3 l) with n 2,n 3 > 1 and 81d 2 > 48r 2 >. Then system (4) has a 2π periodic solution (x(t, ε), y(t, ε)) such that (x(,ε),y(,ε)) (, 2(36r)1/3 3Γ + 1 3 ( ) 6 1/3 Γ), r ( when ε, where Γ = 9d+ ) 1/3. 81d 2 48r 2 Moreover, for r sufficiently small this periodic solution is stable. Theorem 2. Consider ε > sufficiently small, (α,δ,ρ,λ) = (εa,εd,ε n 1 r, ε n 3 l) with n 1,n 3 > 1 and α. Then system (4) has a 2π periodic solution (x(t,ε),y(t,ε)) such that ( (4d ) 1/3 (x(,ε),y(,ε)),), 3a when ε. Theorem 3. Consider ε > sufficiently small and (λ,δ,ρ,α) = (εl,εd,ε n 1 r, ε n 2 a) with n 1,n 2 > 1. Then system (4) has a 2π periodic solution (x(t,ε), y(t,ε)) such that ( (8d ) 1/5 (x(,ε),y(,ε)),), 5l when ε. We remark that the stability of the periodic solutions of Theorems 2 and 3 cannot be decided with the real part of the eigenvalues of their Poincaré map because these real part are zero. Theorem 4. Consider ε > sufficiently small and (x,y,ρ,δ,α,λ) = (εx, εy, εr,ε 2 d,ε n 2 a,ε n 3 l) with ρ. Then system (4) has a 2π periodic solution (x(t,ε),y(t,ε)) such that ( (x(,ε),y(,ε)) = O(ε), δ ) ερ +O(ε), when ε. Moreover this periodic solution is unstable.

PERIODIC SOLUTIONS OF THE φ 6 -VAN DER POL OSCILLATOR 5 Note that in Theorem 4 we have a periodic solution that comes from infinity. As far as we know this kind of behaviour also has not been observed in the papers concerned with system (4). Nevertheless this behaviour is common when we perform a rescaling in the parameters and variables as we made previously and can be observed also in papers [9] and [18]. Theorem 5. Consider ε > sufficiently small, (α, λ, δ, ρ) = (εa, εl, εd, ε n 1 r) with n 1 > 1, αλ < and D = 53747712a5 d 2 78125l 7 + 248d4 l 4. Then system (4) has a 2π periodic solution if D = and three 2π periodic solution if D <. Moreover, there are values of α, λ and δ that realize these number of periodic solutions. Theorem 6. Consider ε > sufficiently small, (ρ, α, δ, λ) = (εr, εa, εd, ε n 3 l) with n 3 > 1 and 324a 4 r 2 9a 2 d 2 r 2 +36a 2 r 4 d 2 r 4. Consider also the values 1 = 324a 4 +d 2 r 2 +9a 2 (d 2 4r 2 ), 2 = 2187a 4 d 4 +27d 4 r 4 16d 2 r 6 +18a 2 (27d 4 r 2 72d 2 r 4 +32r 6 ). The system rx ( 4+x 2 +y2 ) 3ay (x 2 +y2 ) =, 4d ry ( 4+x 2 +y2 ) 3ax (x 2 +y2 ) =, has one solution (x,y ) if 1 2 = or 2 > and three solutions (x i,yi ) if 2 < for i = 1,2,3, all of them satisfying 27α 2( (x i )2 +(y i )2) 2 +ρ 2 ( 4+(x i )2 +(y i )2)( 4+3(x i )2 +3(y i )2). Then if 1 2 = or 2 > system (4) has a 2π periodic solution (x (t,ε), y (t,ε)) such that (x (,ε),y (,ε)) (x,y ) when ε. Additionally if 2 < system (4) has three 2π periodic solution (x i (t,ε),y i (t,ε)) such that (x i (,ε),y i (,ε)) (x i,yi ) when ε for i = 1,2,3. Moreover there are values of ρ, α and δ for which these number of periodic solutions is reached.

6 R.D. EUZÉBIO AND J. LLIBRE Theorem 7. Consider ε > sufficiently small, (α,ρ,λ,δ) = (εa,εr,εl,εd), αρ, (r 2 3a 2 )(r 2 9a 2 ) and C = 4(3a+1l)(54a 3 l 9a 2 (d 2 6l 2 ) 12al(d 2 15l 2 )+ 5l 2 ( 7d 2 +4l 2 ))r 6 +6(a+5l)(6a d+2l)(6a+d+2l)r 8 ), D = 26624268a 14 d 6 3125d 8 r 12 531441a 12 (3125d 8 +96d 6 r 2 + 1536d 4 r 4 124d 2 r 6 )+354294a 1 (3125d 8 r 2 +616d 6 r 4 16 d 4 r 6 )+18a 2 r 1 (9375d 8 +5d 6 r 2 88d 4 r 4 +124d 2 r 6 32768r 8 )+2916a 6 r 6 (15625d 8 +117d 6 r 2 45632d 4 r 4 +23552 d 2 r 6 248r 8 ) 6561a 8 r 4 (46875d 8 +24832d 6 r 2 66816d 4 r 4 + 12288d 2 r 6 +496r 8 )+81a 4 r 8 ( 46875d 8 36d 6 r 2 +2752d 4 r 4 246784d 2 r 6 +6144r 8 ). The system 2rx ( 4+x 2 +y2 ) y (x 2 +y2 )(6a+5l(x2 +y2 ) =, 8d+8ry (x 2 +y2 )(6ax +2ry +5lx (x 2 +y2 )) =, has one solution (x,y ) if D < and three solutions (xi,yi ) if D > for i = 1,2,3, all of them satisfying 4ρ 2( 4+(x i )2 +(y i )2) ( 4+3(x i )2 +3(y i )2 )+((x i )2 +(y i )2 ) 2 (6α+ ρλ((x i )2 +(y i )2 ))(18α+25λ(x i )2 +3(y i )2 )). Then, if D < system (4) has a 2π periodic solution (x (t,ε),y (t,ε)) such that (x (,ε),y (,ε)) (x,y ) when ε. Furthermore, if D > system (4) has three 2π periodic solution (x i (t,ε),y i (t,ε)) such that (x i (,ε), y i (,ε)) (x i,yi ) when ε for i = 1,2,3. Moreover, there are values of α, ρ, λ and δ which realizes the number of these periodic solutions.

PERIODIC SOLUTIONS OF THE φ 6 -VAN DER POL OSCILLATOR 7 Theorem 8. Consider ε > sufficiently small, (ρ, λ, δ, α) = (εr, εl, εd, ε n 2 a) with n 2 > 1. Consider the numbers C = 4d 3 l 3 r 4 +72dl 5 r 4 +6d 3 lr 6 +72dl 3 r 6 +(5d 4 l 2 r 4 22d 2 l 4 r 4 8l 6 r 4 28d 2 l 2 r 6 +16l 4 r 6 6d 2 r 8 + 24l 2 r 8 ), N 2 = 45l2 r 4 +3r 6 3125d 2 l 4, N 3 = 5859375d 4 l 6 +1875d 2 l 6 r 2 +875l 4 ( 7d 2 +12l 2 )r 4 + 25l 2 ( 7d 2 +1556l 2 )r 6 +156l 2 r 8 +9r 1, N 4 = 1171875d 8 l 6 (125l 4 +85l 2 r 2 +r 4 ) 25d 6 l 4 (3125 l 8 +534375l 6 r 2 +74625l 4 r 4 +358375l 2 r 6 +69r 8 ) 3 d 4 l 2 r 2 (21875l 1 21625l 8 r 2 686625l 6 r 4 3715l 4 r 6 65l 2 r 8 3r 1 ) 16l 2 r 6 (4l 1 1449l 8 r 2 2817l 6 r 4 888l 4 r 6 +387l 2 r 8 + 27r 1 )+12d 2 r 4 ( 5l 12 3225l 1 r 2 6624l 8 r 4 283l 6 r 6 +541l 4 r 8 +192l 2 r 1 +9r 12 ) N 5 = 48828125d 8 l 6 46875d 6 l 4 r 4 +64l 2 r 1 (16l 2 +r 2 ) 16d 2 r 8 (2l 4 +29l 2 r 2 +r 4 )+d 4 r 6 (275l 4 +6l 2 r 2 +27r 4 ) M 5 = 25d 4 l2+d 2 (11l4+14l2r 2 +3r 4 )+4(1l6 2l 4 r 2 3l 2 r 4 ), where CM 5. The system 2rx ( 4+x 2 +y2 ) 5ly (x 2 +y2 )2 =, 8d 2ry ( 4+x 2 +y2 ) 5lx (x 2 +y2 )2 =, has one solution (x,y ) if N 2, N 3 or N 4 and N 5 > and three solutions (x i,yi ) if N 5 < for i = 1,2,3, all of them satisfying 125λ 2( (x i )2 +(y i )4) +4ρ 2 ( 4+(x i )2 +(y i )2 )( 4+3(x i )2 +3(y i )2 ). Hence, if N 2, N 3 or N 4 and N 5 >, then system (4) has a 2π periodic solution (x (t,ε),y (t,ε)) such that (x (,ε),y (,ε)) (x,y ) when ε. If N 5 < then system (4) has three 2π periodic solution (x i (t,ε),y i (t,ε)) such that (x i (,ε),y i (,ε)) (x i,yi ) when ε for i = 1,2,3. Moreover, there are values of ρ, λ and δ for which system (4) has one or three periodic solutions. We note that in Theorems 5, 6, 7 and 8 we do not say anything about the kind of stability of the periodic solutions because we do not have the explicit expressions of the real part of the eigenvalues of their Poincaré maps.

8 R.D. EUZÉBIO AND J. LLIBRE We remark that the periodic solutions provided in Theorems 1 until 8 exit when we take small values for the parameters of system (4) obeying some relations among them. For instance Theorem 1 states the existence of one periodic solution for system (4) if each parameter of this system is small and α and λ are much smaller than ρ and δ. This assertion becomes more clear if we observe the replacement done in each theorem and taking into account that ε is sufficiently small. As far as we know periodic solutions of system (4) whose parameters have this characteristic have not been observed in the literature. Besides, it seems that the simultaneously bifurcation of three harmonic periodic solutions in extended Duffing-Van der Pol system s is also new. Furthermore we note that the periodic solutions presented in the present paper are different from those ones stated in [8]. Indeed, in the referred paper the authors ask for V 6 () <, and in our case we have V 6 () =. 2. Proof and discussion of the results 2.1. Proof of the results. In order to apply averaging theory described in section 3 in systems (4) we start performing a rescaling of the variables x and y and of the parameters ρ, α, λ and δ as follows. (5) x = ε m 1 X, y = ε m 2 Y, ρ = ε n 1 r, α = ε n 2 a, λ = ε n 3 l, δ = ε n 4 d, where ε is positive and sufficiently small and m i and n j are non negative integers, for i = 1,2 and j = 1,2,3,4. We recall that since δ >, ρ and λ we have d >, r and l. (6) In the new variables (X,Y) system (4) writes dx dt dy dt = ε m 1+m 2 Y, = ε m 1 m 2 X +ε n 1 ry ε 2m 1+n 1 rx 2 Y ε 3m 1 m 2 +n 2 ax 3 ε 5m 1 m 2 +n 3 lx 5 +ε m 2+n 4 dcos t. In order to have non negative powers of ε we must impose the conditions (7) m 1 = m 2 = m, and n 4 m, where m is a non negative integer. Hence, with conditions (7) system (6) becomes (8) dx dt dy dt = Y, = X +ε n 1 ry ε 2m+n 1 rx 2 Y ε 2m+n 2 ax 3 ε 4m+n 3 lx 5 +ε m+n 4 dcos t.

PERIODIC SOLUTIONS OF THE φ 6 -VAN DER POL OSCILLATOR 9 In this paper we will find periodic solutions of system (8) depending on the parameters r, a, l and d and the powers m, n i of ε, for i = 1,2,3,4. Then we will go back through the rescaling (5) to ensure the existence of periodic solutions in system (4). From now on we assume that the values n 1, m 2 +n 2 2, m2 +n 2 3 and n 4 m are positive and observe that considering these conditions each power of ε in system (8) becomes positive. The reason for a such assumption will be explained later on in section 2.2. Now we shall apply the averaging theory described in section 3. Thus following the notation of the mentioned section and denoting again the variables (X,Y) by (x,y) we have x = (x,y) T and system (13) corresponding to system (8) can be written as (9) ẋ = F (t,x) = (y, x) T. So the solution x(t,z) = (x(t,z),y(t,z)) of system (9) such that x(,z) = z = (x,y ) is x(t,z) = x cos t+y sin t, y(t,z) = y cos t x sin t. Itisclear thattheoriginofcoordinatesofr 2 isaglobal isochronouscenter for system (9) whose circular periodic solutions starting on z = (x,y ) R 2 are parametrized by the above functions x(t,z) and y(t,z). Therefore through every initial condition (x,y ) in R 2 passes a 2π periodic solution of system (9). Physically speaking system (9) models a simple harmonic oscillator and its solutions (x(t, z), y(t, z)) describe the wave behaviour of this oscillator. In this direction the problem of perturbation of the global center (9) is equivalent to the problem of perturbation of a simple harmonic oscillator introducing a damping parameter, as external force and considering a potential V by taking small values for ρ, α, λ and δ as stated in Theorems from 1 to 8. We note that using the notation of section 3 the fundamental matrix Y(t,z) of system (9) satisfying that Y(,z) is the identity of R 2 writes ( ) cos t sin t Y(t,z) =. sin t cos t From the averaging theory we are interested in the simple zeros of the function (1) f(z) = (f 1 (z),f 2 (z)) = Y 1 (t,z)f 1 (t,x(t,z,))dt, where F 1 (t,x) are the terms of order 1 on ε of the vector field associated to system (8) (y, x+ε n 1 ry ε 2m+n 1 rx 2 y ε 2m+n 2 ax 3 ε 4m+n 3 lx 5 +ε m+n 4 dcos t) T. Now we prove our results.

1 R.D. EUZÉBIO AND J. LLIBRE Proof of Theorem 1: Firstly we take n 1 = n 4 = 1, m = and consider n 2,n 3 > 1. So the rescaling (5) becomes (x,y,ρ,δ,α,λ) = (X,Y,εr,εd,ε n 2 a, ε n 3 d) and calling again (X,Y) by (x,y) we verify the hypotheses of Theorem 1. Now the vector field of system (4) becomes (y, x+εry εrx 2 y ε n 2 ax 3 ε n 3 lx 5 +εdcos t) T. Hence as we stated before the terms of order 1 in this vector field are given by F 1 (t,x) = (,ry rx 2 y +dcos t). Then since z = (x,y ) the function f(x,y ) = (f 1 (x,y ),f 2 (x,y )) given in (1) turns into the form f 1 (x,y ) = (sin t)(dcos t+r(y cos t x sin t) r(y cos t x sin t)(x cos t+y sin t) 2 )dt = 1 4 πrx ( 4+x 2 +y 2 ), f 2 (x,y ) = (cos t)(dcos t+r(y cos t x sin t) r(y cos t x sin t)(x cos t+y sin t) 2 )dt = 1 ( ( 4d ry 4+x 2 4 +y 2 )). The zeros of functions f 1 and f 2 are a pair of conjugate complex vectors and a real pair (x,y ) satisfying x = and y = 2 (36r) 1/3 3 ( 9d+ ) 1/3 + 1 ( ) 6 1/3( 9d+ ) 1/3. 81d 81d 2 48r 2 3 r 2 48r 2 Note that y is real because 81d2 > 48r 2 > by assumption. We observe thattheexpressionofy = y (r,d) doesnotchangewhenwegobackthrough the rescaling (5) taking r = ε 1 ρ and d = ε 1 δ according to the hypotheses of Theorem 1. ( Moreover we denote Γ = 9d+ ) 1/3 81d 2 48r 2 and observe that the matrix M = (m ij ) = (f 1,f 2 )/ (x,y ) at (x,y ) is diagonal and its elements are m 11 = πrγ + /36 and m 22 = πrγ, where Γ ± = 12± 61/3 24r 2/3 Γ 2 ± 62/3 Γ 2 r 2/3. The determinant Π of M written in a power series of r around r = is Π = 3 (2π 2 d(4d) 1/3) r 2/3 + (πr)2 +O(r 7/3 ), 8 3 which is positive for r sufficiently small because d >. So the averaging theory described in section 3 guarantee the existence of a 2π periodic solution (x(t,ε),y(t,ε)) such that (x(,ε),y(,ε)) tends to (x,y ) when ε.

PERIODIC SOLUTIONS OF THE φ 6 -VAN DER POL OSCILLATOR 11 On the other hand the trace Σ of M at (x,y ) is Σ = 2 2/3 π(2d) 2/3 r 1/3 2πr 3 +O(r4/3 ). Thus for r sufficiently small the trace Σ of M is negative. So the real part of the eigenvalues of the matrix M are both negative and then the periodic solution (x(t,ε),y(t,ε)) is stable for each ε sufficiently small. So Theorem 1 is proved. Proof of Theorem 2: Now we take n 2 = n 4 = 1 and consider n 1,n 3 > 1. So we are under the hypotheses of Theorem 2. Now, taking m = and doing the rescaling (5) the vector field of system (8) writes (y, x+ε n 1 ry ε n 1 rx 2 y εax 3 ε n 3 lx 5 +εdcos t) T. Thus we have F 1 (t,x) = (, ax 3 + dcos t) and f(x,y ) = (f 1 (x,y ), f 2 (x,y )) becomes f 1 (x,y ) = (sin t)(dcos t a(x cos t+y sin t) 3 dt = 3 4 aπy ( x 2 +y 2 ), f 2 (x,y ) = (cos t)(dcos t a(x cos t+y sin t) 3 dt = 1 4 π( 4d 3ax ( x 2 +y 2 )). The real zero (x,y ) of f 1 and f 2 is ( (4d ) 1/3 (x,y,) ) =. 3a Again the values (x,y ) does not depends on ε when we go back to the original parameters α and δ through the rescaling (5). Besides the matrix M = (f 1,f 2 )/ (x,y ) at (x,y ) now writes ( ) M = 1 ((3/4)αδ 2 ) 1/3 π, ε 3((3/4)αδ 2 ) 1/3 π and then the determinant Π of M is ( Π = 3π2 9α 2 δ 4 ) 1/3 2ε 2 >. 2 Therefore the averaging theory implies the existence of a 2π periodic solution (x(t,ε),y(t,ε)) such that (x(,ε),y(,ε)) tends to (x,y ) when ε.

12 R.D. EUZÉBIO AND J. LLIBRE We observe that since the trace Σ of M is zero and the determinant Π is positive, then the real part of the eigenvalues of M are both zero. So we cannot decide about the stability of the periodic solution of Theorem 2. In what follows we see that the same occurs in Theorem 3. Proof of Theorem 3: We start considering n 3 = n 4 = 1, n 1,n 2 > 1 and m =. Therefore we are under the assumptions of Theorem 3. Doing the rescaling (5) the vector field of system (8) becomes (y, x+ε n 1 ry ε n 1 rx 2 y ε n 2 ax 3 εlx 5 +εdcos t) T. So we conclude that F 1 (t,x) = (, lx 5 +dcos t) and then we get f 1 (x,y ) = f 2 (x,y ) = (sin t)(dcos t l(x cos t+y sin t) 5 dt = 5 8 lπy ( x 2 +y 2 ) 2, (cos t)(dcos t l(x cos t+y sin t) 5 dt = 1 8 π (8d 5lx ( x 2 +y 2 ) 2 ). Now the real zero (x,y ) of f 1 and f 2 is ( (8d ) 1/5 (x,y) =,). 5l We note that as in Theorems 1 and 2 the root (x,y ) do not depend on ε after we perform rescaling (5). However using (5) the matrix M = (f 1,f 2 )/ (x,y ) at (x,y ) is ( ) M = 1 ((5/8)λδ 4 ) 1/5 π, ε 5((5/8)λδ 4 ) 1/5 π whose determinant is Π = 5π2 2ε 2 ( 25λ 2 δ 8 ) 1/5 >. Then the averaging method states the existence of a 2π periodic solution (x(t,ε),y(t,ε)) such that (x(,ε),y(,ε)) tends to (x,y ) when ε. The stability of this periodic solution is unknown because the trace Σ of M is zero and the determinant Π > is positive just as in Theorem 2. 2 Proof of Theorem 4: In what follows we assume that n 1 = m = 1 and n 4 = 2. Now the expression of the vector field of system (8) is (y, x+εry ε 3 rx 2 y ε n 2+2 ax 3 ε n 3+4 lx 5 +εdcos t) T,

PERIODIC SOLUTIONS OF THE φ 6 -VAN DER POL OSCILLATOR 13 and then F 1 (t,x) = (,ry +dcos t). In this case we obtain f 1 (x,y ) = f 2 (x,y ) = = πrx, = π(d+ry ). (sin t)(dcos t+r(y cos t x sin t)dt (cos t)(dcos t+r(y cos t x sin t)dt It is immediate that the only zero (x,y ) of f 1 and f 2 in this case is (x,y ) = (,d/r). So after going back through the rescaling (5) taking r = ε 1 ρ and d = ε 2 δ we have ( (x,y) =, δ ). ερ In addition the matrix M = (f 1,f 2 )/ (x,y ) at (x,y ) is diagonal and its elements are both given by ρπ/ε. Consequently M has determinant Π = (ρπ/ε) 2 > and trace Σ = 2ρπ/ε >. Then averaging theory described in section 3 assure the existence of a 2π periodic solution (x(t,ε),y(t,ε)) such that (x(,ε),y(,ε)) tends to (x,y ) when ε. Moreover since ρ is positive this periodic solution is unstable. Then Theorem 4 is proved. Proof of Theorem 5: Now we take m =, n 2 = n 3 = n 4 = 1 and n 1 > 1. The expression of the vector field of system (8) in this case is (y, x+ε n 1 ry ε n 1 rx 2 y εax 3 εlx 5 +εdcos t) T, and we get F 1 (t,x) = (, ax 3 lx 5 +dcos t). Hence function f = (f 1,f 2 ) writes f 1 (x,y ) = (sin t)(dcos t a(x cos t+y sin t) 3 l(x cos t+y sin t) 5 )dt = 1 8 πy (x 2 +y2 )(6a+5l(x2 +y2 )), f 2 (x,y ) = (cos t)(dcos t a(x cos t+y sin t) 3 l(x cos t+y sin t) 5 )dt = 1 8 π(8d x (x 2 +y2 )(6a+5l(x2 +y2 ))). In order to find the zeros of f(x,y ) we compute a Gröbner basis {b k (x, y ), k = 1,2,3} in the variables x and y for the set of polynomials {f 1 (x,y ), f 2 (x,y )}, where f 1,2 = (8/π)f 1,2. Then we look for the zeros of each b k, k = 1,2,3. It is a known fact that the zeros of a Gröbner basis of {f 1 (x,y )),f 2 (x,y )} are the zeros of f 1 and f 2, and consequently zeros

14 R.D. EUZÉBIO AND J. LLIBRE of f 1 and f 2 too. For more information about a Gröbner basis see [1] and [15]. The Gröbner is b 1 (x,y ) = dy, b 2 (x,y ) = (6ax 2 +5lx4 )y +(6a+1lx 2 )y3 +5ly5, b 3 (x,y ) = 8d+6ax 3 +5lx5 +(6ax +1lx 3 )y2 +5lx y 4. The only zero of b 1 (x,y ) is y = since d = ε 1 δ >. So replacing y into b 3 (x,y ) and simplifying the new expressions we reduce our problem to find a zero of the polynomial (11) p(x ) = 1 5l b 3(x,) = 8d rl + 6a 5l x3 +x5. As we know from algebra there is no general formula to provide the roots of a quintic polynomial as polynomial (11). Nevertheless some techniques may become this task easier. An elegant one can be found in [29]. In this paper it is provide a way to study the number of roots (complex and real, including multiplicity) of a polynomial of any order only performing some calculations on the coefficients of the considered polynomial. We will present a brief summary of the algorithmic of this method in the Appendix A at the end of this paper. By the Fundamental Theorem of Algebra polynomial (11) has 5 roots taking into account their multiplicities and once we want to apply the averaging theory we are only interested in the real simple roots of polynomial (11). Indeed, applying the method due to [29] and summarized in Appendix A we obtain D 4 = 384ad2 l 3, D 5 = 53747712 a 5 d 2 78125 l 7 +248 d4 l 4. Therefore as described in Appendix A, since αλ = ε 2 al < by hypothesis we have that D 4 is negative. So polynomial (11) has a unique simple real root if D 5 = D = according to statement (5) of Appendix A. On the other hand statement (3) says that polynomial (11) has three simple real roots if D 5 = D <. Moreover D and D 4 does not change when we replace (a,l,d) by (ε 1 α,ε 1 λ,ε 1 δ) using rescaling (5). Now using the property on zeros of the Gröbner basis cited previously the function f(x,y ) has a zero (x,) if D = and three zeros (xi,) if D <, with i = 1,2,3. Besides the matrix M = (f 1,f 2 )/ (x,y ) evaluated at = after rescaling (5) is y M = 1 ε π 8 (6αx2 +5λx4 ) π 8 (18αx2 +25λx4 ),

PERIODIC SOLUTIONS OF THE φ 6 -VAN DER POL OSCILLATOR 15 whose determinant is Π(x ) = 1 ( 27 ε 2 16 (απx2 )2 + 15 4 αλ(πx3 )2 + 125 ) 64 (λπx4 )2. Thus by computing a Gröbner basis for polynomials {Π(x ),f 1 (x,),f 2 (x, )} in the variables α, δ, λ and x we get the polynomials b 1 (α,δ,λ,x ) = δ and b 2 (α,δ,λ,x ) = 6α + 5λx 2. So b 1(α,δ,λ,x ) > because δ > and it means that we cannot have f 1, f 2 and Π equal to zero simultaneously. Consequently Π(x i ) for all i =,1,2,3. Then using the averaging theory we conclude the existence of a 2π periodic solution if D = and three 2π periodic solution if D <. Now we exhibit values of α, λ and δ for which system (4) has one or three 2π periodic solutions. Despite, if consider a = 25/18, l = 1 and d = 5/12 we obtain al = 25/18 < and D =. On the other hand, taking d = 18a/31 and l = 42a/155 we get al = 42a 2 /155 < and D = 142533982548/823543 <. So Theorem 5 is proved. We note that by assumption in Theorem 5 we have al <. So considering the notation of Appendix A the coefficients D 2 = 6a/5l and D 3 = 384ad 2 /l 3 corresponding to polynomial (11) are positive and D 4 is negative. Thus the only possible configurations of roots of the polynomial (11) are those one listed in statements (3) and (5) of the Appendix. This fact means that we cannot have 2 or 5 periodic solutions when we consider m =, n 2 = n 3 = n 4 = 1 and n 1 > 1. Proof of Theorem 6: We start fixing the values m =, n 1 = n 2 = n 4 = 1 and n 3 > 1. Now the vector field of system (8) writes (y, x+εry εrx 2 y εax 3 ε n 3 lx 5 +εdcos t) T, and then F 1 (t,x) = (,ry rx 2 y ax 3 +dcos t). With this expression of F 1 function (1) is f 1 (x,y ) = (sin t)(dcos t+r(y cos t x sin t) r(y cos t x sin t)(x cos t+y sin t) 2 a(x cos t+y sin t) 3 )dt = 1 4 π(rx ( 4+x 2 +y2 ) 3ay (x 2 +y2 )), f 2 (x,y ) = (cos t)(dcos t+r(y cos t x sin t) r(y cos t x sin t)(x cos t+y sin t) 2 a(x cos t+y sin t) 3 )dt = 1 4 π(4d ry ( 4+x 2 +y2 ) 3ax (x 2 +y2 )).

16 R.D. EUZÉBIO AND J. LLIBRE In order to find zeros (x,y ) of f(x,y ) as before we compute a Gröbner basis {b k (x,y ), k = 1,...,14} in the variables x and y now for the set of polynomials {f 1 (x,y ),f 2 (x,y )} with f 1,2 = (4/π)f 1,2. We will look for zeros of two elements of the Gröbner basis for the polynomials {f 1 (x,y ),f 2 (x,y )} in the variables x and y. These polynomials are b 1 (x,y ) = 144a 2 dr 3 4d 3 r 3 +(18a 2 d 2 r 2 +144a 2 r 4 4d 2 r 4 )y + 144a 2 dr 3 y 2 +(81a4 d 2 +18a 2 d 2 r 2 +d 2 r 4 )y 3, b 2 (x,y ) = 216a 3 dr 2 +( 324a 4 r 2 9a 2 d 2 r 2 +36a 2 r 4 d 2 r 4 )x + (27a 3 d 2 r+216a 3 r 3 +3ad 2 r 3 )y +(243a 5 d+54a 3 dr 2 + 3adr 4 )y 2. We must observe that b 1 (x,y ) depends only on y. Then for each zero y of b 1(y ) the second polynomial b 2 (x,y ) provides a zero x associated to y because the coefficient of x in b 2 (x,y ) is not zero by hypothesis. Now we will look for zeros of b 1 (y ). Indeed the discriminant of the cubic polynomial b 1 (y ) writes = 16d 2 r 6( 324a 4 +d 2 r 2 +9a 2( d 2 4r 2)) 2 ( 2187a 4 d 4 +27d 4 r 4 16d 2 r 6 +18a 2( 27d 4 r 2 72d 2 r 4 +32r 6)) = 16d 2 r 6 2 1 2. Soif 1 or 2 iszerothen iszeroandconsequentlyb 1 (y )hasasimplereal root y because this polynomial has no root with multiplicity 3. Actually b 1 (y ) = 6(81a 4 d 2 + 18a 2 d 2 r 2 + d 2 r 4 ) >. The same occurs if 2 > since this condition implies <. On the other hand, if 2 < then > and consequently polynomial b 1 (y ) has three simple real roots y i, i = 1,2,3. Additionally, replacing each value y i into b 2(x,y ) we obtain the respective values x i, for =,1,2,3. We note that coming back through the rescaling (5) the signs of, 1 and 2 does not change because each monomial composing has the same degree. Now we will verify the condition M = det(( f/ z)(x i,yi )) for i =,1,2,3. In fact M is π 4ε ( 6αx y +ρ ( 4+3x 2 +y2 ) 2ρx y +3α ( 3x 2 ) +y2 whose determinant Π now is 2ρx y 3α ( x 2 ) +3y2 6αx y +ρ ( 4+x 2 ) +3y2 Π(x,y ) = π2 16ε 2 ( 27α 2( x 2 +y 2 ) 2 +ρ 2 ( 4+x 2 +y 2 )( 4+3x 2 +3y 2 ) ). However for each i =,1,2,3 the determinant Π(x i,yi ) by hypothesis. Therefore using averaging theory described in section 3 system (4) has a 2π periodic solution if 1 2 = or 2 > and three 2π periodic solution if 2 <. ),

PERIODIC SOLUTIONS OF THE φ 6 -VAN DER POL OSCILLATOR 17 Now we present concrete values of α, ρ and δ for which we have one or three periodic solutions. Indeed taking a = r = 1 and d = 6 we have 2 = 3452544 >. Meantime if we consider d = 6a and 9a r = 9+4 6, we obtain 2 = 3412224a 8 /(9 4 6) 2 <. We also observe that the 15473 values (r, a, d) = ((555/4) 9622,185,22) and (r,a,d) = (585 3,195, 585 2) make 1 and 2 equal to zero, respectively. This ends the proof of Theorem 6. Proof of Theorem 7: Forthiscaseweassumem = andn i = 1, i = 1,2,3,4. Then the vector field of system (8) writes (y, x+εry εrx 2 y εax 3 εlx 5 +εdcos t) T, and we have F 1 (t,x) = (,ry rx 2 y ax 3 lx 5 +dcos t). So it follows that f 1 (x,y ) = (sin t)(dcos t+r(y cos t x sin t) r(y cos t x sin t)(x cos t+y sin t) 2 a(x cos t+y sin t) 3 l(x cos t+y sin t) 5 )dt = 1 8 π(2rx ( 4+x 2 +y2 ) y (x 2 +y2 )(6a+5l(x2 +y2 )), f 2 (x,y ) = (cos t)(dcos t+r(y cos t x sin t) r(y cos t x sin t)(x cos t+y sin t) 2 a(x cos t+y sin t) 3 l(x cos t+y sin t) 5 )dt = 1 8 π(8d+8ry (x 2 +y2 )(6ax +2ry +5lx (x 2 +y2 ))). We will find zeros of f(x,y ) through the zeros of a Gröbner basis {b k (x,y )} of it. As before instead of f(x,y ) the Gröbner basis will be related with the functions f 1,2 = (8/π)f 1,2. One of the elements of the Gröbner basis in the variables x and y is b 1 (x,y ) = 144a 2 dr 5 4d 3 r 5 +96adlr 5 +16dl 2 r 5 +(18a 2 d 2 r 4 + 96ad 2 lr 4 +2d 2 l 2 r 4 +144a 2 r 6 4d 2 r 6 +96alr 6 + 16l 2 r 6 )y +(12ad 3 lr 3 +5d 3 l 2 r 3 +144a 2 dr 5 +12 adlr 5 +24dl 2 r 5 )y 2 +81a4 d 2 r 2 +54a 3 d 2 lr 2 +9a 2 d 2 l 2 r 2 +18a 2 d 2 r 4 +3ad 2 lr 4 +9d 2 l 2 r 4 +d 2 r 6 )y 3 +( 45 a 2 d 3 l 2 r 15ad 3 l 3 r +5d 3 l 2 r 3 )y 4 +625d4 l 4 y 5.

18 R.D. EUZÉBIO AND J. LLIBRE We note that again b 1 (x,y ) depends only on y. Moreover, the Gröbner basis has other element b 2 (x,y ) which is linear on x and the coefficients depend on y, where the coefficient C of x is not zero by hypotheses. Its expression is too large and we will omit it here. It can be easily obtained by analgebraicmanipulator. Soforeachzeroy ofb 1(y )thesecondpolynomial b 2 (x,y ) provides a zero x associated to y. Now we will look for zeros of b 1 (y ). Indeed using again Appendix A the conditions on a, r, d and l that provide three zeros, (x i,yi ), where i = 1,2,3 for b 1 (y ) are D 5 = D5 1D2 5 D <, where D is given in Theorem 7 and D1 5 and D5 2 are the values D 1 5 = 1175314176a8 r 2 d 2 ( 9a 2 +r 2 ) 28, D 2 5 = (5949a1 d 4 +16r 14 +36a 2 r 1 (11d 2 4r 2 ) 81a 4 r 6 (d 4 12d 2 r 2 + 16r 4 ) 6561a 8 (3d 4 r 2 4d 2 r 4 )+729a 6 (3d 4 r 4 28d 2 r 6 +16r 8 )) 2. By hypothesis ar = ε 2 αρ and D5 2, then D1 5 and D2 5 are positive. Therefore, in order to have D 5 negative we must have D positive. On the other hand, if D 5 is positive then according to statement (2) of Appendix A we have a unique zero (x,y ) for b 1(y ) because the value D 2 = 1296a2 r 4 ( 3a 2 +r 2 ) 2 d 2 ( 9a 2 +r 2 ) 4, is negative if (r 2 3a 2 )(r 2 9a 2 ). In addition by considering a = ε 1 α, r = ε 1 ρ and d = ε 1 δ we get M = π ( 2ρ( 4+3x 2 +y 2) Γ 1 Γ 2 Γ 4 5λΓ 3 (x 2 +5y2 ) ) 8ε Γ 2 +Γ 4 +5λΓ 3 (5x 2 +y2 ) 2ρ( 4+x2 +3y2 ) Γ 1 where Γ 1 = 2x y (3α + 5l(x 2 + y2 )), Γ 2 = 4ρx y, Γ 3 = x 2 + y2 and Γ 4 = 6α(x 2 +3y2 ). Moreover the determinant Π of M writes Π(x,y ) = 1 ε 2 ( 4ρ 2 ( 4+x 2 +y2 ) ( 4+3x 2 +3y 2 )+(x2 +y2 )2 (6α+ ρλ(x 2 +y2 )(18α+25λx2 +3y2 ))). However, by hypothesis Π(x i,yi ) for i =,1,2,3 and hence we conclude the first part of Theorem 7. Now we provide values for α, ρ, λ and δ for which we have one or three 2π periodic solutions. First we observe that taking (a,r,d) = (1,2,8/3) we obtain D = 2773582112 <. 6561 However if we consider (a,r,d,l) = (1/ 3,1,1, 1/(5 3)), then D = 2752. So Theorem 7 is proved.

PERIODIC SOLUTIONS OF THE φ 6 -VAN DER POL OSCILLATOR 19 Now we will prove that under the conditions m = and n i = 1, i = 1,2,3,4 and the assumptions of Theorem 7 we cannot have 2 or 5 periodic solutions by using the averaging theory and the rescaling (5). Indeed we observe that the only possibility to have 5 periodic solutions that the statement (1) of Appendix A holds. However this condition does not hold since D 2 is negative for any a, r and d satisfying the hypotheses of Theorem 7. In addition we note that the only condition that provides two periodic solutions according with Appendix A is statement (7), and it needs D 5 = D 4 = E 2 =. Nevertheless we start considering this condition and we will find a Gröbner basis for the set of polynomials {D 4,D 5,E 2 } in the variables a, d and r. First we will consider the factor D 2 5 of D 5 because D 1 5 by hypothesis. In this case we obtain 55 polynomials in the Gröbner basis, where the first one is g(d,r) = r 62 g(d,r) and g(d,r) = 91125d 12 4388d 1 r 2 +91936d 8 r 4 138464d 6 r 6 + 66992d 4 r 8 22528d 2 r 1 +32768r 12. Since r we must find a zero of g(d,r). So solving g(d,r) = in the variable r we obtain 6 pairs of values r i ± = ±d q i, where i = 1,...,6 and for each i the value q i is the ith root of the polynomial q(x) = 91125 4388x +91936x 2 138464x 3 +66992x 4 22528x 5 +32768x 6. This polynomial has only complex roots. Indeed, if we study the function q (x) which is given by q (x) = 4388+183872x 3115392x 2 +2643968x 3 11264x 4 +19668x 5, we can apply the method described in Appendix A to show that the respective values D 2 and D 5 related to q (x) are negative and positive, respectively. Therefore q (x) has only one real root q whose approximate value is q = 1673/1. Additionally the approximate value of q(x) evaluated in q is q(q ) = 16558687/1 >. So if we observe that the coefficient of x 6 is positive it follows that the minimum value of q(x) is positive. Then each r± i, i = 1,...,6 is a complex valueand thecorrespondent Gröbnerbasis has no zeros. Then we cannot have a common real zero of D 4, D5 2 and E 2 and so statement (7) in Appendix A does not hold. Consequently it is not possible to obtain 2 periodic solutions. The proof considering the factor D of D 5 instead of D5 2 leads to the same polynomial g(d,r) and then we have the same conclusion.

2 R.D. EUZÉBIO AND J. LLIBRE Proof of Theorem 8: In what follows we take m =, n 1 = n 3 = n 4 = 1 and n 2 > 1. So the vector field of system (8) becomes (y, x+εry εrx 2 y ε n 2 ax 3 εlx 5 +εdcos t) T. and we obtain F 1 (t,x) = (,ry rx 2 y lx 5 +dcos t). So f(x,y ) writes f 1 (x,y ) = (sin t)(dcos t+r(y cos t x sin t) r(y cos t x sin t)(x cos t+y sin t) 2 l(x cos t+y sin t) 5 )dt = 1 8 π(2rx ( 4+x 2 +y2 ) 5ly (x 2 +y2 )2 ), f 2 (x,y ) = (cos t)(dcos t+r(y cos t x sin t) r(y cos t x sin t)(x cos t+y sin t) 2 l(x cos t+y sin t) 5 )dt = 1 8 π(8d 2ry ( 4+x 2 +y2 ) 5lx (x 2 +y2 )2 ). Again we will find zeros of f(x,y ) through the roots of a Gröbner basis {b k (x,y )} of f(x,y ). We will find a Gröbner basis of the functions f 1,2 = (8/π)f 1,2. An element of the Gröbner basis in the variables x and y is b 1 (x,y ) = 4d 3 r 5 +16dL 2 r 5 +(2d 2 l 2 r 4 4d 2 r 6 +16L 2 R 6 )y +(5d 3 l 2 r 3 +24dl 2 r 5 )y 2 +(9d2 l 2 r 4 +d 2 r 6 )y 3 + 5d 3 l 2 r 3 y 4 +625d4 l 4 y 5. As before b 1 (x,y ) depends only on y and another element b 2 (x,y ) of the Gröbner basis is linear on x with coefficients depending on y, where the coefficient C of x in b 2 (x,y ) is not zero by hypotheses. We will not present the expression of b 2 (x,y ) in order to avoid large expressions. Then for each zero y of b 1(y ) we have a second zero x through b 2(x,y ) related to y. In order to apply the method described in [29] and Appendix A for b 1 (y ) we will perform the translation y = ϕ (2r 3 /d(125l) 2 ). With this translation we obtain a new polynomial b 1 into the form b 1(ϕ) = ϕ 5 + 45l2 r 4 3r 6 3125d 2 l 4 ϕ 3 + 2k 1 k 1 d 3 l 6ϕ2 + 4r4 k 2 k 2 d 4 l 8ϕ 4r5 k 3 k 3 d 5 l 1,

PERIODIC SOLUTIONS OF THE φ 6 -VAN DER POL OSCILLATOR 21 where k 1 = 15625d 2 l 4 r 3 +75l 4 r 5 135l 2 r 7 +r 9, k 1 = 39625, k 2 = 39625d 2 l 6 39625l 4 (d 2 8l 2 )r 2 15l 4 r 4 +135l 2 r 6 + 3r 8, k 2 = 48828125, k 3 = 48828125d 4 l 6 78125d 2 (25l 8 5l 6 r 2 +3l 4 r 4 )+2r 4 (15625l 6 375l 4 r 2 +225l 2 r 4 +9r 6 ), k 3 = 3517578125. We observe that the translation performed do not change the number or kind of the zeros of the original polynomial b 1. Now we will apply the method of Appendix A for b 1. Indeed we have D 2 = N 2, 48r 6 D 3 = 9765625d 6 l 1N 3, 16r 12 D 4 = 3814697265625d 12 l 18N 4, 256r 2 D 5 = 149116119384765625d 2 l 26N 5M 5. Therefore if D 5 is negative then b 1 and consequently b 1 has exactly three zeros and hence function f = (f 1,f 2 ) has exactly three zeros (x i,yi ), i = 1,2,3. On the other hand if D 5 is positive and one of the values D 2, D 3 or D 4 is non positive, then f has exactly one zero (x,y ). Moreover using the rescaling (5) the matrix M is the same than the one of the proof of Theorem 7 taking a =. In addition the determinant Π of M writes Π(x,y ) = 1 ε 2 ( 125λ 2 ( x 2 +y4 ) +4ρ 2 ( 4+x 2 +y2 )( 4+3x2 +3y2 )). By hypothesis Π(x i,yi ) for i =,1,2,3, then we have the first part of Theorem 8 proved. Now we exhibit values of ρ, λ and δ for which we have either one, or three periodic solutions. Firstly we consider (r,l,d) = ( 31,1, 1 125 1547875 226851 35 With these values we obtain N 3 = 1 625 (26837718179 141999532263 35), N 5 = 2( 64949147851728715458244 +1963983255479427558525 35), 48828125 ).

22 R.D. EUZÉBIO AND J. LLIBRE which are positive. On the other hand, considering (r,l,d) = ( 31/4,1, 1) we get N 5 = 258261575 /148576 <. This ends the proof of Theorem 8. Now we will show that if m =, n 1 = n 3 = n 4 = 1 and n 2 > 1, then we cannot have 2 periodic solutions by using averaging theory described in section 3 and rescaling (5). Indeed in order that the function b 1 (ϕ), given in the proof of Theorem 8, has two periodic solutions we need statement (7) of Appendix A. We start considering the conditions N 5 = N 4 = E 2 = and we will see that these equalities implies that N 3 cannot be negative. We claim that N 5 and M 5 are factors of D 5 as we saw in the last proof. Also we will not present the expression of E 2 in order to avoid its large expression. It can be obtained by an algebraic manipulator as Mathematica. Additionally we note that N 4 = N 5 = implies D 4 = D 5 =, which are necessary conditions for having exactly 2 real roots. We start replacing the condition N 4 = N 5 = E 2 = for other ones easier when r >. It is easy to see that if r = we have N 3 = and consequently we cannot get 2 roots for b 1. Now we obtain a Gröbner basis for the set of polynomials {N 4,N 5,E 2 } in the variables d, r and l. This basis has twelve elements, where two of them are g 1 (r,l) = 125 1 8 l 8 +2 1 8 l 6 r 2 1l 4 r 4 14l 2 r 6 3r 8, g 2 (d,r,l) = 25d 4 l 2 11d 2 l 4 4l 6 14d 2 l 2 r 2 +8l 4 r 2 3d 2 r 4 +12l 2 r 4. In addition the resultant of these two polynomials with respect to the variable r has the factor g 3 (d,l) = d 8 124d 6 l 2 26944d 4 l 4 +6144 1 6 d 2 l 6 +28672 1 8 l 8. Consequently we will replace the problem of finding a zero of N 4, N 5 and E 2 by the problem of finding a zero of g 1, g 2 and g 3. Moreover we suppose that we have a common zero between each g i and N 3, i = 1,2,3. With this supposition and computing the resultant between g 1 and N 3 we obtain a new polynomial g 4 in the variables r and d given by g 4 (r,d) = 1423828125d 16 +25974d 14 r 2 47933388d 12 r 4 4788124512d 1 r 6 +17118322432d 8 r 8 125591686188d 6 r 1 133461671398648d 4 r 12 + 51918655746653936d 2 r 14 445633368973152r 16. However the resultant between g 3 and g 4 writes Kr 128, wherek is a positive constant. This polynomial cannot be zero since r. It means that when we consider each g i zero, N 3 must be positive or negative, exclusively. But taking the values (r,l,d) = (r,l,1) we have g i, for each i = 1,2,3 and N 3 is approximately 8738483/25 >, where r and l are roots of the

polynomials and PERIODIC SOLUTIONS OF THE φ 6 -VAN DER POL OSCILLATOR 23 r(x) = 125 444x 2 +18464x 4 237568x 6 +8616x 8, l(x) = 1 124x 2 26944x 4 +6144x 6 +28672x 8. respectively. So N 3 is positive when g i, for each i = 1,2,3 and then D 3 is negative. Hence using the factor N 5 of D 5 we cannot obtain 2 periodic solutions. The proof that cannot exists 2 periodic solutions for the case that we consider the factor M 5 of D 5 instead of N 5 is similar and we will omit it here. We remark that we cannot prove analytically the non existence of 5 periodic solutions for system (8) considering m =, n 1 = n 3 = n 4 = 1 and n 2 > 1. Actually, using the algebraic manipulator Mathematica we obtained evidences that this number of periodic solutions cannot happen, but an analytical treatment is not trivial. 2.2. Discussion of the results. Now we discuss some aspects of the averaging method presented in section 3 in order to find periodic solutions in system (4). In fact by using the averaging theory and the rescaling (5) we shall see that we cannot prove that system (4) has periodic solutions except in the cases presented in Theorems from 1 to 8. We start studying the possible values to the powers of ε in system (8). This is important because these powers play an important role in averaging theory described in section 3 because they determine the terms of order 1 that we are interested in. Indeed we note that each one of the 5 different powers of ε in system (8) must be non negative. Therefore considering these powers as zero or positive we have 32 possible combinations of these 5 different powers of ε. Actually some of them are not algebraically possible. Table 1 exhibit only the possible case. Observing Table 1 we see that in fact we have only 18 cases. Moreover in order to apply averaging theory, for which one of the 18 cases we must integrate the equations of the non perturberd part of vector field of system (8). This is not a simple task because the major part of these equations are nonlinear and non autonomous. Indeed in each case from 1 to 8 and from 11 to 16 we could not integrate the equations even using the algebraic manipulators Mathematica or Maple. The non integrable cases from Table 1 are listed above. Caso 1: F (t,x) = (y, x+ry ax 3 rx 2 y lx 5 +dcos t), Caso 2: F (t,x) = (y, x+ry ax 3 rx 2 y lx 5 ), Caso 3: F (t,x) = (y, x+ry ax 3 rx 2 y +dcos t), Caso 4: F (t,x) = (y, x+ry ax 3 rx 2 y), Caso 5: F (t,x) = (y, x+ry rx 2 y lx 5 +dcos t), Caso 6: F (t,x) = (y, x+ry rx 2 y lx 5 ), Caso 7: F (t,x) = (y, x+ry rx 2 y +dcos t),

24 R.D. EUZÉBIO AND J. LLIBRE n 1 = Possible combination of powers of ε for the system (9) 2m+n 1 = 2m+n 2 = 2m+n 2 > 4m+n 3 = m+n 4 = C1 m+n 4 > C2 4m+n 3 > m+n 4 = C3 m+n 4 > C4 4m+n 3 = m+n 4 = C5 m+n 4 > C6 4m+n 3 > m+n 4 = C7 m+n 4 > C8 2m+n 2 = 2m+n 1 > 2m+n 2 > 4m+n 3 = 4m+n 3 > m+n 4 = C9 m+n 4 > C1 2m+n 1 = n 1 > 4m+n 3 = m+n 4 = C11 m+n 4 > C12 2m+n 2 = 4m+n 3 > m+n 4 = C13 m+n 4 > C14 2m+n 1 > 4m+n 3 = m+n 4 = C15 m+n 4 > C16 2m+n 2 > 4m+n 3 > m+n 4 = C17 m+n 4 > C18 Table 1. Here, when appear, this means that the corresponding cases cannot occur. The notation Ci, i = 1,...,18 enumerate each possible case.

PERIODIC SOLUTIONS OF THE φ 6 -VAN DER POL OSCILLATOR 25 Caso 8: F (t,x) = (y, x+ry rx 2 y), Caso 11: F (t,x) = (y, x ax 3 lx 5 +dcos t), Caso 12: F (t,x) = (y, x ax 3 lx 5 ), Caso 13: F (t,x) = (y, x ax 3 +dcos t), Caso 14: F (t,x) = (y, x ax 3 ), Caso 15: F (t,x) = (y, x lx 5 +dcos t), Caso 16: F (t,x) = (y, x lx 5 ). In particular we observe that cases 12, 14 and 16 turn the non perturbed part of system (8) into a Hamiltonians system. In [4] the authors provide a method to apply averaging theory in planar systems when the system is Hamiltonian but the expressions become also to much complicate and again we could not integrate the equations in cases 12, 14 and 16. In cases 9, 1 and 17 we can integrate the expressions but in these cases the hypotheses of averaging method do not apply. In these cases we have the following expressions for F (t,x). Caso 9: F (t,x) = (y, x+ry +dcos t), Caso 1: F (t,x) = (y, x+ry), Caso 17: F (t,x) = (y, x+dcos t). Solution x(t) in case 9 is where 1 1 C e2 t(r 4+r 2 ) +C + e2 t(r+ 4+r 2 ) (d/r)sin t, C = 2d+r ((±r+ ) 4+r )x 2 2y 2r. 4+r 2 In order that x(t) be periodic we must have C and C + equal to zero. But these conditions are verified only if we take x = and y = d/r. So x and y are fixed and then the unperturbed system ẋ = F (t,x) has no submanifold of periodic solutions. Consequently we can not apply averaging theory in this case. On the other hand, in case 1 solution x(t) writes where 1 1 C e2 t(r 4+r 2 ) +C + e2 t(r+ 4+r 2 ) (d/r)sin t, C = (±r+ 4+r 2 ) x 2y 2r 4+r 2. Again we need to chose C and C + equal to zero but this happen only if (x,y ) = (,). So we cannot apply averaging theory because (,) is the equilibrium point of the non perturbed part of system (8) in case 1.

26 R.D. EUZÉBIO AND J. LLIBRE Finally in case 17 solutions x(t) and y(t) have the form x(t) = x cos t+y sin t+ d tcos t, 2 y(t) = y cos t x sin t+ d (tcos t+sin t). 2 It is immediate to note that these solutions are non periodic because d is positive. Case 18 provide positive results and Theorems from 1 to 8 are based in this case. In fact we observe that in case 18 the vector field of system (8) is (y, x+ε n 1 ry ε 2m+n 1 rx 2 y ε 2m+n 2 ax 3 ε 4m+n 3 lx 5 +ε m+n 4 dcos t) T. As we commented before the periodic solutions of our problem correspond to the simple zeros of function (1). So in order to calculate the zeros of the function f(z) given in (1) we must determine F 1 (t,x), where F 1 (t,x) is determinated by the terms of order 1 on ε of the vector field of system (8). However we see that the expression of F 1 (t,x) depends on the values m and n i, i = 1,2,3,4. In fact, if m is positive and observing the conditions (7) the only terms of the vector field of order 1 on ε can be generate by powers n 1 and m+n 4, where these two values can be 1 or greater than 1. This implies four possibilities. On the other hand, if m is zero the powers of ε in the vector field of system (8) depend on n 1, n 2, n 3 and n 4. Again, each one of these powers can be 1 one greater than 1. So, if m is zero we have 16 possibilities for F 1. Then case 18 has 2 subcases corresponding to the different possibilities of F 1 (t,x). These subcases are presented in Table 2. Theorems from 1 to 8 correspond to subcases 8, 1, 11, 17, 5, 3, 1 and 4, respectively. These are the only subcases of case 18 that the averaging method provide positive results. Indeed, in cases 15, 16, 19 and 2 function (1) do not have zeros. On the other hand, in cases 2, 6, 7, 13, 14, 18 the only zero of function (1) is (x,y ) = (,), that correspond to the equilibrium point of the system, and consequently in these subcases the system does not have periodic solutions. In cases 9 and 12 function (1) has real zeros different from (,) but they are non isolated and then the Jacobian of function (1) at the zero is zero, and consequently the averaging theory cannot be applied. 3. The Averaging Theory for Periodic Solutions Now we present the basic results on the averaging theory of first order that we need to proving our results. Consider the problem of bifurcation of T periodic solutions from differential systems of the form (12) ẋ = F (t,x)+εf 1 (t,x)+ε 2 R(t,x,ε),