Effect of Excess Food Nutrient on Producer-Grazer Model under Stoichiometric and Toxicological Constraints

1.   Introduction

Periodic behavior is essential to our life. The mathematical modeling of real-life problems mainly arises in functional equations, such as partial differential equations, ordinary differential equations, and integro-differential equations. It is a matter of the fact that periodic solutions can be created or destroyed only at infinity, or their stability can change. Such qualitative changes in the dynamics are called bifurcation. It is also possible for pairs of limit cycles (one stable and one unstable) to coalesce and vanish in a codimension two generalized Hopf bifurcation (see Guckenheimer and Holmes ^[1]). Many of these functions have the ability to adopt bifurcation, with which we can check the stability. If all the periodic solutions come out of the origin, then these are unstable, otherwise stable (isolated closed trajectory in periodic behavior). Stable limit cycles are an example of attractors. By the Jordan curve theorem, every closed trajectory divides the plane into two regions, an interior region bounded by the curve and an exterior containing all of the nearby and far away exterior points, so that every continuous path connecting a point of one region to a point of the other intersects with that loop somewhere. In the present article, we are mainly concerned with the nonlinear system because trajectories do not simply need to approach or leave a single point. They might even approach a large set, like a circle or another closed curve. It is a fact that limit cycles are a nonlinear phenomenon; because the linear system $x^{^{\prime }} = Ax$ may have closed orbits, they would not be isolated, such that if they have a periodic solution $x(t)$, then all constant multiples should be, see also ^[2].

In engineering, power engineering, we can exploit oscillators for timing and sequencing. We can use them in locomotion control in robots. Van der Pol oscillator comes up in electrical circuits; it is an example of relaxation oscillation. Its equation is $x^{^{\prime \prime }}-\mu \left(1-x^{2}\right) x^{^{\prime }}+x = 0,$ where $\mu$ is a positive constant. If $\mu = 1,$ then trajectories fastly settle on a closed curve, while for $\mu = 0,$ all trajectories become a circle. In analyzing nonlinear systems in the xy-plane, we have so far concentrated on analyzing how the system's trajectories look in the neighborhood of each critical point.

We also observed periodic behavior in our body like; heartbeats, breathing, chewing, locomotion, diverse rhythms inside the brain, etc. More generally in walking, which is also a periodic, a passive stable limit cycle exists. as shown in Figure 1, see ^[3] for more details.

Consider the system of the form:

where $a_{w}$ and $b_{w}$ are homogenous polynomials having degree $w$. The polar form of $\left(1.1\right)$ is written as following:

Where $x$ and $y$ are polynomials in $\cos \theta$ and $\sin \theta$ with degree $\left(w+1\right).$ In Lloyd ^[4], it is presented that if:

Then the Eq $\left(1.3\right)$ can be transformed to non-autonomous first order differential equation as:

where

and

Here $\kappa$ & $\rho$ are homogeneous polynomials in $\sin \theta$ and $\cos \theta.$

We structured the paper as follows. In Section 2, we have discussed the transformation and some essential formulae to calculate periodic multiplicity. In Section 3, we recall some lemma and theorems for the origin. The significant results and conclusions are discussed in Sections 4 and 5 accordingly.

2.   First order cubic system

In the present article, we are considering the differential equation given as:

Where the coefficients $\kappa, \rho, \upsilon$ and involved variable are real-valued functions and $U\in \mathbb{C}$. This paper's central result is the determination of possible maximum periodic solutions with perturbation of the coefficients on the plane. This equation is the part of the following equation, described in ^[5]:

With $\rho _{0}(\tau) = 1$. For $n = 3,$ the Eq$\ \left(\text{2.2}\right)$ is known as Abel's differential equation, we focused on it because of its connection with the second part of Hilbert 16$^{th}$ problem (maximal number of limit cycles and their relative locations of planar polynomial real vector fields of given degree); it is related to ODEs and dynamical systems. The fascination of the problem comes from the fact that it sits at the confluence of analysis, algebra, geometry and even logic. It is known, for instance, that when $\rho _{3}(t)$ does not change sign, the upper bound for the number of limit cycles is 3, see ^[6,7]. When $\rho _{3}(t)\equiv 1$ this upper bound also holds taking into account complex limit cycles, see ^[5,8]. In ^[6], we see that when $\rho _{0}(t)\equiv 0$ and $\rho _{2}(t)$ does not change sign, the maximum number of limit cycles of Abel's equation is also 3. We likewise refer reader to the articles ^[9,10,11]$,$ for additional data with respect to this issue. For the Eq $\left(\text{2.1}\right),$ we are considering a complex dependent variable, so that the number of zeros of a function in a bounded region of the complex plane can't be changed by any small perturbations. We substitute $\upsilon \left(\tau \right)$ $\widetilde{ = }$ $0$ in $\left(\text{2.1}\right),$ as was in ^[12]. Consequently, the Eq $\left(\text{2.1}\right)$ becomes as:

Here $\kappa$ and $\rho$ may be polynomials in $\left(\text{i}\right)$ $\tau$ $\left(\text{ii}\right)$ $cos\tau$ and $sin\tau$, for more detail see ^[12,13]. Also suppose that $\exists$ $\beta \in \mathbb{R}$ such that:

For $U = 0$, the method for computing multiplicity $"\mu "$ is explained in ^{[9,10,12,14,15]}$.$ For $t$ in $[0, \beta]$ and $c$ small, let

Since $U\left(0, 0, c\right) = c,$ we have $\xi _{1}\left(0\right) = 1$ and $\xi _{i}\left(0\right) = 0$ for $i > 1$. The solution $U = 0$ is a centre for the Eq $\left(2.3\right)$ if all the solutions are periodic for $c$ in neighborhood of $0$. The solution $U = 0$ is a centre if$\ \xi _{1}\left(\beta \right) = 1$ and $\xi _{i}\left(\beta \right) = 0$ for $i > 1$. Substituting the Eq $\left(2.4\right)$ in the Eq $\left(2.3\right)$ and equating coefficients of $c$ yields $\overset{\cdot }{\xi }_{1}\left(\tau \right) = 0.$ Hence, $\xi _{1}\left(\tau \right) = 1.$ Moreover, the functions $\overset{\cdot }{\xi } _{i}\left(\tau \right),$ for $i > 1$ are obtained with the help of following equation:

With $\xi _{1}\left(\tau \right) = 1$. Expect to be that $\mathcal{\varkappa }_{i} =$ $\xi _{i}\left(\beta \right),$ at that point $\mathcal{\varkappa } = i$ if $\mathcal{\varkappa }_{1} = 1$ and $\mathcal{\varkappa }_{k} = 0$ for $2\leq k\leq i-2$ but $\mathcal{\varkappa }_{i}\neq 0$, shown in Theorem 2.1. Alwash in ^[12]$,$ presented $\xi _{i}\left(\tau \right)$ and $\mathcal{\varkappa }_{i}$ for $i\leq 8,$ for $i = 9$ are in ^[16]$,$ for $i = 10$ we calculated $\xi _{10}\left(\tau \right)$ and $\mathcal{\varkappa }_{10}$ in ^[9], also presented in Theorem 2.1$.$

The following Theorem is the modification of Theorem 2 in ^[12], with the help of this Theorem periodic multiplicity is calculated. Here, in integral; $\int \kappa \left(\tau \right) \overline{\rho \left(\tau \right) }d\tau$, bar $"-"$ function is like $\overline{\rho \left(\tau \right) } = \int \rho \left(\tau \right) d\tau$.

Theorem 2.1. The solution $U = 0$ of the Eq $\left(\mathit{\text{2.3}}\right)$ has a multiplicity $k,$ wherever $2\leq k\leq 10$ if $\mathcal{ \varkappa }_{n} = 0$ for $2\leq n\leq k-1$ and $\mathcal{\varkappa } _{n}\neq 0$ where

$\mathcal{\varkappa }_{2} = \int_{0}^{\beta }\rho,$

$\mathcal{\varkappa }_{3} = \int_{0}^{\beta }\kappa,$

$\mathcal{\varkappa }_{4} = \int_{0}^{\beta }\kappa \overline{\rho },$

$\mathcal{\varkappa }_{5} = \int_{0}^{\beta }\kappa \overline{\rho }^{2},$

$\mathcal{\varkappa }_{6} = \int_{0}^{\beta }\left(\kappa \overline{\rho } ^{3}-\frac{1}{2}\overline{\kappa }^{2}\rho \right),$

$\mathcal{\varkappa }_{7} = \int_{0}^{\beta }\left(\kappa \overline{\rho } ^{4}+2\kappa \overline{\rho }^{2}\overline{\kappa }\right),$

$\mathcal{\varkappa }_{8} = \int_{0}^{\beta }\left(\kappa \overline{\rho } ^{5}+3\kappa \overline{\rho }^{3}\overline{\kappa }+\kappa \overline{\rho } ^{2}\overline{\overline{\rho }\kappa }-\frac{1}{2}\overline{\kappa }^{3}\rho \right),$

$\mathcal{\varkappa }_{9} = \int_{0}^{\beta }\left(\kappa \overline{\rho } ^{6}-5\kappa \overline{\rho }^{4}\overline{\kappa }-2\overline{\rho }^{3} \overline{\rho \overline{\kappa }}+20\overline{\rho \overline{\kappa }^{2}}+2 \overline{\rho \overline{\kappa }}\rho \overline{\kappa }^{2}\right),$

$and$

$\mathcal{\varkappa }_{10} = \int_{0}^{\beta }\left(\begin{array}{c} \kappa \overline{\rho }^{7}-\frac{1235}{6}\overline{\kappa }\kappa \overline{ \rho }^{5}-\frac{970}{3}\kappa \overline{\kappa }^{2}\overline{\rho } ^{3}-237\rho \overline{\rho }^{2}\overline{\kappa }^{3}-24\kappa \overline{ \kappa }^{2}\rho \overline{\rho }^{2}-70\overline{\rho }^{3}\kappa ^{2} \overline{\kappa }-21\overline{\kappa }^{4}\rho -74\kappa \overline{\kappa } ^{3}\overline{\rho }+ \\ \frac{5}{2}\overline{\kappa }^{2}\rho \overline{\rho }^{4}+32\overline{\rho } ^{4}\kappa \overline{\rho \overline{\kappa }}-16\overline{\rho }^{4}\rho \overline{\kappa }-15\overline{\rho }^{5}\kappa ^{2}-36\overline{\rho }\rho \overline{\kappa }^{2}\overline{\rho \overline{\kappa }}-8\rho \overline{ \rho }^{4}\kappa \overline{\kappa } \end{array} \right).$

3.   Conditions for centre

For $U = 0$ as a centre, conditions that are useful for calculating maximum multiplicity $\mathcal{\varkappa }_{k}, \ 2\leq k\leq 10$ are from ^[12] and are defined below.

Theorem 3.1. Consider that $f,$ $g$ are continuous functions defined on interval $I = \delta \left(\left[0, \beta \right] \right)$ and a differentiable function $\delta$ with$\ \delta \left(\beta \right) = \delta \left(0\right)$ such that:

then origin is the centre for the Eq (2.3).

Corollary 3.1. If $\kappa$ is a constant multiple of $\rho$ and $\int_{0}^{\beta }\rho \left(\tau \right) d\tau = 0,$ then the origin is a centre for the Eq (2.3).

Corollary 3.2. If any $\rho$ or $\kappa$ is identically zero and other has mean value zero then the origin is a centre.

Remark 1. In ^[12], bifurcation method is described that when the coefficients of $\rho (\tau)$ and $\kappa (\tau)$ are slightly perturbed, two periodic solutions bifurcate out of the origin, when bifurcation method is applied. For the number of real periodic solutions, we conclude that if multiplicity $\mu$ is even, the origin is stable $\varkappa _{\mu } < 0$ and unstable if $\varkappa _{\mu } > 0.$ If $\mu$ is odd, then the origin is stable on the right and unstable on the left if $\varkappa _{\mu } < 0;$ however, if $\varkappa _{\mu } > 0,$ the origin is stable on the left and unstable on the right.

4.   Results

4.1.   Polynomial coefficients

This section describes the method for computing the maximum number of limit cycles of a polynomial differential equation in a plane for various classes of different degrees. Suppose $C_{r, q}$ indicates the class for the Eq $\left(\text{2.3}\right),$ with degree$\ r$, $q$ for $\kappa \left(\tau \right)$ and $\rho \left(\tau \right)$ accordingly, for more examples see ^{[9,10,13,14,17]}. The confirmation of the accompanying theorems, stems from papers in ^[12,18]$.$ We use Theorem 2.1 with $\beta = 1,$ as is done in Lloyd et al. If we use all coefficients for polynomials $\kappa \left(\tau \right)$ and $\rho \left(\tau \right)$ then we can easily see that the periodic solutions greater than 4 can't obtained. So, some possible suitable coefficients are restricted in the following classes to find as many periodic solutions as possible. All calculations regarding the different classes are carried out using Maple 18.

Theorem 4.1. Suppose the class $C_{9, 4}$ for the Eq (2.3), if

Then we come to conclusions $\mu _{\max }\left(C_{9, 4}\right) \geq 10.$

Proof. From Theorem $2.1,$ we extract that:

Thus multiplicity of $U = 0$ is $\mu = 2,$ if $\mathcal{\varkappa }_{2}\neq 0.$ And is $\mu = 3,$ if $\mathcal{\varkappa }_{2} = 0$ but $\mathcal{ \varkappa }_{3}\neq 0$. If $\mathcal{\varkappa }_{2} = \mathcal{\varkappa } _{3} = 0,$ then $\kappa (\tau)$ and $\rho (\tau)$ are as below:

And also we compute $\mathcal{\varkappa }_{4}$ as given below:

If $\mathcal{\varkappa }_{4} = 0,$ then either $q = 0$ or:

If $q = 0,$ then $\rho (\tau) = 0$ and for $\mathcal{\varkappa }_{3} = 0,$ origin is the centre derived from Corollary $3.2.$ So, assuming $\left(4.3 \right)$ holds and $q\neq 0$, $\mathcal{\varkappa }_{5}$ is given as:

If $\mathcal{\varkappa }_{5} = 0,$ as $q\neq 0$ implies:

And by using $\left(4.4\right)$ we take $\mathcal{\varkappa }_{6}$ as:

If $\mathcal{\varkappa }_{6} = 0,$ then, as $q\neq 0$ either $f = -6b$ or

If $f = -6b,$ then the Eqs $\ \left(4.2\right)$ and $\left(4.1\right)$ become:

Where, $\overset{\cdot }{\delta }\left(\tau \right) = \tau ^{4}-\frac{1}{5}$, also $\delta \left(0\right) = \delta \left(1\right).$ As a result of Theorem $3.1,$ the origin is the centre, as follows:

So, $f\neq -6b.$ If $\left(4.5\right)$ holds then $\mathcal{\varkappa }_{7}$ is computed as:

If $\mathcal{\varkappa }_{7} = 0,$ since $q\neq 0$, either $b = -\frac{6}{25} q^{2}$ or,

If $b = -\frac{6}{25}q^{2}$ then,

The origin is the centre with $f\left(\delta \right) =$ $\left[e+q^{2}\left(-\frac{6}{5}\tau ^{5}+6\tau \right) \right]$ and $g\left(\delta \right) = q$ according to Theorem $3.1$. So, consider $b\neq -\frac{6}{25}q^{2}.$ Using $\left(4.6\right)$ we calculate $\mathcal{\varkappa }_{8}$ as:

Where,

$\phi = -3325005947550744260462549014528q^{4}-345198526986114785414664480eq^{2}+21339337823946954208255625e^{2}.$

Now, if $\mathcal{\varkappa }_{8} = 0$ then either $\phi = 0$ or

As $q\neq 0$. If $\left(4.7\right)$ holds but $\phi \neq 0, q\neq 0$, $\mathcal{\varkappa }_{9}$ is calculated as follows:

If $\mathcal{\varkappa }_{9} = 0$ then, as $q\neq 0,$ it results $q^{5}\neq 0$, we examine $q$ as:

If Eq $\left(4.7\right) \neq 0,$ $q\neq 0,$ but $\phi = 0$ holds then $e = y_{i}m^{2}$ for $i = 1, 2$ with $y_{1} = 133.1766309,$ $y_{2} = -116.9999242.$ If $\left(4.8\right)$ holds, we can calculate $\mathcal{\varkappa }_{10}$ as:

Here $\mathcal{\varkappa }_{10}$ is not zero. As a result, we can deduce that $\mu _{\max }\left(C_{9, 4}\right) \geq 10$.

For non zero different but the same (either positive or negative) values of the constants in $\kappa$ and $\rho,$ we can see that only one root is real and rest all the zeros are complex and are also in conjugate pairs. It gives stable limit cycles from Remark 1. The stability analysis is shown in Figure 2.

Theorem 4.2. For given below equation:

Consider:

$\kappa (\tau) = \frac{167936}{2352625}(-\frac{286469068411235}{2477977842432} +\epsilon _{1})^{2}+\frac{1989561074}{412868500465}\epsilon _{2}-\frac{641861 }{103320}\epsilon _{3}-\frac{301}{2976}\epsilon _{4}-\frac{17}{945}\epsilon _{5}-\frac{1}{10}\epsilon _{6}+\epsilon _{7}+\\ (-\frac{6}{25}(-\frac{ 286469068411235}{2477977842432}+\epsilon _{1})^{2}-\frac{2722420008}{ 82573700093}\epsilon _{2}+\epsilon _{3})\tau +(-\frac{286469068411235}{ 2477977842432}+ \epsilon _{1})^{2}+\epsilon _{2})\tau ^{4}+(\frac{36}{25}(- \frac{286469068411235}{2477977842432}+\\ \epsilon _{1})^{2}-\frac{754408015152}{ 412868500465}\epsilon _{2}+\frac{1670138}{30135}\epsilon _{3}+\epsilon _{4})\tau ^{5}+(\frac{240204612285}{82573700093}\epsilon _{2}-\frac{1418065}{ 16072}\epsilon _{3}-\frac{1425}{992}\epsilon _{4}+\epsilon _{5})\tau ^{8}+(- \frac{170939097560}{82573700093}\epsilon _{2}+\\ \frac{162185}{2583}\epsilon _{3}+\frac{175}{186}\epsilon _{4}-\frac{176}{189}\epsilon _{5}-\frac{6}{5}(- \frac{286469068411235}{2477977842432}+\epsilon _{1})^{2}+\epsilon _{6})\tau ^{9},$

Choose $\epsilon _{p}\neq 0$ for $1\leq p\leq 8$ and small as compared to $\epsilon _{p-1}$. Then there are eight non-trivial real periodic solutions to Eq $\left(\mathit{\text{4.9}}\right)$.

Proof. See ^[14].

Theorem 4.3. For class $C_{9, 5}$ with

Then we conclude $\mu _{\max }\left(C_{9, 5}\right) \geq 8,$ where $A_{1} = B_{1} = \left(\tau -1\right).$

Proof. By utilizing Theorem $2.1,$ we calculate:

If $\mathcal{\varkappa }_{2} = \mathcal{\varkappa }_{3} = 0$, we calculated $\mathcal{\varkappa }_{4}$ as:

If $\mathcal{\varkappa }_{4} = 0$ then, either $r = 0$ or:

If $r = 0$ then, $\mathcal{\varkappa }_{2} = 0$ gives $m = 0,$ hence, $\rho (\tau) = 0.$ If $\mathcal{\varkappa }_{2} = 0,$ using Corollary $3.2,$ origin is the centre. As a result, take $r\neq 0$. By using $\left(4.12 \right)$ we compute:

If $\mathcal{\varkappa }_{5} = 0$ then:

Because we have already suppose $r\neq 0.$ If $\left(4.13\right)$ holds then:

If $\mathcal{\varkappa }_{6} = 0$, either $b = 0$ or:

as $r\neq 0.$ If $b = 0$, the Eqs $\left(4.10\right)$ and $\left(4.11\right)$ becomes:

Let $\delta \left(\tau \right) = \frac{\left(\tau -1\right) ^{6}}{6}+\frac{ \tau }{6}$ then, $\overset{\cdot }{\delta }\left(\tau \right) = \left(\tau -1\right) ^{5}+\frac{1}{6}$. Also $\delta \left(0\right) = \delta \left(1\right).$ Using these results we write as below:

From Theorem $3.1,$ origin is the centre with $f\left(\delta \right) =$ $f$ and $g\left(\delta \right) = r.$ So, suppose $b\neq 0.$ If $\left(4.14 \right)$ holds then $\mathcal{\varkappa }_{7}$ is:

Now if $\mathcal{\varkappa }_{7} = 0,$ as $r\neq 0$ then:

Using $\left(4.15\right)$ we obtain:

Thus, $\mu _{\max }\left(C_{9, 5}\right) \geq 8.$

For non zero different but the same (either positive or negative) values of the constants in $\kappa$ and $\rho,$ we can see that only two roots are real and rest all the zeros are complex and are also in conjugate pairs. It gives stable limit cycles from Remark 1. The stability analysis is as shown in Figure 2.

Theorem 4.4. Consider class $C_{9, 6}$ for the Eq $\left(\mathit{\text{2.3}} \right)$, if

Then we conclude $\mu _{\max }\left(C_{9, 6}\right) \geq 8$ with $A_{2} = B_{2} = \left(2\tau -1\right).$

Proof. By utilizing Theorem $2.1,$ we obtain:

Thus, multiplicity of $U = 0$ is $\mu = 2$ if $\mathcal{\varkappa }_{2}\neq 0,$ And is $\mu = 3$ if $\mathcal{\varkappa }_{2} = 0$ but $\mathcal{\varkappa }_{3}\neq 0$. If $\mathcal{\varkappa }_{2} = \mathcal{\varkappa }_{3} = 0$ then:

By using Eqs $\left(4.16\right)$ and $\left(4.17\right),$ we have:

Also we calculate $\mathcal{\varkappa }_{4}$ as:

If $\mathcal{\varkappa }_{4} = 0$ then, either $s = 0$ or:

If $s = 0$ then, $\left(4.17\right)$ gives $m = 0$ so $\rho (\tau) = 0,$ and also for $\mathcal{\varkappa }_{2} = 0,$ origin is the centre. So, $s\neq 0$ is taken. If $\left(4.20\right)$ holds $\mathcal{\varkappa }_{5}$ is given as:

If $\mathcal{\varkappa }_{5} = 0,$ then, $e = 0$ because we had already seen that $s\neq 0$. Thus, by substituting $e = 0,$ we have:

Now if $\mathcal{\varkappa }_{6} = 0,$ as $s\neq 0$ either $b = 0$ or:

If we take $b = 0$ gives $j = 0$ and $e = 0$ from $\mathcal{\varkappa }_{5}$ then, Eqs $\left(4.18\right)$ and $\left(4.19\right)$ take given below form:

and:

Let $\delta \left(\tau \right) = \frac{\left(2\tau -1\right) ^{7}}{14}- \frac{\tau }{7},$ Also $\delta \left(0\right) = \delta \left(1\right).$ So:

From Theorem $3.1$, origin is the centre, with $f\left(\delta \right) =$ $g$ and $g\left(\delta \right) = s$. By using $\left(4.21\right)$, we have $\mathcal{\varkappa }_{7}$ as:

If $\mathcal{\varkappa }_{7} = 0,$ recalling that $s\neq 0$ (considered above) then by substituting $g = 0$, we found:

As $s\neq 0$ considered above, we can't proceed further. So, $\mu _{\max }\left(C_{9, 6}\right) \geq 8.$

For non zero different but the same (either positive or negative) values of the constants in $\kappa$ and $\rho,$ we can see that only one root is real and rest all the zeros are complex conjugate pairs. It gives stable limit cycles from Remark 1, as shown in Figure 2.

Theorem 4.5. If for the Eq $\left(\mathit{\text{2.3}}\right),$

With $A_{1} = \tau -1.$

and

For $\epsilon _{l}\neq 0, \ \left(1\leq l\leq 6\right)$ if we take $\left\vert \epsilon _{6}\right\vert \ll \left\vert \epsilon _{5}\right\vert \ll...\ll \left\vert \epsilon _{1}\right\vert.$ Then there are six different real periodic solutions.

Proof. See ^[9].

Theorem 4.6. For the class $C_{11, 3}$ consider that

Then $\mu _{\max }\left(C_{11, 3}\right) \geq 10.$

Proof. We write $w = a+\frac{1}{2}b+\frac{1}{4}d+\frac{1}{5}e+\frac{1}{8}h+\frac{1}{12 }k,$ $w_{1} = m+\frac{1}{4}p;$ using Theorem $2.1,$ we see $\mathcal{ \varkappa }_{2} = w_{1},$ $\mathcal{\varkappa }_{3} = w.$ If $w, w_{1} = 0,$ then by putting 'a' and 'm', $\mathcal{\varkappa }_{4}$ is as shown below:

Now, for $\mathcal{\varkappa }_{4} = 0,$ we put

and calculate

From $\mathcal{\varkappa }_{5} = 0,$ by substituting value of h, $\mathcal{ \varkappa }_{6}$ is

Here $\varpi = e+5b.$ If $\mathcal{\varkappa }_{6} = 0$; as $p\neq 0,$ using value of e from $2574821080b-154952704e+1256094225p^{2} = 0$, $\mathcal{\varkappa }_{7}$ is

with $\varpi _{1} = 3p^{2}+8b.$ For $\mathcal{\varkappa }_{6}, \mathcal{ \varkappa }_{7} = 0,$ suppose that $\varpi = \varpi _{1} = 0,$ (as possible) then for $e = -5b, \ b = -\frac{3}{8}p^{2};$ $\kappa \left(t\right)$ and $\rho \left(t\right)$ can be written as below:

Let $\delta \left(\tau \right) = \frac{\tau ^{4}}{4}-\frac{\tau }{4}$ then, $\overset{\cdot }{\delta }\left(\tau \right) = \tau ^{3}-\frac{1}{4}$, also $\delta \left(0\right) = \delta \left(1\right).$ So, above equations may written as

As a result of Theorem $3.1,$ the origin is the centre with $f\left(\delta \right) = \left[d+b\left(4\tau ^{4}-4\tau \right) \right],$ $g\left(\delta \right) = p;$ and for $b = -\frac{3}{8}p^{2}$ having $f\left(\delta \right) = \left[d+p^{2}\left(\frac{3}{2}\tau ^{4}-\frac{3}{2}\tau \right) \right],$ $g\left(\delta \right) = p.$ Thus, suppose $e\neq -5b, \ b\neq -\frac{3}{8} p^{2}$. Now, we put

and obtain

where $\varpi _{2} = -13428577611822818269376896810070582196773931p^{4}-17734094726830067851160472587240248614dp^{2}+2205020838128046275374234026162794502d^{2}.$

If $\mathcal{\varkappa }_{8} = 0,$ then we substitute $d = -\frac{15853669}{ 40748730}p^{2}$ and obtain $\mathcal{\varkappa }_{9}$ as follows:

If $\mathcal{\varkappa }_{9} = 0,$ recalling that $p\neq 0,$ we put $p = -\frac{ 2120781350311700544}{21005802551120299},$ and calculate

It is a non-zero constant number. Therefore, it is concluded that $\mu _{\max }\left(C_{11, 3}\right) \geq 10.$

For non zero different but the same (either positive or negative) values of the constants in $\kappa$, and $\rho,$ we can see that one root is real and others are complex. It gives stable limit cycles from Remark 1, as shown in Figure 2.

Corollary 4.1. For the equation:

If the polynomials $\kappa (\tau),$ $\rho (\tau)$ are as used in Theorem 4.2 and Theorem 4.6$.$ If $\upsilon,$ $\upsilon _{1}$ are enough small, then $\left(\mathit{\text{4.22}}\right)$ has ten real periodic solutions.

Theorem 4.7. For the Eq (2.3) consider that

With

$\upsilon _{1} = (-\frac{2120781350311700544}{21005802551120299}+\epsilon _{1})^{2}+\frac{4709898985434859}{434993315743423814}~~\epsilon _{2}-\frac{ 4224512475}{1491419776}~~\epsilon _{3}-\frac{8667}{66143}\epsilon _{4}-\frac{ 115}{2376}\epsilon _{5}-\frac{1}{12}\epsilon _{6}+\epsilon _{7}$

$b = -(-\frac{2120781350311700544}{21005802551120299}+\epsilon _{1})^{2}-\frac{ 20027618082927840}{217496657871711907}~~\epsilon _{2}+\epsilon _{3},$

$d = -(-\frac{2120781350311700544}{21005802551120299}+\epsilon _{1})^{2}+\epsilon _{2},$

$h = (-\frac{2120781350311700544}{21005802551120299}+\epsilon _{1})^{2}+\frac{ 1011560654078040}{984147773175167}~~\epsilon _{2}-\frac{47294663}{4236988} \epsilon _{3}-\frac{21088}{30065}\epsilon _{4}+\epsilon _{5},$

$e = -(-\frac{2120781350311700544}{21005802551120299}+\epsilon _{1})^{2}-\frac{ 102398557040420025}{66922048575911356}\epsilon _{2}+\frac{321852635}{19369088 }\epsilon _{3}+\epsilon _{4},$

and

$k = -\frac{439851904713900}{984147773175167}~~\epsilon _{2}+\frac{64632555}{ 13316248}\epsilon _{3}+\frac{10608}{47245}\epsilon _{4}-\frac{91}{99} \epsilon _{5}+\epsilon _{6}.$

If $\epsilon _{l}$ for $\left(1\leq l\leq 8\right),$ have the property that

Then there exists eight distinct non-trivial periodic solutions.

4.2.   Trigonometric coefficients

Now, we consider the Eq $\left(2.3\right)$, with polynomials $\kappa (\tau)$, $\rho (\tau)$ in $\sin \tau$ and $\cos \tau$; here upper limit $\beta$ is $2\pi$ (period of these trigonometric functions) for Theorem 2.1.

4.2.1.   Non-homogeneous coefficients

Theorem 4.8. Consider the class $C_{14, 7},$ If the coefficients are:

Then $\varkappa _{\max }\left(C_{14, 7}\right) \geq 9.$

Proof. Using Theorem $2.1,$ we calculate

If $\varkappa _{5} = 0$, either $a, b = 0$ or:

If $a, b = 0$ then $\rho (\tau) = 0,$ and $\varkappa _{3} = 0,$ gives origin is the centre$.$ Using $\left(4.23\right)$ we calculate $\varkappa _{6} = 0$ and $\varkappa _{7}$ as:

If $\varkappa _{7} = 0$ as $abc\neq 0,$ we put $a = -b,$ and calculate $\varkappa _{8} = 0,$ and $\varkappa _{9}$ as shown below:

If $\varkappa _{9} = 0,$ for $c, b = 0$ the origin is the centre. Hence, concluded that $\mu _{\max }\left(C_{14, 7}\right) \geq 9.$

Theorem 4.9. For the Eq $\left(2.3\right),$ consider coefficients as:

Then $\mu _{\max }\left(C_{16, 8}\right) \geq 9$ is presented$.$

Proof. With Theorem $2.1,$ we calculate $\varkappa _{2} = 0$ and:

If $\varkappa _{2} = \varkappa _{3} = 0,$ as $\pi \neq 0,$ we calculate $\varkappa _{4}$ as:

For $\varkappa _{4} = 0,$ as $\pi \neq 0,$ either $e = 0$ or $a = -b$. If we substitute $e = 0$ then $\kappa (\tau) = 0,$ and $\varkappa _{2} = 0$, origin is centre. Hence suppose that $e\neq 0.$ We substitute for $a = -b$ and take $\varkappa _{5} =$ $\varkappa _{6} = \varkappa _{7} = \varkappa _{8} = 0$ with $\varkappa _{9}$ as:

If $\varkappa _{9} = 0,$ then $c, b = 0$ it can be easily seen that origin is the centre. Hence, concluded that $\mu _{\max }\left(C_{16, 8}\right) \geq 9.$

4.2.2.   Homogeneous coefficients

Theorem 4.10. Let the class$\ C_{5, 5}$ with:

Then $\mu _{\max }\left(C_{5, 5}\right) \geq 10$ is calculated as follows$.$

Proof. Using Theorem $2.1,$ $\varkappa _{2} = \varkappa _{3} = 0$ and by proceeding further we calculate as:

If $\varkappa _{4} = 0,$ either $m = 0$ or:

If $m = 0,$ then $\rho \left(\tau \right) = 0$ and for $\varkappa _{3} = 0,$ $\kappa \left(\tau \right)$ has mean value as zero. So, origin is a centre, by Corollary 3.2. Also with holding the Eq $\left(4.24 \right),$ we calculate$\ \varkappa _{5} = 0$ and $\varkappa _{6}$ as:

For $\varkappa _{6} = 0,$ we substitute value of $c$ as:

and obtained:

If $\varkappa _{7} = 0,$ then we substitute $a = 0,$ because $fm^{2}\pi \neq 0$ and get:

For $\varkappa _{8} = 0,$ we put $f^{2} = -\frac{ 46125464931}{3146842112}m^{4}$ and obtained $\varkappa _{9} = 0$ whereas $\varkappa _{10}$ comes out as:

$\varkappa _{10} = -\frac{ 1840750191m^{9}\pi(59976855709863203909055069m^{6}+2195549986673047538458624m^{2}-1043525219711433352151040)}{102111532729282706230029677127270400}.$

Hence we conclude that $\mu _{\max }\left(C_{5, 5}\right) \geq 10.$

Theorem 4.11. For the class $C_{6, 6},$ if:

Then $\mu _{\max }\left(C_{6, 6}\right) \geq 10$ is given below$.$

Proof. By using Theorem $2.1,$ it is calculated that $\varkappa _{2} = 0,$ $\varkappa _{3} = \frac{5\pi (e+a)}{8}$ and by substituting $a = -e$ from $\varkappa _{3} = 0,$ we calculate $\varkappa _{4}$ as:

If $\varkappa _{4} = 0$ then, as $\pi$ is nonzero, either $e = 0$ or:

If $e = 0,$ then $e, a = 0$, gives $\kappa \left(\tau \right)$ and$\ \rho (\tau)$ are same with only different coefficients; i.e. $\kappa \left(\tau \right)$ is constant multiple of $\rho (\tau).$ As a result of Corollary 3.1, the origin is a centre. So, suppose $e\neq 0.$ Also with holding Eq $\left(4.25\right),$ we calculate$\ \varkappa _{5} = 0$ and $\varkappa _{6}$ as:

If $\varkappa _{6} = 0,$ consider that $g, e\neq 0,$ because for $g = 0$, the Eq $\left(4.25\right)$ gives $i = 0$; so $\rho \left(\tau \right) = 0.$ For $\varkappa _{3} = 0,$ origin is the centre. By substituting $d = -b,$ we calculate $\varkappa _{7} = 0,$ $\varkappa _{8} = 0$ and by proceeding further, we get $\varkappa _{9}$ as:

If $\varkappa _{9} = 0,$ then the only possible substitution is as follows:

With holding the Eq $\left(4.26\right),$ $\varkappa _{10}$ comes out as:

$\varkappa _{10} = -\frac{ \begin{array}{c}g^{3}b\pi (361424202582272bg^{3}-20937295726323g^{5}+ 341930657292484608b^{2}- \\ 18978246855069696bg^{2}) \end{array} }{30247865046792142848}.$

Now, due to lack of formula for $\varkappa _{11},$ we can't proceed further to calculate focal value greater than 10. Hence we conclude that $\mu _{\max }\left(C_{6, 6}\right) \geq 10.$

5.   Further examples

In this section, we have discussed various examples. Here, we are considering the series expansion of some functions. For finding the maximum possible periodic solutions, we shall restrict these expansions up to some extent, for $n = 5$ (say). All the coefficients used below in series expansion are the equal to 1, in general.

Example 5.1. For the Eq (2.3), suppose the series expansion upto power five of the functions; $\kappa \left(t\right) = \frac{1 }{1+t} = c-dt+et^{2}-ft^{3}+gt^{4}-ht^{5}$ and $\rho \left(t\right) = \tanh ^{-1}\left(t\right) = it+\frac{1}{3}kt^{3}+\frac{1}{5}mt^{5}.$ Then we calculate the periodic multiplicity.

Solution 5.1. For finding maximum periodic solutions, we make suitable restrictions of the coefficients and put $e, k = 0$. We write $w = c-\frac{1}{2}d-\frac{1}{4}f+\frac{ 1}{5}g-\frac{1}{6}h,$ $w_{1} = \frac{1}{2}i+\frac{1}{30}m;$ using Theorem 2.1$,$ we see $\mathcal{\varkappa }_{2} = w_{1},$ $\mathcal{\varkappa } _{3} = w.$ If $w, w_{1} = 0,$ then by putting 'c' and 'i', $\mathcal{\varkappa }_{4}$ is as shown below:

Now, for $\mathcal{\varkappa }_{4} = 0,$ as $m = 0$ gives origin is the centre. We put

and calculate

From $\mathcal{\varkappa }_{5} = 0,$ by substituting $g = \frac{7837}{4096}f,$ $\mathcal{\varkappa }_{6}$ is

As for $f, m = 0,$ origin is centre, using value of f, $\mathcal{\varkappa } _{7}$ is

For $\mathcal{\varkappa }_{7} = 0,$ by substituting 'd' in terms of 'm$^{2}$ ' we obtain:

Thus, the multiplicity is 8.

Example 5.2. For the Eq (2.3), suppose the series expansion upto power five of the functions; $\kappa \left(t\right) = \ln \left(1+t\right) = ct-\frac{e}{2}t^{2}+\frac{f}{3}t^{3}-\frac{k}{4}t^{4}+ \frac{h}{5}t^{5}$ and $\rho \left(t\right) = \sinh \left(t\right) = it+\frac{ 1}{3!}lt^{3}+\frac{1}{5!}mt^{5}.$ Then we calculate the periodic multiplicity.

Solution 5.2. For finding maximum periodic solutions, we make suitable restrictions of the coefficients and put $l = 0$. We write $w = \frac{1}{2}c-\frac{1}{3}e-\frac{1}{4 }f+\frac{1}{5}k-\frac{1}{6}h,$ $w_{1} = \frac{1}{2}i+\frac{1}{30}m;$ using Theorem 2.1$,$ we see $\mathcal{\varkappa }_{2} = w_{1},$ $\mathcal{ \varkappa }_{3} = w.$ If $w = 0 = w_{1},$ then by putting values of 'c' and 'i', $\mathcal{\varkappa }_{4}$ is as shown below:

If $\mathcal{\varkappa }_{4} = 0,$ we put $k = \ \frac{77}{27}f-\frac{154}{27} e,$ and calculate

Now, if $\mathcal{\varkappa }_{5} = 0,$ we substitute $f = \frac{96256}{18293} e$ and obtained

From $\mathcal{\varkappa }_{6} = 0,$ if $e, m = 0$, origin is the centre. So, we put value of 'e' and calculate

and

Thus, the multiplicity is 8.

6.   Conclusions

We obtained periodic solutions for two-dimensional non-autonomous differential equations. We discussed two types of coefficients called algebraic and trigonometric (non-homogeneous and homogeneous) coefficients for various classes. Maximum possible upper bound is executed for classes $C_{14, 7}$, $C_{16, 8}$ with non-homogeneous and $C_{5, 5}$ with homogeneous trigonometric coefficients, while $C_{9, 4}$, $C_{9, 5},$ $C_{9, 6}$ and $C_{11, 3}$ with the polynomial coefficients. We attained the highest multiplicity as $10$ for classes $C_{9, 4},$ $C_{11, 3}$ $C_{5, 5},$ and $C_{6, 6}$ which is the highest one. All the result has been done using computer algebra package Maple 18. We have concluded from this extensive analysis that there is no connection between the degree of polynomial and the number of limit cycles. The higher degree classes like $C_{14, 7}$, $C_{16, 8}$ ends up with maximum multiplicity eight in contrast to lower degree polynomial classes $C_{5, 5},$ $C_{6, 6}$ and $C_{9, 4}, C_{11, 3}$ having ten limit cycles. This conclusion is in accordance with the 2nd part of Hilbert's 16th problem. We have calculated the formulae in Section 2; which extends previous works of Alwash ^[12] and Yasmin ^[16] and attempt to step forward in the comprehension of the case n > 9, in the Eq $\left(2.5\right)$. In future, one can generalize the same concept of multiplicity and can calculate multiplicity greater than ten by extending Theorem 2.1. The same results can be validated by the use of simulation against time domain.

Acknowledgments

This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors. All the authors contributed equally from their own sources to fund this manuscript.

Conflict of interest

The authors declare that they have no known financial interest or personal relationship that could have appeared to influence the work reported in this paper.