RETRACTED
This article has been retracted as it was previously published as “Classification of bounded travelling wave solutions of the generalized Zakharov equation,” by H. R. Z. Zangeneh, R. Kazemi and M. Mosaddeghi, in Iranian Journal of Science & Technology.

Journal of Applied Mathematics

Volume 2014 (2014), Article ID 170946, 11 pages

http://dx.doi.org/10.1155/2014/170946

Research Article

Bifurcation of Travelling Wave Solutions of the Generalized Zakharov Equation

Masoud Mosaddeghi

Département de Mathématiques et de Statistique, Université de Montréal, CP 6128, Succursale Centre-Ville, Montréal, QC, Canada, H3C 3J7

Received 11 April 2014; Accepted 26 July 2014; Published 24 August 2014

Academic Editor: Tiecheng Xia

Copyright © 2014 Masoud Mosaddeghi. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

By using bifurcation theory of planar ordinary differential equations all different bounded travelling wave solutions of the generalized Zakharov equation are classified in to different parametric regions. In each of these parametric regions the exact explicit parametric representation of all solitary, kink (antikink), and periodic wave solutions as well as their numerical simulation and their corresponding phase portraits are obtained.

1. Introduction

Many phenomena in physics, engineering, and science are described by nonlinear partial differential equations (NPDEs). Exact travelling wave solution of nonlinear evolution equation is one of the fundamental objects of study in mathematical physic. When these exact solutions exist, they can help one to understand the mechanism of the complicated physical phenomena and dynamical processes modelled by these nonlinear evolution equations. In the past decades a vast variety of the powerful and direct methods to find the explicit solutions of NPDE have been developed, such as Hirota bilinear method [1, 2], inverse scattering transform method [3], Backlund and Darboux transforms method [4], Lie group method [5] F-expansion method [6], sine-cosine method [7], homotopy perurbation method [8], homogenous balance method [9], algebraic method [10], and Jacobi elliptic function expansion method [11]. Certainly, the bifurcation theory of planar dynamical systems is an efficient method too [12, 13]. In this paper we consider the generalized Zakharov equation by using the bifurcation theory. The Zakharov equations [14] are used extensively in considerations of the evolution of Langmuir turbulence when strong turbulence effects are considered. The Zakharov equation has various applications in physics in a theory of deep-water waves [15], communication [16], and nonlinear pulse propagation in fibers [17]. The original derivation of these equations was based on a simplified model involving fluid concepts. The model leads to two equations: one of these describes the evolution of the envelope of the Langmuir waves with the nonlinearity included through a term involving a density fluctuation, and the other describes the evolution of the density fluctuation due to the ponderomotive force exerted by the Langmuir waves. Now we consider the generalized Zakharov equations which have the forms

In (1) represents the coefficient of dispersion and , , represent the coefficients of nonlinearity. When , (1) reduce to the famous Zakharov equations (see [18, 19] for details) which describe the propagation of Langmuir waves in plasmas. The complex number represents the envelope of the electric field, and is equilibrium value from the fluctuation of the ion density. The parameter is proportional to the ion acoustic speed. We find exact solutions of (1) via the bifurcation theory of planar dynamical systems. The purpose of this paper is to give the bifurcation sets of the bounded travelling wave solutions, that is, solitary wave solutions, kink (anti kink) wave solutions, and periodic wave solutions. Also we obtain the explicit representation for some of these solutions in different parametric region determined by the bifurcation set. To find the travelling wave solutions of (1) we consider the travelling wave solutions of the form

where , , and are arbitrary constants; denotes the wave speed, is the frequency, is the soliton wave number, and is a real function that represents the shape of the pulse. By substituting for and from (2) into (1) we get the following ordinary differential equations:

By integration of the first equation in (3) twice with respect to and taking the integration constant to be zero, we obtain

now we substitute (4) in the second equation of (3) to obtain

where

Now let . Then we drive the following travelling wave system which is a planar Hamiltonian system:

Because the phase portraits of the Hamiltonian system (7) determine travelling wave solutions of (1), we find the bifurcation set for which the qualitative behavior of phase portraits of (7) changes. Here we consider only bounded travelling waves because in physical model only bounded travelling waves are meaningful. Suppose that is a continuous solution of system (7) for and , . We recall that

(i)if , then is called a solitary or impulse wave solution, and(ii)if then is called kink (antikink) wave solution.

Usually a solitary wave solution, a kink (antikink) wave and periodic travelling wave solutions of (1) correspond to a homoclinic orbits or cuspidal loop, heteroclinic orbit or eye-figure loop, and periodic orbit of (7), respectively. Thus we need to find all periodic orbits, homoclinic orbits, and heteroclinic orbits of system (7) which depend on the systems’ parameters.

The rest of this paper is organized as follows. In Section 2, we give the bifurcation set and corresponding phase portrait of system (7). In Section 3, using the information obtained about the phase portraits of bounded solutions of (7), we obtain the numerical simulation for corresponding bounded travelling wave solutions of system (1). In Section 4, we give exact explicit parametric representation of different possible solitary wave solutions, periodic travelling wave solutions, and kink (antikink) wave solutions of (1).

2. Bifurcation of Phase Portrait of (7)

In this section, we consider bifurcation set and phase portraits of (7). First we consider the generic case . With some time scaling and without loss of generality, we can assume , so that we can have Hamiltonian system:

with the Hamiltonian , where correspond to and is the corresponding potential function. It is clear that critical points of are zeros of . It is known that isolated minimum, maximum, and inflection points of correspond to center, saddle point, and cusp point of system (8), respectively (e.g., see [20]). Also it is known that the global structure of phase portraits of system (8) will not change qualitatively unless one of the conditions listed below is violated [20].

(i)There are only finitely many critical points of .(ii)Each critical point of is nondegenerate; that is, for critical point .(iii)No two maximum values of are equal.(iv) as ; that is, is unbounded for both and .

Potential functions satisfying the above four conditions are called the generic potential functions. In our case it is clear that conditions and are satisfied for all values of and . Therefore to find the bifurcation set, we first need to find conditions where critical points of become degenerate. So we set

By solving (9) we find the bifurcation set to be and which correspond to , respectively. These bifurcation sets divide the parametric plane into 7 distinct regions (see Figures 1 and 2). In each parametric region, the number and type of critical points remain unchanged. To see the type and number of critical points it is sufficient to consider only a typical equation for a particular value of and in each region. Critical points of are and where which correspond to , respectively. It is easy to verify that number of critical points of will change from one nondegenerate critical point in region to five nondegenerate critical points in region and three nondegenerate critical points in region . On the boundary of these regions, that is, regions , , , and , critical points are degenerate. To classify the critical points and determine the phase portraits of system (7) we need to consider the two cases and separately.

Figure 1: Bifurcation sets and phase portraits of (8) with in different parametric regions.

Figure 2: Bifurcation sets and phase portraits of (8) with in different parametric regions.

Case I  . Since is a sixth order polynomial and the coefficient of is positive, it is easy to verify that the only nondegenerate critical point of in region or is which is a minimum point. Therefore in this region system (7) has a global center (see Figure 1).

In region , has three critical points where is a minimum and are inflection points. These points correspond to a center and cusp points of (7), respectively. Also because of symmetry, . Therefore system (7) will have an eye-figure loop connecting the cusp points and an oval of periodic orbits encircling the origin and also a band of periodic orbits outside the eye-figure loop (see Figure 1).

In region the potential function always has three nondegenerate minimum points and two nondegenerate maximum points between them. Therefore in this region system (7) has three center and two saddle points. Furthermore, since is even and the critical points are symmetric with respect to -axis, we have

There will be a cycle, heteroclinic to these two saddle points, and since

this heteroclinc cycle includes a band of periodic orbits encircling the origin. Also since and > , there will be two orbits homoclinic to these saddle points. These homoclinc orbits include a band of periodic orbits encircling centers at points of system (7) (see Figure 1).

In regions    and , the only critical point of is which is a degenerate minimum point; therefore system (7) has a global degenerate center at origin (see Figure 1).

In region (), has three nondegenerate critical points at and , where is a maximum and are minimum points. Also . Therefore system (7) has orbits double homoclinc to saddle point at the origin. Also there will be ovals of periodic orbits inside each of the homoclinic orbits encircling centers at and a band of periodic orbits outside the double homoclinic (figure-eight loop) orbits (see Figure 1).

In region (, ), the situation is similar to that of region with the difference that the maximum point of and saddle point of system (7) are degenerate in this case (see Figure 1).

Therefore we have proved the following lemma.

Lemma 1. Phase portrait of system (7) corresponding to is classified as follows (see Figure 1):

(i)In region   (, , ), phase portrait of system (8) consists of two saddle points, three centers, a cycle heteroclinic to saddle points, two homoclinic orbits to saddle points, and bands of periodic orbits inside heteroclinic cycle and homoclinic orbits and outside figure-eight loop.(ii)In region   (, ) phase portrait of (7) consists of a center, two cusp points, an eye-figure loop, and two bands of periodic orbits inside and outside the eye-figure loop.(iii)In regions   ( or , , ) phase portrait consists of a global nondegenerate center.(iv)In regions () and (, ), phase portrait consists of a degenerate global center.(v)In region (), phase portrait consists of a nondegenerate saddle point at the origin, two centers at , a double homoclinic orbit to the saddle point, and bands of periodic orbits inside and outside the double homoclinic orbit.(vi)In region   , phase portrait consists of a degenerate saddle point at the origin, two centers, a double homoclinic orbit to the saddle point, and bands of periodic orbits inside and outside the double homoclinic orbit.

Case II  . In this case parametric region is divided into eight locally topologically equivalent regions by bifurcation set (see Figure 2). Using a similar analysis to Case I, we can derive the phase portraits of system (7). In region (, , ) potential function has three nondegenerate maximum points at , and two minimum points at . Also because of symmetry, . Therefore saddle points always lie on the same level curve of the Hamiltonian . Furthermore along parabola , , all three saddle points lie on the same potential level; that is, ; therefore region will be divided into three subregions which we will denote by , , and . In region bounded orbits of phase portraits of system (7) consist of two cycles heteroclinic to saddle points at origin and and bands of periodic orbit inside the heteroclinic cycles. In region , ; therefore the phase portrait of bounded orbits of (7) consists of a double homoclinic orbit to origin, a cycle heteroclinic to saddle points , bands of periodic orbits inside homoclinic orbits, and a band of periodic orbits outside double homoclinic orbit and inside the heteroclinic cycles. In region , . Therefore the phase portraits of bounded periodic orbits consists of homoclinic orbits to saddle points and band of periodic orbits inside each of the homoclinic orbits (see Figure 2). Phase portraits of system (7) in other regions were derived in a similar manner, omitted here for the sake of brevity. Therefore we proved the following.

Lemma 2. Phase portrait of system (7) corresponding to is classified as follows (see Figure 2).

(i)In region () phase portrait of bounded orbits consists of a center at the origin, two saddle points at , cycle heteroclinic to the saddle points, and a band of periodic orbits encircling the origin inside the heteroclinic cycle.(ii)In region (, ), phase portrait of bounded orbits consists of a degenerate center at the origin, two saddle points at , cycle heteroclinic to the saddle points, and a band of periodic orbits encircling the origin inside the heteroclinic cycle.(iii)In region (, , ), it has three saddle points at , , and two centers at . In regions there are a double homoclinic orbit to origin, a cycle heteroclinic to saddle points , bands of periodic orbits inside homoclinic orbits, and a band of periodic orbits outside double homoclinic orbit and inside the heteroclinic cycles. In region there are two cycles heteroclinic to saddle points at origin and and bands of periodic orbit inside the heteroclinic cycles. In region there are homoclinic orbits to saddle points and band of periodic orbits inside each of the homoclinic orbits.(iv)In region (, , ), there are a saddle point at and two cusp points at .(v)In region ( or , , ) there is only a saddle point at .(vi)In regions (, ) and (), there is a nonhyperbolic saddle point at the origin .

Case III  . Now we consider the degenerate case . If , (7) becomes

Then without loss of generality we can assume , so that we can have the Hamiltonian system

with Hamiltonian , where correspond to and is the corresponding potential function. Similar to the above discussions of the Cases I and II, we see that bifurcation occurs in . If , then system (5) becomes a linear differential equation. Then we have the following lemma.

Lemma 3. Phase portrait of system (12) corresponding to is classified as follows.

(i)If , then the phase portrait of system (12) consists of a global center at the origin for , a degenerate global center at the origin for and a nondegenerate saddle point at the origin, two centers at , a double homoclinic orbit to the saddle point, and bands of periodic orbits inside and outside the double homoclinic orbit for .(ii)If , then the phase portrait of system (12) consists of two nondegenerate saddle points at , cycle heteroclinic to the saddle points, and a band of periodic orbits encircling the origin inside the heteroclinic cycle for , a degenerate saddle at the origin for , and a saddle at the origin for .(iii)If , then for we have a global center at the origin and for we have a nondegenerate saddle point at the origin.

3. The Numerical Simulation of Bounded Travelling Waves

It is well known that the bounded travelling waves of system (1) correspond to the bounded integral curves of system (7). In Lemmas 1, 2, and 3, we have classified all bounded integral curves of system (7). In this section we give numerical simulation for a typical member of bounded travelling waves of system (7) in the form of as follows.

Case I  (homoclinic loops). These orbits only exist in regions , when and in regions , when and for , . Homoclinic orbits of system (7) correspond to solitary travelling waves of (5). Let , , and which correspond to a point in region of Figure 1. Now we consider system (7) and choose initial conditions , . So, they lie on the left branch of double homoclinic orbit (figure-eight loop). In physics this type of travelling waves is called solitary wave with valley form (see Figure 3(a)). Now let , , and which correspond to a point in parametric region in Figure 1. Again we use initial conditions to be on the homoclinic orbit of system (7). Let , , so that they lie on the right branch of double homoclinic orbit (figure-eight loop). This type of travelling wave in physics is called solitary wave with peak form (see Figure 3(b)).

(a)

(b)

Figure 3: The simulation of solitary waves corresponding to the homoclinic orbits of (7). (a) Solitary wave of valley form. (b) Solitary wave of peak type.

Case II  (nilpotent homoclinic loops). These orbits only exist in region when . Nilpotent Homoclinic loop of system (7) corresponds to solitary wave of system (5). As in the previous part we choose two sets of parameters , , and which correspond to a point in region in Figure 1. Now we consider system (7) and choose initial conditions , , so that they lie on the left branch of degenerate double homoclinic (figure-eight loop) orbits which correspond to solitary wave with valley form (see Figure 4(a)). For the same parameters and initial conditions , , so that they lie on the right branch of degenerate double homoclinic (figure eight loop) orbits which correspond to solitary wave with peak form (see Figure 4(b)).

(a)

(b)

Figure 4: The simulation of the solitary waves corresponding to nilpotent homoclinic orbits of (7). (a) Solitary wave of valley form. (b) Solitary wave of peak form.

Case III  (heteroclinic orbits). These orbits exist only in region when and in regions , , and when and for , . Upper and lower heteroclinic orbits of system (7) correspond to kink and antikink travelling waves of system (5), respectively. Again we consider system (7) and choose , , and which correspond to a point in parametric region in Figure 1. As we have mentioned above in this case we have orbits heteroclinic to saddle points and . Now we use initial conditions , , and , , on upper and lower heteroclinic orbits, respectively, and we get Figures 5(a) and 5(b).

(a)

(b)

Figure 5: The simulation of the kink and antikink waves corresponding to the heteroclinic orbits of (7). (a)Antikink waves. (b) Kink waves.

Case IV  (eye-figure loop). This loop exists only in region when . Upper and lower orbits of eye-figure loop of system (7) which again correspond to kink and antikink travelling waves of system (5), respectively. As in the previous part, we choose a set of parameters , , and which correspond to a point in region in Figure 1. Now we consider system (7) and choose initial conditions , , and , , on upper and lower orbits of eye-figure loop, respectively, and we get Figures 6(a) and 6(b).

(a)

(b)

Figure 6: The simulation of the kink and antikink waves corresponding to the eye-figure loop of (7). (a) Antikink waves. (b) Kink waves.

Case V  (periodic orbits). These periodic orbits are global (regions , , ) when and for , , or local which lie inside homoclinic orbits (regions , ) when and (regions , ) when , inside and outside of eye-figure loop (regions ) when , inside heteroclinic cycles (region ) when and (region , , ) when and for , , , inside and outside of double figure-eight loop (regions ) when , inside and outside of nilpotent double homoclinic loop, (regions ) when and inside and outside of figure eight loop (regions ) when and for , , . Periodic orbits of system (7) correspond to periodic travelling waves of system (5). Here we choose a periodic orbit inside the heteroclinc orbits in region in Figure 1. Of course we could choose a global center or period orbits inside homoclinic orbits and period orbits inside and outside of eye-figure loop as well, but the figures are qualitatively the same. Let , , and . Heteroclinic orbit corresponding to these sets of parameters passes through saddle points and and includes the center of system (7). We choose three sets of initial conditions , , , , and , , close to center , somewhere in middle and very close to heteroclinc orbit, respectively (see Figures 7(a), 7(b), and 7(c)). We notice that period of these periodic orbits increases as we move away from the center toward the heteroclinic orbits. When we choose two sets of initial conditions , , and , , close and far from the figure-eight loop in region , respectively (see Figures 8(a) and 8(b)). Finally we choose two sets of initial conditions , , and , , close and far from the double figure-eight loop, respectively (see Figures 9(a) and 9(b)).

(a)

(b)

(c)

Figure 7: Simulation of periodic waves corresponding to periodic orbits inside heteroclinc cycle of (7). (a) Medium period. (b) Short period. (c) Long period.

(a)

(b)

Figure 8: Simulation of periodic waves corresponding to periodic orbits outside figure-eight loop of (7). (a) Close to figure-eight loop. (b) Far from figure-eight loop.

(a)

(b)

Figure 9: Simulation of periodic waves corresponding to periodic orbits outside double figure-eight loop of (7). (a) Close to double figure-eight loop. (b) Far from double figure-eight loop.

Remark 4. We notice the difference between solitary waves in Figures 3 and 4 and kink and antikink in Figures 5 and 6 which shows their asymptotic behavior as . In Figure 3 stable and unstable manifolds of equilibrium point intersect transversally but in Figure 4 they intersect tangentially. A similar condition holds in Figures 5 and 6, respectively.

We notice that the number of inflection points on periodic waves during one-half period in each of Figures 7, 8, and 9 is zero, one, and two which correspond to the shape periodic orbits of system (7).

4. Explicit Formulas for Bounded Integral Curves of (7)

In this section we give explicit formulas for bounded travelling waves of system (7). In the first step we consider the bounded graphics of system (7).

Case 1  

Region . There are two homoclinic orbits and two heteroclinic orbits of system (7) defined by connecting saddle points which are passing through points . Therefore

where and . On the right and down branch of the homoclinic orbit we have . Since , by integration along the right homoclinic orbit for negative , we get:

Then after some algebraic calculations we get solitary wave solutions:

Along the above heteroclinic orbit we have . As above, we can obtain kink and antikink wave solutions:

Region . There is eye-figure loop of system (7), defined by , connecting the cuspidal points . Therefore . Along the above heteroclinic orbit we have , which have kink and antikink solutions:

Region . There are two orbits homoclinic to origin for system (7) defined by that are passing through points . Therefore , where and . The solitary wave solutions are

Region . There are two nilpotent homoclinic orbits of system (7), defined by connecting to degenerate saddle point and passing through points . Then we have where . The solitary wave solutions are

Case 2  

Region . There are two heteroclinic orbits of system (7), defined by connecting the two saddle points . Therefore , where , which have kink and antikink solutions:

Region . There are two heteroclinic orbits of system (7), defined by connecting with tow saddle points . Therefore . The kink and antikink solutions are given by

Region . There are two homoclinic orbits of system (7) defined by connecting with saddle point and passing through points . Therefore

where and , of which we have two solitary wave solutions:

In this region there are also two heteroclinic orbits of system (7), defined by connecting the saddle points . Therefore , where . Kink and antikink solutions are given by

Region . There are four heteroclinic orbits of system (8) defined by . Two of these orbits connect and and the other two connect and . Therefore , where . Kink and antikink solutions are of the form

where .

Region . There are two homoclinic orbits of system (7) defined by connecting with saddle points and passing through points . Therefore , where and . There are solitary wave solutions which are given by

Case 3  

Region . There are two orbits homoclinic to origin for system (7) defined by which are passing through points . Therefore . Two solitary wave solutions are given by

Region . There are two heteroclinic orbits of system (7), defined by connecting saddle points defined by . Kink and antikink solutions are

Case 4 (periodic orbits). Now we calculate explicit formulas for bounded periodic travelling waves. Here we only consider periodic orbits of (7) which are located inside the right homoclinic loops of Figure 2 () in regions III(a) and III(c) and inside the right heteroclinic loop of Figure 2 () in region III(b). Suppose that the periodic orbits pass through and so that . Therefore this periodic orbit lies on level curve where is the Hamiltonian function. Define . is a sixth order polynomial with respect to , where , , and are their roots with . Therefore

Since , we get

With change of variable of the form on the right hand side above, we can derive

Now by using integral tables for elliptic integrals [21] we have

where is inverse Jacobian elliptic function with modulus (see [21]), , and

Then after some algebraic calculations we get solitary wave solution:

Formulas for other periodic orbits of system (7) can be derived in a similar manner which we omitted here for the sake of brevity.

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

References

R. Hirota, “Exact solution of the korteweg—de vries equation for multiple Collisions of solitons,” Physical Review Letters, vol. 27, no. 18, pp. 1192–1194, 1971.

R. Hirota, The Direct Method in Soliton Theory, vol. 155 of Cambridge Tracts in Mathematics, Cambridge University Press, Cambridge, UK, 2004.

M. J. Ablowitz and P. A. Clarkson, Soliton, Nonlinear Evolution Equations and Inverse Scattering, vol. 149 of London Mathematical Society Lecture Note Series, Cambridge University Press, Cambridge, Mass, USA, 1991.

C. Rogers and W. K. Schief, Bäcklund and Darboux Transformations, Geometry and Modern Applications in Soliton Theory, Cambridge Texts in Applied Mathematics, Cambridge University Press, Cambridge, UK, 2002.

P. J. Olver, Applications of Lie Groups to Differential Equations, Springer, New York, NY, USA, 1993.

E. G. Fan, Integrable Systems and Computer Algebra, Science Press, 2004.

C. Yan, “A simple transformation for nonlinear waves,” Physics Letters A, vol. 224, no. 1-2, pp. 77–84, 1996.

J. He, “Application of homotopy perturbation method to nonlinear wave equations,” Chaos, Solitons and Fractals, vol. 26, no. 3, pp. 695–700, 2005.

M. L. Wang, “Exact solutions for a compound KdV-Burgers equation,” Physics Letters A: General, Atomic and Solid State Physics, vol. 213, no. 5-6, pp. 279–287, 1996.

H. Junqi, “An algebraic method exactly solving two high-dimensional nonlinear evolution equations,” Chaos, Solitons and Fractals, vol. 23, no. 2, pp. 391–398, 2005.

S. Liu, Z. Fu, and Q. Zhao, “Jacobi elliptic function expansion method and periodic wave solutions of nonlinear wave equations,” Physics Letters A, vol. 285, no. 1-2, pp. 69–74, 2001.

J. Li and L. Zhang, “Bifurcations of traveling wave solutions in generalized Pochhammer-Chree equation,” Chaos, Solitons and Fractals, vol. 14, no. 4, pp. 581–593, 2002.

J. Li and T. L. He, “Exact traveling wave solutions and bifurcations in a nonlinear elastic rod equation,” Acta Mathematicae Applicatae Sinica, vol. 26, no. 2, pp. 283–306, 2010.

V. E. Zakharov and A. B. Shabat, “Exact theory of two-dimensional self-focusing and one dimensionsl self-modulation of waves in nonlinear media,” Journal of Experimental and Theoretical Physics, vol. 34, no. 1, pp. 62–69, 1972.

E. Kit and L. Shemer, “Spatial versions of the Zakharov and Dysthe evolution equations for deep-water gravity waves,” Journal of Fluid Mechanics, vol. 450, pp. 201–205, 2002.

M. V. Goldman, “Strong turbulence of plasma waves,” Reviews of Modern Physics, vol. 56, no. 4, pp. 709–735, 1984.

D. Anderson, “Variational approach to nonlinear pulse propagation in optical fibers,” Physical Review A, vol. 27, no. 6, pp. 3135–3145, 1983.

V. E. Zakharov, “Collapse of Langmuir waves,” Soviet Physics JETP, vol. 35, no. 5, pp. 908–914, 1972.

V. E. Zakharov and V. S. Synakh, “The nature of the self-focusing singularity,” Soviet Physics-JETP, vol. 41, pp. 465–468, 1976.

J. K. Hale and H. Koçak, Dynamics and Bifurcation, Springer, New York, NY, USA, 1991.

P. F. Byrd and M. D. Friedman, Handbook of Elliptic Integrals for Engineers and Scientists, Springer, Berlin, Germany, 1971.