Some Latest Families of Exact Solutions to Date–Jimbo–Kashiwara–Miwa Equation and Its Stability Analysis

: The present study demonstrates the derivation of new analytical solutions for the Date– Jimbo–Kashiwara–Miwa equation utilizing two distinct methodologies, speciﬁcally the modiﬁed Kudryashov technique and the ( g (cid:48) )-expansion procedure. These innovative concepts employ symbolic computations to provide a dynamic and robust mathematical procedure for addressing a range of nonlinear wave situations. Additionally, a comprehensive stability analysis is performed, and the acquired results are visually represented through graphical representations. A comparison between the discovered solutions and those already found in the literature has also been performed. It is anticipated that the solutions will contribute to the existing literature related to mathematical physics and soliton theory.


Introduction
Many intricate natural phenomena have been described using nonlinear partial differential equations (NPDEs).Many researchers concentrate on this topic, since finding the exact solutions to these equations have enhanced our comprehension of how they function, how they are applied, and how they are created [1,2].For the generalized Schrödinger equation, Hosseini et al. used a modified Jacobi elliptic expansion approach to discover exact solutions [3].Fadhal et al. utilized exponential rational function (ERF) and modified simple equation (MSE) procedures to verify traveling wave solutions of the Sasa-Satsuma equation with beta derivative [4].Zhang et al. obtained new solutions of the (2+1)-dimensional fractional KMM system consisting of Jacobi elliptic functions [5].In the presence of cubicquintic nonlinearity and fourth-order dispersion, Raza et al. were able to obtain novel optical solitons [6].Ismael et al. worked on the modulation instability analysis of the coupled Schrödinger-Boussinesq system with the beta-derivative [7].Zafar et al. handled some solutions to the DNA Peyrard-Bishop equation with the beta-derivative by use of the Kudryashov approach [8].Martinez applied the subequation procedure [9].Hosseini et al. obtained dark soliton solutions of some equations with the beta-derivative [10].Pandir et al. employed the modified exponential function technique [11].Yazgan et al. handled the sine-Gordon expansion method [12,13].Ghanbari and Gomez Aguilar utilized the generalized exponential function procedure [14], Kudryashov employed the simplest equation method to the Chavy-Waddy-Kolokolnikov model [15], Sebogodi et al. applied the symmetry reduction method to (2+1)-dimensional combined potential Kadomtsev-Petviashvili-B-type Kadomtsev-Petviashvili [16], Sebogadi et al. used the ansatz method to obtain the traveling wave solutions of the generalized Chaffee-Infante equation in (1+3) dimensions [17], Podile et al. applied the multiple exp-function technique to the e (2+1)-dimensional Hirota-Satsuma-Ito equation [18], and so on [19][20][21][22][23].
We investigate the Date-Jimbo-Kashiwara-Miwa equation in this work with [24]: This equation is an integrable extension of the Kadomtsev-Petviashvili (KP) hierarchy [25].The existence and uniqueness of the Kadomtsev-Petviashvili (KP) hierarchy solutions were proved in [26].Numerical study of the Kadomtsev-Petviashvili (KP) equation were accurate, and the initial conditions are given in [27].The efficient and well-established two-variable (G'/G, 1/G)-expansion procedure is employed to construct exact dynamical wave solutions for the Date-Jimbo-Kashiwara-Miwa equation.Through this approach, precise wave solutions are systematically generated with high efficiency and accuracy [28].The (2+1)-dimensional Date-Jimbo-Kashiwara-Miwa equation has successfully shown N-soliton waves, fusion solutions, multiple M-lump solutions, and the collision phenomenon between one-M-lump and one or two soliton solutions [29].Entirely novel complex analytical solutions for the governing model have been discovered for the Date-Jimbo-Kashiwara-Miwa equation, employing two potent approaches, namely SGEM and IBSEFM [30].The multisoliton solutions for the Date-Jimbo-Kashiwara-Miwa equation are uncovered through the utilization of the Hirota simple technique, followed by the derivation of M-lump solutions for two types of time-dependent scenarios using the long-wave procedure [31].Analytical soliton solutions for the (2+1)-dimensional Date-Jimbo-Kashiwara-Miwa equation are obtained by employing the extended tanh function procedure and the simple equation technique [32], Adem et al. derived the complexiton solutions with the help of the extended transformed rational function algorithm [33], Wazwaz furnished Painleve analysis to show that the equation is completely integrable in the Painleve sense [34].
The current work is arranged as follows: In Section 2, the main steps of the modified Kudryashov technique and (g )expansion method are provided.In Section 3, we present the construction of precise solutions for the Date-Jimbo-Kashiwara-Miwa equation.Also, a graphical representation by taking parameters as special values and a stability analysis of the obtained solutions are given in Section 4. We present the discussion in Section 5. Some conclusions are shown in Section 6.

Techniques
In the current section, we present the modified Kudryashov technique and (g )-expansion algorithms.To this end, we first discuss the necessary auxiliary information as follows.

Auxiliary Information
For the purpose at hand, we postulate a system of NPDEs with the subsequent configuration: wherein P incorporates u and its diverse partial derivatives with respect to t and classical derivatives with respect to x and y.
The wave transformation can be expressed in the following manner: where k, γ, and ω are constants to be determined.Subsequently, Equation ( 3) is substituted into Equation ( 2), leading to the discovery of a system of equations.This system is then solved to establish the conditions for the parameters and utilize the obtained results.
Consequently, we obtain the following ordinary differential equation (ODE), which will be integrated the possible amount of times with respect to ξ.

The Modified Kudryashov (MK) Technique
In accordance with the methodology, the solutions to Equation ( 4) are represented in the following manner [35][36][37][38]: where i (i = 0, 1, . . ., N) are constants that are determined later, N is calculated by the homogeneous balance principle, and the function ψ(ξ) is given by: where Equation ( 6) satisfies the following ODE: By surrogating the Equation ( 5) into Equation ( 4) while taking into account Equation (7), a system of algebraic equations is derived involving m , a, k, γ, χ, and ω.Ultimately, by solving this system of equations, the exact solutions for Equation ( 2) can be computed.

The (g )-Expansion Procedure
We hypothesize that the solution to Equation (4) can be expressed in the following manner: Here, the constants a i (i = 0, . . ., N) are yet to be determined and are found at a later stage.The value of N is determined using the procedure of homogeneous balance.Here, the function g fulfills the following ODE: The solutions of Equation ( 9) can be expressed as follows: 1: When the discriminant ∆ = 4ac − b 2 is less than zero, the solution is defined as follows: 2: When the discriminant ∆ = 4ac − b 2 is equal to zero, the solution is defined as follows: and 3: When the discriminant ∆ = 4ac − b 2 is greater than zero, the solution is defined as follows: By substituting Equation (8) into Equation ( 4) along with Equation ( 9) and the coefficients of all powers of g i to zero, we attain a system of equations.This obtained system encompasses the variables a i (i = 0, . . ., N), a, b, c, k, γ, and ω.Solving these determining equations enables us to obtain the values of a i and ω.Substituting these derived values into (8), along with Equations ( 10)-( 12), yields all feasible solutions [39-41].

Application of the Methods
In this particular section of the study, our objective is to seek new traveling wave solutions for Equation (1).To accomplish this, we commence by conducting a mathematical analysis of the problem at hand.

Mathematical Analysis
If we substitute Equation (3) into Equation ( 1), the following equation can be derived: Integrating Equation ( 13) with respect to ξ twice and assuming the integral constants are zero, we find setting u = U in Equation ( 14), we acquire the ensuing equation: We obtain the balancing number as N = 2 if we balance U 2 with U .

Application of the MK Procedure
In this subsection, we demonstrate the utilization of the MK procedure.If we take an auxiliary solution for Equation (15) as follows: then substitute the Equation ( 16) in Equation ( 15) without ignoring Equation ( 7), and then collect the coefficients of the ψ i (ξ), we obtain the determining equation system as follows: If we solve the obtained system, we obtain solution families as follows: Case 1: If we substitute the obtained values of the coefficients (17) in Equation ( 16) without ignoring Equation ( 6) and relation u = U, we obtain the following solution: Case 2: If we use values (19) with the same procedure above, we obtain the following exact solution:

Application of the (g )-Expansion Procedure
We organize this subsection to demonstrate the application of the (g )-expansion technique to the given equation.If we take an auxiliary solution for Equation ( 15) as follows: then substitute the Equation ( 21) in Equation ( 15) without ignoring Equation ( 9), and then collect the coefficients of the (g ) i , we obtain the determining equation system as follows: If we solve the obtained system, we obtain solution families as follows: Case 1: Case 2:

Stability Analysis and Figures
The stability property of the solutions are closely related to the momentum in the Hamilton system.Kao and Pasumarthy studied the relation of the stability analysis and Hamiltonian systems in [42].From this point of view, the following formula is given for the Hamiltonian system of the solution: where u(ξ) is the solution of the model; then, we calculate the momentum of the Hamilton system as follows: where σ is the optional constant [43,44].
If we substitute k = 0.1, γ = 0.6, = 0.3, υ = 0.5, y = 1, a = 0.02, b = 0.9, c = 0.2 in Equation ( 23), we obtain ω = 6.475236000.When we consider the solution in the square area of [−2, 2] and perform the necessary operations, we find the condition as follows: According to the result, we can say that our solution is stable for the assumed conditions.
Also, we give the figures of some of the results in this section.Firstly, we give the plots for Equation (18). Figure 1 represents the kinky periodic solitary wave solution for the considered equation when we set the parameters as the following special values: k = 0.1, γ = 0.2, = 0.5, υ = 0.2, y = 1, a = 3, χ = 2 in Equation (18).

Discussion
This paper introduces two different solution types obtained through the use of the modified Kudryashov technique and (g )-expansion technique.These solutions differ from those obtained in previous studies [28,32,45].The simulations and analysis were conducted using the Maple software program.It is important to highlight that the solutions' precision was confirmed by substituting them into the original equation.

Conclusions
Within this publication, we constructed some exact traveling wave solutions for the nonlinear fractional Date-Jimbo-Kashiwara-Miwa equation by using the modified Kudryashov technique and (g )-expansion technique.It is demonstrated that the con- sidered procedures offer practical methods for solving NPDEs in mathematical physics.Comparing the results obtained in this paper, we conclude that our solutions are new.Also, we investigated the shape of the figures for different solutions.Additionally, the exact solutions found in this paper may be highly helpful in many applications of applied mathematics for interpreting certain physical phenomena, such as plasma physics, fusion energy, astrophysics, space studies, etc.