Research article Special Issues

Exact and numerical solutions of the generalized breaking soliton system: Insights into non-linear wave dynamics

  • In this paper, we examined the (2+1)-dimensional generalized breaking soliton system (GBSS), an adaptable framework that accurately describes the three-dimensional, wave-dominated interactions occurring in many non-linear media, i.e., fluids, plasmas, and optical fibers. We used an improved F-expansion technique to generate new families of exact solitonic and periodic wave solutions, significantly enlarging the well-studied solution space and providing insight into the complicated interaction between multi-pulse solitons. Validation of these results and an assessment of their stability were carried out by developing a numerical scheme based on finite difference and undertaking a detailed error and stability analysis, demonstrating unconditional stability across a range of parameter values. The results provide new insights into the interplay of dispersion, non-linearity, and cross-wave coupling in governing soliton formation and energy transport in multidimensional systems. In addition to its theoretical importance, this work can provide valuable practical information on engineering applications such as soliton-based communications and wave control applications in fluid systems. This study offers a new methodology to investigate more complex non-linear wave phenomena by integrating the power of symbolic computation with that of robust numerical verification, opening new opportunities for further developments in soliton-driven technologies.

    Citation: Amer Ahmed, Abdulghani R. Alharbi, Ishak Hashim. Exact and numerical solutions of the generalized breaking soliton system: Insights into non-linear wave dynamics[J]. AIMS Mathematics, 2025, 10(3): 5124-5142. doi: 10.3934/math.2025235

    Related Papers:

    [1] Jun Hu, Jie Wu, Mengzhe Wang . Research on VIKOR group decision making using WOWA operator based on interval Pythagorean triangular fuzzy numbers. AIMS Mathematics, 2023, 8(11): 26237-26259. doi: 10.3934/math.20231338
    [2] Manar A. Alqudah, Artion Kashuri, Pshtiwan Othman Mohammed, Muhammad Raees, Thabet Abdeljawad, Matloob Anwar, Y. S. Hamed . On modified convex interval valued functions and related inclusions via the interval valued generalized fractional integrals in extended interval space. AIMS Mathematics, 2021, 6(5): 4638-4663. doi: 10.3934/math.2021273
    [3] Yanhong Su, Zengtai Gong, Na Qin . Complex interval-value intuitionistic fuzzy sets: Quaternion number representation, correlation coefficient and applications. AIMS Mathematics, 2024, 9(8): 19943-19966. doi: 10.3934/math.2024973
    [4] Tatjana Grbić, Slavica Medić, Nataša Duraković, Sandra Buhmiler, Slaviša Dumnić, Janja Jerebic . Liapounoff type inequality for pseudo-integral of interval-valued function. AIMS Mathematics, 2022, 7(4): 5444-5462. doi: 10.3934/math.2022302
    [5] Mustafa Ekici . On an axiomatization of the grey Banzhaf value. AIMS Mathematics, 2023, 8(12): 30405-30418. doi: 10.3934/math.20231552
    [6] Muhammad Bilal Khan, Muhammad Aslam Noor, Thabet Abdeljawad, Bahaaeldin Abdalla, Ali Althobaiti . Some fuzzy-interval integral inequalities for harmonically convex fuzzy-interval-valued functions. AIMS Mathematics, 2022, 7(1): 349-370. doi: 10.3934/math.2022024
    [7] Waqar Afzal, Mujahid Abbas, Jongsuk Ro, Khalil Hadi Hakami, Hamad Zogan . An analysis of fractional integral calculus and inequalities by means of coordinated center-radius order relations. AIMS Mathematics, 2024, 9(11): 31087-31118. doi: 10.3934/math.20241499
    [8] Le Fu, Jingxuan Chen, Xuanchen Li, Chunfeng Suo . Novel information measures considering the closest crisp set on fuzzy multi-attribute decision making. AIMS Mathematics, 2025, 10(2): 2974-2997. doi: 10.3934/math.2025138
    [9] Li Li, Mengjing Hao . Interval-valued Pythagorean fuzzy entropy and its application to multi-criterion group decision-making. AIMS Mathematics, 2024, 9(5): 12511-12528. doi: 10.3934/math.2024612
    [10] Scala Riccardo, Schimperna Giulio . On the viscous Cahn-Hilliard equation with singular potential and inertial term. AIMS Mathematics, 2016, 1(1): 64-76. doi: 10.3934/Math.2016.1.64
  • In this paper, we examined the (2+1)-dimensional generalized breaking soliton system (GBSS), an adaptable framework that accurately describes the three-dimensional, wave-dominated interactions occurring in many non-linear media, i.e., fluids, plasmas, and optical fibers. We used an improved F-expansion technique to generate new families of exact solitonic and periodic wave solutions, significantly enlarging the well-studied solution space and providing insight into the complicated interaction between multi-pulse solitons. Validation of these results and an assessment of their stability were carried out by developing a numerical scheme based on finite difference and undertaking a detailed error and stability analysis, demonstrating unconditional stability across a range of parameter values. The results provide new insights into the interplay of dispersion, non-linearity, and cross-wave coupling in governing soliton formation and energy transport in multidimensional systems. In addition to its theoretical importance, this work can provide valuable practical information on engineering applications such as soliton-based communications and wave control applications in fluid systems. This study offers a new methodology to investigate more complex non-linear wave phenomena by integrating the power of symbolic computation with that of robust numerical verification, opening new opportunities for further developments in soliton-driven technologies.



    Random numbers are very important in statistical, probability theory, and mathematical analysis in such complex cases, where the real numbers are difficult to record. The random numbers are generated from the uniform distribution when an interval is defined for their selection of random numbers. The random numbers are generated in sequence and depict the behavior of the real data. In addition, random data can be used for estimation and forecasting purposes. According to [1], "The method is based on running the model many times as in random sampling. For each sample, random variates are generated on each input variable; computations are run through the model yielding random outcomes on each output variable. Since each input is random, the outcomes are random. In the same way, they generated thousands of such samples and achieved thousands of outcomes for each output variable. In order to carry out this method, a large stream of random numbers was needed". To generate random numbers, a random generator is applied. The random numbers have no specific pattern and are generated from the chance process. Nowadays, the latest computer can be used to generate random numbers using a well-defined algorithm; see [1]. Bang et al. [2] investigated normality using random-number generating. Schulz et al. [3] presented a pattern-based approach. Tanyer [4] generated random numbers from uniform sampling. Kaya and Tuncer [5] proposed a method to generate biological random numbers. Tanackov et al. [6] presented a method to generate random numbers from the exponential distribution. Jacak et al. [7] presented the methods to generate pseudorandom numbers. More methods can be seen in [8,9,10].

    The neutrosophic statistical distributions were found to be more efficient than the distributions under classical statistics. The neutrosophic distributions can be applied to analyze the data that is given in neutrosophic numbers. Sherwani et al. [11] proposed neutrosophic normal distribution. Duan et al. [12] worked on neutrosophic exponential distribution. Aliev et al. [13] generated Z-random numbers from linear programming. Gao and Ralescu [14] studied the convergence of random numbers generated under an uncertain environment. More information on random numbers generators can be seen in [15,16,17,18]. In recent works, Aslam [19] introduced a truncated variable algorithm for generating random variates from the neutrosophic DUS-Weibull distribution. Additionally, in another study [20], novel methods incorporating sine-cosine and convolution techniques were introduced to generate random numbers within the framework of neutrosophy. Albassam et al. [21] showcased probability/cumulative density function plots and elucidated the characteristics of the neutrosophic Weibull distribution as introduced by [22]. The estimation and application of the neutrosophic Weibull distribution was also presented by [21].

    In [22], the Weibull distribution was introduced within the realm of neutrosophic statistics, offering a more inclusive perspective compared to its traditional counterpart in classical statistics. [21] further examined the properties of the neutrosophic Weibull distribution introduced by [22]. Despite an extensive review of existing literature, no prior research has been identified regarding the development of algorithms for generating random numbers using both the neutrosophic uniform and Weibull distributions. This paper aims to bridge this gap by presenting innovative random number generators tailored specifically for the neutrosophic uniform distribution and the neutrosophic Weibull distribution. The subsequent sections will provide detailed explanations of the algorithms devised to generate random numbers for these distributions. Additionally, the paper will feature multiple tables showcasing sets of random numbers across various degrees of indeterminacy. Upon thorough analysis, the results reveal a noticeable decline in random numbers as the degree of indeterminacy increases.

    Let xNU=xNL+xNUIxNU;IxNUϵ[IxLU,IxUU] be a neutrosophic random variable that follows the neutrosophic uniform distribution. Note that the first part xNL denotes the determinate part, xNUIxNU the indeterminate part, and IxNUϵ[IxLU,IxUU] the degree of indeterminacy. Suppose f(xNU)=f(xLU)+f(xUU)INU;INUϵ[ILU,IUU] presents the neutrosophic probability density function (npdf) of neutrosophic uniform distribution (NUD). Note that the npdf of NUD is based on two parts. The first part xNL, f(xLU) denotes the determinate part and presents the probability density function (pdf) of uniform distribution under classical statistics. The second part xNUIxNU, f(xUU)INU denotes the indeterminate part and IxNUϵ[IxLU,IxUU], INUϵ[ILU,IUU] are the measures of indeterminacy associated with neutrosophic random variable and the uniform distribution. The npdf of the uniform distribution by following [22] is given as

    f(xNU)=(1(bLaL))+(1(bUaU))IxNU;IxNUϵ[IxLU,IxUU],aNxNUbN, (1)

    where bNϵ[bL,bU] and aNϵ[aL,aU] are neutrosophic parameters of the NUD. The simplified form when L=U=SU of Eq (1) can be written as

    f(xNSU)=(1(bNSaNS))(1+IxNS);IxNSϵ[IxLS,IxUS],aNxNUbN. (2)

    Note here that the npdf of uniform distribution is a generalization of pdf of the uniform distribution. The neutrosophic uniform distribution reduces to the classical uniform distribution when IxUU = 0. The neutrosophic cumulative distribution function (ncdf) of the neutrosophic uniform distribution is given by

    F(xNU)=(xNLaL(bLaL))+(xNUaU(bUaU))IxNU;IxNUϵ[IxLU,IxUU],aNxNUbN. (3)

    Note that the first part presents the cumulative distribution function (cdf) of the uniform distribution under classical statistics, and the second part is the indeterminate part associated with ncdf. The ncdf reduces to cdf when IxUU = 0. The simplified form of ncdf of the Uniform distribution when L=U=S can be written as

    F(xNSU)=(xNSaNS(bNSaNS))(1+INS);INSϵ[ILS,IUS],aNxNUbN. (4)

    Aslam [22] introduced the neutrosophic Weibull distribution (NWD) originally. The neutrosophic form of the Weibull distribution is expressed by

    f(xNW)=f(xLW)+f(xUW)INW;INWϵ[ILW,IUW]. (5)

    The following npdf of the Weibull distribution is taken from [22] and reported as

    f(xNW)={(βα)(xLα)β1e(xLα)β}+{(βα)(xUα)β1e(xUα)β}INW;INWϵ[ILW,IUW]. (6)

    The simplified form of the npdf of the Weibull distribution when L=U=SW is expressed by

    f(xNSW)={(βα)(xSα)β1e(xSα)β}(1+INS);INSϵ[ILS,IUS], (7)

    where α and β are the scale and shape parameters of the Weibull distribution. The npdf of the Weibull distribution reduces to pdf of the Weibull distribution when INS=0. The ncdf of the Weibull distribution is expressed by

    F(xNSW)=1{e(xNSWα)β(1+INW)}+INW;INWϵ[ILW,IUW]. (8)

    The ncdf of the Weibull distribution reduces to cdf of the Weibull distribution under classical statistics when INW = 0. The neutrosophic mean of the Weibull distribution is given as [22]

    μNW=αΓ(1+1/β)(1+INW);INWϵ[ILW,IUW]. (9)

    The neutrosophic median of the Weibull distribution is given by

    ˜μNW=α(ln(2))1/β(1+INW);INWϵ[ILW,IUW]. (10)

    This section presents the methodology to generate random variates from the proposed neutrosophic uniform distribution and the neutrosophic Weibull distribution. Let uNϵ[uL,uU] be a neutrosophic random uniform from uNUN([0,0],[1,1]). The neutrosophic random numbers from NUD and NWD will be obtained as follows:

    Let

    uN=F(xNU)=(xNLaL(bLaL))+(xNUaU(bUaU))IxNU;IxNUϵ[IxLU,IxUU],aNxNUbN,

    or

    uN=F(xNU)=(xNSaNS(bNSaNS))(1+INS);INSϵ[ILS,IUS],aNxNUbN.

    The neutrosophic random numbers xNSU from NWD can be obtained using the following Eq (11)

    xNSU=aNS+(uN(1+INS))(bNSaNS);uNϵ[uL,uU],INSϵ[ILS,IUS]. (11)

    The random number from the Weibull distribution using classical statistics can be obtained when INS = 0 using the following Eq (12)

    x=a+u(ba);axb. (12)

    The neutrosophic random numbers from the NWD will be obtained using the following methodology.

    Let

    uN=F(xNSW)=1{e(xNSWα)β(1+INW)}+INW;INWϵ[ILW,IUW],uNϵ[uL,uU].

    The neutrosophic random numbers from NWD can be obtained through the following expression

    xNSW=α[ln(1(uNINW)1+INW)]1β;INWϵ[ILW,IUW],uNϵ[uL,uU]. (13)

    The NWD reduces to neutrosophic exponential distribution (NED) when β=1. The neutrosophic random numbers from the NED can be obtained as follows:

    xNSE=αln(1(uNINW)1+INW);INWϵ[ILW,IUW],uNϵ[uL,uU]. (14)

    The random numbers from the Weibull distribution using classical statistics can be obtained as

    xNSW=αln(1u)1β. (15)

    The random numbers from the exponential distribution using classical statistics can be obtained as

    xNSW=αln(1u). (16)

    The following routine can be run to generate n random numbers from the NUD.

    Step-1: Generate a uniform random number uN from uNUN([0,0],[1,1]).

    Step-2: Fix the values of INS.

    Step-3: Generate values of xNSU using the expression

    xNSU=aNS+(uN(1+INS))(bNSaNS);uNϵ[uL,uU],INSϵ[ILS,IUS].

    Step-4: From the routine, the first value of xNSU will be generated.

    Step-5: Repeat the routine k times to generate k random numbers from NUD.

    The following routine can be run to generate n random numbers from the NUD.

    Step-1: Generate a uniform random number uN from uNUN([0,0],[1,1]).

    Step-2: Fix the values of INS, α and β.

    Step-3: Generate values of xNSW using the expression

    xNSW=α[ln(1(uNINW)1+INW)]1β;INWϵ[ILW,IUW],uNϵ[uL,uU].

    Step-4: From the routine, the first value of xNSW will be generated.

    Step-5: Repeat the routine k times to generate k random numbers from NWD.

    To illustrate the proposed simulation methods, two examples will be discussed in this section.

    Suppose that xNSU is a neutrosophic uniform random variable with parameters ([20,20],[30,30]) and a random variate xNSU under indeterminacy is needed. To generate a random number from NUD, the following steps have been carried out.

    Step-1: Generate a uniform random number uN=0.05 from uNUN([0,0],[1,1]).

    Step-2: Fix the values of INS=0.1.

    Step-3: Generate values of xNSU using the expression xNSU=20+(0.05(1+0.1))(3020)=20.5.

    Step-4: From the routine, the first value of xNSU=20.5 will be generated.

    Step-5: Repeat the routine k times to generate k random numbers from NUD.

    Step-1: Generate a uniform random number uN=0.30 from uNUN([0,0],[1,1]).

    Step-2: Fix the values of INS=0.20, α=5, and β=0.5.

    Step-3: Generate values of xNSW using the expression xNSW=5[ln(1(uNINW)1+INW)]10.5=0.04.

    Step-4: From the routine, the first value of xNSW=0.04 will be generated.

    Step-5: Repeat the routine k times to generate k random numbers from NWD.

    In this section, random numbers are generated by simulation using the above-mentioned algorithms for NUD and NWD. To generate random numbers from NUD, several uniform numbers are generated from uNUN([0,0],[1,1]) and placed in Tables 1 and 2. In Tables 1 and 2, several values of INS are considered to generate random numbers from the NUD. Table 1 is depicted by assuming that NUD has the parameters aNS = 10 and bNS = 20 and Table 2 is shown by assuming that NUD has the parameters aNS = 20 and bNS = 30. From Tables 1 and 2, the following trends can be noted in random numbers generated from NUD.

    Table 1.  Random numbers from NUD when aNS = 10 and bNS = 20.
    u INS
    0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1
    0.05 10.5 10.45 10.42 10.38 10.36 10.33 10.31 10.29 10.28 10.26 10.3 10.2
    0.1 11 10.91 10.83 10.77 10.71 10.67 10.63 10.59 10.56 10.53 10.5 10.5
    0.15 11.5 11.36 11.25 11.15 11.07 11.00 10.94 10.88 10.83 10.79 10.8 10.7
    0.2 12 11.82 11.67 11.54 11.43 11.33 11.25 11.18 11.11 11.05 11.0 11.0
    0.25 12.5 12.27 12.08 11.92 11.79 11.67 11.56 11.47 11.39 11.32 11.3 11.2
    0.3 13 12.73 12.50 12.31 12.14 12.00 11.88 11.76 11.67 11.58 11.5 11.4
    0.35 13.5 13.18 12.92 12.69 12.50 12.33 12.19 12.06 11.94 11.84 11.8 11.7
    0.4 14 13.64 13.33 13.08 12.86 12.67 12.50 12.35 12.22 12.11 12.0 11.9
    0.45 14.5 14.09 13.75 13.46 13.21 13.00 12.81 12.65 12.50 12.37 12.3 12.1
    0.5 15 14.55 14.17 13.85 13.57 13.33 13.13 12.94 12.78 12.63 12.5 12.4
    0.55 15.5 15.00 14.58 14.23 13.93 13.67 13.44 13.24 13.06 12.89 12.8 12.6
    0.6 16 15.45 15.00 14.62 14.29 14.00 13.75 13.53 13.33 13.16 13.0 12.9
    0.65 16.5 15.91 15.42 15.00 14.64 14.33 14.06 13.82 13.61 13.42 13.3 13.1
    0.7 17 16.36 15.83 15.38 15.00 14.67 14.38 14.12 13.89 13.68 13.5 13.3
    0.75 17.5 16.82 16.25 15.77 15.36 15.00 14.69 14.41 14.17 13.95 13.8 13.6
    0.8 18 17.27 16.67 16.15 15.71 15.33 15.00 14.71 14.44 14.21 14.0 13.8
    0.9 19 18.18 17.50 16.92 16.43 16.00 15.63 15.29 15.00 14.74 14.5 14.3
    0.95 19.5 18.64 17.92 17.31 16.79 16.33 15.94 15.59 15.28 15.00 14.8 14.5

     | Show Table
    DownLoad: CSV
    Table 2.  Random numbers from NUD when aNS = 20 and bNS = 30.
    u INS
    0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1
    0.05 20.5 20.5 20.4 20.4 20.4 20.3 20.3 20.3 20.3 20.3 20.3 20.2
    0.1 21.0 20.9 20.8 20.8 20.7 20.7 20.6 20.6 20.6 20.5 20.5 20.5
    0.15 21.5 21.4 21.3 21.2 21.1 21.0 20.9 20.9 20.8 20.8 20.8 20.7
    0.2 22.0 21.8 21.7 21.5 21.4 21.3 21.3 21.2 21.1 21.1 21.0 21.0
    0.25 22.5 22.3 22.1 21.9 21.8 21.7 21.6 21.5 21.4 21.3 21.3 21.2
    0.3 23.0 22.7 22.5 22.3 22.1 22.0 21.9 21.8 21.7 21.6 21.5 21.4
    0.35 23.5 23.2 22.9 22.7 22.5 22.3 22.2 22.1 21.9 21.8 21.8 21.7
    0.4 24.0 23.6 23.3 23.1 22.9 22.7 22.5 22.4 22.2 22.1 22.0 21.9
    0.45 24.5 24.1 23.8 23.5 23.2 23.0 22.8 22.6 22.5 22.4 22.3 22.1
    0.5 25.0 24.5 24.2 23.8 23.6 23.3 23.1 22.9 22.8 22.6 22.5 22.4
    0.55 25.5 25.0 24.6 24.2 23.9 23.7 23.4 23.2 23.1 22.9 22.8 22.6
    0.6 26.0 25.5 25.0 24.6 24.3 24.0 23.8 23.5 23.3 23.2 23.0 22.9
    0.65 26.5 25.9 25.4 25.0 24.6 24.3 24.1 23.8 23.6 23.4 23.3 23.1
    0.7 27.0 26.4 25.8 25.4 25.0 24.7 24.4 24.1 23.9 23.7 23.5 23.3
    0.75 27.5 26.8 26.3 25.8 25.4 25.0 24.7 24.4 24.2 23.9 23.8 23.6
    0.8 28.0 27.3 26.7 26.2 25.7 25.3 25.0 24.7 24.4 24.2 24.0 23.8
    0.9 29.0 28.2 27.5 26.9 26.4 26.0 25.6 25.3 25.0 24.7 24.5 24.3
    0.95 29.5 28.6 27.9 27.3 26.8 26.3 25.9 25.6 25.3 25.0 24.8 24.5

     | Show Table
    DownLoad: CSV

    1) For fixed INS, aNS = 10 and bNS = 20, as the values of u increase from 0.05 to 0.95, there is an increasing trend in random numbers.

    2) For fixed u, aNS = 10 and bNS = 20, as the values of INS increase from 0 to 1.1, there is a decreasing trend in random numbers.

    3) For fixed values of u and INS, as the values of aNS and bNS increases, there is an increasing trend in random numbers.

    The random numbers for NWD are generated using the algorithm discussed in the last section. The random numbers for various values of u, INS, α, and β are considered. The random numbers when α=5 and β=0 are shown in Table 3. The random numbers when α=5 and β=1 are shown in Table 4. The random numbers when α=5 and β=2 are shown in Table 5.

    Table 3.  Random numbers from NUD when α=5 and β=0.5.
    u INS
    0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1
    0.05 0.01 0.01 0.07 0.15 0.23 0.31 0.38 0.43 0.48 0.52 0.56 0.58
    0.1 0.06 0.00 0.03 0.10 0.18 0.25 0.32 0.38 0.43 0.48 0.51 0.54
    0.15 0.13 0.01 0.01 0.06 0.13 0.20 0.27 0.33 0.39 0.43 0.47 0.51
    0.2 0.25 0.05 0.00 0.03 0.08 0.15 0.22 0.28 0.34 0.39 0.43 0.47
    0.25 0.41 0.11 0.01 0.01 0.05 0.11 0.18 0.24 0.30 0.35 0.39 0.43
    0.3 0.64 0.21 0.04 0.00 0.02 0.07 0.13 0.20 0.25 0.31 0.35 0.39
    0.35 0.93 0.34 0.09 0.01 0.01 0.04 0.10 0.16 0.21 0.27 0.31 0.36
    0.4 1.30 0.53 0.17 0.03 0.00 0.02 0.06 0.12 0.17 0.23 0.28 0.32
    0.45 1.79 0.77 0.29 0.08 0.01 0.01 0.04 0.09 0.14 0.19 0.24 0.28
    0.5 2.40 1.08 0.44 0.15 0.03 0.00 0.02 0.06 0.11 0.16 0.21 0.25
    0.55 3.19 1.48 0.64 0.24 0.07 0.01 0.00 0.03 0.08 0.12 0.17 0.22
    0.6 4.20 1.99 0.91 0.38 0.13 0.02 0.00 0.02 0.05 0.10 0.14 0.19
    0.65 5.51 2.63 1.24 0.55 0.21 0.06 0.01 0.00 0.03 0.07 0.11 0.16
    0.7 7.25 3.47 1.67 0.77 0.32 0.11 0.02 0.00 0.01 0.05 0.09 0.13
    0.75 9.61 4.55 2.21 1.06 0.47 0.18 0.05 0.00 0.00 0.03 0.06 0.10
    0.8 12.95 5.99 2.92 1.42 0.67 0.28 0.10 0.02 0.00 0.01 0.04 0.08
    0.9 26.51 10.70 5.03 2.48 1.23 0.58 0.25 0.09 0.02 0.00 0.01 0.04
    0.95 44.87 14.87 6.67 3.26 1.63 0.79 0.36 0.14 0.04 0.00 0.00 0.02

     | Show Table
    DownLoad: CSV
    Table 4.  Random numbers from NUD when α=5 and β=1.
    u INS
    0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1
    0.05 0.26 - - - - - - - - - - -
    0.1 0.53 0.00 - - - - - - - - - -
    0.15 0.81 0.23 - - - - - - - - - -
    0.2 1.12 0.48 0.00 - - - - - - - - -
    0.25 1.44 0.74 0.21 - - - - - - - - -
    0.3 1.78 1.01 0.44 0.00 - - - - - - - -
    0.35 2.15 1.31 0.68 0.20 - - - - - - - -
    0.4 2.55 1.62 0.93 0.41 0.00 - - - - - - -
    0.45 2.99 1.96 1.20 0.63 0.18 - - - - - - -
    0.5 3.47 2.32 1.49 0.86 0.38 0.00 - - - - - -
    0.55 3.99 2.72 1.79 1.11 0.58 0.17 - - - - - -
    0.6 4.58 3.15 2.13 1.37 0.80 0.35 0.00 - - - - -
    0.65 5.25 3.63 2.49 1.66 1.03 0.54 0.16 - - - - -
    0.7 6.02 4.16 2.89 1.96 1.27 0.74 0.33 0.00 - - - -
    0.75 6.93 4.77 3.33 2.30 1.54 0.96 0.51 0.15 - - - -
    0.8 8.05 5.47 3.82 2.67 1.82 1.19 0.70 0.31 0.00 - - -
    0.9 11.51 7.32 5.02 3.52 2.48 1.70 1.11 0.66 0.29 0.00 - -
    0.95 14.98 8.62 5.78 4.04 2.85 1.99 1.35 0.85 0.45 0.13 - -

     | Show Table
    DownLoad: CSV
    Table 5.  Random numbers from NUD when α=5 and β=2.
    u INS
    0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1
    0.05 1.13 - - - - - - - - - - -
    0.1 1.62 0.00 - - - - - - - - - -
    0.15 2.02 1.08 - - - - - - - - - -
    0.2 2.36 1.55 0.00 - - - - - - - - -
    0.25 2.68 1.92 1.03 - - - - - - - - -
    0.3 2.99 2.25 1.48 0.00 - - - - - - - -
    0.35 3.28 2.56 1.84 0.99 - - - - - - - -
    0.4 3.57 2.85 2.16 1.42 0.00 - - - - - - -
    0.45 3.87 3.13 2.45 1.77 0.96 - - - - - - -
    0.5 4.16 3.41 2.73 2.07 1.37 0.00 - - - - - -
    0.55 4.47 3.69 3.00 2.35 1.70 0.92 - - - - - -
    0.6 4.79 3.97 3.26 2.62 2.00 1.33 0.00 - - - - -
    0.65 5.12 4.26 3.53 2.88 2.27 1.65 0.90 - - - - -
    0.7 5.49 4.56 3.80 3.13 2.52 1.93 1.28 0.00 - - - -
    0.75 5.89 4.88 4.08 3.39 2.77 2.19 1.59 0.87 - - - -
    0.8 6.34 5.23 4.37 3.65 3.02 2.44 1.87 1.24 0.00 - - -
    0.9 7.59 6.05 5.01 4.20 3.52 2.92 2.36 1.81 1.21 0.00 - -
    0.95 8.65 6.57 5.37 4.49 3.78 3.16 2.59 2.06 1.50 0.82 - -

     | Show Table
    DownLoad: CSV

    From Tables 35, the following trends can be noted in random numbers generated from NUD.

    1) For fixed INS, α=5 and β=0.5, as the values of u increase from 0.05 to 0.95, there is an increasing trend in random numbers generated from NWD.

    2) For fixed u, α=5, and β=0.5, as the values of INS increase from 0 to 1.1, there is an increasing trend in random numbers.

    3) For fixed values of INS and α, as the values of β increase, there is an increasing trend in random numbers.

    The algorithms to generate the random variables from NUD and NWD are depicted in Figures 1 and 2.

    Figure 1.  Algorithm to generate random numbers from NUD.
    Figure 2.  Algorithm to generate random numbers from NWD.

    In this section, the performance of simulations using classical simulation and neutrosophic simulation will be discussed using the random numbers from the NUD and the NWD distribution. As explained earlier, the proposed simulation method under neutrosophy will be reduced to the classical simulation method under classical statistics when no uncertainty is found in the data. To study the behavior of random numbers, random numbers from NUD when INS=1.1, aNS = 20, and bNS = 30 are considered and depicted in Figure 3. In Figure 3, it can be seen that the curve of random numbers from the classical simulation is higher than the curve of random numbers from the neutrosophic simulation. From Figure 3, it is clear that the proposed neutrosophic simulation method gives smaller values of random numbers than the random numbers generated by the neutrosophic simulation method. The random numbers from NWD when INS=0.9, α=5, and β=0.5 are considered and their curves are shown in Figure 4. From Figure 4, it can be seen that random numbers generated by neutrosophic simulation are smaller than the random numbers generated by the classical simulation method under classical statistics. The random numbers generated by the neutrosophic simulation are close to zero. The random numbers from NWD when INS=0.1, α=5, and β=1 (exponential distribution) are considered and their curves are shown in Figure 5. From Figure 5, it can be seen that the curve of random numbers generated by neutrosophic simulation is lower than the curve of random numbers generated by the classical simulation method under classical statistics. The random numbers from NWD when INS=0.1, α=5, and β=2 are considered and their curves are shown in Figure 6. From Figure 6, it can be seen that the curve of random numbers generated by neutrosophic simulation is lower than the curve of random numbers generated by the classical simulation method under classical statistics. From Figures 46, it can be concluded that the proposed simulation gives smaller values of random numbers as compared to the classical simulation method under classical statistics.

    Figure 3.  Random numbers behavior from NUD when INS=1.1, aNS = 20, and bNS = 30.
    Figure 4.  Random numbers behavior from NWD when INS=0.9, and when α=5, and β=0.5.
    Figure 5.  Random numbers behavior from NWD when INS=0.1, and when α=5, and β=1.
    Figure 6.  Random numbers behavior from NWD when INS=0.1, and when α=5, and β=2.

    The simulation method under neutrosophic statistics and classical methods was discussed in the last sections. From Tables 1 and 2, it can be seen that random numbers from the NUD can be generated when INS<1, INS=1 and INS>1. On the other hand, the random numbers from the NWD can be generated for INS<1, INS=1, and INS>1 when the shape parameter β<1. From Table 4 and 5, it can be noted that for several cases, the NWD generates negative results or random numbers do not exist. Based on the simulation studies, it can be concluded that the NWD generates random numbers INS<1, INS=1, and INS>1 only when β<1. To generate random numbers from NWD when β1, the following expression will be used

    xNSW=α[ln(1uN+INW1+INW)]1β;1uN+INW0.

    In this paper, we initially introduced the NUD and presented a novel method for generating random numbers from both NUD and the NWD. We also introduced algorithms for generating random numbers within the context of neutrosophy. These algorithms were applied to generate random numbers from both distributions using various parameters. We conducted an extensive discussion on the behavior of these random numbers, observing that random numbers generated under neutrosophy tend to be smaller than those generated under uncertain environments. It is worth noting that generating random numbers from computers is a common practice. Tables 15 within this paper offer valuable insights into how the degree of determinacy influences random number generation. Additionally, these tables can be utilized for simulation purposes in fields marked by uncertainty, such as reliability, environmental studies, and medical science. From our study, we conclude that the proposed method for generating random numbers from NUD and NWD can be effectively applied in complex scenarios. In future research, exploring the statistical properties of the proposed NUD would be advantageous. Additionally, investigating the proposed algorithm utilizing the accept-reject method could be pursued as a future research avenue. Moreover, there is potential to develop algorithms using other statistical distributions for further investigation.

    The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

    The authors are deeply thankful to the editor and reviewers for their valuable suggestions to improve the quality and presentation of the paper.

    The authors declare no conflicts of interest.



    [1] A. Seadawy, I. Mujahid, Optical soliton solutions for non-linear complex Ginzburg–Landau dynamical equation with laws of nonlinearity Kerr law media, Int. J. Mod. Phys. B, 34 (2020), 2050179. https://doi.org/10.1142/S0217979220501799 doi: 10.1142/S0217979220501799
    [2] A. Wazwaz, Gaussian solitary waves for the logarithmic Boussinesq equation and the logarithmic regularised Boussinesq equation, Ocean Eng., 94 (2015), 111–115. https://doi.org/10.1016/j.oceaneng.2014.11.024 doi: 10.1016/j.oceaneng.2014.11.024
    [3] B. Li, Y. Ma, T. Yang, Stable optical soliton in the ring-cavity fiber system with carbon nanotube as saturable absorber, Superlattice. Microst., 113 (2018), 366–372. https://doi.org/10.1016/j.spmi.2017.11.016 doi: 10.1016/j.spmi.2017.11.016
    [4] J. Yang, Y. Gao, C. Su, D Zuo, Y. Feng, Solitons and quasi-periodic behaviors in an inhomogeneous optical fiber, Commun. Nonlinear Sci., 42 (2017), 477–490. https://doi.org/10.1016/j.cnsns.2016.05.013 doi: 10.1016/j.cnsns.2016.05.013
    [5] H. Ismael, T. Sulaiman, H. Nabi, W. Mahmoud, M. Osman, Geometrical patterns of time variable Kadomtsev–Petviashvili (Ⅰ) equation that models dynamics of waves in thin films with high surface tension, Nonlinear Dyn., 111 (2023), 9457–9466. https://doi.org/10.1007/s11071-023-08319-8 doi: 10.1007/s11071-023-08319-8
    [6] M. Zhang, Y. Ma, B. Li, Novel loop-like solitons for the generalized Vakhnenko equation, Chinese Phys. B, 22 (2013), 030511. https://dx.doi.org/10.1088/1674-1056/22/3/030511 doi: 10.1088/1674-1056/22/3/030511
    [7] A. R. Seadawy, Stability analysis for Zakharov–Kuznetsov equation of weakly non-linear ion-acoustic waves in a plasma, Comput. Math. Appl., 67 (2014), 172–180. https://doi.org/10.1016/j.camwa.2013.11.001 doi: 10.1016/j.camwa.2013.11.001
    [8] Y. Ma, B. Li, The wrinkle-like N-solitons for the thermophoretic motion equation through graphene sheets, Physica A, 494 (2018), 169–174. https://doi.org/10.1016/j.physa.2017.12.014 doi: 10.1016/j.physa.2017.12.014
    [9] B. Ren, P. Chu, Dynamics of D'Alembert wave and soliton molecule for a (2+1)-dimensional generalized breaking soliton equation, Chinese J. Phys., 74 (2021), 296–301. https://doi.org/10.1016/j.cjph.2021.07.025 doi: 10.1016/j.cjph.2021.07.025
    [10] B. Ren, J. Lin, Z. Lou, Consistent Riccati expansion and rational solutions of the Drinfel'd–Sokolov–Wilson equation, Appl. Math. Lett., 105 (2020), 106326. https://doi.org/10.1016/j.aml.2020.106326 doi: 10.1016/j.aml.2020.106326
    [11] M. Remoissenet, Basic concepts and the discovery of solitons. Waves Called Solitons: Concepts and Experiments, Heidelberg: Springer Berlin Heidelberg, 1999. https://doi.org/10.1007/978-3-662-03790-4_1
    [12] C. Becker, S. Stellmer, P. Soltan-Panahi, S. Dörscher, M. Baumert, E. Richter, et al., Oscillations and interactions of dark and dark–bright solitons in Bose–Einstein condensates, Nat. Phys., 4 (2008), 496–501. https://doi.org/10.1038/nphys962 doi: 10.1038/nphys962
    [13] A. Slunyaev, Persistence of hydrodynamic envelope solitons: Detection and rogue wave occurrence, Phys. Fluids, 33 (2021), 036606. https://doi.org/10.1063/5.0042232 doi: 10.1063/5.0042232
    [14] G. Dematteis, T. Grafke, E. Eijnden, Rogue waves and large deviations in deep sea, P. Natl. Acad. Sci., 115 (2018), 855–860. https://www.pnas.org/doi/abs/10.1073/pnas.1710670115
    [15] A. Alharbi, Traveling-wave and numerical solutions to non-linear evolution equations via modern computational techniques, AIMS Math., 9 (2024), 1323–1345. https://doi.org/10.3934/math.2024065 doi: 10.3934/math.2024065
    [16] S. Residori, M. Onorato, U. Bortolozzo, F. Arecchi, Rogue waves: A unique approach to multidisciplinary physics, Contemp. Phys., 58 (2017), 53-–69. https://doi.org/10.1080/00107514.2016.1243351 doi: 10.1080/00107514.2016.1243351
    [17] A. Demircan, S. Amiranashvili, C. Brée, C. Mahnke, F. Mitschke, G. Steinmeyer, Rogue events in the group velocity horizon, Sci. Rep., 2 (2012), 850. https://doi.org/10.1038/srep00850 doi: 10.1038/srep00850
    [18] G. Xu, Integrability of a (2+1)-dimensional generalized breaking soliton equation, Appl. Math. Lett., 50 (2015), 16–22. https://doi.org/10.1016/j.aml.2015.05.015 doi: 10.1016/j.aml.2015.05.015
    [19] L. Ling, L. Zhao, B. Guo, Darboux transformation and classification of solution for mixed coupled non-linear Schrödinger equations, Commun. Nonlinear Sci., 32 (2016), 285–304. https://doi.org/10.1016/j.cnsns.2015.08.023 doi: 10.1016/j.cnsns.2015.08.023
    [20] S. Jia, Y. Gao, C. Zhao, Z. L. Feng, Solitons, breathers and rogue waves for a sixth-order variable-coefficient non-linear Schrödinger equation in an ocean or optical fiber, Eur. Phys. J. Plus, 132 (2017), 34. https://doi.org/10.1140/epjp/i2017-11318-y doi: 10.1140/epjp/i2017-11318-y
    [21] R. Hirota, Exact solution of the modified Korteweg-de Vries equation for multiple collisions of solitons, J. Phys. Soc. Japan, 33 (1972), 1456–1458. https://doi.org/10.1143/JPSJ.33.1456 doi: 10.1143/JPSJ.33.1456
    [22] B. Ghanbari, K. Nisar, Determining new soliton solutions for a generalized non-linear evolution equation using an effective analytical method, Alex. Eng. J., 59 (2020), 3171–3179. https://doi.org/10.1016/j.aej.2020.07.032 doi: 10.1016/j.aej.2020.07.032
    [23] A. Alharbi, Numerical solutions to two-dimensional fourth order parabolic thin film equations using the Parabolic Monge-Ampere method, AIMS Math., 8 (2023), 16463–16478. https://doi.org/10.3934/math.2023841 doi: 10.3934/math.2023841
    [24] A. Alharbi, Traveling-wave and numerical solutions to a Novikov-Veselov system via the modified mathematical methods, AIMS Math., 8 (2023), 1230–1250. https://doi.org/10.3934/math.2023062 doi: 10.3934/math.2023062
    [25] M. Iqbal, A. Seadawy, S. Althobaiti, Mixed soliton solutions for the (2+1)-dimensional generalized breaking soliton system via new analytical mathematical method, Results Phys., 32 (2022), 105030. https://doi.org/10.1016/j.rinp.2021.105030 doi: 10.1016/j.rinp.2021.105030
    [26] M. Jawad, A. Ja'afar, Soliton solutions for nonlinear systems (2+1)-dimensional equations, IOSR J. Math., 1 (2012), 27–34. http://dx.doi.org/10.9790/5728-0162734 doi: 10.9790/5728-0162734
    [27] W. Hereman, M. Takaoka, Solitary wave solutions of non-linear evolution and wave equations using a direct method and MACSYMA, J. Phys. A-Math. Gen., 23 (1990), 4805. https://dx.doi.org/10.1088/0305-4470/23/21/021 doi: 10.1088/0305-4470/23/21/021
    [28] T. Alharbi, A. Alharbi, Traveling-wave and numerical investigations to non-linear equations via modern computational techniques, AIMS Math., 9 (2024), 12188–12210. https://doi.org/10.3934/math.2024595 doi: 10.3934/math.2024595
    [29] T. Alharbi, A. Alharbi, A study of traveling wave structures and numerical investigations into the coupled nonlinear schrödinger equation using advanced mathematical techniques, Mathematics, 11 (2023), 4597. https://doi.org/10.3390/math11224597 doi: 10.3390/math11224597
    [30] M. Ahmed, A. Zaghrout, H. Ahmed, Construction of solitons and other solutions for NLSE with Kudryashov's generalized non-linear refractive index, Alex. Eng. J., 64 (2023), 391–397. https://doi.org/10.1016/j.aej.2022.09.015 doi: 10.1016/j.aej.2022.09.015
    [31] I. Samir, E. Salah, E. A. El-Dahab, H. M. Ahmed, M. Ammar, W. Alexan, et al., Unraveling solitons dynamics in system of dispersive NLSE with Kudryashov's law of nonlinearity using improved modified extended tanh function method, Alex. Eng. J., 91 (2024), 419–428. https://doi.org/10.1016/j.aej.2024.02.020 doi: 10.1016/j.aej.2024.02.020
    [32] R. Attia, H. Qazi, A. Kamran, A. Jamshad, Z. Aniqa, Soliton solutions of nonlinear evolution equations by basic (G'/G)-expansion method, Math. Model. Eng. Probl., 7 (2020), 242–-250. http://dx.doi.org/10.18280/mmep.070210 doi: 10.18280/mmep.070210
    [33] F. Calogero, A. Degasperis, Non-linear evolution equations solvable by the inverse spectral transform.—Ⅰ, Nuovo Cimento B, 32 (1976), 201–242. https://doi.org/10.1007/BF02727634 doi: 10.1007/BF02727634
    [34] A. Wazwaz, A new integrable (2+1)-dimensional generalized breaking soliton equation: N-soliton solutions and traveling wave solutions, Commun. Theor. Phys., 66 (2016), 385. https://dx.doi.org/10.1088/0253-6102/66/4/385 doi: 10.1088/0253-6102/66/4/385
    [35] Y. Ma, B. Li, Interactions between soliton and rogue wave for a (2+1)-dimensional generalized breaking soliton system: Hidden rogue wave and hidden soliton, Comput. Math. Appl., 78 (2019), 827–839. https://doi.org/10.1016/j.camwa.2019.03.002 doi: 10.1016/j.camwa.2019.03.002
    [36] T. Taklo, W. Choi, Group resonant interactions between surface and internal gravity waves in a two-layer system, J. Fluid Mech., 892 (2020), A14. https://doi.org/10.1017/jfm.2020.180 doi: 10.1017/jfm.2020.180
    [37] V. Kruglov, J. Harvey, Solitary waves in optical fibers governed by higher-order dispersion, Phys. Rev. A, 98 (2018), 063811. https://link.aps.org/doi/10.1103/PhysRevA.98.063811
    [38] M. Abdou, The extended F-expansion method and its application for a class of non-linear evolution equations, Chaos Soliton. Fract., 31 (2007), 95–104. https://doi.org/10.1016/j.chaos.2005.09.030 doi: 10.1016/j.chaos.2005.09.030
    [39] L. Shampine, M. Reichelt, The MATLAB ODE suite, SIAM J. Sci. Comput., 18 (1997), 1–22. https://doi.org/10.1137/S1064827594276424 doi: 10.1137/S1064827594276424
  • Reader Comments
  • © 2025 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(431) PDF downloads(60) Cited by(1)

Figures and Tables

Figures(8)  /  Tables(2)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog