Welcome to visit Communications in Theoretical Physics,
Mathematical Physics

Study for obtaining the exact solutions for the propagation of pulses in optical fibers

  • Arzu Akbulut , 1 ,
  • Waseem Razzaq 2 ,
  • Filiz Tascan 3
Expand
  • 1Department of Mathematics, Faculty of Arts and Science, Bursa Uludağ University, Bursa, Türkiye
  • 2 Math Center, House no 87 Rahmanyia Colony, Vehari, Pakistan
  • 3Department of Mathematics and Computer Sciences, Faculty of Arts and Science, Eskisehir Osmangazi University, Eskisehir, Türkiye

Received date: 2023-11-04

  Revised date: 2024-03-30

  Accepted date: 2024-05-30

  Online published: 2024-08-16

Copyright

© 2024 Institute of Theoretical Physics CAS, Chinese Physical Society and IOP Publishing

Abstract

In this study, the exact solutions for the propagation of pulses in optical fibers are obtained. Special values are given in the model used, and two nonlinear differential equations are obtained. Nonlinear equations are reduced to ordinary differential equations with the help of wave transformations. Then, the obtained differential equations are solved by two different methods, namely the modified simplest equation and the modified Kudryashov procedures. The solutions are given by hyperbolic, trigonometric and rational functions and the results are useful for optics, engineering and other related areas. Finally three-dimensional, contour and two-dimensional shapes are given for some solutions. These figures are important for understanding the motion of the wave. The given methods are applied to the equations for the first time. To the best of the authors' knowledge, these results are new and have not been obtained in the literature. The results are useful for applied mathematics, physics and other related areas.

Cite this article

Arzu Akbulut , Waseem Razzaq , Filiz Tascan . Study for obtaining the exact solutions for the propagation of pulses in optical fibers[J]. Communications in Theoretical Physics, 2024 , 76(10) : 105003 . DOI: 10.1088/1572-9494/ad526c

1. Introduction

In the literature, nonlinear partial differential equations (NLPDEs) have an important role; therefore there has been much research into obtaining exact solutions for NLPDEs. As a result of these studies, many new methods for finding exact solutions of NLPDEs have been discovered. For example, Wazwaz applied the tanh technique to the Dodd–Bullough–Mikhailov equation [1], Hashemi et al utilized the Nucci reduction method with the (2+1)-dimensional Korteweg–de Vries equation [2], Biswas et al handled the method of undetermined coefficients [3], Zhou et al studied the Riccati equation mapping approach [4], Biswas et al applied the trial equation procedure to the Radhakrishnan–Kundu–Lakshmanan equation with a couple of integration schemes [5], Karakoc et al discussed the modified regularized long wave equation [6], Raza et al studied the extended trial approach [7], Hosseini et al applied the unified method to the generalized third-order nonlinear Schrödinger equation [8], Zafar et al utilized the modified extended tanh expansion method in the Biswas–Arshed model with full nonlinear form [9], Kaplan et al applied the exponential rational function procedure to the wave packet envelope [10], Rabie et al applied the extended F-expansion technique to the generalized derivative nonlinear Schrödinger equation with quintic nonlinearity [11], Hashemi et al handled the Nucci reduction method [12], Kudryashov and Lavrova employed the simplest equation method with the Chavy–Waddy–Kolokolnikov model [13] and Niwas and Kumar used the inverse (G'/G)-expansion method [14]. Finally, Kumar and colleagues used a generalized exponential rational function method [15], similarity reduction [16], the extended Jacobian elliptic function expansion method [17] and Lie symmetry analysis [18, 19].
The propagation pulse in the equation describing the propagation of pulses in optical fibers serves as a cornerstone for understanding, predicting and optimizing various phenomena in optical communication, nonlinear optics and other related fields, driving advances in diverse technological areas. This equation is important for the understanding signal transmission and exploring nonlinear phenomena. The equation describing the propagation of pulses in optical fibers, often referred to as the nonlinear Schrödinger equation or its variants, holds immense importance in several fields, particularly optical communications and nonlinear optics. Therefore, we have concluded that it is crucial to work on this model. In this paper, the propagation pulse in an optical fiber is explained using the following new PDE:
$\begin{eqnarray}\begin{array}{l}{\rm{i}}{W}_{t}+{W}_{{xx}}+{\alpha }_{1}\,| W{| }^{2m-2n}W\\ \,+\,{\beta }_{1}\,| W{| }^{2m-n}W+{\gamma }_{1}\,| W{| }^{2m}W+{\delta }_{1}\,| W{| }^{2m+n}W\\ \,+\,{\lambda }_{1}\,| W{| }^{2n+2m}W=0,\end{array}\end{eqnarray}$
where W(x, t) is the complex function representing an optical wave and α1, β1, γ1, δ1 and λ1 are the parameters. If we substitute m = 0 in equation (1) , the Schrödinger equation is obtained as follows:
$\begin{eqnarray}\begin{array}{l}{\rm{i}}{W}_{t}+{W}_{{xx}}+{\alpha }_{1}\,| W{| }^{-2n}W+{\beta }_{1}\,| W{| }^{-n}W\\ \,+\,{\gamma }_{1}\,W+{\delta }_{1}\,| W{| }^{n}W+{\lambda }_{1}\,| W{| }^{2n}W=0.\end{array}\end{eqnarray}$
As a result, we can say that equation (1) is a generalization of the Schrödinger equation.
The first equation
If we substitute m = n in equation (1), the following NLPDE is obtained:
$\begin{eqnarray}\begin{array}{l}{\rm{i}}{W}_{t}+{W}_{{xx}}+{\alpha }_{1}\,W+{\beta }_{1}\,| W{| }^{n}W\\ \,+\,{\gamma }_{1}\,| W{| }^{2n}W+{\delta }_{1}\,| W{| }^{3n}W+{\lambda }_{1}\,| W{| }^{4n}W=0.\end{array}\end{eqnarray}$
The second equation
If we substitute m = 2n, in equation (1), the following NLPDE is obtained:
$\begin{eqnarray}\begin{array}{l}{\rm{i}}{W}_{t}+{W}_{{xx}}+{\alpha }_{1}\,| W{| }^{2n}W+{\beta }_{1}\,| W{| }^{3n}W\\ \,+\,{\gamma }_{1}\,| W{| }^{4n}W+{\delta }_{1}\,| W{| }^{5n}W+{\lambda }_{1}\,| W{| }^{6n}W=0.\end{array}\end{eqnarray}$
The given models are handled in some studies in the literature, for example Kudryashov obtained the first integrals of the given models [20], Zayed and Alngar applied three different methods to obtain the optical soliton solutions [21] and Raza et al used the sine-Gordon expansion and modified auxiliary equation procedures [22].
In this paper we will obtain the exact solutions of equations (3) and (4) via the modified simplest equation and the modified Kudryashov procedures. For this purpose, a description of the procedures will be given in section 2. Then, the procedures will be applied in section 3. Some plots of the results will be given in section 4. Finally, conclusions will be given in section 5.

2. Description of the procedures

We take a PDE of the form
$\begin{eqnarray}P(W,{W}_{x},{W}_{t},{W}_{{xx}},{W}_{{xt}},{W}_{{tt}},\ldots )=0,\end{eqnarray}$
where P is a polynomial of W and W is a complex and unknown function.
Using the traveling wave transform
$\begin{eqnarray}W(x,t)=V(\varepsilon ){{\rm{e}}}^{{\rm{i}}{\rm{\Phi }}},\varepsilon =x-{c}_{1}t,{\rm{\Phi }}={k}_{1}x-{w}_{1}t,\end{eqnarray}$
we can modify equation (5), resulting in a nonlinear ordinary differential equation (ODE)
$\begin{eqnarray}Q(V,V^{\prime} ,V^{\prime\prime} ,V\prime\prime\prime ,\ldots )=0,\end{eqnarray}$
where a prime denotes the derivative respect to ϵ.

2.1. The modified simplest equation procedure

We extend the solution of equation (7) in a finite series form
$\begin{eqnarray}F(\varepsilon )=\displaystyle \sum _{i=1}^{m}{b}_{i}{\sigma }^{i}(\varepsilon ),\end{eqnarray}$
where bi (i = 1, 2, ..., m) are non-variables to be found and bm should be different from zero. m is a balance term and is calculated according to the principle of homogeneous balance. Let us determine the integer m by balancing the highest-order linear term with the highest-order nonlinear term in equation (7). The function σ(ϵ) fulfills some ODEs. In this research, we use the Riccati equations as the simplest equation
$\begin{eqnarray}\sigma ^{\prime} (\varepsilon )={\sigma }^{2}(\varepsilon )+\delta ,\end{eqnarray}$
where δ is a non-variable. After that we get family of solutions to equation (9) due to variations of δ.
If δ < 0, solutions are given as follows:
$\begin{eqnarray}\sigma (\varepsilon )=-\sqrt{-\delta }\tanh \sqrt{-\delta }\varepsilon ,\end{eqnarray}$
$\begin{eqnarray}\sigma (\varepsilon )=-\sqrt{-\delta }\coth \sqrt{-\delta }\varepsilon ,\end{eqnarray}$
$\begin{eqnarray}\begin{array}{l}\sigma (\varepsilon )=\sqrt{-\delta }\left(-\tanh \left(2\sqrt{-\delta }\,\varepsilon \right)\right.\\ \,\left.\pm \,\mathrm{isech}\left(2\sqrt{-\delta }\,\varepsilon \right)\right),\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{l}\sigma (\varepsilon )=\sqrt{-\delta }\left(-\coth \left(2\sqrt{-\delta }\,\varepsilon \right)\right.\\ \,\left.\pm \,\mathrm{csch}\left(2\sqrt{-\delta }\,\varepsilon \right)\right),\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{l}\sigma (\varepsilon )=-\displaystyle \frac{\sqrt{-\delta }}{2}\left(\tanh \left(\displaystyle \frac{\sqrt{-\delta }}{2}\,\varepsilon \right)\right.\\ \,\left.+\,\coth \left(\displaystyle \frac{\sqrt{-\delta }}{2}\,\varepsilon \right)\right).\end{array}\end{eqnarray}$
If δ > 0, the solutions are as follows:
$\begin{eqnarray}\sigma (\varepsilon )=\sqrt{\delta }\tan \sqrt{\delta }\varepsilon ,\end{eqnarray}$
$\begin{eqnarray}\sigma (\varepsilon )=-\sqrt{\delta }\cot \sqrt{\delta }\varepsilon ,\end{eqnarray}$
$\begin{eqnarray}\begin{array}{l}\sigma (\varepsilon )=\sqrt{\delta }\left(\tan \left(2\sqrt{\delta }\,\varepsilon \right)\right.\\ \,\left.\pm \,\sec \left(2\sqrt{\delta }\,\varepsilon \right)\right),\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{l}\sigma (\varepsilon )=\sqrt{\delta }\left(-\cot \left(2\sqrt{\delta }\,\varepsilon \right)\right.\\ \,\left.\pm \,\csc \left(2\sqrt{\delta }\,\varepsilon \right)\right),\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{l}\sigma (\varepsilon )=\displaystyle \frac{\sqrt{\delta }}{2}\left(\tan \left(\displaystyle \frac{\sqrt{\delta }}{2}\,\varepsilon \right)\right.\\ \,\left.-\,\cot \left(\displaystyle \frac{\sqrt{\delta }}{2}\,\varepsilon \right)\right).\end{array}\end{eqnarray}$
If δ = 0, the solution is as follows:
$\begin{eqnarray}\sigma (\varepsilon )=-\displaystyle \frac{1}{\varepsilon }.\end{eqnarray}$
Substituting equations (8) into (7) with (9) and setting all the coefficients of function σi to zero, the system of algebraic equations in the terms bi, c1, k1 and w1 is obtained; then we solve these algebraic equations. Finally, the obtained results are substituted into the auxiliary solutions [23, 24].

2.2. The modified Kudryashov procedure

We transform equation (7) using the wave transformation below
$\begin{eqnarray}F\left(\varepsilon \right)=\mathop{\displaystyle \sum _{i=0}}\limits^{m}{b}_{i}{\left(\sigma \left(\varepsilon \right)\right)}^{i},{b}_{m}\ne 0,\end{eqnarray}$
where ${b}_{i}\left(\,{,}(,i=0,1,\,\ldots ,\,m\right)$ are constants that will be determined later, m is calculated by the homogeneous balance principle and the function $\sigma \left(\varepsilon \right)$ is given by
$\begin{eqnarray}\sigma \left(\varepsilon \right)=\displaystyle \frac{1}{1+{\rm{\Gamma }}{a}^{\varepsilon }},\end{eqnarray}$
where (22) satisfies the following ODE:
$\begin{eqnarray}{\sigma }^{{}^{{\prime} }}\left(\varepsilon \right)=\left({\sigma }^{2}\left(\varepsilon \right)-\sigma \left(\varepsilon \right)\right)\mathrm{ln}a.\end{eqnarray}$
Without ignoring equation (23), equation (21) is substituted into equation (7) to produce a set of algebraic equations for bi, a, Γ, c1, k1 and w1. Eventually, after solving the resulting system, the exact solutions to equation (5) are determined [25, 26].

3. Application of the procedures

In this section, the given models will be reduced to an ODE by transformation (6). Then, two different procedures will be applied to the obtained reduced equations. Equations (3) and (4) can be reduced to the following nonlinear ODEs:
$\begin{eqnarray}\begin{array}{l}V^{\prime\prime} -({w}_{1}-{k}_{1}^{2}+{\alpha }_{1})V+{\beta }_{1}\,{V}^{n+1}+{\gamma }_{1}\,{V}^{2n+1}\\ \quad +\,{\delta }_{1}\,{V}^{3n+1}+2\,{\lambda }_{1}\,{V}^{4n+1}=0,\end{array}\end{eqnarray}$
and
$\begin{eqnarray}\begin{array}{l}V^{\prime\prime} -({w}_{1}-{k}_{1}^{2})V+{\alpha }_{1}\,{V}^{2n+1}+{\beta }_{1}\,{V}^{3n+1}\\ \quad +\,{\gamma }_{1}\,{V}^{4n+1}+{\delta }_{1}\,{V}^{5n+1}+2\,{\lambda }_{1}\,{V}^{6n+1}=0.\end{array}\end{eqnarray}$
The imaginary parts of equations (3) and (4) give the speed of the traveling wave as
$\begin{eqnarray}{c}_{1}=2{k}_{1}.\end{eqnarray}$
In order to obtain the closed form solutions for equation (24), we apply the following transformation:
$\begin{eqnarray}V={F}^{\displaystyle \frac{1}{2n}},\end{eqnarray}$
$\begin{eqnarray}\begin{array}{l}\displaystyle \frac{1}{2n}\left(\displaystyle \frac{1}{2n}-1\right){\left(F^{\prime} \right)}^{2}+\displaystyle \frac{1}{2n}\,F\,F^{\prime\prime} +({w}_{1}-{k}_{1}^{2}+{\alpha }_{1}){F}^{2}+{\beta }_{1}\,{F}^{\tfrac{5}{2}}\\ \quad +{\gamma }_{1}\,{F}^{3}+{\delta }_{1}\,{F}^{\tfrac{7}{2}}+2\,{\lambda }_{1}\,{F}^{4}=0.\end{array}\end{eqnarray}$
To obtain closed form solutions for equation (28), the constraint condition β1 = δ1 = 0 is imposed. Equation (28) becomes
$\begin{eqnarray}\begin{array}{l}\displaystyle \frac{1}{2n}\left(\displaystyle \frac{1}{2n}-1\right){\left(F^{\prime} \right)}^{2}+\displaystyle \frac{1}{2n}\,F\,F^{\prime\prime} +({w}_{1}-{k}_{1}^{2}\\ \quad +\,{\alpha }_{1}){F}^{2}+{\gamma }_{1}\,{F}^{3}+2\,{\lambda }_{1}\,{F}^{4}=0.\end{array}\end{eqnarray}$
To obtain the closed form solutions for equation (25), we apply the following transformation:
$\begin{eqnarray}V={F}^{\displaystyle \frac{1}{3n}}.\end{eqnarray}$
Equation (25) takes the following form:
$\begin{eqnarray}\begin{array}{l}\displaystyle \frac{1}{3n}\left(\displaystyle \frac{1}{3n}-1\right){\left(F^{\prime} \right)}^{2}+\displaystyle \frac{1}{3n}\,F\,F^{\prime\prime} +({w}_{1}-{k}_{1}^{2}){F}^{2}\\ \quad +\,{\alpha }_{1}\,{F}^{\tfrac{8}{3}}+{\beta }_{1}\,{F}^{3}+{\gamma }_{1}\,{F}^{\tfrac{10}{3}}\\ \quad +\,{\delta }_{1}\,{F}^{\tfrac{11}{3}}+{\lambda }_{1}\,{F}^{4}=0.\end{array}\end{eqnarray}$
To obtain closed form solutions for equation (31), the constraint condition α1 = γ1 = δ1 = 0 is been imposed. Equation (31) becomes
$\begin{eqnarray}\begin{array}{l}\displaystyle \frac{1}{3n}\left(\displaystyle \frac{1}{3n}-1\right){\left(F^{\prime} \right)}^{2}+\displaystyle \frac{1}{3n}\,F\,F^{\prime\prime} +({w}_{1}-{k}_{1}^{2}){F}^{2}\\ \quad +\,{\beta }_{1}\,{F}^{3}+{\lambda }_{1}\,{F}^{4}=0.\end{array}\end{eqnarray}$
Balancing ${\left(F^{\prime} \right)}^{2}$ and F4 in equations (29) and (32), we then derive m = 1.

3.1. The modified simplest equation procedure

We choose the auxiliary solution as follows for the equations because the balance number is 1:
$\begin{eqnarray}F(\varepsilon )={b}_{0}+{b}_{1}\sigma (\varepsilon ).\end{eqnarray}$
The first equation
By substituting equation (33) into (29) and using some mathematical operations, we obtain a nonlinear algebraic system. Solving this system for different values of b0, b1, λ1 and w1 results in Cases 1 and 2.
Case 1:
$\begin{eqnarray*}\begin{array}{l}\left\{{b}_{0}\to \displaystyle \frac{-n\delta -\delta }{{\gamma }_{1}{n}^{2}},{b}_{1}\to -\displaystyle \frac{{\rm{i}}(n+1)\sqrt{\delta }}{{\gamma }_{1}{n}^{2}},\right.\\ \left.{\lambda }_{1}\to \displaystyle \frac{{\gamma }_{1}^{2}{n}^{2}(2n+1)}{4{\left(n+1\right)}^{2}\delta },{w}_{1}\to \displaystyle \frac{{k}_{1}^{2}{n}^{2}-{\alpha }_{1}{n}^{2}+\delta }{{n}^{2}}\right\}.\end{array}\end{eqnarray*}$
When δ < 0, the following solutions are obtained:
$\begin{eqnarray}\begin{array}{l}W(x,t)={{\rm{e}}}^{{\rm{i}}({k}_{1}\,x-{w}_{1}\,t)}\\ \,\times \,{\left(-\displaystyle \frac{(n+1)\sqrt{\delta }\left(\sqrt{\delta }+{\rm{i}}\sqrt{-\delta }\tanh \left(\varepsilon \sqrt{-\delta }\right)\right)}{{\gamma }_{1}{n}^{2}}\right)}^{\tfrac{1}{2n}},\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{l}W(x,t)={{\rm{e}}}^{{\rm{i}}({k}_{1}\,x-{w}_{1}\,t)}\\ \,\times \,{\left(-\displaystyle \frac{(n+1)\sqrt{\delta }\left(\sqrt{\delta }+{\rm{i}}\sqrt{-\delta }\coth \left(\varepsilon \sqrt{-\delta }\right)\right)}{{\gamma }_{1}{n}^{2}}\right)}^{\tfrac{1}{2n}},\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{l}W(x,t)={{\rm{e}}}^{{\rm{i}}({k}_{1}x-{w}_{1}t)}\\ \,\times \,{\left(-\displaystyle \frac{(n+1)\sqrt{\delta }\left(\sqrt{\delta }+{\rm{i}}\sqrt{-\delta }\left(-\tanh \left(2\varepsilon \sqrt{-\delta }\right)\pm i{\rm{sech}} \left(2\varepsilon \sqrt{-\delta }\right)\right)\right)}{{\gamma }_{1}{n}^{2}}\right)}^{\tfrac{1}{2n}},\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{l}W(x,t)={{\rm{e}}}^{{\rm{i}}({k}_{1}x-{w}_{1}t)}\\ \,\times \,{\left(-\displaystyle \frac{(n+1)\sqrt{\delta }\left(\sqrt{\delta }+{\rm{i}}\sqrt{-\delta }\left(-\coth \left(2\varepsilon \sqrt{-\delta }\right)\pm \mathrm{csch}\left(2\varepsilon \sqrt{-\delta }\right)\right)\right)}{{\gamma }_{1}{n}^{2}}\right)}^{\tfrac{1}{2n}},\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{l}W(x,t)={{\rm{e}}}^{{\rm{i}}({k}_{1}x-{w}_{1}t)}\\ \,\times \,{\left(-\displaystyle \frac{(n+1)\left(\delta -{\rm{i}}\sqrt{-{\delta }^{2}}\coth \left(\varepsilon \sqrt{-\delta }\right)\right)}{{\gamma }_{1}{n}^{2}}\right)}^{\tfrac{1}{2n}}.\end{array}\end{eqnarray}$
When δ > 0, the following solutions are obtained:
$\begin{eqnarray}\begin{array}{l}W(x,t)={{\rm{e}}}^{{\rm{i}}({k}_{1}x-{w}_{1}t)}\\ \,\times \,{\left(-\displaystyle \frac{{\rm{i}}(n+1)\delta \left(\tan \left(\varepsilon \sqrt{\delta }\right)-{\rm{i}}\right)}{{\gamma }_{1}{n}^{2}}\right)}^{\tfrac{1}{2n}},\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{l}W(x,t)={{\rm{e}}}^{{\rm{i}}({k}_{1}x-{w}_{1}t)}\\ \,\times \,{\left(\displaystyle \frac{{\rm{i}}(n+1)\delta \left(\cot \left(\varepsilon \sqrt{\delta }\right)+{\rm{i}}\right)}{{\gamma }_{1}{n}^{2}}\right)}^{\tfrac{1}{2n}},\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{l}W(x,t)={{\rm{e}}}^{{\rm{i}}({k}_{1}x-{w}_{1}t)}\\ \,\times \,{\left(-\displaystyle \frac{{\rm{i}}(n+1)\delta \left(\left(\tan \left(2\varepsilon \sqrt{\delta }\right)\pm \sec \left(2\varepsilon \sqrt{\delta }\right)\right)-{\rm{i}}\right)}{{\gamma }_{1}{n}^{2}}\right)}^{\tfrac{1}{2n}},\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{l}W(x,t)={{\rm{e}}}^{{\rm{i}}({k}_{1}x-{w}_{1}t)}\\ \,\times \,{\left(-\displaystyle \frac{{\rm{i}}(n+1)\delta \left(\left(-\cot \left(2\varepsilon \sqrt{\delta }\right)\pm \csc \left(2\varepsilon \sqrt{\delta }\right)\right)-{\rm{i}}\right)}{{\gamma }_{1}{n}^{2}}\right)}^{\tfrac{1}{2n}},\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{l}W(x,t)={{\rm{e}}}^{{\rm{i}}({k}_{1}x-{w}_{1}t)}\\ \,\times \,{\left(\displaystyle \frac{{\rm{i}}(n+1)\delta \left(\cot \left(\varepsilon \sqrt{\delta }\right)+{\rm{i}}\right)}{{\gamma }_{1}{n}^{2}}\right)}^{\tfrac{1}{2n}},\end{array}\end{eqnarray}$
where ϵ = xc1 t.
Case 2:
$\begin{eqnarray*}\begin{array}{l}\left\{{b}_{0}\to \displaystyle \frac{-n\delta -\delta }{{\gamma }_{1}{n}^{2}},{b}_{1}\to \displaystyle \frac{{\rm{i}}(n+1)\sqrt{\delta }}{{\gamma }_{1}{n}^{2}},\right.\\ \left.{\lambda }_{1}\to \displaystyle \frac{{\gamma }_{1}^{2}{n}^{2}(2n+1)}{4{\left(n+1\right)}^{2}\delta },{w}_{1}\to \displaystyle \frac{{k}_{1}^{2}{n}^{2}-{\alpha }_{1}{n}^{2}+\delta }{{n}^{2}}\right\}.\end{array}\end{eqnarray*}$
When δ < 0, the following solutions are obtained:
$\begin{eqnarray}\begin{array}{l}W(x,t)={{\rm{e}}}^{{\rm{i}}({k}_{1}x-{w}_{1}t)}\\ \,\times \,{\left(-\displaystyle \frac{(n+1)\sqrt{\delta }\left(\sqrt{\delta }-{\rm{i}}\sqrt{-\delta }\tanh \left(\varepsilon \sqrt{-\delta }\right)\right)}{{\gamma }_{1}{n}^{2}}\right)}^{\tfrac{1}{2n}},\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{l}W(x,t)={{\rm{e}}}^{{\rm{i}}({k}_{1}x-{w}_{1}t)}\\ \,\times \,{\left(-\displaystyle \frac{(n+1)\sqrt{\delta }\left(\sqrt{\delta }-{\rm{i}}\sqrt{-\delta }\coth \left(\varepsilon \sqrt{-\delta }\right)\right)}{{\gamma }_{1}{n}^{2}}\right)}^{\tfrac{1}{2n}},\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{l}W(x,t)={{\rm{e}}}^{{\rm{i}}({k}_{1}x-{w}_{1}t)}\\ \,\times \,{\left(-\displaystyle \frac{(n+1)\sqrt{\delta }\left(\sqrt{\delta }-{\rm{i}}\sqrt{-\delta }\left(-\tanh \left(2\varepsilon \sqrt{-\delta }\right)\pm \mathrm{isech}\left(2\varepsilon \sqrt{-\delta }\right)\right)\right)}{{\gamma }_{1}{n}^{2}}\right)}^{\tfrac{1}{2n}},\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{l}W(x,t)={{\rm{e}}}^{{\rm{i}}({k}_{1}x-{w}_{1}t)}\\ \,\times \,{\left(-\displaystyle \frac{(n+1)\sqrt{\delta }\left(\sqrt{\delta }-{\rm{i}}\sqrt{-\delta }\left(-\coth \left(2\varepsilon \sqrt{-\delta }\right)\pm \mathrm{csch}\left(2\varepsilon \sqrt{-\delta }\right)\right)\right)}{{\gamma }_{1}{n}^{2}}\right)}^{\tfrac{1}{2n}},\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{l}W(x,t)={{\rm{e}}}^{{\rm{i}}({k}_{1}x-{w}_{1}t)}\\ \,\times \,{\left(-\displaystyle \frac{(n+1)\left(\delta +{\rm{i}}\sqrt{-{\delta }^{2}}\coth \left(\varepsilon \sqrt{-\delta }\right)\right)}{{\gamma }_{1}{n}^{2}}\right)}^{\tfrac{1}{2n}}.\end{array}\end{eqnarray}$
When δ > 0, the following solutions are obtained:
$\begin{eqnarray}\begin{array}{l}W(x,t)={{\rm{e}}}^{{\rm{i}}({k}_{1}x-{w}_{1}t)}\\ \,\times {\left(\displaystyle \frac{{\rm{i}}(n+1)\delta \left(\tan \left(\varepsilon \sqrt{\delta }\right)+{\rm{i}}\right)}{{\gamma }_{1}{n}^{2}}\right)}^{\tfrac{1}{2n}},\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{l}W(x,t)={{\rm{e}}}^{{\rm{i}}({k}_{1}x-{w}_{1}t)}\\ \,\times \,{\left(-\displaystyle \frac{{\rm{i}}(n+1)\delta \left(\cot \left(\varepsilon \sqrt{\delta }\right)-{\rm{i}}\right)}{{\gamma }_{1}{n}^{2}}\right)}^{\tfrac{1}{2n}},\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{l}W(x,t)={{\rm{e}}}^{{\rm{i}}({k}_{1}x-{w}_{1}t)}\\ \,\times \,{\left(\displaystyle \frac{{\rm{i}}(n+1)\delta \left(\left(\tan \left(2\varepsilon \sqrt{\delta }\right)\pm \sec \left(2\varepsilon \sqrt{\delta }\right)\right)+{\rm{i}}\right)}{{\gamma }_{1}{n}^{2}}\right)}^{\tfrac{1}{2n}},\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{l}W(x,t)={{\rm{e}}}^{{\rm{i}}({k}_{1}x-{w}_{1}t)}\\ \,\times \,{\left(\displaystyle \frac{{\rm{i}}(n+1)\delta \left(\left(-\cot \left(2\varepsilon \sqrt{\delta }\right)\pm \csc \left(2\varepsilon \sqrt{\delta }\right)\right)+{\rm{i}}\right)}{{\gamma }_{1}{n}^{2}}\right)}^{\tfrac{1}{2n}},\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{l}W(x,t)={{\rm{e}}}^{{\rm{i}}({k}_{1}x-{w}_{1}t)}\\ \,\times \,{\left(-\displaystyle \frac{{\rm{i}}(n+1)\delta \left(\cot \left(\varepsilon \sqrt{\delta }\right)-{\rm{i}}\right)}{{\gamma }_{1}{n}^{2}}\right)}^{\tfrac{1}{2n}},\end{array}\end{eqnarray}$
where ϵ = xc1 t.
The second equation
By substituting equation (33) into equation (32) and using some mathematical operations, we obtain a nonlinear algebraic system. Solving this system for different values of b0, b1, λ1 and ω1 results in Cases 1 and 2.
Case 1:
$\begin{eqnarray*}\begin{array}{l}\left\{{b}_{0}\to -\displaystyle \frac{2(3n\delta +2\delta )}{9{\beta }_{1}{n}^{2}},{b}_{1}\to -\displaystyle \frac{2{\rm{i}}(3n+2)\sqrt{\delta }}{9{\beta }_{1}{n}^{2}},\right.\\ \left.{\lambda }_{1}\to \displaystyle \frac{9{\beta }_{1}^{2}{n}^{2}(3n+1)}{4{\left(3n+2\right)}^{2}\delta },{\omega }_{1}\to \displaystyle \frac{9{k}^{2}{n}^{2}+4\delta }{9{n}^{2}}\right\}.\end{array}\end{eqnarray*}$
When δ < 0, the following solutions are obtained:
$\begin{eqnarray}\begin{array}{l}W(x,t)={{\rm{e}}}^{{\rm{i}}({k}_{1}x-{\omega }_{1}t)}\\ \,\times \,{\left(-\displaystyle \frac{2(3n+2)\sqrt{\delta }\left(\sqrt{\delta }+{\rm{i}}\sqrt{-\delta }\tanh \left(\varepsilon \sqrt{-\delta }\right)\right)}{9{\beta }_{1}{n}^{2}}\right)}^{\tfrac{1}{3n}},\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{l}W(x,t)={{\rm{e}}}^{{\rm{i}}({k}_{1}x-{\omega }_{1}t)}\\ \,\times \,{\left(-\displaystyle \frac{2(3n+2)\sqrt{\delta }\left(\sqrt{\delta }+{\rm{i}}\sqrt{-\delta }\coth \left(\varepsilon \sqrt{-\delta }\right)\right)}{9{\beta }_{1}{n}^{2}}\right)}^{\tfrac{1}{3n}},\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{l}W(x,t)={{\rm{e}}}^{{\rm{i}}({k}_{1}x-{\omega }_{1}t)}\\ \,\times \,{\left(-\displaystyle \frac{2(3n+2)\sqrt{\delta }\left(\sqrt{\delta }+{\rm{i}}\sqrt{-\delta }\left(-\tanh \left(2\varepsilon \sqrt{-\delta }\right)\pm \mathrm{isech}\left(2\varepsilon \sqrt{-\delta }\right)\right)\right)}{9{\beta }_{1}{n}^{2}}\right)}^{\tfrac{1}{3n}},\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{l}W(x,t)={{\rm{e}}}^{{\rm{i}}({k}_{1}x-{\omega }_{1}t)}\\ \,\times \,{\left(-\displaystyle \frac{2(3n+2)\sqrt{\delta }\left(\sqrt{\delta }+{\rm{i}}\sqrt{-\delta }\left(-\coth \left(2\varepsilon \sqrt{-\delta }\right)\pm \mathrm{csch}\left(2\varepsilon \sqrt{-\delta }\right)\right)\right)}{9{\beta }_{1}{n}^{2}}\right)}^{\tfrac{1}{3n}},\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{l}W(x,t)={{\rm{e}}}^{{\rm{i}}({k}_{1}x-{\omega }_{1}t)}\\ \,\times \,{\left(-\displaystyle \frac{2(3n+2)\sqrt{\delta }\left(\sqrt{\delta }-{\rm{i}}\sqrt{-\delta }\coth \left(\varepsilon \sqrt{-\delta }\right)\right)}{9{\beta }_{1}{n}^{2}}\right)}^{\tfrac{1}{3n}}.\end{array}\end{eqnarray}$
When δ > 0, the following solutions are obtained:
$\begin{eqnarray}\begin{array}{l}W(x,t)={{\rm{e}}}^{{\rm{i}}({k}_{1}x-{\omega }_{1}t)}\\ \,\times \,{\left(-\displaystyle \frac{2{\rm{i}}(3n+2)\delta \left(\tan \left(\varepsilon \sqrt{\delta }\right)-{\rm{i}}\right)}{9{\beta }_{1}{n}^{2}}\right)}^{\tfrac{1}{3n}},\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{l}W(x,t)={{\rm{e}}}^{{\rm{i}}({k}_{1}x-{\omega }_{1}t)}\\ \,\times \,{\left(\displaystyle \frac{2{\rm{i}}(3n+2)\delta \left(\cot \left(\varepsilon \sqrt{\delta }\right)+{\rm{i}}\right)}{9{\beta }_{1}{n}^{2}}\right)}^{\tfrac{1}{3n}},\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{l}W(x,t)={{\rm{e}}}^{{\rm{i}}({k}_{1}x-{\omega }_{1}t)}\\ \,\times \,{\left(-\displaystyle \frac{2{\rm{i}}(3n+2)\delta \left(\left(\tan \left(2\varepsilon \sqrt{\delta }\right)\pm \sec \left(2\varepsilon \sqrt{\delta }\right)\right)-{\rm{i}}\right)}{9{\beta }_{1}{n}^{2}}\right)}^{\tfrac{1}{3n}},\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{l}W(x,t)={{\rm{e}}}^{{\rm{i}}({k}_{1}x-{\omega }_{1}t)}\\ \,\times \,{\left(-\displaystyle \frac{2{\rm{i}}(3n+2)\delta \left(\left(-\cot \left(2\varepsilon \sqrt{\delta }\right)\pm \csc \left(2\varepsilon \sqrt{\delta }\right)\right)-{\rm{i}}\right)}{9{\beta }_{1}{n}^{2}}\right)}^{\tfrac{1}{3n}},\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{l}W(x,t)={{\rm{e}}}^{{\rm{i}}({k}_{1}x-{\omega }_{1}t)}\\ \,\times \,{\left(\displaystyle \frac{2i(3n+2)\delta \left(\cot \left(\varepsilon \sqrt{\delta }\right)+{\rm{i}}\right)}{9{\beta }_{1}{n}^{2}}\right)}^{\tfrac{1}{3n}},\end{array}\end{eqnarray}$
where ϵ = xc1 t.
Case 2:
$\begin{eqnarray*}\begin{array}{l}\left\{{b}_{0}\to -\displaystyle \frac{2(3n\delta +2\delta )}{9{\beta }_{1}{n}^{2}},{b}_{1}\to \displaystyle \frac{2{\rm{i}}(3n+2)\sqrt{\delta }}{9{\beta }_{1}{n}^{2}},\right.\\ \left.{\lambda }_{1}\to \displaystyle \frac{9{\beta }_{1}^{2}{n}^{2}(3n+1)}{4{\left(3n+2\right)}^{2}\delta },{\omega }_{1}\to \displaystyle \frac{9{k}^{2}{n}^{2}+4\delta }{9{n}^{2}}\right\}.\end{array}\end{eqnarray*}$
When δ < 0, the following solutions are obtained:
$\begin{eqnarray}\begin{array}{l}W(x,t)={{\rm{e}}}^{{\rm{i}}({k}_{1}x-{\omega }_{1}t)}\\ \,\times \,{\left(-\displaystyle \frac{2(3n+2)\sqrt{\delta }\left(\sqrt{\delta }-{\rm{i}}\sqrt{-\delta }\tanh \left(\varepsilon \sqrt{-\delta }\right)\right)}{9{\beta }_{1}{n}^{2}}\right)}^{\tfrac{1}{3n}},\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{l}W(x,t)={{\rm{e}}}^{{\rm{i}}({k}_{1}x-{\omega }_{1}t)}\\ \,\times \,{\left(-\displaystyle \frac{2(3n+2)\sqrt{\delta }\left(\sqrt{\delta }-{\rm{i}}\sqrt{-\delta }\coth \left(\varepsilon \sqrt{-\delta }\right)\right)}{9{\beta }_{1}{n}^{2}}\right)}^{\tfrac{1}{3n}},\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{l}W(x,t)={{\rm{e}}}^{{\rm{i}}({k}_{1}x-{\omega }_{1}t)}\\ \,\times \,{\left(-\displaystyle \frac{2(3n+2)\sqrt{\delta }\left(\sqrt{\delta }-{\rm{i}}\sqrt{-\delta }\left(-\tanh \left(2\varepsilon \sqrt{-\delta }\right)\pm \mathrm{isech}\left(2\varepsilon \sqrt{-\delta }\right)\right)\right)}{9{\beta }_{1}{n}^{2}}\right)}^{\tfrac{1}{3n}},\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{l}W(x,t)={{\rm{e}}}^{{\rm{i}}({k}_{1}x-{\omega }_{1}t)}\\ \,\times \,{\left(-\displaystyle \frac{2(3n+2)\sqrt{\delta }\left(\sqrt{\delta }-{\rm{i}}\sqrt{-\delta }\left(-\coth \left(2\varepsilon \sqrt{-\delta }\right)\pm \mathrm{csch}\left(2\varepsilon \sqrt{-\delta }\right)\right)\right)}{9{\beta }_{1}{n}^{2}}\right)}^{\tfrac{1}{3n}},\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{l}W(x,t)={{\rm{e}}}^{{\rm{i}}({k}_{1}x-{\omega }_{1}t)}\\ \,\times \,{\left(-\displaystyle \frac{2(3n+2)\sqrt{\delta }\left(\sqrt{\delta }+{\rm{i}}\sqrt{-\delta }\coth \left(\varepsilon \sqrt{-\delta }\right)\right)}{9{\beta }_{1}{n}^{2}}\right)}^{\tfrac{1}{3n}}.\end{array}\end{eqnarray}$
When δ > 0, the following solutions are obtained:
$\begin{eqnarray}\begin{array}{l}W(x,t)={{\rm{e}}}^{{\rm{i}}({k}_{1}x-{\omega }_{1}t)}\\ \,\times \,{\left(\displaystyle \frac{2{\rm{i}}(3n+2)\delta \left(\tan \left(\varepsilon \sqrt{\delta }\right)+i\right)}{9{\beta }_{1}{n}^{2}}\right)}^{\tfrac{1}{3n}},\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{l}W(x,t)={{\rm{e}}}^{{\rm{i}}({k}_{1}x-{\omega }_{1}t)}\\ \,\times \,{\left(-\displaystyle \frac{2{\rm{i}}(3n+2)\delta \left(\cot \left(\varepsilon \sqrt{\delta }\right)-{\rm{i}}\right)}{9{\beta }_{1}{n}^{2}}\right)}^{\tfrac{1}{3n}},\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{l}W(x,t)={{\rm{e}}}^{{\rm{i}}({k}_{1}x-{\omega }_{1}t)}\\ \,\times \,{\left(\displaystyle \frac{2{\rm{i}}(3n+2)\delta \left(\left(\tan \left(2\varepsilon \sqrt{\delta }\right)\pm \sec \left(2\varepsilon \sqrt{\delta }\right)\right)+{\rm{i}}\right)}{9{\beta }_{1}{n}^{2}}\right)}^{\tfrac{1}{3n}},\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{l}W(x,t)={{\rm{e}}}^{{\rm{i}}({k}_{1}x-{\omega }_{1}t)}\\ \,\times \,{\left(\displaystyle \frac{2{\rm{i}}(3n+2)\delta \left(\left(-\cot \left(2\varepsilon \sqrt{\delta }\right)\pm \csc \left(2\varepsilon \sqrt{\delta }\right)\right)+{\rm{i}}\right)}{9{\beta }_{1}{n}^{2}}\right)}^{\tfrac{1}{3n}},\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{l}W(x,t)={{\rm{e}}}^{{\rm{i}}({k}_{1}x-{\omega }_{1}t)}\\ \,\times \,{\left(-\displaystyle \frac{2{\rm{i}}(3n+2)\delta \left(\cot \left(\varepsilon \sqrt{\delta }\right)-{\rm{i}}\right)}{9{\beta }_{1}{n}^{2}}\right)}^{\tfrac{1}{3n}},\end{array}\end{eqnarray}$
where ϵ = xc1 t.

3.2. The modified Kudryashov procedure

The first equation
If we take the following auxiliary solution for equation (29):
$\begin{eqnarray}F\left(\varepsilon \right)={b}_{0}+{b}_{1}\sigma \left(\varepsilon \right),\end{eqnarray}$
then substitute equation (74) into equation (29) and collect the coefficients of ${\sigma }^{i}\left(\varepsilon \right)$, we obtain the determining equation system. If we solve the obtained system, we get the following solution family (Cases 1 and 2).
Case 1:
$\begin{eqnarray}\begin{array}{l}{b}_{0}\to 0,{b}_{1}\to \displaystyle \frac{{\left(\mathrm{ln}a\right)}^{2}\left(n+1\right)}{2{\gamma }_{1}{n}^{2}},\\ {\lambda }_{1}\to -\displaystyle \frac{{\gamma }_{1}^{2}{n}^{2}\left(2n+1\right)}{{\left(\mathrm{ln}a\right)}^{2}\left({n}^{2}+2n+1\right)},\\ {\omega }_{1}\to -\displaystyle \frac{-4{k}^{2}{n}^{2}+4{\alpha }_{1}{n}^{2}+{\left(\mathrm{ln}a\right)}^{2}}{4{n}^{2}}.\end{array}\end{eqnarray}$
If we substitute the obtained values of the coefficients from equations (75) into equation (3) without ignoring equations (26), (27) and (74), we get the following solution:
$\begin{eqnarray}\begin{array}{l}W\left(x,t\right)={{\rm{e}}}^{{\rm{i}}\left({k}_{1}x-{\omega }_{1}t\right)}\\ \,\times \,{\left(\displaystyle \frac{{\left(\mathrm{ln}a\right)}^{2}\left(n+1\right)}{2{\gamma }_{1}{n}^{2}\left(1+{\rm{\Gamma }}{a}^{x-2{k}_{1}t}\right)}\right)}^{\tfrac{1}{2n}}.\end{array}\end{eqnarray}$
Case 2:
$\begin{eqnarray}\begin{array}{l}{b}_{0}\to \displaystyle \frac{{\left(\mathrm{ln}a\right)}^{2}\left(n+1\right)}{2{\gamma }_{1}{n}^{2}},\\ {b}_{1}\to -\displaystyle \frac{{\left(\mathrm{ln}a\right)}^{2}\left(n+1\right)}{2{\gamma }_{1}{n}^{2}},\\ {\lambda }_{1}\to -\displaystyle \frac{{\gamma }_{1}^{2}{n}^{2}\left(2n+1\right)}{{\left(\mathrm{ln}a\right)}^{2}\left({n}^{2}+2n+1\right)},\end{array}\end{eqnarray}$
$\begin{eqnarray*}{\omega }_{1}\to -\displaystyle \frac{-4{k}^{2}{n}^{2}+4{\alpha }_{1}{n}^{2}+{\left(\mathrm{ln}a\right)}^{2}}{4{n}^{2}}.\end{eqnarray*}$
If we use the values from (77), we obtain the following exact solution:
$\begin{eqnarray}\begin{array}{l}W\left(x,t\right)={{\rm{e}}}^{{\rm{i}}\left({k}_{1}x-{\omega }_{1}t\right)}\\ \quad \times \,{\left(\displaystyle \frac{{\left(\mathrm{ln}a\right)}^{2}{a}^{x-2{k}_{1}t}\left(n+1\right){\rm{\Gamma }}}{2{\gamma }_{1}{n}^{2}\left(1+{\rm{\Gamma }}{a}^{x-2{k}_{1}t}\right)}\right)}^{\tfrac{1}{2n}}.\end{array}\end{eqnarray}$
The second equation
If we take the auxiliary solution for equation (32) as (74), then substitute equation (74) in (32), then collect the coefficients of ${\sigma }^{i}\left(\varepsilon \right),$ we obtain the determining equation system. If we solve the obtained system, we get the following solution family (Cases 1 and 2):
Case 1:
$\begin{eqnarray}\begin{array}{l}{b}_{0}\to 0,{b}_{1}\to \displaystyle \frac{{\left(\mathrm{ln}a\right)}^{2}\left(3n+2\right)}{9{\beta }_{1}{n}^{2}},\\ {\lambda }_{1}\to -\displaystyle \frac{9{\beta }_{1}^{2}{n}^{2}\left(3n+1\right)}{{\left(\mathrm{ln}a\right)}^{2}\left(9{n}^{2}+12n+4\right)},\\ {\omega }_{1}\to -\displaystyle \frac{-9{k}^{2}{n}^{2}+{\left(\mathrm{ln}a\right)}^{2}}{9{n}^{2}}.\end{array}\end{eqnarray}$
If we substitute the obtained values of the coefficients equations (79) in (4) without ignoring equations (26), (27) and (74), we get the following solution:
$\begin{eqnarray}\begin{array}{l}W\left(x,t\right)={{\rm{e}}}^{{\rm{i}}\left({k}_{1}x-{\omega }_{1}t\right)}\\ \quad \times \,{\left(\displaystyle \frac{{\left(\mathrm{ln}a\right)}^{2}\left(3n+2\right)}{9{\beta }_{1}{n}^{2}\left(1+{\rm{\Gamma }}{a}^{x-2{k}_{1}t}\right)}\right)}^{\tfrac{1}{3n}}.\end{array}\end{eqnarray}$
Case 2:
$\begin{eqnarray}\begin{array}{l}{b}_{0}\to \displaystyle \frac{{\left(\mathrm{ln}a\right)}^{2}\left(3n+2\right)}{9{\beta }_{1}{n}^{2}},\\ {b}_{1}\to -\displaystyle \frac{{\left(\mathrm{ln}a\right)}^{2}\left(3n+2\right)}{9{\beta }_{1}{n}^{2}},\\ {\lambda }_{1}\to -\displaystyle \frac{9{\beta }_{1}^{2}{n}^{2}\left(3n+1\right)}{{\left(\mathrm{ln}a\right)}^{2}\left(9{n}^{2}+12n+4\right)},\end{array}\end{eqnarray}$
$\begin{eqnarray*}{\omega }_{1}\to -\displaystyle \frac{-9{k}^{2}{n}^{2}+{\left(\mathrm{ln}a\right)}^{2}}{9{n}^{2}}.\end{eqnarray*}$
If we use values from (81), the exact solution is given by
$\begin{eqnarray}\begin{array}{l}W\left(x,t\right)={{\rm{e}}}^{{\rm{i}}\left({k}_{1}x-{\omega }_{1}t\right)}\\ \,\times \,{\left(\displaystyle \frac{{\left(\mathrm{ln}a\right)}^{2}{a}^{x-2{k}_{1}t}\left(3n+2\right){\rm{\Gamma }}}{9{\beta }_{1}{n}^{2}\left(1+{\rm{\Gamma }}{a}^{x-2{k}_{1}t}\right)}\right)}^{\tfrac{1}{3n}}.\end{array}\end{eqnarray}$

4. Plots of the results

Plots of some of the results will be given in this section. First, we give the plots for equation (34) in figure 1. When we set the parameters ω1 = − 0.8, n = 2, k1 = 0.5, α1 = 0.2 and γ1 = 0.7 in equation (34), the three-dimensional (3D) and contour plots shown in figure 1 are obtained. When we set t = 0, t = 0.5 and t = 1, red, green and blue lines, respectively, are obtained in the two-dimensional (2D) plot in figure 1; this figure represents the kink shape solution.
Figure 1. Plots of equation (34), when ω1 = − 0.8, n = 2, k1 = 0.5, α1 = 0.2 and γ1 = 0.7.
We show the plots for equation (59) in figure 2. When we set the parameters ω1 = 0.8, n = 2, k1 = 0.5 and β1 = 0.2 in equation (59), the 3D and contour plots obtained are shown in figure 2. When we set t = 0, t = 0.5 and t = 1, red, green and blue lines, respectively, are obtained in the 2D plot in figure 2; this figure also represents the periodic shape solution.
Figure 2. Plots of equation (59), when ω1 = 0.8, n = 2, k1 = 0.5 and β1 = 0.2.
We give the plots for equation (76) in figure 3. When we set the parameters in equation (76) as a = 0.2, Γ = 0.1, n = 2, k1 = 0.8, γ1 = 1 and α1 = 0.3, the 3D and contour plots are obtained are shown in figure 3. When we set t = 0, t = 0.5 and t = 1, red, green and blue lines, respectively, are obtained in the 2D plot in figure 3; this figure represents the kink shape solution.
Figure 3. Plots of equation (76), when a = 0.2, Γ = 0.1, n = 2, k1 = 0.8, γ1 = 1 and α1 = 0.3.
We show the plots for equation (82) in figure 4. We set the parameters in equation (82) as a = 0.2, Γ = 0.1, n = 2, k1 = 0.8 and β1 = 1, and the 3D and contour plots obtained are shown in figure 4. When we set t = 0, t = 0 and t = 1, red, green and blue lines, respectively, are obtained in 2D plot shown in figure 4; this figure represents kink shape solution.
Figure 4. Plots of equation (82) when a = 0.2, Γ = 0.1, n = 2, k1 = 0.8 and β1 = 1.

5. Conclusion

In this study, a new form of propagation pulse in optical fibers is studied to obtain the exact solutions. We assumed different values for the arbitrary constants of the given model and obtained two different nonlinear differential equations. Then, we applied suitable wave transformations to the models. Thus, the given models were reduced to nonlinear ordinary differential equations. Then, we applied the two different procedures, the modified simplest equation and the modified Kudryashov procedures, to the reduced equations. The solutions are given by hyperbolic, trigonometric and rational functions. When we compared our methods with the other methods in the literature our methods gave simpler results: our solutions are expressed as different types of functions and our methods are easier to apply to equations. We concluded that the methods discussed in our study are effective and powerful. However, the methods have some disadvantages: while successful for specific types of equations, the modified Kudryashov method may not be applicable to all nonlinear equations. Its effectiveness depends heavily on the specific structure of the equation. We give the 3D, contour and 2D plots for some results. These figures are important for understanding the motion of the wave. The results can be used in various fields such as physics, engineering and optics. These methods can apply to the different equations or researchers could apply different methods to the given models.
1
Wazwaz A M 2005 The tan h method: solitons and periodic solutions for the Dodd-Bullough-Mikhailov and the Tzitzeica-Dodd-Bullough equations Chaos Soliton Fract. 25 55

DOI

2
Hashemi M S, Mirzazadeh M, Ahmad H 2023 A reduction technique to solve the (2+1)-dimensional KdV equations with time local fractional derivatives Opt. Quant. Electron 55 721

DOI

3
Biswas A et al 2017 Cubic–quartic optical solitons in Kerr and power law media Optik 144 357

DOI

4
Zhou Q et al 2016 Dark and singular optical solitons with competing nonlocal nonlinearities Opt. Appl. 46 79

5
Biswas A et al 2018 Optical soliton perturbation for Radhakrishnan-Kundu-Lakshmanan equation with a couple of integration schemes Optik 163 126

DOI

6
Karakoc S B G, Geyikli T, Bashan A 2013 A numerical solution of the modified regularized Long Wave MRLW equation using quartic B splines Mathematics 3 231

7
Raza N, Aslam M R, Rezazadeh H 2019 Analytical study of resonant optical solitons with variable coefficients in Kerr and non-Kerr law media Opt. Quant. Electron 51 59

DOI

8
Hosseini K et al 2020 Investigation of different wave structures to the generalized third-order nonlinear Schrödinger equation Optik 206 164259

DOI

9
Zafar A, Raheel M, Bekir A 2020 Exploring the dark and singular soliton solutions of Biswas-Arshed model with full nonlinear form Optik 204 164133

DOI

10
Kaplan M, Akbulut A 2022 A mathematical analysis of a model involving an integrable equation for wave packet envelope J. Math. 2022 3486780

DOI

11
Rabie W B et al 2023 Investigation of solitons and conservation laws in an inhomogeneous optical BER through a generalized derivative nonlinear Schrödinger equation with quintic nonlinearity Opt. Quant. Electron 55 825

DOI

12
Hashemi M S et al 2024 Analytical treatment with the Nucci reduction technique on the p-forced nonlinear Klein-Gordon equation Opt. Quant. Electron 56 34

DOI

13
Kudryashov N A, Lavrova S F 2023 Painleve test, phase plane analysis and analytical solutions of the Chavy-Waddy-Kolokolnikov Model for the description of bacterial colonies Mathematics 11 3203

DOI

14
Niwas M, Kumar S 2023 Multi-peakons, lumps, and other solitons solutions for the (2+1)-dimensional generalized Benjamin-Ono equation: an inverse (G'/G)-expansion method and real-world applications Nonlinear Dyn. 111 22499

DOI

15
Kumar S et al 2021 Abundant different types of exact soliton solution to the (4+1)-dimensional Fokas and (2+1)-dimensional breaking soliton equations Commun. Theor. Phys. 73 105007

DOI

16
Kumar S, Kumar D, Wazwaz A M 2019 Group invariant solutions of (3+1)-dimensional generalized B-type Kadomstsev Petviashvili equation using optimal system of Lie subalgebra Phys. Scr. 94 065204

DOI

17
Kumar S, Rani S, Mann N 2022 Diverse analytical wave solutions and dynamical behaviors of the new (2+1)-dimensional Sakovich equation emerging in fluid dynamics Eur. Phys. J. Plus 137 1226

DOI

18
Kumar S, Dhiman S K 2024 Exploring cone-shaped solitons, breather, and lump-forms solutions using the lie symmetry method and unifed approach to a coupled breaking soliton model Phys. Scr. 99 025243

DOI

19
Kumar S, Kumar D, Kumar A 2021 Lie symmetry analysis for obtaining the abundant exact solutions, optimal system and dynamics of solitons for a higher-dimensional Fokas equation Chaos Soliton Fract. 142 110507

DOI

20
Kudryashov N A 2019 A generalized model for description of propagation pulses in optical fiber Optik 189 42

DOI

21
Zayed E M E, Alngar M E M 2021 Optical soliton solutions for the generalized Kudryashov equation of propagation pulse in optical BER with power nonlinearities by three integration algorithms Math. Models Methods Appl. Sci. 44 315

DOI

22
Raza N et al 2022 An exploration of novel soliton solutions for propagation of pulses in an optical fiber Opt. Quant. Electron 54 462

DOI

23
Razzaq W, Zafar A, Akbulut A 2022 The Modified simplest equation procedure for conformable time-fractional Boussinesq qquations Int. J. Mod. Phys. B 36 2250095

DOI

24
Chen C, Jiang Y L 2018 Simplest equation method for some time-fractional partial differential equations with conformable derivative Comput. Math. Appl. 75 2978

DOI

25
Hosseini K, Mayeli P, Ansari R 2017 Modified Kudryashov method for solving the conformable time-fractional Klein-Gordon equations with quadratic and cubic nonlinearities Optik 130 737

DOI

26
Akbulut A 2023 Obtaining the soliton type solutions of the conformable time-fractional complex Ginzburg-Landau Equation with Kerr law nonlinearity by using two kinds of Kudryashov methods J. Math. 2023 4741219

DOI

Outlines

/