All articles published by MDPI are made immediately available worldwide under an open access license. No special
permission is required to reuse all or part of the article published by MDPI, including figures and tables. For
articles published under an open access Creative Common CC BY license, any part of the article may be reused without
permission provided that the original article is clearly cited. For more information, please refer to
https://www.mdpi.com/openaccess.
Feature papers represent the most advanced research with significant potential for high impact in the field. A Feature
Paper should be a substantial original Article that involves several techniques or approaches, provides an outlook for
future research directions and describes possible research applications.
Feature papers are submitted upon individual invitation or recommendation by the scientific editors and must receive
positive feedback from the reviewers.
Editor’s Choice articles are based on recommendations by the scientific editors of MDPI journals from around the world.
Editors select a small number of articles recently published in the journal that they believe will be particularly
interesting to readers, or important in the respective research area. The aim is to provide a snapshot of some of the
most exciting work published in the various research areas of the journal.
The integral representation of the two-parameter Mittag-Leffler function is considered in the paper that expresses its value in terms of the contour integral. For this integral representation, the transition is made from integration over a complex variable to integration over real variables. It is shown that as a result of such a transition, the integral representation of the function has two forms: the representation “A” and “B”. Each of these representations has its advantages and drawbacks. In the paper, the corresponding theorems are formulated and proved, and the advantages and disadvantages of each of the obtained representations are discussed.
The Mittag-Leffler function is an entire function defined by a power series
where is the Gamma-function. This function was introduced by Mittag-Leffler in a number of works published from 1902 to 1905 in connection with his development of a method of summing divergent series. For more detailed information on the content of these works and on the history of the introduction of the Mittag-Leffler function, we refer the reader to the book [1] (see chapter 2 in [1]). The function itself was introduced in the work [2]. In the paper [3], the integral representation for this function was obtained that expresses its value in terms of the contour integral.
In this paper the two parameter Mittag-Leffler function
is studied. This function was first introduced by A. Wiman in 1905 [4,5]. Later in 1953 this function was rediscovered in the works of Humbert and Agarval [6,7,8]. A new function was introduced by replacing the additive unit in the Gamma function argument in for an arbitrary complex parameter . At the same time, irrespective of Humbert and Agarval, the function (1) was studied by M.M. Djrbashian in the papers [9,10] (see also [11] Chapter 3, §2, 4). As we can see, the two parameter Mittag-Leffler function is connected with the classic Mittag-Leffler function by a simple relation . For more detailed information on the properties of the function we refer the reader to the book [1], as well as to other review works [12,13,14]. In this paper, integral representations of the function will be obtained and studied.
The integral representation of the Mittag-Leffler function is important from the point of view of its practical use, as well as for studying the asymptotic properties and zeros of this function. The integral representation expressed through the contour integral is used for these purposes. Several such representations are known for the Mittag-Leffler function. One of the earliest integral representations of the function was given in the book [15] (see §18.1, formula (20)). Further development of the issue of the integral representation of the Mittag-Leffler function and the study of its asymptotic properties was carried out in the works of M.M. Djrbashian. In the work [9] the integral representation was obtained that expressed the Mittag-Leffler function through the contour integral. Later it was included in his monograph [11] (see chapter 3, §2, Lemma 3.2.1). Using this integral representation, asymptotic formulas and the distribution of the zeros of the Mittag-Leffler function were obtained. Further, the results of the work [9,11] were used in the books of [1,16], as well as in the works of [17,18,19,20,21] to develop the methods and algorithms of calculating the Mittag-Leffler function. However, despite a wide use of the integral representation for the Mittag-Leffler function that was obtained in the works of [9,11], it turned out that there was a mistake in that representation. This fact was pointed out in the work of [22]. In this regard, there is an issue of obtaining the correct integral representation for the Mittag-Leffler function.
One of the possible solutions to this issue was given in the works of [22,23]. In these works, the representation of the Mittag-Leffler function was obtained that expresses its value through the contour integral. As it was noted earlier, this representation is used to study the asymptotic properties of the Mittag-Leffler function. However, for practical use and for calculating the value of the function, it is convenient to have integral representations expressing the function in terms of the integrals of real variables. This paper is devoted to obtaining such integral representations for the Mittag-Leffler function. The starting point of the solution to this problem is the integral representation of the function obtained in the work [22]. Running a little ahead, we will say that the transition from integration over a complex variable to integration over real variables leads to the appearance of two forms of the integral representation of the Mittag-Leffler function. The first form will be abbreviated as the representation “A”, the second as the representation “B”. The representation “A” is a direct consequence of the integral representation obtained in the work [22]. It is obtained as a result of the transition from the contour integral to integration over real variables. To obtain the representation “B” in addition to performing such a transition, it is necessary to carry out a terminal transition . This leads to the fact that the integral representation “B” is valid only for parameter values satisfying the condition . As a result, both the representation “A” and the representation “B” have its advantages and drawbacks which will be discussed in detail in the paper.
It should be pointed out that in this paper the letters are used to denote the parameters of the Mittag-Leffler function (1). These notations were introduced by M.M. Djrbashian in his works [9,10,11]. Parameters and of the function (1) are connected with the commonly used notations and of the parameters of the Mittag-Leffler function by simple relations , .
2. Integral Representation “A”
The purpose of this paper is to obtain integral representations for the function expressing this function in terms of integrals over real variables. It is convenient to use the representations of such a kind in practical problems as well as for calculating the values of the Mittag-Leffler function. As an example of the use of such integral representations, we can mention the work [24]. In that paper, the inverse Fourier transform of the characteristic function of the fractionally stable law was performed which is expressed in terms of the Mittag-Leffler function. To perform the inverse Fourier transform, we used the integral representation of the Mittag-Leffler function which in the current paper is formulated in Corollary 3 item 2 (see (84)). Using this integral representation in the paper [24] expressions of density and distribution function of a fractionally stable law were obtained. Articles [25,26] are another example of usage of the Mittag-Leffler function. In these articles are shown that solutions of a master equation for the fractional Poisson process [25] and a fractional relaxation equation for dielectrics [26] are expressed thorough the Mittag-Leffler function. A Monte Carlo method was used for calculation the obtained solutions. However, the use of the integral representation of the Mittag-Leffler function for calculation of the solutions would significantly increase the accuracy of the results.
The starting point of this paper is the integral representation of the function , obtained in the work [22]. The following theorem was formulated in this work
Theorem1.
For any real satisfying the conditions , , , , any and any satisfying the condition
the Mittage-Leffler function can be represented in the from
where the contour of integration has the form (see Figure 1)
The proof of this theorem can be found in the work [22].
The form of the contour of integration on the complex plane is given in Figure 1. The cut of the complex plane goes along the positive part of a real axis. The contour of integration consists of the half-line , the arc of the circle radius and the half-line . The contour is traversed in a positive direction. The parameters and have the meaning of inclination angles of the half-lines and in relation to the contour axis . In Figure 1 the contour axis coincides with a real axis. The values of the angles and are measured from the negative part of the real axis. The reference directions of these angles in Figure 1 are shown by arrows.
We will obtain the integral representations that of interest to us by going from integration over the complex variable to integration over variables r and interconnected by the relation . As a result, the following theorem is true for the function .
Theorem2
(The representation “A”). For any real , and , satisfying the conditions
for any complex μ and any complex z satisfying the condition
the Mittag-Leffler function can be represented in the form
Here
where
and
where
Proof.
For convenience we will introduce the notation
As a result, the representation (3) will take the form
where the contour of integration is defined by (4).
Substituting in (14) the variable of integration and directly calculating this integral we obtain
We consider the integral . Getting rid of the complexity in the denominator we have
where the notation was introduced.
Similarly, for the integral we get
Summing now the integrals and we obtain
where
We consider now the integral . At the beginning, we will get rid of the complexity in the denominator. To do this, we will multiply and divide the integrand by the complex conjugate of the denominator and open the brackets in the denominator. Then, in the resulting expression, we substitute the definition of the function (see (13)) and in the indices of the exponents we use the Euler formula . As a result, we get
where
Using (16) and (17) in (15) we obtain the representation (7). It is important to pay attention that in the process of proving no additional limitations on the values of parameters , , , and the argument z were imposed and it means that the ranges of admissible values pass from Theorem 1 without change. Thus, the representation (7) is valid for any real , any real , satisfying the conditions (5), any complex and any complex z satisfying the condition (6). ☐
The proved theorem formulates an integral representation for the Mittag-Leffler function that expresses this function in terms of the sum of improper and definite integrals. To be definite, we will call this integral representation of the Mittag-Leffler function the representation “A”. As we can see, the representation “A” is a direct consequence of the representation (3). It is obtained by passing from the contour integral to integrals over the real variable. Moreover, the improper integral in (7) corresponds to the sum of integrals along the half-lines and of the contour and the definite integral corresponds to the integral along the arc of a circle . It should be noted that this integral is taken along the arc of a circle of radius , where . The representation (7) is valid for arbitrary values , any and any and that satisfy the condition (5).
However, in general case, for arbitrary values of parameters and , satisfying the condition (5), the kernel function turns to be lengthy. The representation (7) takes the more compact form in case when the half-lines and of the contour of integration run symmetrically relative to the real axis i.e., when . We will formulate the obtained result in the form of a corollary.
Corollary1.
For any real , any complex μ the following integral representations of the Mittag-Leffler function are true:
1.
at any real , any real satisfying the condition and any complex z satisfying the condition
(1) According to Theorem 2, the Mittag-Leffler function can be represented in the form (7). This representation is true for arbitrary and satisfying the conditions (5). In case if the conditions (5) take the form
Using the definition (see (9)) and the fact that the function defined by (10) is an odd function according to the variable
we have
Now substituting this result in (25) we obtain (19). Since in the proof process no additional restrictions on the values of parameters , and on the value were imposed, then the conditions for these parameters go from Theorem 2 without change. Thus, we come to the conditions of the corollary.
(2) We consider the case . As we can see, the case considered is a particular case of the previous one. It follows from (5) that this case can be implemented if . For the range of values from (23) we get
Now we consider the representation (18). In the case under consideration it will be written in the form
We denote
From (10) it follows that . Now using this result in (19) we get (21). ☐
From the proved corollary it follows that if the parameter values and coincide, then in this the kernel function is significantly simplified. Recall that the parameters and are inclination angles of half-lines and in the contour relative to the axis of this contour (see Figure 1). Since in the Theorem 1 the axis of the contour coincides with the real axis, then the selection of means that half-lines and run symmetrically in relation to the real axis. The kernel function takes the simplest form in the case .
Further, it is necessary for us to know the position of singular points of the integrand of the representation (3). This issue was studied in the work [27]. For completeness of the statement here we give the result obtained in the work [27] and formulate it in the form of a lemma.
Lemma1.
For any real and any complex values of the parameter the integrand of the representation (3)
relative to the variable ζ has two singular points and . The point is a pole of the first order. The point is:
1.
the regular point of the function , with the values of parameters , where (the positive integer), and , where (the non-negative integer);
2.
a pole of the order , if , where (the positive integer), and , where (the positive integer);
3.
the branch point, for any other values of parameters .
The proof of this lemma can be found in the work [27].
We will make the following remark to Corollary 1.
Remark1.
In Corollary 1 the special case under consideration . We will assume that and study the behavior of the Formula (20) in this case. As a result, we obtain
From here it is clear, if , where , then . Consequently,
Thus, with integer values of μ the first summand in (28) becomes zero and to calculate the value of it remains to calculate the second integral. To calculate this integral, numerical methods can be used. However, in this case, this integral can be calculated analytically using the residue theory.
In fact, we return to the integral representation formulated in Theorem 1. Recall that we consider the case . Using the notation (27), the representation (3) takes the form
where the contour of integration , defined by (4), is written in the form
We represent the complex parameter μ in the form and make use of Lemma 1. According to this lemma, the function at values and , where has one singular point which is a pole of the first order. The point , in this case, is the regular point. In case, if and , where the function has two singular points: the point is a pole of the first order and the point is a pole of the order . As we can see, in both cases the point is not a branch point. As a result, in these two cases, the function is the entire function of a complex variable ζ. From here it follows that when , and , where the arc of the circle that enters the contour (31) is the closed circle of radius . The half-lines and pass along the positive part of a real axis in mutually opposite directions. With all other values of the parameter μ (when or ), according to Lemma 1, the point is a branch point of the function . In this case, the circle of the contour (31) will not close up and half-lines and will go along the upper and lower banks of the cut of the complex plane which runs along the positive part of a real axis.
It is clear from here that the result (29) is a consequence of Lemma 1 In fact, in the case when the parameter μ takes integer real values, the first and second items of Lemma 1 turn out to be true. As we have already pointed out, in this case the arc of the circle of the contour (31) is a closed circle and the half-lines and run along the positive part of a real axis in mutually opposite directions. Consequently, the sum of the integrals along the half-lines and will be equal to zero. Next, it is necessary to recall that the improper integral in the expression (28) just corresponds to the sum of the integrals along the half-lines and . Therefore, with integer real values μ it should be equal to zero which has been obtained. A definite integral in (28) corresponds to integration along the closed circle. Therefore, one can use the theory of residues to calculate it.
The calculation of the integral in (30) using the theory of residues with integer real values of the parameter was carried out in the work [27]. For completeness of the statement, we will give the results obtained in this paper and formulate them in the form of a corollary to Lemma 1.
Corollary2.
For the values of the parameters , , any complex z, satisfying the condition and for integer real values of the parameter the Mittag-Leffler function has the form:
1.
if (i.e., ), then
2.
if (i.e., ), then
The proof of this corollary can be found in the paper [27].
3. Integral Representation “B”
The integral representation “A” consists of the sum of two integrals. As it has been found earlier, the improper integral in (7) corresponds to the sum of integrals along the half-lines and of the contour in the representation (3), a definite integral corresponds to the integral along the arc of the circle . As a result, in analytical studies of the function , as well as in the solution of problems where it is encountered, one should conduct studies of these two integrals. This causes certain difficulties. As it will be shown below in the representation (7) one can get rid of an integral on the arc of the circle and write the integral representation for the function in the form of an improper integral. This representation will be much easier to use. However, as a result of such a transition, some restrictions are imposed on the parameter values . The integral representation of the Mittag-Leffler function represented in the following theorem will be called the representation “B”.
Theorem3
(Representation “B”).For any real and any complex μ satisfying the condition for the function , the following integral representations are valid:
1.
for any real satisfying the conditions
and any complex z satisfying the condition the Mittag-Leffler function can be represented in the form
if and , , then for any complex z satisfying the condition , the Mittag-Leffler function can be represented in the form
where is an arbitrary real number satisfying the condition ;
4.
if and , then for any complex z satisfying the condition , the Mittag-Leffler function can be represented in the form
where is an arbitrary real number satisfying the condition and has the form (19).
Proof.
The starting point of the proof is Theorem 1 and the integral representation (3) which is defined in it. In view of the notation (13), the representation (3) will be written in the form
The problem consists in calculating this contour integral.
We consider an auxiliary integral
where the contour (see Figure 2) consists of the arc of the circle of radius with the center at the origin of coordinates, the segment , the arc of the circle the radius with the center at the origin of coordinates and the segment , which are defined in the following way:
The contour is traversed in the positive direction. The cut of the complex plane goes along the positive part of a real axis.
As it follows from Lemma 1, depending on a value of the parameter , the integrand (38) has one or two poles. If , then there are two poles in the points and . If , then there is one pole in the point . In addition, at non-integer values and the point is the branch point. In the case when and the segment will pass through a singular point (Figure 2). The similar situation will be at and . In this case, the segment will pass through the point . These two cases will be considered separately. The case , also requires a separate consideration. In this case each of the segments and passes through a singular point . In view of the foregoing, it is necessary to consider four cases:
, , ,
, , ,
, , .
Case 1. We consider the first case at the beginning. Directly calculating the integral I and substituting the variable of integration we obtain
Now we let in this expression and we will study the behavior of , and . We consider the integral at the beginning. For this integral the relation is true
We consider the integrand. To get rid of the complexity in the denominator we multiply and divide the integrand by . We also represent and , and for a power function we will use the representation . As a result, we obtain
where
For the numerator we have the estimate
Thus,
From this it follows that
Now we consider the behavior of the integral at . We have the estimate
Consequently, it is necessary to study the behavior of the integrand at . For convenience, we introduce the notation . Further, similarly to the previous case, we get rid of the complexity in the denominator. To do this, we multiply and divide the integrand by the complex conjugate value of the denominator, i.e., by and represent , . For a power function we will make use of the representation . Using (13) we get
where
For the numerator we have the estimate
Consequently, we obtain
We get the similar result for the integral
Thus, at the limits and will converge to the corresponding definite integrals
We will pay attention to the fact that in the case under consideration when the parameters satisfy the condition (40), the integrand of the integral I (see (38)) inside the region limited by the contour , is an analytical function. Consequently, from (41) we have
Letting in this expression and taking into consideration (42) we obtain
Substituting here the expression for and using (43) we obtain
Here it should be pointed out that at values the integrals (43) have a singularity in the point and at values the singularity disappears.
Now let us return to the Mittag-Leffler function. Directly calculating (37) and replacing the variable of integration we get
We consider the first summand in this expression. Getting rid of the complexity in the denominator and using (13) we get
where has the form (10). Similarly, for the second summand in (46) we have
Now substituting (47) and (48) in (46) and making simple transformations we finally obtain
where is defined by the expression (8). The values and satisfy the conditions (40). The condition for the value did not change and passes from Theorem 1 without any changes. Thus, the first item of the theorem is proved.
Case 2. Now we consider the case when , and . As it was mentioned earlier, in this case the segment of the auxiliary contour (see Figure 2) will go through the singular point . That is why, it is necessary to change the contour of integration in such a way that one could bypass this point leaving it outside the contour. The auxiliary contour that we get will consist of the arc of the circle with the center at the origin of coordinates and radius (see Figure 3), the segment , the circle with the center at the origin of coordinates and radius satisfying the condition , the segment , the arc of the circle with the center in the point () and radius (here and ) and the segment .
For further study we need to parametrize this contour. As a result, in the plane the arc of the circle can be written in the form . The segment is written in the form . The arc of the circle has the form . The segments and are written in the form and .
To parametrize the arc of the circle , we consider the mapping of the complex plane on the complex plane u. This mapping is a conformal mapping and is a left shift of the entire complex plane by the value 1. As a result of such a shift, the point in the plane is mapped onto the point of the plane u. Thus, the circle with the center in the point and the radius of the plane will be mapped on the circle with the center in the point and the radius on the plane u. Thus, to parametrize the arc of the circle in the plane it is enough to parametrize the arc of the circle in the complex plane u and then to map the complex plane u on to the plane with the help of an inverse conformal mapping .
We will represent the complex number u in the form . Then, an arc of a circle in the complex plane u has the form , where . Now we map the complex plane u onto the plane . We have . From this we obtain
Here it is necessary to note that is a multivalued function. The principal branch of this function takes values in the interval . However, as one can see from the definition of the auxiliary contour , the center of the arc of the circle lies in the point in which . Therefore, it is necessary to choose the required branch in in such a way that the mapping (49) could map an arc of a circle of the plane u on to an arc of a circle of the plane with the center in the point and . As a result, the mapping (49) will take the form
For the arc of the circle we obtain . Here it should be pointed out that these formulas produce the mapping of the circle on to the circle only in the case . If , then these formulas will not be true.
Thus, in view of the remarks made the auxiliary contour can be represented in the form
The contour is traversed in a positive direction.
Next, we consider an auxiliary integral
By calculating this integral we have
One should pay attention to the fact that inside the region limited by the contour the integrand of the integral is an analytical function of the complex variable . Consequently, according to the Cauchy integral theorem
From here we get that
Thus, it is possible to substitute integration on the arc of the circle by integration on the remaining parts of the contour .
We consider each integral in the right part of this expression. The integral has already been considered by us earlier in this lemma. Using (42) we obtain that in the case considered
Consider the integral . Representing the complex number in the form and using (51) we obtain
We introduce the notation
Now we transform the function . Using the definition for (see (13)) and getting rid of the complexity in the denominator we obtain
By introducing the notation , that coincides with (10) we obtain that can be represented in the form
Thus, for the integral we get
Now consider the integral . Representing the complex number in the form and using (51) we get
Comparing this expression with (55), we notice that the integrand of the integral obtained is similar to the integrand of the integral (55). Therefore, using (56) we get that the integrand can be represented in the form (57). Thus, for the integral we obtain
Using now (56) and (57) the integral takes the form
In the sum (52) it remains to consider the integral . This integral is taken along the arc of the circle with the center in the point (at that ) and radius . The contour is traversed in the direction from the point D to the point E (see Figure 3). As it was shown earlier, the arc of the circle in the plane can be given in the form (see (51))
As one can see, the arc traverse in the direction from the point D to the point E corresponds to a change of the parameter from to . Thus, in the case under consideration the contour on the complex plane can be represented in the form
where , and
Representing now in the integral the complex number in the form (62) we find
One should pay attention that traversing along the arc of the circle in the direction from the point D to the point E (see Figure 3) corresponds to the negative traversing direction along the contour. Here the value of the parameter corresponds to the point D and the value of the parameter to the point E. That is why in the integral derived one must transpose the limits of integration.
We now consider the integrand of this integral. We introduce the following notation for it
We transform this expression using (13) and getting rid of the complexity in the denominator we get
Now we get back to the expression (53) and let in this expression. It is necessary to point out that the integrals , and do not depend on and, consequently, they will not change with such a passage to the limit. In view of this, we have
Taking into account (54) and using the expressions (58), (59), (60) and (67) we get
From this we obtain
Now we get back to the Mittag-Leffler function. Assuming that the contour of integration in (37) takes the form
By calculating (37) directly and representing a complex number in the form we obtain
using the notation (56), this expression can be written in the form
where the form of the functions is defined by (57). Now we will make use of the representation here (68) for the integral . As a result, we obtain
where . It remains to consider how to change the condition (2) in this case. Since, in the case considered , then the condition (2) takes the form . Thus, we have obtained the statement of the theorem for the second case.
Case 3. Consider now the case , , . In this case the segment of the auxiliary contour (see Figure 2) will run along the positive part of a real axis and, thus, will run through a singular point . It should be pointed out that in this case an argument of points of the segment is equal to . Consequently, in this case . That is why, as in the previous case, we will change the contour in such a way that the contour bypasses this singular point leaving it outside the contour. The contour formed (Figure 4) consists of the arc of the circle with the center in the point and radius , the segment , the arc of the circle with the center in the point () and radius , the segment , the arc of the circle with the center at the origin of coordinates and radius and and the segment .
As one can see, only the arc of the circle remains non-parametric. For its parametrization we will fulfill the procedure as in the previous case. Considering the conformal mapping and representing the complex number u in the form we get that the arc of the circle can be represented in the form (50), where the parameter varies within the limits from to 0. Note that the value of the parameter corresponds to the point C and the value to the point B. It should be pointed out that the center of a circle lies in the point with . Consequently, in the formulas (50) one should choose the principal branch , i.e., to take . As a result, an arc of a circle can be represented in the form
In view of the foregoing, the contour is written in the form
The contour is traversed in a positive direction.
We consider now an auxiliary integral
In a similar way to how it was done in the previous case, we calculate this integral. Using the definition of the contour we have
It is clear from this expression that inside the region limited by the contour the integrand of the integral is an analytical and continuous function of a variable . Consequently, according to the Cauchy theorem
From this we obtain that
We consider now each of the integrals in the right part separately.
The integral was considered by us earlier. Using (42), we obtain
Now consider the integral . Representing in the form we find
From this it is clear that this integral is similar to the integral (55). That is why one can use the study results of this integral. Using (57), we obtain
In the same way, using (70), (56) and (57) for the integrals and we find
It remains to consider the integral . This integral is taken along the arc of the with the center in the point , where . This means that the contour is traversed from the point B of the complex plane in which to the point C in which . In other words, the starting and ending points of this contour have the same value of the argument. As it was shown above, the arc of a circle can be represented in the form (69). Traversing this arc in the direction from the point B to the point C corresponds to a change of the parameter from 0 to . Thus, the points of the arc of the circle can be represented in the form
Using now the representation (75), in the integral we obtain
Comparing this expression with (64) it is clear that these two integrals are similar. Therefore, using (65) and (66) for the integral we get
We return to the expression (71) and let in this expression. Note that the integrals , , do not depend on and, consequently, with such a passage to the limit they will not change. As a result we obtain
Using here the expressions (72)–(74), (76) we obtain
As a result, performing the passage to the limit we have
We return now to the Mittag-Leffler function. Assuming that the contour of integration in (37) takes the form
and the condition (2) will be written in the form .
By direct calculating (37) and representing the complex number in the form we get
using the notation (56) this expression can be written in the form
where the form of functions is defined as (57). We substitute the integral in this expression by the expression (77). As a result, we obtain
Thus, we have obtained the statement of the theorem for the third case.
Case 4. Consider now the case , . In this case the segments and of the auxiliary contour (see Figure 2) run along the positive part of a real axis. However, the argument of these segments differs by the angle . Thus, the segment passes through a singular point , and the segment passes through a singular point (). Consequently, the contour must be changed so as to bypass these two points leaving them outside the contour.
We consider the auxiliary contour (Figure 5) consisting of the arc of the circle with the center at the origin of coordinates and radius , , the segment , the arc of the circle with the center in the point and radius defined by the expression (69), the segment , the arc of the circle with the center at the origin of coordinates and radius , the segment , the arc of the circle with the center in the point and radius having the form (61) and the segment . The contour is traversed in a positive direction. As a result, the contour is written in the form
Next, we consider an auxiliary integral
where is defined by (13). By calculating this integral directly we obtain
We immediately note that inside the region bounded by the contour , the function is an analytical function of the complex variable . Consequently, according to the Cauchy integral theorem,
From this it directly follows that
Next, we calculate each integral on the right-hand side of this expression. However, these integrals were already calculated by us earlier. The integral corresponds to the integral . Therefore, using (73) we obtain
The integral corresponds to the integral . Using (74), we obtain
The integral corresponds to the integral . Using (76)
The integral was considered by us when dealing with case 1. Using (42) for the case considered, we obtain
The integral corresponds to the integral . Using (59), we obtain
The integral corresponds to the integral . Therefore, using (60) we obtain
The integral corresponds to the integral . Using (67), we obtain
We get back to the expression (78) and let in it. Note that the integrals , , , do not depend on . Therefore, with such a passage to the limit, they do not change. Taking into account the above expressions for the integrals, the sum (78) takes the form
Performing the passage to the limit in this expression we obtain
Now we consider the Mittag-Leffler (37). Taking into account that in the studied case we get
Directly calculating this integral and representing the complex number in the form we have
Now using here (56) we obtain that the expression may be represented in the form
Now here we make use of the representation (79) for the integral . As a result, we obtain
In the expression derived there is the difference of . We transform this difference. For this we use the representation (57) for the function . As a result, we obtain
From the definition for the function (10) it follows that it is an odd function relative to the parameter
In the total, we have
Here the function was obtained by us earlier and has the form (19). As a result, the expression (81) takes the form
The derived expression proves the theorem completely. ☐
The proved theorem shows that in the representation (3) at the values of parameters satisfying the conditions (32) the integral over the arc of the circle of the contour can be replaced with integration with respect to the segments and of the contour (39). Such a transition makes it possible to replace the contour integral with one improper integral over the real variable of the complex-valued function, which simplifies the use and study of the Mittag-Leffler function. The integral representation of the Mittag-Leffler function obtained in this theorem will be called the integral representation “B”.
Consider in detail the differences between two derived representations. The integral representation “A” (see (7)) of the Mittag-Leffler function consists of two summands. The first summand corresponds to the sum of integrals along the half-lines and and the second summand is the integral along the arc of a circle . Such a representation is not convenient in analytical studies of the Mittag-Leffler function since one has to deal with two integrals. The representation “B” turns out to be more convenient and it consists of one improper integral (see (33)). However, this convenience is fraught with the appearance of constraints imposed on values of the parameters of the function . As it was shown in Theorem 3, the representation “B” is true only at values . If , then it is necessary to use the representation “A” which is the main integral representation for the Mittag-Leffler function. It should also be noted that in the representation “B” with the parameter values and values or the segments or of the auxiliary contour (see Figure 2) pass through the pole . This leads to the need to deform the contour of integration so as to bypass the pole. As a result, in integral representations (34), (35), (36) there are summands describing the pole bypass along the arcs of the circle and . In addition, one has to split the integrals along the half-lines and into parts to exclude the sections corresponding to the integrals along the arcs of the circle and . All this leads to the appearance of additional terms in the corresponding integral representations. It should also be noted that at parameter values in the representation (36) the integral corresponds to the sum of integrals along the segments and the contour of integration (see Figure 5) the integral corresponds to the sum of integrals along and . As one can see, the specified sections of the integration contour go along the positive part of a real axis. As it was shown in the Lemma 1 at values of parameters or and or , and is not an integer, the point is a branch point of the integrand (37). Consequently, the indicated sections of the integration contour will go along different sides of the cut of the complex plane . As a result, and . In case, if , and is an integer, Lemma 1 shows that in this case the point is a regular point. Consequently, the segments and , as well as the half-lines and will go along one straight line of the complex plane, but in different directions, and arcs of the circles and will close. As a result, we obtain and , and the sum of integrals corresponds to the integral over the closed contour and will be equal to the residue in the point . Thus, we come to the condition of Corollary 2.
The formulas of the integral representation obtained in Theorem 3 are rather lengthy. They can be simplified and reduced to a simpler form. The representation “B” takes the simplest form in the case when and coincide, i.e., . We will formulate the result obtained in the form of a corollary.
Corollary3.
For any real , any complex μ and z satisfying the conditions and
the Mittag-Leffler function can be represented in the form:
(1) According to Theorem 3 the representation (33) is true at values and satisfying the conditions (32). In this corollary it is assumed that . In view of this assumption, the conditions (32) take the form if and if . The representation (33) will be written in the form
From this it is clear that it is necessary to consider how the kernel function is transformed at the parameter values specified. This issue was already considered when proving item 1 of Corollary 1 (see (25)). It was shown there that
The first part of the corollary is proved.
(2) Consider the case . Note immediately that the value of an angle cannot exceed . Consequently, this case is implemented only at values . However, the value should be excluded from consideration. In fact, at the segments and of the auxiliary contour (see Figure 2) pass through the pole . Consequently, it is necessary to deform the contour so as to bypass this pole. However, this has already been done by us in item 4 of Theorem 3. By putting the value in the representation (83) we obtain
Using the definition of the function (see (10)) we obtain . Now using this result when calculating the kernel function and introducing the notation we obtain
It has been shown in the paper that when passing from the integral representation formulated in Theorem 1 to integration over real variables, the integral representation of the Mittag-Leffler function can be written in two forms: the representation “A” and “B”. The integral representation “A” was given in Theorem 2 and the representation “B” in Theorem 3. Each of these representations has its own advantages and drawbacks. The representation “A” is true for any complex and at values there is no need to bypass a pole in the point . This excludes the necessity to consider particular cases and greatly simplifies the presentation itself. These facts are the advantages of the representation “A”. The disadvantages of this representation include the fact that it consists of the sum of two integrals: improper and definite. This leads to certain difficulties when working with this representation, since one has to study the behavior of these two integrals. The representation “B” is valid only for the parameter values satisfying the condition . Particular cases also arise at values in which one has to bypass a singular point . That is why, in the problems in which there is a necessity to investigate these particular cases this form of representation is not very convenient. These facts are the disadvantage of the representation “B”. In all other cases the representation “B” consists of one improper integral and turns out to be more convenient in use than the representation “A”. This is the advantage of the representation “B”.
The forms of the representations “A” and “B” obtained in Theorems 2 and 3 are given for the case of four parameters . The parameters and are referred to the Mittag-Leffler function directly, and parameters and describe the integration contour. As a result, in the general case, the integral representations “A” and “B” of the Mittag-Leffler function turn out to be four-parameter. For definiteness, we call this four-parameter description parametrization 1. However, in parametrization 1 though the representations “A” and “B” of the function have a more general form but they are awkward enough. The representations “A” and “B” take a simpler form if . In this case, the integral representations “A” and “B” can be described in three parameters: . We call this parametrization 2. An even simpler view of the form “A” and “B” is adopted if . In this case integral representations of the Mittag-Leffler function can be described in two parameters and . We call this case of parametrization—parametrization 3. Thus, the integral representations of the Mittag-Leffler function in the form “A” and parameterizations 2 and 3 are given in Corollary 1, and the representation “B” in parametrizations 2 and 3 in Corollary 3.
Taking into account of the geometric meaning of the parameters and one can give the geometric interpretation of three introduced parametrizations of the representations “A” and “B”. In fact, the parameters and describe an inclination angle of half-lines and in the contour (see Figure 1). Thus, in case if the half-lines and independently lie in the range of angles , then we have parametrization 1. If these half-lines lie symmetrically in relation to the positive part of a real axis, then we obtain parametrizations 2 and 3. With this, parametrization 3 corresponds to the angle of inclination .
In conclusion it should be pointed out that the representations for the function , formulated in Theorems 2 and 3 are valid for the vaues satisfying the condition . This constraint appears as a result of the use of the proof of Theorem 1 of the integral representation for the gamma function obtained in the work [22] (see Appendix in [22]). The presence of this constraint somewhat narrows the possibilities of using the obtained integral representations of the function . Nevertheless, it is possible to get rid of this constraint and expand the range of admissible values to the entire complex plane. However, this requires additional studies that are beyond the scope of this paper.
Funding
This work was supported by the Russian Foundation for Basic Research (projects No 19-44-730005, 20-07-00655, 18-44-730004) and the Ministry of Science and Higher Education of the Russian Federation (projects No. 0830-2020-0008, RFMEFI60719X0301).
Acknowledgments
The author thanks to M. Yu. Dudikov for translation the article into English.
Conflicts of Interest
The author declares no conflict of interest.
References
Gorenflo, R.; Kilbas, A.A.; Mainardi, F.; Rogosin, S.V. Mittag-Leffler Functions, Related Topics and Applications; Springer Monographs in Mathematics; Springer: Berlin/Heidelberg, Germany, 2014. [Google Scholar] [CrossRef]
Mittag-Leffler, M.G. Une généralization de l’intégrale de Laplace–Abel. Comp. Rend. Acad. Sci.1903, 136, 537–539. [Google Scholar]
Mittag-Leffler, G. Sur la représentation analytique d’une branche uniforme d’une fonction monogène: Cinquième note. Acta Math.1905, 29, 101–181. [Google Scholar] [CrossRef]
Wiman, A. Über den Fundamentalsatz in der Teorie der Funktionen Ea(x). Acta Math.1905, 29, 191–201. [Google Scholar] [CrossRef]
Wiman, A. Über die Nullstellen der Funktionen Ea(x). Acta Math.1905, 29, 217–234. [Google Scholar] [CrossRef]
Humbert, P.; Agarwal, R.P. Sur la fonction de Mittag-Leffler et quelques-unes de ses généralisations. Bull. Sci. Math. II Sér.1953, 77, 180–185. [Google Scholar]
Humbert, P. Quelques résultats rélatifs à la fonction de Mittag-Leffler. C. R. Acad. Sci.1953, 236, 1467–1468. [Google Scholar]
Agarwal, R.P. A propos d’une note de M. Pierre Humbert. C. R. Acad. Sci.1953, 236, 2031–2032. [Google Scholar]
Djrbashian, M.M. On the integral representation of functions continuous on several rays (generalization of the Fourier integral). Izv. Akad. Nauk SSSR Ser. Mat.1954, 18, 427–448. (In Russian) [Google Scholar]
Djrbashian, M.M. On the asymptotic expansion of a function of Mittag-effler type. Akad. Nauk. Armjan. SSR Dokl.1954, 19, 65–72. (In Russian) [Google Scholar]
Djrbashian, M.M. Integral Transfororms and Representation of the functions in the Complex Domain; Nauka. Glav. Red. Fiz.-Mat. Lit: Moscow, Russia, 1966; p. 672. (In Russian) [Google Scholar]
Popov, A.Y.; Sedletskii, A.M. Distribution of roots of Mittag-Leffler functions. J. Math. Sci.2013, 190, 209–409. [Google Scholar] [CrossRef]
Rogosin, S. The role of the Mittag-Leffler function in fractional modeling. Mathematics2015, 3, 368–381. [Google Scholar] [CrossRef]
Gorenflo, R.; Mainardi, F.; Rogosin, S. Mittag-Leffler function: Properties and applications. In Handbook of Fractional Calculus with Applications. Volume 1: Basic Theory; Kochubei, A., Luchko, Y., Eds.; De Gruyter: Berlin, Germany; Boston, MA, USA, 2019; pp. 269–296. [Google Scholar] [CrossRef]
Bateman, H. Higher Transcendental Functions; McGraw-Hill Book Company: New York, NY, USA, 1955; Volume 3, p. 292. [Google Scholar]
Podlubny, I. Fractional Differential Equations; Volume 198, Mathematics in Science and Engineering; Academic Press: Amsterdam, The Netherlands, 1999; p. 340. [Google Scholar]
Gorenflo, R.; Loutchko, J.; Luchko, Y. Computation of the Mittag-Leffler function Eα,β(z) and its derivative. Fract. Calc. Appl. Anal.2002, 5, 491–518. [Google Scholar]
Seybold, H.J.; Hilfer, R. Numerical results for the generalized Mittag-Leffler function. Fract. Calc. Appl. Anal.2005, 8, 129–139. [Google Scholar]
Hilfer, R.; Seybold, H.J. Computation of the generalized Mittag-Leffler function and its inverse in the complex plane. Integral Transform. Spec. Funct.2006, 17, 637–652. [Google Scholar] [CrossRef]
Seybold, H.; Hilfer, R. Numerical Algorithm for Calculating the Generalized Mittag-Leffler Function. SIAM J. Numer. Anal.2009, 47, 69–88. [Google Scholar] [CrossRef]
Parovik, R.I. Calculation specific functions of Mittag-Leffler in the computer mathematics «Maple». Bull. KRASEC. Phys. Math. Sci.2012, 2, 51–61. [Google Scholar] [CrossRef]
Saenko, V.V. An integral representation of the Mittag-Leffler function. arXiv2020, arXiv:2001.09606. [Google Scholar]
Saenko, V.V. Integral representation of the Mittag-Leffler function. In Book of Abstracts 2nd International Conference on Mathematical Modeling in Appliend Sciences, BSU Belgorod-Russia (20–24 August 2019); Debbouche, A., Ed.; Belgorod State University: Belgorod, Russia, 2019; pp. 157–158. [Google Scholar]
Saenko, V.V. Integral Representation of the Fractional Stable Density. J. Math. Sci.2020, 248, 51–66. [Google Scholar] [CrossRef]
Uchaikin, V.V.; Cahoy, D.O.; Sibatov, R.T. Fractional processes: From Poisson to branching one. Int. J. Bifurc. Chaos2008, 18, 2717–2725. [Google Scholar] [CrossRef] [Green Version]
Uchaikin, V.; Sibatov, R.; Uchaikin, D. Memory regeneration phenomenon in dielectrics: The fractional derivative approach. Phys. Scr.2009, T136, 014002. [Google Scholar] [CrossRef]
Saenko, V.V. Singular points of the integral representation of the Mittag-Leffler function. arXiv2020, arXiv:2004.08164. [Google Scholar]
Figure 1.
The contour of integration . The region that corresponds to the condition (2) is shaded in grey.
Figure 1.
The contour of integration . The region that corresponds to the condition (2) is shaded in grey.
Saenko, V.
Two Forms of the Integral Representations of the Mittag-Leffler Function. Mathematics2020, 8, 1101.
https://doi.org/10.3390/math8071101
AMA Style
Saenko V.
Two Forms of the Integral Representations of the Mittag-Leffler Function. Mathematics. 2020; 8(7):1101.
https://doi.org/10.3390/math8071101
Chicago/Turabian Style
Saenko, Viacheslav.
2020. "Two Forms of the Integral Representations of the Mittag-Leffler Function" Mathematics 8, no. 7: 1101.
https://doi.org/10.3390/math8071101
APA Style
Saenko, V.
(2020). Two Forms of the Integral Representations of the Mittag-Leffler Function. Mathematics, 8(7), 1101.
https://doi.org/10.3390/math8071101
Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.
Article Metrics
No
No
Article Access Statistics
For more information on the journal statistics, click here.
Multiple requests from the same IP address are counted as one view.
Saenko, V.
Two Forms of the Integral Representations of the Mittag-Leffler Function. Mathematics2020, 8, 1101.
https://doi.org/10.3390/math8071101
AMA Style
Saenko V.
Two Forms of the Integral Representations of the Mittag-Leffler Function. Mathematics. 2020; 8(7):1101.
https://doi.org/10.3390/math8071101
Chicago/Turabian Style
Saenko, Viacheslav.
2020. "Two Forms of the Integral Representations of the Mittag-Leffler Function" Mathematics 8, no. 7: 1101.
https://doi.org/10.3390/math8071101
APA Style
Saenko, V.
(2020). Two Forms of the Integral Representations of the Mittag-Leffler Function. Mathematics, 8(7), 1101.
https://doi.org/10.3390/math8071101
Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.