Jump to ContentJump to Main Navigation
Show Summary Details
More options …

Open Physics

formerly Central European Journal of Physics

Editor-in-Chief: Seidel, Sally

Managing Editor: Lesna-Szreter, Paulina

IMPACT FACTOR 2018: 1.005

CiteScore 2018: 1.01

SCImago Journal Rank (SJR) 2018: 0.237
Source Normalized Impact per Paper (SNIP) 2018: 0.541

ICV 2017: 162.45

Open Access
See all formats and pricing
More options …
Volume 16, Issue 1


Volume 13 (2015)

Complementary wave solutions for the long-short wave resonance model via the extended trial equation method and the generalized Kudryashov method

Rodica Cimpoiasu / Alina Streche Pauna
Published Online: 2018-07-19 | DOI: https://doi.org/10.1515/phys-2018-0057


In this paper the nonlinear long-short (LS) wave resonance model is analyzed through a new perspective. We obtain the classification of exact solutions by making use of the complete discrimination system for the trial equation method and through the generalized Kudryashov method. These methods do generate complementary wave solutions such as bright and dark solitons, rational functions, Jacobi elliptic functions as well as singular and periodic wave solutions. Some among them extend the already reported solutions through other techniques. For some types of solutions adequately graphical representations are displayed. The concerned methods could also be used in order to study other interesting nonlinear evolution processes in n dimensions.

Keywords: Long-short wave resonance model; solitary wave solutions; extended trial equation method; generalized Kudryashov method; symbolic computation

PACS: 02.30.Jr; 04.20.Jb; 05.45.Yv

1 Introduction

It is well known that nonlinear phenomena occur in various areas of science and engineering, such as fluid mechanics, plasma, solid-state physics, biophysics, optical fibers, technology of space, control engineering problems, hydrodynamics etc. They could be modeled through nonlinear partial differential equations (NPDEs). Thus to search for various solutions of the NPDEs does become one of the most exciting and extremely active areas of investigation of these complex physical phenomena. During the recent years, a lot of powerful methods such as the inverse scattering method [1, 2], the Hirota bilinear transformation [3], the generalized Riccati equation method [4], the GG -expansion method [12, 13], the Lie symmetry method [7, 8, 9, 10, 11], the soliton ansatz method [12, 13], the generalized conditional symmetry approach [14, 15] and other techniques were established in order to obtain exact travelling wave solutions of nonlinear physical problems [16, 17, 18, 19].

A general theory able to explain the interactions between short and long waves has been developed (see Ref. [20]). It is known that the long wave-short wave resonance is a special case of the three wave resonances which does appear when the phase velocity of a long wave does match with the group velocity of a short wave and it could also appear when the second order nonlinearity would occur in the process. The physical importance of the so called (LS) type equations [21] consists in the fact that the dispersion of the shortwave is balanced by the nonlinear interaction between the long wave and the short wave, while the evolution of the long wave is driven by the self-interaction of the short wave. The (LS) wave resonance model is described by the following equations:


where L represents the amplitude of the long wave and it is a real function, S is the envelope of the short wave and it is a complex function, while α and β are real constants. Some works such as [22, 23], have been published about the qualitative research of global solutions for the long–short wave resonance equations. In [24] Lax’s formulation for Eqs. (1) have been provided and Cauchy problem has been solved through the inverse scattering method. The orbital stability of the solitary waves associated to (LS) type equations has been proven in [25, 26]. The first integral method [27] has been recently applied [28] to the governing system. The obtained solutions have been expressed in terms of trigonometric, hyperbolic and exponential functions. Multiple exact travelling wave solutions have been also derived in [29] by making use of the extended hyperbolic functions method for nonlinearwave equations. This approach is explicitly related to the solutions of the projective Riccati equations. Another perspective is offered in [30], where the solutions of the (LS) system are assumed as polynomial expansions in two variables which do satisfy to the projective Riccati equations yet their explicit solutions are not sought for.

However, in our contribution, we will classify the types of solutions related to Eq. (1) according to the values of some parameters through two distinct approaches. In Sections 2 and 3 we will provide the fundamental usefulness of the extended trial function method [31] and respectively of the generalized Kudryashov method [32]. We will study in Section 4 how some types of travelling wave solutions are successfully pointed out by making use of the concerned methods. We will outline some solutions such as (a) the solitary wave solution of bell-type, (b) the solitary wave solution of kink-type for S and bell-type for L, (c) the compound solitary wave solution of the bell-type and the kink-type for S and L, (d) the singular travelling wave solutions and (e) the solitary wave solution of Jacobi elliptic function type. Some of them are more general than the ones derived through other algorithms. Section 5 is dedicated to some conclusions and final remarks.

2 Basics on the extended trial equation method

In this section we will describe the extended trial equation method for finding the travelling wave solutions of NPDEs. The main steps of the method are the the following:

  • Step 1

    For a given NPDE:


    the wave transformation is applied:


    Here k, w and ξ0 are constants. The travelling wave variable (3) allows us to reduce Eq. (2) into an ODE for u = u(η):


  • Step 2

    Suppose that the solution of Eq. (4) can be expressed as follows:


    where Γ(η) satisfies the trial equation:


    Here Φ(Γ) and Ψ(Γ) are polynomials.

    The coefficients ai, i = 0,δ, ξk, k = 0,p, νr, r = 0,n are constants to be determined later. A relationship between p, n, δ could be defined by taking into consideration the homogeneous balance between the highest order derivatives and the nonlinear terms that do appear in the ODE (4).

  • Step 3

    When the solution (5) together with (6) is substituted into (4) an equation for the polynomial Ω(Γ) is generated:


    Equating the coefficients of each power Γj, j = 0,s to zero we generate a set of algebraic equations for ai, ξk and νr.

  • Step 4

    Simplify Eq. (6) to the elementary integral shape:


    By making use of a complete discrimination system [33] for polynomials in order to classify the roots of Φ(Γ), we can solve (8) and we can as well derive the travelling wave solutions of underlying Eq. (2).

3 Basics on the generalized Kudryashov method

In this section we will describe a generalized version of the Kudryashov method for finding solitary wave solutions of NPDEs. When taking into consideration the NPDE:


the main steps of the method can be summarized as follows:

  • Step 1

    We are looking for travelling wave solutions of (9) under the form:


    where μk, ρk, k = 1,2 are arbitrary constants. The wave variables (10) allow us to convert (9) into an ODE:


  • Step 2

    We assume that the exact solutions of Eq.(11) can be expressed as follows:


    where ai, i =0,N, bj, j =0,M are arbitrary constants and Q(ρ) is the solution of Kudryashov equation [34]:


    The relationship between the integers M and N can be found by taking into consideration the homogeneous balance between the highest order derivatives and the nonlinear terms which appear in the travelling wave Eq. (11).

  • Step 3

    For each couple of values respectively taken by M and N, by substituting Eqs. (12) and (13) into Eq. (11) yields to a polynomial R(Q). By setting all the coefficients of R(Q) to zero, the parameters ai, bi,μk,ρk can be explicitly determined by solving the equations of their algebraic relations.

4 Travelling wave solutions related to (LS)-type equations

In this section, we will apply the methods which we have summarized in the previous section in order to construct the travelling wave solutions for the system (1). Let us present in the next two subsections some preliminary results.

4.1 Relationship between the real functions in the case of (LS)-type equations

Suppose that the governing system (1) admits travelling wave solutions under the form:


where μ = μ1 x + μ2 t, ρ = ρ1 x + ρ2 t are wave variables and S(ρ), L(ρ) are real-valued functions.

By substituting (14) into the system (1) the following ODEs are generated:


where the prime does indicate the derivative with respect to ρ.

From the first equation of (15) we derive the parameter relation:


When integrating the last equation from (15) with respect to ρ, then taking into consideration (16), we derive the relationship between functions S and L under the form:


where the integration constant is chosen to be zero.

The previous result together with the second equation from system (15) generate the Liénard ODE:


where we use the denotations:


Consequently the search for travelling wave solutions in the case of (LS)-type equations is reduced to a discussion about travelling wave solutions related to the ODE (18). This discussion will be the object of the next subsection.

4.2 Travelling wave solutions of (1) via the extended trial equation method

Let us suppose that the solutions of Eq. (18) are expressed in accordance with (5)-(6). Taking into consideration the homogeneous balance between S and S3, we should require that:


In order to obtain the travelling wave solutions for (1) let us choose n = 0, δ = 1 and p = 4. This choice corresponds to solutions under the form:


where the constants a0 and a1 do satisfy to a02+a120 and have to be determined later. Relying upon the above description of the method, we can obtain:


By substituting (21) and (22) into (18), the left-hand side of the underlying ODE becomes a third order polynomial in Γ. Setting the coefficients of power Γm,m = 0,3 to zero yields an algebraic system for parameters ai, i = 0, 1ξk, k =0,4, μj, j = 1, 2, ρ1, ν0. The solution given by the Maple program is:


In the above mentioned particular case, the integral (8) becomes:


where K=ν0/ξ4. By substituting the solution (23) into (24), we get:


where K=4ρ12μ1a12 with μ1 < 0.

If we denote by αn, n = 1,4 the roots of the polynomial involved in (25), we reach five distinct cases which allow us to find the unknown Γ function. They are:

  • 1)

    When α1 is a one and only multiple root of the polynom, we get


  • 2)–3)

    When the polynom admits two distinct roots αn, n = 1, 2, we obtain:


  • 4)

    When the polynom admits three distinct roots with α1 > α2 > α3, we have:


  • 5)

    When all the roots of the polynomial are distinct with α1 > α2 > α3 > α4, the result is:




Let us choose for simplicity ρ0 = 0, a0 =–a1 α1. After introducing Γ(ρ) from (26)-(30) into (21) and (17), we can classify the solutions of the (LS) system (1) as follows:

  1. solitary wave solutions of bell-type given by



    where K1 = a1K.

  2. singular wave solutions


    or bright soliton solutions


    where K2=2(α1α2)(α1α3)a1α3α2,P=2α1α2α3α3α20,R=(α1α2)(α1α3)K. Here K2 is the amplitude of the soliton, and R is the inverse width of the soliton. Thus, we can assert that solitons exist for a1 < 0. In our situation P cannot be zero, because R does not admit real values. However, the solutions S(ρ) proportional with sech function could be derived in [30] through another technique, namely an improved GG-expansion method.

  3. Jacobi elliptic function type solutions:


    where K3 = 2a1(α1-α3)(α4α2),M = α4α2, N = α1α4,


For Jacobi elliptic wave solution S5(t,x) with appropriate parametric choices, the surface plot of imaginary values is represented in Figure 1. The real values of S5 do present a similar graph except for a phase difference.

The graphical representation of imaginary values of the short wave’ s envelope S5(t,x) from (35) is shown at K3 = –8, μ1 = –1, μ2 = 1,125, M = –2, N = 3, l = –0,866, φ = 0,25(2x + 4t). The real values of S5(t,x) do present a similar graph except for a phase difference.
Figure 1

The graphical representation of imaginary values of the short wave’ s envelope S5(t,x) from (35) is shown at K3 = –8, μ1 = –1, μ2 = 1,125, M = –2, N = 3, l = –0,866, φ = 0,25(2x + 4t). The real values of S5(t,x) do present a similar graph except for a phase difference.

Remark 1

When the modulus l → 1, the solution (35) can be converted into the hyperbolic function solution of kink type for the (LS)-system:


where α3 = α4.

Remark 2

By taking the modulus l → 0 into (35) it yields to the periodic wave solutions of the (LS)-system:


where α2 = α3.

4.3 Travelling wave solutions of (1) via the generalized Kudryashov method

Let us now expand S(ρ) from (18) in accordance with (12) and (13). By balancing the terms S and S3, we require that:


If we choose N = 2, M = 1, the solution of (18) will be sought under the form:


where ai, i = 0, 2, bj, j = 0, 1 are constants to be determined later and Q(ρ)= 11±eρ are the solutions of the auxiliary Eq. (13).

When substituting (39) together with (13) into the (18) and gathering all the terms with the same power of Q, the left-hand side of (18) is converted into a polynomial with Qr, r = 0, 6. When equating each coefficient of this polynomial to zero, a set of 7 algebraic equations for ai, i = 0,2, bi, j = 0, 1, μ1, μ2, ρ1 is derived. By solving this algebraic system with the help of Maple program, we get solitary wave solutions for (LS)-system as follows:

The graphical representations of the real values of the short wave’s envelope S9 and respectively of the long wave’s amplitude L9, are shown at parametric choices a2 = 1, ρ1 = 2, ρ2 = –0,1, μ1=0,025, μ2 = 5,001. The imaginary values of S9 do present similar graph except for a phase difference.
Figure 2

The graphical representations of the real values of the short wave’s envelope S9 and respectively of the long wave’s amplitude L9, are shown at parametric choices a2 = 1, ρ1 = 2, ρ2 = –0,1, μ1=0,025, μ2 = 5,001. The imaginary values of S9 do present similar graph except for a phase difference.

  1. The solitary waves solution of kink-type:


    generated through the choice of Q(ρ)=11eρ and through taking under consideration the parameters:


  2. Dark soliton solution of kink-type for S and bell-type for L expressed under the form:


    derived by assuming that Q(ρ)=11+eρ as well as the relationships among parameters:


    The dynamics of dark soliton solution (42) of kink-type for S9 and of bell-type for L9 with suitable parametric choices is shown in Figure 2. The imaginary values of S9 present a similar graph except for a phase difference.

    By choosing Q(ρ)=11eρ, singular travelling wave solutions expressed in terms of coth function can be found.

  3. Compound wave solution of the bell-type and the kink-type for S and L with the explicit expressions:


    which are related to Q(ρ)=11eρ and to the parameter relationships:


5 Concluding remarks

In this paper we have discussed about the possibilities offered by the extended trial equation method and by the generalized Kudryashov method in order to find new wave solutions for the long–short (LS) wave resonance Eq. (1). At least for the considered cases (21) respectively (39), these two methods lead us to complementary classes of solutions.

The key idea of the trial equation method is to reduce an ODE into a solvable integral Eq. (8) which involves polynomial functions. For the (LS)-type equations we have chosen to take advantage of the solutions to the trial function (forth degree polynomial) by making use of a mathematical tool denominated as the complete discrimination system. We have pointed out solutions such as: solitary wave solution of the bell-type, bright soliton solution, Jacobi elliptic function type solution, singular wave solution, periodic function type solution. To our knowledge, such classification involving the complete discrimination system for polynomials has not been reported. According to the specific parameters, some wave solutions (32), (34)-(37) do extend the ones obtained by the means of other techniques [27, 35, 36].

The generalized Kudryashov method is useful in order to reach the analytical solutions for the (LS) wave resonance model. We have been able to derive wave solution of kink-type (40), dark soliton solution of the kink-type for S and the bell-type for L (42) as well as a compound wave solution of the bell-type and the kink-type for S and L (44).

When we have taken into consideration some other values of parameters which should verify (20) and respectively (38), some new solitary wave solutions could be identified. The proposed methods have their own advantages. They offer both a concise and an efficient way for solving NPDEs, in two or more dimensions. Other methods, such as multi-symplectic method [37, 38], or generalized multi-symplectic method [39, 40] will be employed in future works.


The authors are grateful for the financial support offered by the project “Computational Methods in Astrophysics and Space Sciences”, no. 181/20.07.2017, of the Romanian National Authority for Scientific Research, Program for Research –Space Technology and Advanced Research.


  • [1]

    Ablowitz M.J., Clarkson P.A., Solitons, 1990, Cambridge University Press, Cambridge Google Scholar

  • [2]

    Vakhnenko V.O., Parkes E.J., Morrison A.J., A Bäcklund transformation and the inverse scattering transform method for the generalised Vakhnenko equation, Chaos Soliton Fract., 2003, 17, 683–692. CrossrefGoogle Scholar

  • [3]

    Hirota R., Backlund Transformations, Springer, Berlin, 1980 Google Scholar

  • [4]

    Lu D., Hong B., Tian L., New explicit exact solutions for the generalized coupled Hirota–Satsuma KdV system, Comput. Math. Appl., 2007, 53, 1181–1190. Web of ScienceCrossrefGoogle Scholar

  • [5]

    Bian C., Pang J., Jin L., Ying X., Solving two fifth order strong nonlinear evolution equations by using the G/G-expansion method, Commun. Nonlinear Sci. Numer. Simul., 2010, 15, 2337–2343.Web of ScienceCrossrefGoogle Scholar

  • [6]

    Kudryashov N.A., A note on the G/G-expansion method, Appl. Math. Comput., 2010, 217(4), 1755–1758.Web of ScienceGoogle Scholar

  • [7]

    Olver P.J., Applications of Symmetry Methods to Partial Differential Equations, Springer-Verlag, New-York, 1986 Google Scholar

  • [8]

    Cimpoiasu R., Constantinescu R., Cimpoiasu V.M., Integrability of dynamical systems with polynomial Hamiltonians, Rom. J. Phys., 2005, 50(3-4), 317–324. Google Scholar

  • [9]

    Cimpoiasu R., Constantinescu R., Lie Symmetries for Hamiltonian Systems Methodological Approach, Int. J.Theor. Phys., 2006, 45(9), 1769–1782. CrossrefGoogle Scholar

  • [10]

    Xu M.J., Tian S.F., Tu J.M., Zhang T.T., Lie symmetry analysis, conservation laws, solitary and periodic waves for a coupled Burger equation, Superlattice Microst., 2017, 101, 415-428. Web of ScienceCrossrefGoogle Scholar

  • [11]

    Tu J.M., Tian S.F., Xu M.J., Zhang T.T., On Lie symmetries, optimal systems and explicit solutions to the Kudryashov–Sinelshchikov equation, Appl. Math. Comput., 2016, 275, 345-352. Web of ScienceGoogle Scholar

  • [12]

    Bhrawy A.H., Abdelkawy M.A., Kumar S., Biswas A., Solitons and other solutions to Kadomtsev-Petviashvili equation of B-type, Rom. J. Phys., 2013, 58, 729–748. Google Scholar

  • [13]

    Biswas A., Solitary waves for power-law regularized long-wave equation and R (m, n) equation, Nonlinear Dyn., 2010, 59, 423–426. Web of ScienceCrossrefGoogle Scholar

  • [14]

    Cimpoiasu R., Generalized conditional symmetries and related solutions of the Grad-Shafranov equation, Phys. Plasmas, 2014, 21(4), 042118-1–042118-7. Web of ScienceGoogle Scholar

  • [15]

    Su K., Cao J., Third-Order Conditional Lie–Bäcklund Symmetries of Nonlinear Reaction-Diffusion, Adv. Math. Phys., 2017, 2017, 2825416-1–2825416-9. Google Scholar

  • [16]

    Tian S.F., Zhang H.Q., On the Integrability of a Generalized Variable-Coefficient Forced Korteweg-de Vries Equation in Fluids, Stud. Appl. Math., 2014, 132(3), 212-246. Web of ScienceCrossrefGoogle Scholar

  • [17]

    Wang X.B., Tian S.F., Qin C.Y., Zhang T.T., Characteristics of the breathers, rogue waves and solitary waves in a generalized (2+1)-dimensional Boussinesq equation, Europhys. Lett., 2016, 115, 10002-1-10002-5. Web of ScienceGoogle Scholar

  • [18]

    Wang X.B., Tian S.F., Xu M.J., Zhang T.T., On integrability and quasi-periodic wave solutions to a (3+1)-dimensional generalized KdV-like model equation, Appl. Math. Comput., 2016, 283, 216-233. Web of ScienceGoogle Scholar

  • [19]

    Tian S.F., Initial-boundary value problems for the general coupled nonlinear Schrodinger equation on the interval via the Fokas method, J. Differ. Equations, 2017, 262, 506-558. Web of ScienceCrossrefGoogle Scholar

  • [20]

    Benny D.J., A general theory for interactions between short and long waves, Stud. Appl. Math., 1977, 56(1), 81–94. CrossrefGoogle Scholar

  • [21]

    Djordjevic V.D., Redekopp L.G., On two-dimensional packets of capillary-gravity waves, J. Fluid Mech., 1977, 79(4), 703–714. CrossrefGoogle Scholar

  • [22]

    Xin J., Guo B., Han Y., Huang D., The global solution of the (2+ 1)-dimensional long wave–short wave resonance interaction equation, J. Math. Phys., 2008, 49(7), 073504-1-073504-13. Web of ScienceGoogle Scholar

  • [23]

    Guo B.L., Wang B.X., The global solution and its long time behavior for a class of generalized LS type equations, Prog Nat. Sci., 1996, 6(5), 533–546. Google Scholar

  • [24]

    Ma Y.C., The Complete Solution of the Long-Wave–Short-Wave Resonance Equations, Stud. Appl. Math., 1978, 59(3), 201–221. CrossrefGoogle Scholar

  • [25]

    Laurencot P.H., On a nonlinear Schrödinger equation arising in the theory of water waves, Nonlinear Anal.-Theor., 1995, 24(4), 509–523. CrossrefGoogle Scholar

  • [26]

    Guo B., Chen I., Orbital stability of solitary waves of the long–short resonance equations, Math. Meth. Appl. Sci., 1998, 21, 883–894. CrossrefGoogle Scholar

  • [27]

    Feng Z.S., The first-integral method to study the Burgers–Korteweg–de Vries equation, J. Phys. A-Math. Gen., 2002, 35, 343–349. CrossrefGoogle Scholar

  • [28]

    Jafari H., Soltani R., Khalique C.M., Baleanu D., On the exact solutions of nonlinear long-short wave resonance equations, Rom. Rep. Phys., 2015, 67(3), 762–772. Google Scholar

  • [29]

    Shang Y., The extended hyperbolic function method and exact solutions of the long–short wave resonance equations, Chaos Soliton Fract., 2008, 36(3), 762–772. Web of ScienceCrossrefGoogle Scholar

  • [30]

    Cimpoiasu R., Travelling wave solutions for the long-short wave resonance model by an improved G/G-expansion method, Rom. J. Phys., 2018, (accepted paper), http://www.nipne.ro/rjp/accpaps/016-Cimpoi_3A41FA.pdf 

  • [31]

    Gurefe Y., Misirli E., Sonmezoglu A., Ekici M., Extended trial equation method to generalized nonlinear partial differential equations, App. Math. Comput., 2013, 219(10), 5253–5260. CrossrefGoogle Scholar

  • [32]

    Zayed E.M., Moatimid G.M., Al-Nowehy A.G., The generalized Kudryashov method and its applications for solving nonlinear PDEs in mathematical physics, Sci. J. Math. Res., 2015, 5, 19–39. Google Scholar

  • [33]

    Liu C.S., Applications of complete discrimination system for polynomial for classifications of traveling wave solutions to nonlinear differential equations, Comput. Phys. Commun., 2010, 181(2), 317–324. Web of ScienceCrossrefGoogle Scholar

  • [34]

    Kudryashov N.A., One method forfinding exact solutions of nonlinear differential equations. Commun. Nonlinear Sci. Numer. Simul., 2012, 17, 2248–2253. CrossrefGoogle Scholar

  • [35]

    Shang Y., Huang Y., Yuan W., The extended hyperbolic functions method and new exact solutions to the Zakharov equations, Appl. Math. Comput., 2008, 200(1), 110–122. Web of ScienceGoogle Scholar

  • [36]

    Zuo J.M., Application of the extended G/G-expansion method to solve the Pochhammer–Chree equations, Appl. Math. Comput., 2010, 217(1), 376–383. Web of ScienceGoogle Scholar

  • [37]

    Bridges T.J., Multi-symplectic structures and wave propagation, Math. Proc. Camb. Phil. Soc., 1997, 121, 147-190. CrossrefGoogle Scholar

  • [38]

    Hu W., Deng Z., Zhang Y., Multi-symplectic method for peakon-antipeakon collision of quasi-Degasperis-Procesi equation, Comput. Phys. Commun., 2014, 185, 2020-2028. Web of ScienceCrossrefGoogle Scholar

  • [39]

    Hu W., Deng Z., Yin T., Almost structure-preserving analysis for weakly linear damping nonlinear Schrodinger equation with periodic perturbation, Commun Nonlinear Sci. Numer Simulat., 2017, 42, 298–312. CrossrefGoogle Scholar

  • [40]

    Hu W., Deng Z., Competition between geometric dispersion and viscous dissipation in wave propagation of KdV-Burgers equation, J. Vib. Control, 2015, 21, 2937-2945. Web of ScienceCrossrefGoogle Scholar

About the article

Received: 2018-03-21

Accepted: 2018-04-26

Published Online: 2018-07-19

Citation Information: Open Physics, Volume 16, Issue 1, Pages 419–426, ISSN (Online) 2391-5471, DOI: https://doi.org/10.1515/phys-2018-0057.

Export Citation

© 2018 Rodica Cimpoiasu and Alina Streche Pauna, published by De Gruyter. This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 License. BY-NC-ND 4.0

Comments (0)

Please log in or register to comment.
Log in