Intelligent Automation & Soft Computing DOI:10.32604/iasc.2022.019767
Article

Periodic Solutions for Two Dimensional Quartic Non-Autonomous Differential Equation

Saima Akram1,*, Allah Nawaz1, Muhammad Bilal Riaz2 and Mariam Rehman3

1Centre for Advanced Studies in Pure and Applied Mathematics, Bahauddin Zakariya University, Multan, 60000, Pakistan
2Department of Mathematics, University of Management and Technology, Lahore, 54000, Pakistan
3Department of Information Technology, Government College University, 38000, Faisalabad
*Corresponding Author: Saima Akram. Email: saimaakram@bzu.edu.pk
Received: 25 April 2021; Accepted: 15 June 2021

Abstract: In this article, the maximum possible numbers of periodic solutions for the quartic differential equation are calculated. In this regard, for the first time in the literature, we developed new formulae to determine the maximum number of periodic solutions greater than eight for the quartic equation. To obtain the maximum number of periodic solutions, we used a systematic procedure of bifurcation analysis. We used computer algebra Maple 18 to solve lengthy calculations that appeared in the formulae of focal values as integrations. The newly developed formulae were applied to a variety of polynomials with algebraic and homogeneous trigonometric coefficients of various degrees. We were able to validate our newly developed formulae by obtaining maximum multiplicity nine in the class C4,1 using algebraic coefficients. Whereas the maximum number of periodic solutions for the classes C4,4; C5,1; C5,5; C6,1; C6:6; C7,1 is eight. Additionally, the stability of limit cycles belonging to the aforementioned classes with algebraic coefficients is briefly discussed. Hence, we conclude from the above-stated facts that our new results are a credible, authentic and pleasant addition to the literature.

Keywords: Limit cycle; nonlinear equation; quartic differential equation; algebraic and trigonometric coefficients; focal values

1  Introduction

Most of the real-world problems in nature are multidimensional and when modeled, they arise as higher-order ordinary differential equations. We are interested in those models which are periodic and depend on time and are usually known as non-autonomous. This article contains several recent developments and advances in calculations of periodic solutions and their applications in various areas of the mathematical, physical and engineering sciences. We have investigated upper bounds for the non-autonomous ordinary differential equation (ODE) of the cubic degree [1–3]. The primary question striking in our minds is to investigate the maximum number of periodic solutions when the degree of non-autonomous (ODE) is increased from three to four; so, we started working for the quartic system. The analysis of periodic solutions is vital because they frequently arise as real-world problems from financial matters such as modelling economic processes to complex space robotics, from galaxies to weather forecasting models. Many physics problems are related to nonlinear analysis like relativity, elasticity, chaotic dynamics and Navier stokes equation in fluid dynamics. Bendixson’s negative criterion, based on the connection between compound matrices and (ODE’s), is applied to prove the nonexistence of periodic orbits and then nonexistence of Hopf bifurcation. Almost every bit of life needs this analysis, even inside the body as cardiac rhythms as blood flow, please see the example, [4–6]. The aesthetic appeal of periodic solutions may explain why so many people have become intrigued by these ideas. Alwash et al. [7,8] developed the theoretical foundations of periodic solutions of differential equations.

Alwash et al. [8] examine equation of the form:

y˙=g(z)y3+d(z)y2+u(z)y,(1)

to find the upper bound for the number of limit cycles for such system (1) in accordance with the second part of Hilbert’s sixteenth problem. The coefficients g; d; and u are real-valued continuous functions, but the independent variable is complex. We are principally concerned about the multiplicity of y=0 as a periodic solution, i.e., solutions satisfying

y(0)=y(σ),(2)

for σ ∈ R and look for information about the number of periodic solutions. We refer the reader to [8] for more extra pieces of information, however right here we do not forget that the multiplicity of a solution ϕ(t) of Eq. (1) satisfying Eq. (2) is the multiplicity of ϕ(0) as a root of the holomorphic function q:c→y(t,0,c)−c;y(t,t1,c) is the solution satisfies y(t,t1,c)=c . Here we take a multiplicity of the periodic solution as ‘k’ of the function ϕ(t) . Rouche’s theorem is used for the function q, it is assumed that we have at most k periodic solutions in the regions nearby ϕ (counting multiplicity): see, for example, ([2], Theorem 2.4) when sufficiently small perturbation of the equation is in an account for the following quartic equation:

y˙=ϵy4+g(z)y3+d(z)y2+u(z)y,(3)

For sufficiently small perturbation ϵ, the quartic equation preserves the same number of periodic solutions as Eq. (1). We used scaling y→ϵ1/3y to get the leading coefficient of Eq. (3) as one and this Eq. (3) becomes as:

y˙=y4+ϵ−2/3(z)y3+ϵ−1/3(z)y2+u(z)y,(4)

To make the above equation simple, which is helpful for lengthy calculations that arise in numerical integration, we rewrite Eq. (4) as:

y˙=y4+g(z)y3+d(z)y2+u(z),(5)

We are mainly focused on finding the multiplicity of periodic solutions of y = 0 greater than one for Eq. (1) as was done by researchers in [7], without loss of generality; we could take u(z) ≃ 0; as was taken in [1]. In this paper, we are considering the corresponding equation:

y˙=y4+g(z)y3+d(z)y2,(6)

There is another reason for our interest in Eq. (5). We say this system to be non-autonomous, nonlinear because all the ‘ y ’ on the right-hand side appear to a power greater than one. Because of the non-autonomous equation, equilibrium states are not usually associated with Eq. (6). The general form of the above Eq. (6) is as follows:

y˙=yn+p1(z)yn−1+p2(z)yn−2+,…,+p0(z);here p0(z)=1.(7)

Moreover, the coefficients were considered to be periodic functions, please see the example, [8–11]. However, the questions related to Hilbert’s sixteenth problem are reduced to polynomial equations in which p0 has zeros. They showed that for n = 3, Eq. (7) has precisely three periodic solutions. Here, the results of [12] no longer hold; they have presented some examples in [13] which showed that there is no idea about the upper bound for the number of periodic solutions when n≥4 and p0(z) = 1; until some coefficients are restricted. Earlier, Shahshahani [14]; determined the multiplicity of Eq. (6) when g(z),d(z) and u(z) are polynomial functions of z. We also present in Section 3 that the conjecture made in [14] is false.

In Section 2, some results from the literature are employed and formulae for calculating the maximum number of periodic solutions, greater than eight are formulated and presented. In Section 3, we calculated the periodic multiplicity for various classes of coefficients. In the last section, we present conclusions and discussions.

2  Calculation of Multiplicity

In this section, we recall some important concepts necessary to understand the presented method; for instance, see [1]. For the sake of multiplicity of the zero solution of Eq. (5), we observe that in a neighborhood of y=0 , we can write

y(z,0,c)=∑i=0∞ui(z)ci(8)

for 0≤z≤σ ; the functions ui(z) are continuous and u1(0)=1 and ui(0)=0 for i>1 . The multiplicity (μ) is “ μ=k ” if

u1(σ)=1,u2(σ)=u3(σ)=…=ui−1(σ)=0.}(9)

However, ui(σ)≠0 . To calculate the functions ui(z) , we placed the expression (8) into (5) and achieved a set of linear differential equations which can be solved recursively. We will study from Eq. (9) that u1(z)=u1(z)υ(z) , where u1(z) is set as:

u_1(z)=e∫oσv(s)ds.

Hence, y=0 is a multiple solution iff ∫oσv(s)ds . Since we were right here interested in equations, wherein the multiplicity of the origin is greater than one, we recall that ∫oσv=0 . At that point, we apply the transformation y→ypn−1 , wherein p(z)=e∫ozv and obtain

y=(p(z))3y4+(p(z))2y3+p(z)υ(z)y2.(10)

Since p is periodic, the function q is unchanged through the transformation and therefore the periodic solution of Eq. (5) and Eq. (10) have identical preliminary points and multiplicities. Additionally, the transformation of the independent prompts to the Eq. (6). For Eq. (6), the functions ui(z) , for i>1 are calculated by the following relation:

u˙i=(∑j+k+l+m=ij,k,l,m≥1)(ujukulum)+γ(z)(∑j+k+l=ij,k,l≥1)(ujukul)+δ(z)(∑j+k=ij,k≥1)(ujuk).(11)

With u1(z)=1 . For i≤8 functions ui(z) and ϰi are provided in [1], for i=9 , we have calculated u9(z) and ϰ9 ; and present them below in Theorems 2.1 and 2.2. Some tremendously complicated computations are involved; hence, we omit some details and we write results.

In this article, we present a novel approach for finding periodic solutions. By using Theorem 2.2, we can find the highest periodic solutions as 9. We base this method on the construction of Akram et al. [1]. The approach can apply to other planar systems as well.

When we started working for quartic non-autonomous (ODE), for this investigation, formulae for finding upper bounds greater than 8 are not present in the literature up to now. We came up with many challenges while working with this problem of quartic type. The first task is to compute these new formulas, which are previously unavailable in the literature. By putting in many efforts, we succeeded in constructing the new formula u9 and ϰ9 , which are given in theorems 2.1 and 2.2. In the below theorem, we use the notation “ (.)¯ ” for indefinite integral such as δ¯=∫δ(s)ds .

Theorem 2.1 For Eq. (6); the functions u2, u3, … , u9 in the expansion (11) are as follows:

u2=δ¯,

u3=δ¯2+γ¯,

u4=δ¯3+2δ¯γ¯+δ¯γ¯+z,

u5=δ¯4+3δ¯2γ¯+δ¯2γ¯+3δ¯δ¯γ¯+2δ?+32γ¯2+2zδ¯,

u6=δ¯5+4δ¯3γ¯+δ¯3γ¯+3δ¯2δ¯γ¯+3δ¯2¯+4δ¯δ?+2δ¯δ¯2γ¯−12δγ¯2¯+92δ¯γ¯2+3δ¯γ¯γ¯+γ?+3z(δ¯2+γ¯),

u7=δ¯6+5δ¯4γ¯+δ¯4γ¯+4δ¯3δ¯γ¯+2δ¯δ¯3γ¯+4δ¯3¯+3δ¯2δ¯2γ¯+6δ¯δ¯2¯+6δ¯2δ?+272δ¯2γ¯2+3δ¯2γ¯γ¯−δγ¯2¯δ¯−2δδ¯γ¯2¯+6(δ¯γ¯)2+2δ¯γ?+4δ¯γ¯¯−4(δγ¯¯)2+6δ?⁡γ¯+52γ¯3+4z(δ¯3+3δ¯γ¯+2δγ¯¯)+2z,

u8=δ¯7+6δ¯5γ¯+δ¯5γ¯+2δ¯δ¯4γ¯+5δ¯4δ¯γ¯+3δ¯2δ¯3γ¯+4δ¯3δ¯2γ¯+8δ¯δ¯3¯+9δ¯2δ¯2¯−8δ¯3δ?+5δ¯4¯+752δ¯3γ¯2+3δ¯3γ¯γ¯+12δ¯2γ¯δ¯γ¯−δ¯2δγ¯¯γ¯−6δ¯2δγ¯2¯−25δ¯2δγ¯¯γ¯−32δ¯2δγ¯2¯−4δ¯2γ¯δγ¯¯+5δ¯(δγ¯¯)2−4δ¯δδ¯γ¯2¯+9δ¯2γ¯¯−δ¯2γ?−2δ?⁡δ¯γ¯+9δ¯2¯γ¯+3δ¯2γ¯+16δ¯δ?⁡γ¯+8δ¯δ¯γ¯¯+10δ?⁡δ¯γ¯−2δ¯?−12δγ¯3¯−32δγ¯2¯γ¯+10δ¯γ¯3+152δ¯γ¯γ¯2+32γ¯2¯+3γ¯γ?+z(5δ¯$4$+15δ¯2γ¯+5δ¯2γ¯+10δ¯δ¯γ¯+10γ?+152γ¯2)+5z2δ¯,

u9=δ¯8−δ¯3γ¯+4δ¯3γ+58δ¯γ¯δ¯2¯+12δ¯2+20(δ)2+10δ¯2γ¯δ¯¯−2δγ¯2+2δ¯γ¯2¯+24δ¯γ¯γ¯−2γ+38zδ¯2¯−38δ¯2+24zγ+35δ¯6γ¯−1612δ¯4γ¯2+3δ¯4γ¯γ¯+30δ¯γ¯δ¯3γ¯+90δδ¯2δ¯γ¯¯+4δ¯3γ¯γ¯δ¯−4δ¯3γ¯γ¯δ¯+12δ¯3¯γ¯+11γ¯δ¯2δ¯2γ¯+363γ¯δ¯δ¯2¯+34γ¯δ¯2δ+392δ¯2γ¯3+152γ¯2δ¯2γ¯+21γ¯(δ¯γ¯)2−2γ¯δ¯δγ¯2¯+4γ¯γ⁡δ¯+358γ¯4+12γ¯δ¯γ¯¯+15γ¯2δ+583γ¯δ¯3z−58γ¯δδ¯2z¯+24γ¯δγ¯¯−24γ¯2δ¯+9γ¯z2−36γγ¯δ¯δ¯γ¯¯−66γ¯δδ¯2z¯+42γ¯δ¯γ¯z−42γ¯δ¯γ¯¯−42γ¯2δ¯z+42γ¯2δz¯+9γ¯δ¯3δ¯2γ¯−9γ¯δ¯5γ¯+27γ¯δδ¯2δ¯2γ¯¯−12γ¯(δγ¯¯)2+24γ¯2δδγ¯¯¯−6γδδ¯γ¯2¯¯+2δ¯δ¯5γ¯+7γ¯δ¯6+3δ¯2δ¯4γ¯+325δ¯5δ¯γ¯+4δ¯3δ¯3γ¯+5δ¯4δ¯2γ¯+12δ¯3δ¯2¯+12δ¯2δ¯3¯+10δ¯4δ+10δ¯δ¯4¯+1013δ¯4γ¯+8γ¯δ¯δ¯3γ¯+19γ¯δ¯2δ¯2γ¯−2δ¯δ¯2δγ¯¯γ¯+823γδ¯3δγ¯¯¯−12δ¯γ¯2δδ¯2¯−523γ¯δ¯3δγ¯¯−2δ¯3δγ¯2¯−8δδ¯2γ¯δγ¯¯¯+δ¯2(δγ¯¯)2−6δ¯2δδ¯γ¯2¯+18δ¯δ¯2γ¯¯−2δ¯δ¯2γ−4δ¯δ⁡δ¯γ¯+12γ¯δ¯δ¯2¯+4γ⁡δ¯3−12γ¯δ⁡δ¯2+12δ¯2δ¯γ¯+10δ¯2δ¯γ¯−10δ¯3γ¯−4δ¯δ−δγ¯3¯δ¯+2γ¯δ¯δγ¯2¯+492γ¯3δ¯2+24γ¯2δ¯δ¯γ¯+3γ¯2¯δ¯+8γ¯δγ¯¯−20γ⁡δ¯¯+10δ¯3γ¯δ¯γ¯+4δ⁡δ¯δ¯γ¯−4δ⁡δ¯2γ¯+2δ¯γ¯δ¯δ¯2γ¯−2δ¯2γδ¯2γ¯¯+δ¯5z+563γ¯δ¯3z24γ¯2δ¯z−9δ¯2δγ¯¯z+9δ¯2δγ¯¯¯z+9δ¯2(δ¯γ¯)2z+4δ⁡δ¯z+9δ¯2z2+12δδ¯2γ¯z+10δ¯2δ¯γ¯z+20γ⁡δ¯z−20γ¯δ¯¯z+24γ¯2δ¯z−9δ¯2δγ¯¯z+9δ¯2δγ¯¯¯z+9δ¯2(δ¯γ¯)2z+4δ⁡δ¯z+9δ¯2z2

By using the above mentioned functions of Theorem 2.1, we can obtain Theorem 2.2, under some suitable conditions.

Theorem 2.2 The solution y=0 of (6) has a multiplicity k, wherever 2≤k≤9 if and only if ϰ2=ϰ3=,…,ϰk−1=0 and ϰk≠0 where

ϰ2=∫oσδ,

ϰ3=∫oσγ,

ϰ4=∫oσ(γδ¯+1),

ϰ5=∫oσ(γδ¯2+2δ¯),

ϰ6=∫oσ(2γδ¯3+6δ¯2−δγ¯2+2γ¯),

ϰ7=∫oσ(γδ¯4+4δ¯3−2δδ¯γ¯2+4δ¯γ¯),

ϰ8=∫oσ(γδ¯5+5δ¯4−6δ¯2−δ¯2δγ¯¯γ−δ¯2γ¯+9δ¯2γ¯−2γδ¯δ?−12γ¯3δ−2δ?+32γ¯2),

and

ϰ9=∫oσ(γδ¯6+10δ¯5+959γ¯δ¯3+2903δ¯γ¯δ¯2+20δ¯2¯+503δ¯δ¯2γ¯−103δ¯γ¯2+103δ¯γ¯2+40γ¯δ¯γ¯−1003γ?+150δδ¯2δ¯γ¯−203γ¯δδ¯3γ¯−60γγ¯δ¯δ¯γ¯−70γ¯δ¯γ¯−15δ¯5γ¯γ+40δγ¯δγ¯¯+45δδ¯2γ¯δ¯2γ¯−10γδδ¯γ¯2¯+4109γδ¯3δγ¯¯−403δδ¯2γ¯δγ¯¯−1003δ¯γ?+15δ¯2δγ¯¯−203δ¯δ?⁡γ−1003γδ¯2δ¯2γ¯+70δγ¯2z−6203γ¯δ¯2δz−1003γ¯δ¯z).

Assume that ϰi=ui(σ) , at that point μ=i if ϰ1=1 and ϰk=0 for 2≤k≤i−2 but ϰi≠0 . These ϰi ’s are known as focal values.

Proof By definition, the multiplicity of the zero solution is k if an(σ)=0 for 2≤n≤k−1 and uk(σ)≠0 . Write ξk=uk(σ) and let ϰk be the value of ξk when ξn=0 for n<k .

Since u2(z)=∫oσδ, we have ξ2=∫oσδ . With ξ2=0 , we have ξ3=∫oσγ and ξ4=∫oσ(γδ¯+1) ; hence ϰ3 and ϰ4 as stated. Next, suppose that ξ2,ξ3=0 and substitute the relations ∫oσδ=∫oσγ=0 into the expressions for u5(σ) and u6(σ) given in theorem 2.1; we obtain ϰ4=∫oσ(γδ¯+1),ϰ5=∫oσ(γδ¯2+2δ¯) .

Finally, we suppose that ξ4=0 (as well as ξ2=ξ3=0 ). We substitute the relations ∫oσδ=∫oσγ=∫oσ(γδ¯+1)=0 into the expressions for u7(σ) and u8(σ) ; then ϰ7=∫oσ(γδ¯4+4δ¯3−2δδ¯γ¯2+4δ¯γ¯) ,

and ϰ8=∫oσ(γδ¯5+5δ¯4−6δ¯2−δ¯2δγ¯¯γ−δ¯2γ¯+9δ¯2γ¯−2γδ¯δ−−12γ¯3δ−2δ−+32γ¯2) .

Continuing in the same way we get: ϰ9=∫oσ(γδ¯6+10δ¯5+959γ¯δ¯3+2903δ¯γ¯δ¯2+20δ¯2¯+503δ¯δ¯2γ¯−103δ¯γ¯2+103δ¯γ¯2+40γ¯δ¯γ¯−1003γ?+150δδ¯2δ¯γ¯−203γ¯δδ¯3γ¯−60γγ¯δ¯δ¯γ¯−70γ¯δ¯γ¯−15δ¯$5$γ¯γ+40δγ¯δγ¯¯+45δδ¯2γ¯δ¯2γ¯−10γδδ¯γ¯2¯+4109γδ¯3δγ¯¯−403δδ¯2γ¯δγ¯¯−1003δ¯γ?+15δ¯2δγ¯¯−203δ¯δ?⁡γ−1003γδ¯2δ¯2γ¯+70δγ¯2z−6203γ¯δ¯2δz−1003γ¯δ¯z ).

Now we are going to define the center and some necessary conditions for the center.

Definition An equilibrium point surrounded in its immediate neighborhood (not necessarily over the whole plane) by a closed path is called a center.

Conditions for Center

Corollary 2.1 If any g(z) or d(z) is identically zero and the other has a mean value of zero, then the origin is a center.

Corollary 2.2 For continuously odd differentiable functions g(z) and d(z) of period σ: The origin is a center.

3  Main Results

We consider various classes in which two types of coefficients, namely (z or cos (z) and sin (z)) for g(z) and d(z) in Eq. (6), are considered for the calculations of periodic solutions.

3.1 Trigonometric Coefficients

We have to consider Eq. (6), with coefficients in sin(z) and cos(z) ; here, the period is σ=2π . We now go off to various classes of Eq. (6). For this, we make use of theorem 2.2; to verify the calculations in this module. The related computations are complicated and of a bulk extent. This symbolic representation like ∫02πzcosm(z)sinn(z)dz was overcome using computer algebra programming Maple18.

Lemma 3.1 If possibly g(z) and d(z) are polynomials in sin(z) and cos(z) , all of whose terms are of odd degree. Then multiplicity is μ = 4 or μ > 6.

The above lemma is from Alwash et al. [7], but after solving theorems in the next section, we have concluded that this is also valid if g(z) and d(z) are of even degree, given in the following lemma.

Lemma 3.2 If g(z) and d(z) are polynomials in sin(z) and cos(z) ; all of whose terms are of even degree. Then multiplicity is μ = 4 or μ > 6.

Which are illustrated in the following theorems, throughout this paper, we use the symbol “C” to represent the class and we show the maximum multiplicity for classes with μmax.

Homogeneous Polynomial for Classes C4,4, C5,5 and C6,6

Theorem 3.1 Consider the class C4,4 for Eq. (6). If the coefficients are:

γ(z)=acos4(z)+bcos3(z)sin(z)+csin4(z)+fcos2(z)sin2(z);

ϰ8=(1538563419904)π.

Then we calculate μmax(C4,4)≥8.

Proof By using Theorem 2.2, we now calculate ϰ2=34π(d+e) ,

ϰ3=π(3a+3c+f)4.(12)

From ϰ2,ϰ3=0 , we substitute d=−e,f=−3a−3c and calculate ϰ4=π(16+be)8 . From ϰ4=0, by using the value of “ b ” in terms of “ e ” we calculate ϰ5 as ϰ5=−πe2(a+c)16 .

If ϰ5=0 then, either e=0 or:

a+c=0.(13)

If e=0 then δ(z)=0 and for ϰ3=0 , we can say from corollary 2.1, that the origin is a center. Utilizing (11) we put a=−c and get

ϰ6=π(3e3−4)4e.(14)

with ϰ 7 as: ϰ7=−2cπ9. For ϰ7=0 , we substitute values of “c” and calculate ϰ8 as:

ϰ8=(1538563419904)π.

Hence our conclusion is μmax(C4,4)≥8 .

Theorem 3.2 Let the class C5,5 for Eq. (6); with coefficients as follows:

γ(z)=asin5(z)+bcos3(z)sin2(z)+ccos4(z)sin(z)+dcos5(z),

δ(z)=esin5(z).

Then μmax(C5,5)≥8 is presented.

Proof It is easily calculated that ϰ 2 = ϰ 3 = 0 and:

ϰ4=−π(161eb+689ed−3840)1920.(16)

If ϰ4=0 , then, as π≠0 , we substitute d=3840689e−161689b and calculate ϰ5=0 with ϰ6 as:

ϰ6=e2π(7576305eb+2183801728)1693286400.(17)

If ϰ6=0 , then for e=0 , gives that origin is the center. Thus we put b=−(21838017287576305e) and get ϰ7=(64434587363662640)ecπ.

Now, ϰ7=0 gives c=0 , because eπ≠0 and ϰ8 is as follows:

ϰ8=πϖ10687682598305969912832000e2.(18)

where  ϖ=3463030403249893347931717632+1080907356115092660813824e4a+455749311538699264881025e6 .

Hence, μmax(C5,5)≥8 is calculated.

Theorem 3.3 Consider class C6,6 for Eq. (6), with coefficients as follows:

γ(z)=asin6(z)+ccos5(z)sin(z)+dcos6(z),

δ(z)=ecos6(z)+fsin4(z)cos2(z).

Then μmax(C6,6)≥8 .

Proof From theorem 2.2, we calculated ϰ2=π(f+5e)8, ϰ3=5π(a+d)8. From ϰ2,ϰ3=0 , we substitute a=−d,f=−5e and calculate ϰ4=π(1536+49ce)768 and from ϰ4=0 , by using the value of “c” in terms of “e” we calculate ϰ 5 as:

ϰ5=(652048)de2.(19)

If ϰ5=0 then, as e≠0 , so we put d=0 and get ϰ6=π(−5529600+443401e3)1843968e .

Further, for ϰ6=0 we put =(5529600443401)13 , which results ϰ7=0 and ϰ8 as follows:

ϰ8=1032896577437435266601938594527152689134056056318983503541054008197120000π.(20)

Here it is a nonzero constant number. So, μmax(C6,6)≥8 .

Now, for the polynomial z; we will explain the calculations of the maximum possible periodic solution for many classes of the Eq. (6).

3.2 Polynomial Coefficient

Let Ci,j denotes the classes for the Eq. (6), where i,j denotes the degree for γ(z) and δ(z) accordingly. Using polynomial coefficient “z” for various higher-order classes, we have calculated the possible maximum number of periodic solutions 9 for class C4,1 and 8 for the classes C5,1 , C6,1 , C7,1 . The verification of the presented below theorems stems from papers by Alwash et al. [7,8] and by Saima et al. please see the example, [1–3].

Remark 1. In the perturbation strategies defined in [1], the whole complement of (j − 2) real periodic solutions fails to yield; if it occurs, then there is j<k such that ϰj=0 whenever ϰj−1≠0 . This occurs while the multiplicity is necessarily odd. However, for the number of real periodic solutions, we will say from the “exchange of stability” argument that, If multiplicity μ is even, the origin is stable for ϰμ<0 and unstable if ϰμ > 0. If μ is odd, then the origin is stable on the right and unstable on the left if ϰμ < 0; while it is stable on the left and unstable on the right if ϰμ > 0.

Theorem 3.4 Consider Ci,1 be the class of equations of the form (6) having degree i for γ(z) and 1 for δ(z) , where i=4,5,6,7 . Then we conclude the results μmax(Ci,1)≥8 , for i=5,6,7 and μmax(C4,1)≥9 accordingly.