Novel Kinds of Fractional λ –Kinetic Equations Involving the Generalized Degenerate Hypergeometric Functions and Their Solutions Using the Pathway-Type Integral

: In recent years, fractional kinetic equations (FKEs) involving various special functions have been widely used to describe and solve signiﬁcant problems in control theory, biology, physics, image processing, engineering, astrophysics, and many others. This current work proposes a new solution to fractional λ − kinetic equations based on generalized degenerate hypergeometric functions (GDHFs), which has the potential to be applied to calculate changes in the chemical composition of stars such as the sun. Furthermore, this expanded form can also help to solve various problems with phenomena in physics, such as fractional statistical mechanics, anomalous diffusion, and fractional quantum mechanics. Moreover, some of the well-known outcomes are just special cases of this class of pathway-type solutions involving GDHFs, with greater accuracy, while providing an easily calculable solution. Additionally, numerical graphs of these analytical solutions, using MATLAB Software (latest version 2023b), are also considered.


Introduction
The calculus of derivatives and integrals of any order is the subject of fractional calculus.Because of its scientific uses over the previous three decades, its importance and appeal to researchers have grown to in this field [1,2].Furthermore, fractional calculus has several applications in physics, finance, biology, image processing, engineering, and other disciplines.It can, for example, be used to simulate the behavior of viscoelastic materials, which have both elastic and viscous properties.Moreover, it can also solve differential equations that standard calculus cannot.We encourage readers to use the cited sources for more information on fractional calculus and its applications [3][4][5].
FKEs are types of differential equations that describe a system's evolution over time and are based on the concept of fractional calculus [6,7].One of the most essential applications of FKEs is modeling anomalous diffusion.Anomalous diffusion is a type of diffusion in which the mean square displacement of a particle grows sublinearly with time.This type of diffusion is observed in various systems, including turbulent fluids, porous media, and biological tissues (see, e.g., [8][9][10][11]).Furthermore, FKEs have also been used to model a variety of other systems, including, for example, the relaxation of viscoelastic materials, the spread of disease, the evolution of traffic patterns, and the behavior of financial markets [12][13][14].
Kinetic equation (KE) is defined by (cf.[1][2][3]) The kinetic equation of fractional-order (KEFO) is presented by (see, e.g., [6][7][8][9]) where 0 D −α w is the fractional integral operator (cf.[1,2]) given by Recently, several generalizations of the KEFO have been proposed by several researchers using various special functions.For instance, Alqarni et al. [15] introduced solutions involving the generalized incomplete Wright hypergeometric functions.Meanwhile, Khan et al. [16] considered solutions for KEFO associated with the (p, q)−extended τ−hypergeometric and confluent hypergeometric functions, and Abubakar [17] focused on a solution for KEFO using the (p, q; l)-extended τ-Gauss hypergeometric function.In a similar vein, Fuli He et al. [18] derived the solution of KEFO in terms of the Hadamard product of (p, k)−hypergeometric functions.Further, Hidan et al. [19] discussed the solution of FKEs involving extended (k, t)−Gauss hypergeometric matrix functions.Additionally, the general fractional kinetic model involving the hypergeometric super hyperbolic sine function via the Gauss hypergeometric series was recently presented by Geng et al. [20].Building upon this research in this context, we present a new generalized structure of the λ−KEFO involving a degenerate generalized hypergeometric function, which was further supported by graphical presentations derived through pathway-type transform approach methodologies.Our investigation provides detailed and different results from those given in previous works and derives several special cases of known and new outcomes.

Some Definitions Related to the Concept of DEGENERATE
Carlitz [21,22] initiated the theory of degenerate polynomials by introducing the degenerate forms of the conventional Bernoulli and Euler polynomials.Later, mathematicians investigated degenerate versions of several special numbers and polynomials, such as degenerate Stirling numbers and polynomials [23,24], degenerate Bernoulli and Euler polynomials [25], degenerate generalized Bell polynomials [26], degenerate generalized Laguerre polynomials [27], degenerate Gould-Hopper polynomials [28], degenerate Hermite polynomials [29], and so on (see [30][31][32] and the references therein).It is also worth mentioning that the study of degenerate versions of polynomials and numbers extends to special functions, such as the Euler zeta function [33], the gamma, digamma, and polygamma functions [34,35], and degenerate hypergeometric functions [36,37].Degenerate exponentials are especially potent tools that can be used to solve various problems in probability, statistics, and other areas of mathematics, such as differential and integral equations.

Pathway-Type Transform
The pathway-type integral transform (P ς −transform) of a function f (w), of a real variable w, denoted by P ς [ f (w); ϕ], is an important concept in mathematics.It has been defined initially in [15,39,40] under certain conditions on f (w), as well as the condition ς > 1.This transform can be used to convert a given function from its original domain into another domain and can thus provide insight into the properties and behavior of the original functions.Moreover, the pathway-type integral transform is regarded as one of the most helpful mathematics tools and has numerous applications (see, e.g., [41][42][43]).Indeed, the pathway-type integral transform is strongly linked to the Laplace transform, the Mellin transform, and the Fourier transform; see for example [15,43,44].
The pathway-type transform (P ς −transform) is given in [15,39,40] in the form with and lim where (L[., .]) is the Laplace transform (see, [1,2]).Moreover, the P ς -transform of a function f (w) converges under certain conditions as established in [39] by the following theorem.
(ii) for arbitrary positive a, the integral a υ f (w) dw resort to a finite limit as υ → 0 + ,; (iii) f (w) = O(e εw ), ε > 0 as w → ∞ where O(.) is the standard big O notation which means f (w) is of order not exceeding e εw .

Main Results
The purpose of this section is to study the solutions of the new generalized form of the fractional λ−kinetic equations involving GDHFs, as well as some special cases. is where E θ,ϑ (η) is the generalized Mittag-Leffler function defined in [15] as Proof.Projecting the Equation ( 22) to the P ς -transform ( 16) and using the relations ( 20), (21), and (9), we observe that where K(ϕ) = P ς P ς (w) : ϕ .Upon using the following geometric series and simple computation yields Applying the inverting P ς −transform to the last equation and making use of the relation (24), we arrive at the desired result.
Proof.We obtain the required result in the same manner as provided in the proof of Theorem 2; thus, the details have been avoided.
Proof.The proof is similar to that of Theorem 2 with choosing σ = ℘ and α = β.
. The solution of the following equation is given by where E θ,ϑ (η) is defined in (24).
Proof.By taking the P ς -transform of both sides of (30) and then using (15), we can obtain the proof as in Theorem 2.

Special Cases
Some examples of special cases of the theorems mentioned earlier are archived as Example 1.By insert (10) into (22), we have the solution of which is where E θ,ϑ (η) is defined in (24).
Remark 3. The results already established in [6][7][8][9][15][16][17] can be easily obtained as special cases by applying the functions in Remark 2 to Theorems 2-6, which gives us a more comprehensive understanding of the topic and allows for further exploration.

Numerical Representations of the Solutions
In this section, we give graphical representations of the solutions discussed in the previous section for certain values of their parameters.Figure 1 depicts the plots of the solutions to Equation (23) with the chosen values for the parameters as K 0 = 1, m = 2, n = 3, and w = 1, 2, . . ., 5 for different values of λ in the interval (0, 1).In Figure 1A-C, we set the value of α to 0.05 while that for β, σ, and ℘ equal 1.We notice that the shape of the solutions when the values of λ in the interval ( 1 r+γ 1 , 1 r+γ 1 −1 ) tend to be negative because of the negative values of the degenerate gamma function (6) in this interval, as shown in Figure 1A.However, for λ being in (0, 1 r+γ 1 ) and ( 1 r+γ 1 −1 , 1), we observe positive valued solutions, which are large for small values of λ and smaller for its values when approaching one, as shown in Figure 1B and 1C, respectively.
To illustrate the impact of α in the solutions of Equation ( 23), we fix the values of λ to be 0.118, 0.055, 0.75 while choosing 1 to be the values of β, σ, ℘, K 0 , and setting m = 2, n = 3 in Figure 2. A similar pattern to the solutions in Figure 1 can be seen clearly for various values of α in Figure 2A-C, with small-scale deviation in the solutions as α increases.Moreover, the large values of β lead to the large solutions while the small values of β lead to the small solutions in Figure 3, with fixed values of σ = ℘ = K 0 = 1, α = 0.05, and λ = 0.055.Furthermore, we could change the values of α, β, and λ to obtain more accurate solutions that do not appear in previous works.To show the significance of the derivatives of the GDHF in the solutions of Equation ( 33) in Figure 3, we fix the values of β, σ, ℘, K 0 to equal 1, and set m = 2, n = 3 with the different values of λ in Figure 4A of in Figure 4B.

Conclusions
The KE has been widely studied due to its usefulness in astrophysical issues and others.Recently, FKEs have been investigated to describe anomalous reactions in dynamical systems [1][2][3][4].Various researchers have established solutions to these families of FKEs using Laplace transform [6][7][8][9], Millen transform [49], Sumudu transform [10,11], and pathwaytype transform [15,39,40].To expand upon this research base, the authors developed a new and generalized form of the λ−KEFO involving GDHFs.This new generalization can be used for computing changes in chemical composition, such as those found within stars like our sun [20,50].Furthermore, we used the P ς − transform approach, but other authors applied different transform techniques, and existing results are particular cases of these results, implying that the current work is a generalization.Further, our analysis yields detailed and different results from those presented in previous studies and various situations of known and new results.In addition, plotting the solutions numerically was also supplied to show their conduct and examine unique cases for FKEs.