#### Abstract

Mittag-Leffler function is a natural generalization of the exponential function. Recent applications of Mittag-Leffler function have reshaped the scientific literature due to its fractional effects that cannot be obtained by using exponential function. Present motivation is to define a new special function by modification in the original gamma function with Mittag-Leffler function. Properties of this modified function are discussed by investigating a new series representation involving delta function. Hence, the results are also validated with the earlier obtained results for gamma function as special cases. Furthermore, the new function is used to generate a probability density function, and its statistical properties are explored. Similar properties of existing distributions can be deduced.

#### 1. Introduction and Preliminaries

##### 1.1. Motivation

Use of fractional operators has become a remarkable tool to study new technology and researches [1–17]. One attractive aspect of this approach is to model complex systems by taking into account the effects of filtration and memory. How the fractional operators are defined? It might be summarized as follows: “each fractional operator is the action of a particular special function on a general function, which usually represents the solution of a well-defined problem by using a suitable integral transform.” See for example [18–22]. Hence, a basic motivation and a significant contribution of this research are to study a special function containing the Mittag-Leffler function (also known as “queen of fractional calculus” [23]), which is represented in terms of complex delta function. Hence, its action on a suitably chosen function over a specific domain is easily obtainable in view of the standard properties of the delta function. This gives insights for the emergence of a new fractional operator as one of the natural consequences of this research.

The review of literature [1–23] demonstrates that many of the new researches depend on the involvement of Mittag-Leffler function in place of exponential function. More recently, Mainardi and Masina [24] have presented modified exponential integral by using Mittag-Leffler function. Mainardi et al. [25] have explored a new model which depends on . This model represents the Maxwell body at one end point, namely, , and the Becker body at other end point, namely . Referring to these phenomena, creep rule is generalized in the sense that Becker model involves Mittag-Leffler function of order in place of exponential function. Following them, Paris [26] has discussed asymptotic behavior of the modified exponential integral and studied its more general form by using two-parameter Mittag-Leffler functions. Taking motivation from these researches [24–26], the aim of the present paper is to study the modified gamma function by using Mittag-Leffler function. Before going on to our main results, we describe necessary background and preliminaries in the subsequent paragraphs.

##### 1.2. Preliminaries

Mittag-Leffler [27] proposed a function that appears a natural substitute of exponential function in the case of fractional order. It is known as Mittag-Leffler function defined by

For one-parameter Mittag-Leffler function reduces to the exponential function There are many generalizations of Mittag-Leffler function. For example, two- and three-parameter Mittag-Leffler functions are defined byrespectively. Here, denotes the set of complex numbers and denotes the real part of any complex number. Three-parameter Mittag-Leffler function is also known as Prabhakar function [20]. The decades of this century are mainly devoted to the analysis and applications of this function, see for example [1–24], and related references therein. For our purposes, we would be using the following form of Mittag-Leffler function [28]:

Here, denotes the set of real numbers. Similarly, the definition of fractional circular functions is [28, 29]which are obviously related to the Mittag-Leffler function defined in equation (3). Similarly, the hyperbolic functions can be defined as [28, 29]

The focus point of the current investigation is to study a modified (generalized) form of gamma function [30], defined by

It is a basic special function so that Pochhammer symbols have the following relationship with the gamma function:

Here, denotes the set of natural numbers or positive integers. Further detailed analysis of the gamma function can be found in [31–34] and related bibliography therein. It is important to mention that this contribution remained unachievable without reviewing [35–46]. As we will be using delta function for our results and therefore we briefly describe about the spaces of test functions and distributions in the subsequent paragraphs.

A dual space also known as space of generalized functions (distributions) corresponds to a space of test functions. Such phenomenon involves the representation of singular functions and hence the certain functions can be taken into account and various calculus operations can be applied in this way. Let us begin with the following simple identities:and for a nonvanishing parameter

Gelfand and Shilov [47] have provided a comprehensive discussion and explanation of distributions (or generalized functions). According to them, [47], the commonly used test functions that have compact support are the elements of a space of test functions denoted by and another space denoted by is of infinitely differentiable as well as fast decaying functions. These spaces and have their dual spaces denoted by and , respectively. When we talk about the closeness of such spaces under Fourier transform, then these spaces behave differently. Hence, the spaces and are closed but and are not closed under Fourier transform. As a result, the elements of have Fourier transforms that form distributions for entire functions spaces whose Fourier transforms belong to [48]. In addition to it, since the entire function is nonzero only for a specific interval , so the following relation among these spaces is valid

For the benefit of this investigation, space comprises of entire and analytic functions sustaining the subsequent criteriawhere and are numbers that depend on function . The subsequent equations ([46], Vol 1, p. 169, equation (8); [47], p. 159, equation (4); and [48], p. 201, equation (10)) are vital to follow the steps of the proofs in this research

Hence, and for each ,

In addition to this, the other functions whose Fourier transform involves delta function are , and . Further detailed study about similar type of functions is given in [47–50].

Plan of the remaining paper is given as follows: Section 2.1 deals with the definition of the new function by investigating a series form of modified gamma function in terms of complex delta function. Section 2.2 discusses the conditions of persistence along with possible uses of the new series. Authentication of such results is discussed in Section 2.3. Advance properties of the new function as a distribution are analyzed in Section 2.4. Further examples by using new representation are discussed in Section 2.5. Generalized gamma distribution is discussed in Section 2.6. This work is concluded in Section 3. The results presented here are validated with the earlier obtained results in [35–37].

#### 2. Results

##### 2.1. Modified Form of the Gamma Function and Its Series Representation

By following the methodology of [24], we introduceand for , we get the original gamma function (6). According to the particular asymptotic behavior of Mittag-Leffler function ([2], p. 286, equations (2.8-2.9)], this integral converges for when . Next, by following Lebedeve [30], we rewrite (16) as

Here, first integral is an analytic function of and second integral is an entire function Now, we replace by its power series such as

Uniform convergence of the involved integral allows one to reverse the order of integration and summation in equation (18). It shows that has singularities at

Continuing in the same way, we investigate a new series representation of in the form of the following theorem

Theorem 1. *The modified gamma function has the following representation:*

*Proof. *Let us first replace and in equation (16), then we getNext, we replace Mittag-Leffler function by its series form as given in (1),and then by rearranging the involved exponential termswe getNext, the uniform convergence of the involved integral allows us to interchange the operations of integration and summation. Hence, by using equation (12), we can writeThe above expressions (23) and (24) lead to the result (19) as mentioned above.

Corollary 1. *The modified gamma function has the subsequent demonstration*

*Proof. *This can be obtained by using equations (14) and (19) as follows:Next, by using this relation in (19), we obtain the desired result.

Corollary 2. *The modified gamma function has the subsequent demonstration*

*Proof. *This result can be obtained by a small alteration in equation (19) as follows:which leads to the required result in view of (19).

Corollary 3. *The modified Gamma function has the subsequent demonstration:*

*Proof. *This can be deduced in view of the relation (14) as follows:Next, by making use of this result in (27), we obtain the required form.

Corollary 4. *Gamma function has the subsequent representation ([36], p. 2092, equation (2.7))*

*Proof. *This can be deduced by fixing in equation (19).

Corollary 5. *Gamma function has the subsequent series representation [37]*

*Proof. *This can be deduced by fixing in equation (25).

Corollary 6. *Gamma function has the following representation [37]:*

*Proof. *This can be deduced by fixing in equation (27).

Corollary 7. *Gamma function has the subsequent representation [37]:*

*Proof. *This can be deduced by fixing in equation (29).

*Remark 1. *It is to be remarked that the newly obtained series representation is only meaningful if defined as a distribution acting over a space of test functions. We can observe it by simply multiplying the both sides of equation (27) with ,which produces the following:As proved above that has poles at , here, the right hand side sum will disappear because is in the denominatorwhich is obviously false. At the same time, if we consider the following inner product,we obtainNow, the right-hand side will vanish due to the fact that has poles at , hencewhich is false. This discussion demonstrates that we have to be very careful and rigorous in selecting the set of functions for a meaningful definition of this new representation. This is a part of Section 2.2.

##### 2.2. Existence and Applications of New Representation

This exemplification of the generalized gamma function is an infinite summation over delta function, which is certain only if defined in the way as distributions (generalized functions). Therefore, it is curious to show that the obtained series of delta function is a distribution (generalized function) on a particular space of test functions as proved in the subsequent theorem.

Theorem 2. *Prove that the modified gamma function acts as a distribution over the space of slowly increasing complex test functions denoted by .*

*Proof. *Let us first take the following expression for any , and :It shows that is a linear function over Now, we take an arbitrary sequence and use the fact that acts as continuous function over ,Then, we can show that (tends to zero). Hence, the modified gamma function is a distribution (continuous linear functional) over because of the convergence of (27) given byIn the above equation, we have applied the following property for delta function,In this way, the attained functions are test functions of slow growth (bounded by a polynomial); also, we can observe that summation of the coefficients in (44) isHence, we get a finite and definite answer, which shows that is convergent, for Here, we have used the well-known fact that “the product of a rapid decay and slow growth functions is convergent.” However, it is also obvious from the famous Abel’s theorem. Hence, the same is true for its special cases which are given in (25)–(29).

Let us move to a larger class of functions and consider In view of equation (27) and shifting characteristic of delta functions, the product is computed as follows:By putting in (47), we can get the following known result for gamma function:

*Remark 2. *Particular engineering problems involve sum of delta functions. It is remarked that when one considers the product of with sequence of delta functions , then it yields the modified form . One can validate this result for gamma function at .

The above discussion explores the existence of further results. To illustrate, we take in equation (47), which then leads to Laplace transform of . Hence, we consider the validation of such outcomes in the subsequent part.

##### 2.3. Validation of the Results Obtained by New Representation

In this section, we will testify our results with the existing results. Let us take and in (16), then the modified gamma function is expressible asand putting , in the above equation, we obtain

We know that Fourier transform satisfies the following dual property for any function

Hence, by applying this property for equations (49)–(50), we get the following result:or

This can also be accessible as a particular case of the main result (47) when one puts . So, the results obtained by using the new series representation are accurate when compared with well-known results. Furthermore, when one takes in (49), then the following identity can be obtained:which is also a special case of (53). Hence, it is testified that the new representation of creates new results that are unachievable in view of known representations but the special cases of these novel results are reliable with the existing methods. Continuing in this way, by taking in (53) and (54), we obtainand in (55) yields the following well-known result:

##### 2.4. Advanced Distributional Properties of Modified Gamma Function

Here, by taking motivation from ([48], Chapter 7), we state and prove a list of basic properties of the modified gamma function.

Theorem 3. *For an arbitrary test function , do satisfy the subsequent properties as a distribution over the space of test functions:*(1)*.*(2)*.*(3)*.*(4)*.*(5)*.*(6)*For any regular distribution , is a distribution over *(7)*.*(8)*.*(9)*(10)**.*(11)*.*(12)* where .*(13)*.*(14)*.*(15)*.*(16)*,**where are any complex or real numbers.*

*Proof. *First six results (1)–(6) can be proved by using the basic properties of delta function and by following the procedure of Theorem 2. Here, we will prove the next result (7). Hence, by taking th derivative of representation (27) w.r.t and using (14), one can obtainwhich is a well-defined and finite product. Then, one can prove the next result as follows:Again, by using the property of delta function, one can get the next result (9),Similarly, the results (10)–(15) are obvious to hold in view of Fourier transform and the properties of complex delta function. Let us start proving the next identity (10) from the above list,Next, we prove identity (11) from the above list,and identity (12),where represents transpose of the involved test function *f*. Hence, the above identity (13) can be observed by the following:and by applying Parseval’s identity of Fourier transform, the proof of identity (14) isNext, the identity (15) follows from following details:and so on, we getContinuing in this way, one can prove the last identity (16) by making use of the result ([48], p. 201),Therefore, we haveas required.

*Remark 3. *According to the fact stated and proved in ([48], p. 203), “Fourier transform maps onto and the inverse Fourier transform maps onto .” It is remarked that is an element of , and is also in in view of (49) and (52), respectively.

##### 2.5. Further Examples by Using New Representation

Since delta function is a linear mapping which can map each function to its value at zero. By taking into account this fact, we consider the following examples. Let us consider (29) and restrict the variable, , to real numbers, then we have

Hence,

Next, we consider the following examples.

*Example 1. *Let , then ,

*Example 2. *Let , then ,Similarly, , then

*Example 3. *Let , then ,

*Example 4. *Let , then ,Here, we restrict the value of very small such that . These examples show that our new representation of the modified gamma function is meaningful for all those functions that have derivatives of all orders at point 0. Similar results hold for its special case gamma function given by equation (34).

Next, we consider Laplace transform of delta function that yields the Laplace transform of the modified gamma and the original gamma function itself. We know thatTherefore, we getPutting in the above equation giveswhich yields furtherA special case taking in (78) leads to the well-known result given above in (56). Hence, it validates that the results obtained by new representation are consistent with the known results.

*Remark 4. *It is to be remarked that the success of the application of this new representation lies in the fact that the sum over the coefficients in the new representation is finite and well-defined as given in (46).

##### 2.6. Generalized Gamma Distribution

Using a probability density function, Dirac delta function is frequently used to represent certain probability distributions. For example, considering the probability density function u(x) of a discrete distribution with points and corresponding probabilities it can be written as

Moreover, gamma distribution and its generalizations are of interest for a large audience (see [51] and the cited bibliography therein). Let us consider the following new probability density function of a statistical distribution for a random variable as follows:where is defined in (16) and it is suitably assumed that for the involved parameters,

It can also be verified that

Next, we discuss the usual properties of a random variable which is distributed according to the probability density function given in (81).

The general moment of order for the random variable is obtained aswhere denotes the mathematical expectation.

Mean and expected value for the probability distribution function is the first moment of the variable and obtained as a special case of the above for ,and for ,

Hence, the variance of the random variable is computed as

By using (84), the characteristic function of the random variable is given as follows:

Similarly, the generalized moment generating function of the random variable by using Mittag-Leffler function instead of exponential function is given asand the usual moment generating function of the random variable is obtained as a special case of (89),

The cumulative distribution function of the random variable for is computed as

The Lorenz curve as introduced by Lorenz [52] is commonly used to represent the disparity of money and earnings in economics. For our new probability function, it is computed for a random variable as follows:

Hence, it is to be remarked that the results for the other probability functions, for example,can be obtained as special case of the above statistical measures by taking .

#### 3. Conclusion

During this study, we explored a new special function as a modified form of gamma function so that the gamma function can be obtained as a special case of it. We discussed some basic properties of this function. The combination of distribution theory with different integral transforms is well explored for the analysis of partial differential equations (PDE). Numerous practical questions are impossible to be answered by applying the known techniques but became possible by using this combination. Hence, a new fractional operator naturally emerges as a consequence of this research which can act on a suitably chosen function as follows:

A new probability density function involving the new function is introduced, and its properties are studied. We end this discussion by proposing a generalized representation of the familiar beta function,in terms of the new function as follows:It is customary to consider Mittag-Leffler function as given in equation (3) for certain applications. However, it can be noticed that if we consider (1) instead of (3) in (16) for 0 < ℜ (α) < 1, then (18) would be convergent for a greater set of values u ∈ ℂ\{0, −1, −2, …}. Similarly, equation(2) can be used for further extensions.

#### Data Availability

No data were generated or analyzed for this submission.

#### Conflicts of Interest

The authors declare no conflicts of interest.