In various phenomena of physics, mathematics, and engineering the nonlinear partial differential equations are important to study them. In the fluid dynamics, the analytical solutions of the nonlinear evolution equations of shallow water waves and the commonly studied equations (the Korteweg-de Vries (KdV) equation , modified KdV equation, Boussinesq equation  , Green-Naghdi equation, Gardeners equation, and Whitham-Broer-Kaup and Jaulent-Miodek (JM) equations) are usually not available. Therefore, finding a new methods and techniques for dealing these types of equations is still an open area. The theory of solitons has contributed to understanding many experiments and complex phenomena in mathematical physics. Thus, it is of interest to evaluate new solitary wave solutions of these equations. So that, during the past five decades, a lot of method was discovered by a diverse group of scientists to solve the nonlinear partial differential equations. For example, the modified simple equation-method
 , the -expansion method    , the modified -expansion method , the extended Jacobi elliptic function method   , the -expansion method   , the modified extended -expansion method , the Riccati-Bernoulli Sub-ODE method
 , the exp-function method , the extended tanh-function method  , the Modified extended tanh-function method  , Quadratic B-Spline Galer-kin Finite Element Method  and so on. The purpose of this paper is to find the exact solution of Jaulent-Miodek (JM)-equation    (contain some parameters) using the two different methods mention above. If these parameters take definite values the solaitry wave solutions can be derived from it. The objective of this article is to investigate more applications than obtained in   to justify and demonstrate the advantages of extended Jacobian elliptic function expansion method and the Riccati-Bernoulli Sub-ODE method. Here, we apply these two methods to Jaulent-Miodek (JM)-equation; firstly we apply the extended Jacobian elliptic function expansion method.
2. Extended Jacobian Elliptic Function Expansion Method   
The solution according to the extended Jacobian elliptic functions takes the form:
where and are respectively the Jacobian elliptic sine function, the Jacobian elliptic cosine function and the Jacobian elliptic function of the third kind and the other Jacobi functions which is denoted by Glaisher’s symbols and are generated by these three kinds of functions, namely
That has the relations
with modulus m where 0 < m < 1.
In addition we know that
The derivatives of other Jacobi elliptic functions are obtained by using Equation (1.5).
According to the balance rules, we can balance the highest order derivative term and nonlinear term so that N in Equation (1.1) can be determined. In addition we see that when m ⇒ 1, , , and degenerate as tanh , and respectively, while when m = 0, it will be , and , therefore Equation (1.1) degenerate as the following forms
Therefore the extended Jacobian elliptic function expansion method is more general than sine-cosine method, the tan function method and Jacobi elliptic function expansion method.
Here, we will apply extended Jacobian elliptic function expansion method described in Section 2 to find the exact traveling wave solutions and then the solitary wave solutions for the Jaulent-Miodek (JM)-equation.
Consider the Jaulent-Miodek (JM)-equation    
Let us preceding the transformation , substitute at Equation (1.7) we get
Let , then Equation (1.8) became
Integrating (1.9) once, we get
Balancing the nonlinear term with the highest order derivative term, we find m + 2 = 3m Þ m = 1.
Consequently, according to the constructed method by taking the first pair of the elliptic function, the solution is
Substitute about at Equation (1.10) arid equating all the coefficients of to zero, we get this system of algebraic equations
Solving this system of algebraic equations by Maple, we get the following results:
So that the exact solution of Equation (1.10)
Now, if m → 1 we can obtain the hyperbolic solution (Figure 1):
Repeating, these previous work for the remaining pairs of the Jacobian elliptic functions, we can obtained the other six different new solutions.
Figure 1. The plot of Equation (1.16) in three dimensions with values: .
3. The Riccati-Bernoulli Sub-ODE Method  
According to the Riccati-Bernoulli Sub-ODE method the suggested solution is
where and m are constants to be determined later. It is important to knots that when and , Equation (2.1) is a Riccati equation. When and , Equation (2.1) is a Bernoulli equation.
Differentiate (2.1) once we get
Substituting the derivatives of u into Equation (2.1) yields an algebraic equation of u, by consider the symmetry of the right-hand item of Equation (2.1) and setting equivalence for the highest power exponents of u we can determine m Comparing the coefficients of yields a set of algebraic equations for a, b, c and which solving to get .
According to the obtained values of these constants and use the transformation the Riccati-Bernoulli Sub-ODE equation admits the following solutions:
1) When , the solution of Equation (2.1) is
2) When and , the solution of Equation (2.1) is
3) When and , the solution of Equation (2.1) is
4) When and , the solution of Equation (2.1) is
5) When and , the solution of Equation (15) is
6) When and the solution of Equation (15) is
where is an arbitrary constant.
Now, we apply this method for solving the Jaulent-Miodek (JM)-equation (1.10) mention above,
According to the Riccati-Bernolli Sub-ODE method,
Differentiate once we get,
Substitute about at Equation (2.12) and by a suitable choose of m and equating the coefficients of different power of w to zero, we get this system of algebraic equations,
Solving this system by Maple, we get .
According to this solution and the constructed method we take only the two cases (3) and (5).
Case (3): when m ¹ 1, b ¹ 0 and c = 0:
1) When a > 0, the solution is (Figure 2):
2) When a < 0, the solution is (Figure 3):
Figure 2. The plot of Equation (2.15) in three dimensions with values: .
Figure 3. The plot of Equation (2.16) in three dimensions with values: .
Case (5): when m ¹ 1, a ¹ 0 and .
1) When the solution is (Figure 4):
Or (Figure 5):
Figure 4. The plot of Equation (2.17) in three dimensions with values: .
Figure 5. The plot of Equation (2.18) in three dimensions with values: .
The extended Jacobian elliptic function expansion method (which depends on the balance rule) has been used successfully to find the exact traveling wave solutions of the Jaulent-Miodek (JM) model. Also, the non-balanced Riccati-Bernoulli Sub-ODE method which is a powerful technique does not depend on the balance rule has been applied to get the solitary wave solution of this model. The comparison between our obtained results in this article with that obtained previously is included. It can be concluded that these two methods are reliable, effective, reduce the volume of calculations and can be applied to many other nonlinear evolution equations.
5. Captions, Citations and Physical Meaning of the Figures
We can briefly give the captions and citations for the obtained figures instead of the corresponding physical meaning. Let us start from Figure 1 which represents the plot of the solitary solution Equation (1.16) in three dimensions according to the values of the obtained parameters realized for the first pairs of the extended Jacobi elliptic function which is hyperbolic one. Furthermore, six new different solitary solutions are realized using the other remaining pairs of the Jacobian elliptic functions. Figure 2 and Figure 3 denotes to the plot of the solitary solutions Equation (2.15), Equation (2.16) respectively in three dimensions according to the values of the obtained parameters realized by the Riccati-Bernolli Sub-ODE method in the interval (x = −4:4, t = −4:4 and y = −4:4) which are periodic singular dark soliton, while Figure 4 denotes to the plot of the solitary solution Equation (2.17) in the intervale (x = −2:2, t = −2:2 and y = −2:2) which is kink shape. Furthemore, Figure 5 denotes to the plot of the solitary solution Equation (2.18) in the intervale (x = −2:2, t = −2:2 and y = −2:2) which is hyperbolic soliton solution.
All parts contained in the research carried out by the researcher through hard work and a review of the various references and contributions in the field of mathematics and the physical Applied.
 Antonowicz, M. and Rauch-Wojciechowski, S. (1992) Lax Representation for Restricted Flows of the KdV Hierarchy and for the Kepler Problem. Physics Letters A, 171, 303-310.
 Jawad, A.J.M., Petkovic, M.D. and Biswas, A. (2010) Modified Simple Equation Method for Nonlinear Evolution Equations. Applied Mathematics and Computation, 217, 869-877.
 Kamruzzaman, K., Ali Akbar, M. and Hj Mohd Ali, N. (2013) The Modified Simple Equation Method for Exact and Solitary Wave Solutions of Nonlinear Evolution Equation: The GZK-BBM Equation and Right-Handed Non Commutative Burgers Equations. ISRN Mathematical Physics, 2013, Article ID: 146704.
 Wang, M.L., Zhang, J.L. and Li, X.Z. (2008) The (GG)-Expansion Method and Travelling Wave Solutions of Nonlinear Evolutions Equations in Mathematical Physics. Physics Letters A, 372, 417-423.
 Zahran, E.H.M. and Khater, M.M.A. (2015) The Tow Variable (G’/G, 1/G)-Expansion Method for Solving Nonlinear Dynamics of Microtubules—A New Model. Global Journal of Science Frontier Research, 15, 87-94.
 Zahran, E.H.M. and Khater, M.M.A. (2014) Exact Travelling Wave Solutions for the System of Shallow Water Wave Equation Modified Liouvill Equation Using Extended Jacobi Elliptic Function Expansion Method. American Journal of Computational Mathematics, 4, 455-463.
 Shehata, M.S.M. (2015) Extended Jacobi Elliptic Function Expansion Method and its Applications for Solving Some Nonlinear Evolution Equations in Mathematical Physics. International Journal of Computer Applications, 109, 1-4.
 Zahran, E.H.M. (2015) Exact Traveling Wave Solutions for Nano-Solitons of Ionic Waves Propagation along Microtubules in Living Cells & Nano-Ionic Currents of MTs. World Journal of Nano Science and Engineering, 5, 78-87.
 Zahran, E.H.M. (2015) Exact Traveling Wave Solutions for Nano-Ionic Solitons & Nano-Ionic Currents of MTs Using Exp(-φ(ξ))-Expansion Method. Advances in Nano Particles, 4, 25-36.
 Zahran, E.H.M. (2015) Exact Traveling Wave Solutions for Nonlinear Fractional Partial Differential Equations Arising in Soliton Using the Exp(-φ(ξ))-Expansion Method. International Journal of Computer Applications, 109, 12-17.
 Shehata, M.S.M. (2015) The Exp(-φ(ξ))-Method and Its Applications for Solving some Nonlinear Evolution Equations in Mathematical Physics. American Journal of Computational Mathematics, 5, 468-480.
 Zahran, E.H.M. (2015) Travelling Wave Solutions of Non Linear Evolution Equation Via Modified Exp(-φ(ξ))-Expansion Method. Journal of Computational and Theoretical Nano Science, 12, 5716-5724.
 Shehata, M.S.M. (2016) A New Solitary Wave Solution of the Perturbed Nonlinear Schrodinger Equation Using a Riccati Bernoulli Sub-ODE Method. International Journal of Physical Sciences, 11, 80-84.
 Bekir, A. and Boz, A. (2007) Exact Solutions for a Class of Nonlinear Partial Differential Equations Using Exp-Function Method. International Journal of Nonlinear Sciences and Numerical Simulation, 8, 505-512.
 Shehata, M.S.M. (2016) Exact Traveling Wave Solutions for Nonlinear Evolutions Equation. Journal of Computational and Theoretical Nano Science, 13, 534-538.
 Zahran, E.H.M. and Khater, M.M.A. (2016) Modified Extended Tanh-Function Method and Its Applications to the Bogoyavlenskii Equation. Applied Mathematical Modeling, 40, 1769-1775.
 Lu, D.C., Khater Mostafa, M.A. and Zahran, E.H.M. (2017) Solitary Wave Solutions of the Benjamin Bona-Mahoney-Burgers Equation with Dual Power-Law Nonlinearity. Applied Mathematics & Information Sciences, 11, 1-5.
 Ucar, Y., Karaagas, B. and Esen, A. (2015) A New Approach on Numerical Solutions of the Improved Boussinesq Type Equation Using Quadratic B-Spline Galerkin Finite Element Method. Applied Mathematics and Computation, 270, 148-155.
 Wazwaz, A.-M. (2009) Multiple Kink Solutions and Multiple Singular Kink Solutions for (2+1)-Dimensional Nonlinear Models Generated by the Jaulent-Miodek Hierarchy. Physics Letters A, 373, 1844-1846.
 Mohebbi, A., Asgari, Z. and Dehghan, M. (2012) Numerical Solution of Nonlinear Jaulent-Miodek and Whitham-Broer-Kaup Equations. Communications in Nonlinear Science and Numerical Simulation, 17, 4602-4610.
 Liu, H. and Yan, F. (2011) The Bifurcation and Exact Travelling Wave Solutions for (2+1)-Dimensional Nonlinear Models Generated by the Jaulent-Miodek Hierarchy. International Journal of Nonlinear Science, 11, 200-205.