Research Article | Open Access

Khaled A. Gepreel, Taher A. Nofal, Ali A. Al-Thobaiti, "The Modified Rational Jacobi Elliptic Functions Method for Nonlinear Differential Difference Equations", *Journal of Applied Mathematics*, vol. 2012, Article ID 427479, 30 pages, 2012. https://doi.org/10.1155/2012/427479

# The Modified Rational Jacobi Elliptic Functions Method for Nonlinear Differential Difference Equations

**Academic Editor:**Jianke Yang

#### Abstract

We modified the rational Jacobi elliptic functions method to construct some new exact solutions for nonlinear differential difference equations in mathematical physics via the lattice equation, the discrete nonlinear Schrodinger equation with a saturable nonlinearity, the discrete nonlinear Klein-Gordon equation, and the quintic discrete nonlinear Schrodinger equation. Some new types of the Jacobi elliptic solutions are obtained for some nonlinear differential difference equations in mathematical physics. The proposed method is more effective and powerful to obtain the exact solutions for nonlinear differential difference equations.

#### Introduction

In recent years, the study of difference equations has acquired a new significance, due in large part to their use in the formulation and analysis of discrete-time systems, the numerical integration of differential equation by finite-difference schemes and the study of deterministic chaos.

Since the work of Fermi et al. in the 1960s [1], DDEs have been the focus of many nonlinear studies. On the other hand, a considerable number of well-known analytic methods are successfully extended to nonlinear DDEs by researchers [2â€“17]. However, no method obeys the strength and the flexibility for finding all solutions to all types of nonlinear DDEs. Zhang et al. [18] and Aslan [19] used the (-expansion function method to some physically important nonlinear DDEs. Qiong and Bin [12] constructed the Jacobi elliptic solutions for nonlinear DDEs. Recently, Zhang [20] and Gepreel [21] and Gepreel and Shehata [22] have used the Jacobi elliptic function method for constructing new and more general Jacobi elliptic function solutions of the nonlinear differential difference equations. The main objective of this paper is to modify the rational Jacobi elliptic functions method to obtain the exact wave solutions for nonlinear DDEs. We use this method to calculate the exact wave solutions for some nonlinear DDEs in mathematical physics via the lattice equation, the discrete nonlinear Schrodinger equation with a saturable nonlinearity, the discrete nonlinear Klein-Gordon equation, and the quintic discrete nonlinear Schrodinger equation.

#### 2. Description of the Modified Rational Jacobi Elliptic Functions Method

In this section, we would like to outline an algorithm for using the modified rational Jacobi elliptic functions method to solve nonlinear DDEs. For a given nonlinear DDEs,
where , and are integers, , denotes the set of all *r*th order derivatives of , with respect to .

The main steps of the algorithm for the modified rational Jacobi elliptic functions method to nonlinear DDEs are outlined as follows.

*Step 1. *We seek the traveling wave solutions of the following form:
where
where , and the phase are constants to be determined later. The transformations (2.2) are reduced (2.1) to the following ordinary differential difference equations:
where .

*Step 2. *We suppose the modified rational Jacobi elliptic functions solutions of (2.4) in the following form:
where , , , are constants to be determined later while satisfies a discrete Jacobi elliptic differential equation:
and , , are arbitrary constants.

*Step 3. *Since the general solution of the proposed equation (2.6) is difficult to obtain and so the iteration relations corresponding to the general exact solutions we discuss the solutions of the proposed discrete Jacobi elliptic differential equation (2.6) at some special cases to , , and , to cover all the Jacobi elliptic functions as follows.

*Type 1. * If , , . In this case, (2.6) has the solution , where is the Jacobi elliptic sine function, and is the modulus. In this case from using the properties of Jacobi elliptic functions (see [22]), the series expansion solutions (2.5) take the following form:
Further using the properties of Jacobi elliptic functions, the iterative relations can be written in the following form:
where
, is the *j*th component of shift vector .

*Type 2. *If , , and . In this case, (2.6) has the solution . From using the properties of Jacobi elliptic functions, the series expansion solutions (2.5) take the following form:

*Type 3. *If , , and . In this case, (2.6) has the solution . From using the properties of Jacobi elliptic functions, the series expansion solutions (2.5) take the following form:
From the properties of the Jacobi elliptic functions, we can deduce the iterative relation to the above kind of solutions from Types 2 and 3 as we show in Type 1.

*Step 4. *Determining the degree of (2.5) by balancing the nonlinear term(s) and the highest order derivatives of in (2.4). It should be noted that the leading terms will not affect the balance because we are interested in balancing the terms of .

*Step 5. *Substituting in each type form (2.1)â€“(2.4) and the given values of into (2.4). Cleaning the denominator and collecting all terms with the same degree of , , together, the left-hand side of (2.4) is converted into a polynomial in , , . Setting each coefficient of this polynomial to zero, we derive a set of algebraic equations for , , , , , and .

*Step 6. *Solving the over determined system of nonlinear algebraic equations by using Maple or Mathematica. We end up with explicit expressions for , , , , , and .

*Step 7. *Substituting , , , , , and into , in the corresponding type from (2.1)â€“(2.4), we can finally obtain exact solutions for (2.1).

#### 3. Applications

In this section, we use the proposed method to construct the rational Jacobi elliptic wave solutions for the nonlinear DDEs via the lattice equation, the discrete nonlinear Schrodinger equation with a saturable nonlinearity, the discrete nonlinear Klein-Gordon equation, and the quintic discrete nonlinear Schrodinger equation, which are very important in mathematical physics and have been paid attention to by many researchers.

##### 3.1. Example 1. The Lattice Equation

In this section, we study the lattice equation which takes the following form [23â€“26]: where , , and are nonzero constants. This equation contains hybrid lattice equation, mKdV lattice equation, modified Volterra lattice equation, and Langmuir chain equation for some special values of , , . According to the above steps, to seek traveling wave solutions of (3.1), we construct the transformation where , , and are constants. The transformation (3.2) permits us converting equation (3.1) into the following form: where . Considering the homogeneous balance between the highest order derivative and the nonlinear term in (3.3), we get . Thus, the solution of (3.3) has the following form: where , , and are constants to be determined later and satisfies a discrete Jacobi elliptic ordinary differential equation (2.6). When we discuss the solutions of (2.6), we get the following types.

*Type 1. *If , , and . In this case, the series expansion solution of (3.3) has the form:
With the help of Maple, we substitute (3.5) and (2.8) into (3.3). Cleaning the denominator and collecting all terms with the same degree of , , together, the left-hand side of (3.3) is converted into polynomial in , , . Setting each coefficient of this polynomial to be zero, we derive a set of algebraic equations for , , , , , , , and . Solving the set of algebraic equations by using Maple or Mathematica software package, we have the following.

*Family 1.*
where , , , , and are arbitrary constants.

From (3.5) and (3.6), the solution of (3.3) takes the following form:
where . Figure 1 illustrates the behavior of the exact solution (3.7).

**(a)**

**(b)**

*Family 2.*
where , , , , and are arbitrary constants.

From (3.5) and (3.8), the solution of (3.3) has the following form:
where .

*Type 2. *If , , and . In this case, the solution of (3.3) has the form:
With the help of Maple, we substitute (3.10) into (3.3). Cleaning the denominator and collecting all terms with the same degree of , , together, the left-hand side of (3.3) is converted into polynomial in , , . Setting each coefficient of this polynomial to zero, we derive a set of algebraic equations in , , , , , , , and . Solving the set of algebraic equations by using Maple or Mathematica software package, we get the following.

*Family 1.*
where , , , , and are arbitrary constants.

In this case, the solution of (3.3) takes the following form:
where .

The Jacobi elliptic functions could be generated into the hyperbolic functions when tends to one in the other hand, they are generated into trigonometrical functions when tends to zero.

When , the trigonometrical solution (3.12) takes the following form:
Also if , the hyperbolic solution (3.12) takes the following form:

*Family 2.*
where , ,, , and are arbitrary constants.

In this case, the solution of (3.3) takes the following form:
where .

When , the trigonometrical solution (3.16) takes the following form:
When , the hyperbolic solution (3.16) takes the following form:

*Type 3. *If , , and . In this case, the series expansion solution of (3.3) has the form:
Consequently, using the Maple or Mathematica we get the following results.

*Family 1.*
where , ,, , and are arbitrary constants.

In this case, the solution of (3.3) takes the following form:
where .

*Family 2.*
where , ,, , and are arbitrary constants.

In this case, the solution of (3.3) takes the following form:
where .

##### 3.2. Example 2. The Discrete Nonlinear Schrodinger Equation

The discrete nonlinear Schrodinger equation (DNSE) is one of the most fundamental nonlinear lattice models [8]. Its arise in nonlinear optics as a model of infinite wave guide arrays [27] and has been recently implemented to describe Bose-Einstein condensates in optical lattices. The class of DNSE model with saturable nonlinearity is also of particular interest in their own right, due to a feature first unveiled in [28]. In this section, we study the DNSE with a saturable nonlinearity [29, 30] form: which describes optical pulse propagations in various doped fibers, is a complex valued wave function at sites while and are real parameters. We make the transformation: where , , , and are arbitrary constants. The transformation (3.25) permits us converting equation (3.24) into the following nonlinear difference equation: We assume that (3.26) has a solution of the form: where , , and are constants to be determined later and satisfies a discrete Jacobi elliptic differential equation (2.6). When we discuss the solutions of (3.26), we have the following types.

*Type 1. *If , , and . In this case, the series expansion solution of (3.26) has the form:
With the help of Maple, we substitute (3.28) and (2.8) into (3.26), cleaning the denominator and collecting all terms with the same order of , , together, the left-hand side of (3.26) is converted into polynomial in , , . Setting each coefficient of this polynomial to zero, we derive a set of algebraic equations for , , , , , , and . Solving the set of algebraic equations by using Maple or Mathematica software package, we obtain the following.

*Family 1.*

In this case, the solution of (3.24) takes the following form: where denotes and . Figure 2 illustrates the behavior exact solution (3.30).

**(a)**

**(b)**

**(c)**

**(d)**

*Family 2.*

In this case, the solution of (3.24) takes the following form: where denotes and .

*Type 2. * If , , and . In this case, the solution of (3.26) has the form:

Consequently, using Maple or Mathematica we get the following results.

*Family 1.*

In this case, the solution of (3.24) takes the following form:
where .

In the special case, when , the trigonometrical solution (3.35) takes the following form:

Also if , the hyperbolic solution (3.35) takes the following form:

*Family 2.*

In this case, the solution of (3.24) takes the following form:
where denotes and .

If , the trigonometrical solution (3.39) takes the following form:

Also if , the hyperbolic solution (3.39) takes the following form:

*Type 3. *if , , and . In this case, the series expansion solution of (3.26) has the form:

Consequently, using Maple or Mathematica we get the following results.

*Family 1.*

In this case, the solution of (3.24) takes the following form: