Next Article in Journal
Multi-Objective Environmental Economic Dispatch of an Electricity System Considering Integrated Natural Gas Units and Variable Renewable Energy Sources
Next Article in Special Issue
Fractal Generalization of the Scher–Montroll Model for Anomalous Transit-Time Dispersion in Disordered Solids
Previous Article in Journal
On Solutions of an Extended Nonlocal Nonlinear Schrödinger Equation in Plasmas
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Two Forms of the Integral Representations of the Mittag-Leffler Function

by
Viacheslav Saenko
S.P. Kapitsa Research Institute of Technology, Ulyanovsk State University, L. Tolstoy St. 42, Ulyanovsk 432017, Russia
Mathematics 2020, 8(7), 1101; https://doi.org/10.3390/math8071101
Submission received: 4 June 2020 / Revised: 29 June 2020 / Accepted: 30 June 2020 / Published: 5 July 2020
(This article belongs to the Special Issue Fractional Calculus in Anomalous Transport Theory)

Abstract

:
The integral representation of the two-parameter Mittag-Leffler function E ρ , μ ( z ) 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 E ρ , μ ( z ) 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.

1. Introduction

The Mittag-Leffler function is an entire function defined by a power series
E ρ ( z ) = k = 0 z k Γ ( 1 + k / ρ ) , ρ > 0 , z C ,
where Γ ( x ) 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 E ρ ( z ) 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
E ρ , μ ( z ) = k = 0 z k Γ ( μ + k / ρ ) , ρ > 0 , μ C , z C
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 E ρ ( z ) 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 E ρ , μ ( z ) is connected with the classic Mittag-Leffler function E ρ ( z ) by a simple relation E ρ , 1 ( z ) = E ρ ( z ) . For more detailed information on the properties of the function E ρ , μ ( z ) 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 E ρ , μ ( z ) 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 E ρ , μ ( z ) 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 E ρ , μ ( z ) 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 ε 0 . This leads to the fact that the integral representation “B” is valid only for parameter values μ satisfying the condition μ < 1 + 1 / ρ . 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 α = 1 / ρ , β = μ .

2. Integral Representation “A”

The purpose of this paper is to obtain integral representations for the function E ρ , μ ( z ) 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 E ρ , μ ( z ) , obtained in the work [22]. The following theorem was formulated in this work
Theorem 1.
For any real ρ , δ 1 ρ , δ 2 ρ , ϵ satisfying the conditions ρ > 1 / 2 , π 2 ρ < δ 1 ρ min ( π , π / ρ ) , π 2 ρ < δ 2 ρ min ( π , π / ρ ) , ϵ > 0 , any μ C and any z C satisfying the condition
π 2 ρ δ 2 ρ + π < arg z < π 2 ρ + δ 1 ρ + π
the Mittage-Leffler function can be represented in the from
E ρ , μ ( z ) = ρ 2 π i γ ζ exp ( z ζ ) ρ ( z ζ ) ρ ( 1 μ ) ζ 1 d ζ .
where the contour of integration γ ζ has the form (see Figure 1)
γ ζ = S 1 = { ζ : arg ζ = δ 1 ρ π , | ζ | 1 + ϵ } , C ϵ = { ζ : δ 1 ρ π arg ζ δ 2 ρ π , | ζ | = 1 + ϵ } , S 2 = { ζ : arg ζ = δ 2 ρ π , | ζ | 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 S 1 , the arc of the circle C ϵ radius 1 + ϵ and the half-line S 2 . The contour γ ζ is traversed in a positive direction. The parameters δ 1 ρ and δ 2 ρ have the meaning of inclination angles of the half-lines S 1 and S 2 in relation to the contour axis γ ζ . In Figure 1 the contour axis γ ζ coincides with a real axis. The values of the angles δ 1 ρ and δ 2 ρ 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 ζ = r e i φ . As a result, the following theorem is true for the function E ρ , μ ( z ) .
Theorem 2
(The representation “A”). For any real ρ > 1 / 2 , ϵ > 0 and δ 1 ρ , δ 2 ρ satisfying the conditions
π 2 ρ < δ 1 ρ min π , π ρ , π 2 ρ < δ 2 ρ min π , π ρ ,
for any complex μ and any complex z satisfying the condition
π 2 ρ δ 2 ρ + π < arg z < π 2 ρ + δ 1 ρ + π ,
the Mittag-Leffler function can be represented in the form
E ρ , μ ( z ) = 1 + ϵ K ρ , μ ( r , δ 1 ρ , δ 2 ρ , z ) d r + δ 1 ρ π δ 2 ρ π P ρ , μ ( 1 + ϵ , φ , z ) d φ .
Here
K ρ , μ ( r , φ 1 , φ 2 , z ) = ρ 2 π i ( z r ) ρ ( 1 μ ) e i ρ π ( 1 μ ) [ A ρ , μ ( r , φ 2 , φ 1 , z ) A ρ , μ ( r , φ 1 , φ 2 , z ) ] ( r 2 + 2 r cos φ 1 + 1 ) ( r 2 + 2 r cos φ 2 + 1 ) ,
where
A ρ , μ ( r , ω 1 , ω 2 , z ) = exp ( z r ) ρ e i ρ π cos ( ρ ω 1 ) ( r 2 + 2 r cos ω 2 + 1 ) e i η ( r , ω 1 , z ) r + e i ω 1 ,
η ( r , φ , z ) = ( z r ) ρ e i ρ π sin ( ρ φ ) + ρ ( 1 μ ) φ ,
and
P ρ , μ ( r , φ , z ) = ρ 2 π exp ( z r ) ρ cos ( ρ φ ) ( z r ) ρ ( 1 μ ) e i χ ( r , φ , z ) r r e i φ r 2 2 r cos φ + 1 ,
where
χ ( r , φ , z ) = ( z r ) ρ sin ( ρ φ ) + ρ ( 1 μ ) φ .
Proof. 
For convenience we will introduce the notation
ϕ ρ , μ ( ζ , z ) = ρ 2 π i exp ( ζ z ) ρ ( ζ z ) ρ ( 1 μ ) .
As a result, the representation (3) will take the form
E ρ , μ ( z ) = γ ζ ϕ ρ , μ ( ζ , z ) ζ 1 d ζ ,
where the contour of integration γ ζ is defined by (4).
Substituting in (14) the variable of integration ζ = r e i φ and directly calculating this integral we obtain
E ρ , μ ( z ) = S 1 ϕ ρ , μ ( ζ , z ) ζ 1 d ζ + C ϵ ϕ ρ , μ ( ζ , z ) ζ 1 d ζ + S 2 ϕ ρ , μ ( ζ , z ) ζ 1 d ζ = 1 + ϵ ϕ ρ , μ ( r e i φ , z ) r e i φ 1 e i φ d r φ = δ 1 ρ π + i δ 1 ρ π δ 2 ρ π ϕ ρ , μ ( r e i φ , z ) r e i φ 1 r e i φ d φ r = 1 + ϵ + 1 + ϵ ϕ ρ , μ ( r e i φ , z ) r e i φ 1 e i φ d r φ = δ 2 ρ π = 1 + ϵ ϕ ρ , μ ( r e i ( δ 1 ρ π ) , z ) r e i ( δ 1 ρ π ) 1 e i ( δ 1 ρ π ) d r + i δ 1 ρ π δ 2 ρ π ϕ ρ , μ ( ( 1 + ϵ ) e i φ , z ) ( 1 + ϵ ) e i φ 1 ( 1 + ϵ ) e i φ d φ + 1 + ϵ ϕ ρ , μ ( r e i ( δ 2 ρ π ) , z ) r e i ( δ 2 ρ π ) 1 e i ( δ 2 ρ π ) d r = I S 1 + I C ϵ + I S 2 .
We consider the integral I S 1 . Getting rid of the complexity in the denominator we have
I S 1 = 1 + ϵ ϕ ρ , μ r e i ( δ 1 ρ π ) , z e i ( δ 1 ρ π ) r e i ( δ 1 ρ π ) 1 d r = 1 + ϵ ϕ ρ , μ r e i ( δ 1 ρ π ) , z e i ( δ 1 ρ π ) r e i ( δ 1 ρ + π ) 1 r e i ( δ 1 ρ π ) 1 r e i ( δ 1 ρ + π ) 1 d r = 1 + ϵ ϕ ρ , μ r e i ( δ 1 ρ π ) , z r + e i δ 1 ρ r 2 + 2 r cos δ 1 ρ + 1 d r = ρ 2 π i 1 + ϵ exp r z e i ( δ 1 ρ + π ) ρ z r e i ( δ 1 ρ + π ) ρ ( 1 μ ) r + e i δ 1 ρ r 2 + 2 r cos δ 1 ρ + 1 d r = ρ 2 π i 1 + ϵ exp ( z r ) ρ e i ρ π ( cos ( ρ δ 1 ρ ) i sin ( ρ δ 1 ρ ) ) r 2 + 2 r cos δ 1 ρ + 1 ( z r ) ρ ( 1 μ ) e i ρ π ( 1 μ ) e i ρ ( 1 μ ) δ 1 ρ ( r + e i δ 1 ρ ) d r = ρ 2 π i 1 + ϵ exp ( z r ) ρ e i ρ π cos ( ρ δ 1 ρ ) r 2 + 2 r cos δ 1 ρ + 1 ( z r ) ρ ( 1 μ ) e i ρ π ( 1 μ ) e i η ( r , δ 1 ρ , z ) r + e i δ 1 ρ d r ,
where the notation η ( r , φ , z ) = ( z r ) ρ e i ρ π sin ( ρ φ ) + ρ ( 1 μ ) φ was introduced.
Similarly, for the integral I S 2 we get
I S 2 = 1 + ϵ ϕ ρ , μ r e i ( δ 2 ρ π ) , z e i ( δ 2 ρ π ) r e i ( δ 2 ρ π ) 1 d r = 1 + ϵ ϕ ρ , μ r e i ( δ 2 ρ π ) , z r + e i δ 2 ρ r 2 + 2 r cos δ 2 ρ + 1 d r = ρ 2 π i 1 + ϵ exp ( z r ) ρ e i ρ π cos ( ρ δ 2 ρ ) ( z r ) ρ ( 1 μ ) e i ρ π ( 1 μ ) e i η ( r , δ 2 ρ , z ) r + e i δ 2 ρ r 2 + 2 r cos δ 2 ρ + 1 d r .
Summing now the integrals I S 1 and I S 2 we obtain
I S 1 + I S 2 = ρ 2 π i 1 + ϵ exp ( z r ) ρ e i ρ π cos ( ρ δ 1 ρ ) r 2 + 2 r cos δ 1 ρ + 1 ( z r ) ρ ( 1 μ ) e i ρ π ( 1 μ ) e i η ( r , δ 1 ρ , z ) r + e i δ 1 ρ d r + ρ 2 π i 1 + ϵ exp ( z r ) ρ e i ρ π cos ( ρ δ 2 ρ ) ( z r ) ρ ( 1 μ ) e i ρ π ( 1 μ ) e i η ( r , δ 2 ρ , z ) r + e i δ 2 ρ r 2 + 2 r cos δ 2 ρ + 1 d r = 1 + ϵ K ρ , μ ( r , δ 1 ρ , δ 2 ρ , z ) d r ,
where
K ρ , μ ( r , φ 1 , φ 2 , z ) = ρ 2 π i ( z r ) ρ ( 1 μ ) e i ρ π ( 1 μ ) [ A ρ , μ ( r , φ 2 , φ 1 , z ) A ρ , μ ( r , φ 1 , φ 2 , z ) ] ( r 2 + 2 r cos φ 1 + 1 ) ( r 2 + 2 r cos φ 2 + 1 ) , A ρ , μ ( r , ω 1 , ω 2 , z ) = exp ( z r ) ρ e i ρ π cos ( ρ ω 1 ) ( r 2 + 2 r cos ω 2 + 1 ) e i η ( r , ω 1 , z ) r + e i ω 1 .
We consider now the integral I C ϵ . 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 ϕ ρ , μ ( ζ , z ) (see (13)) and in the indices of the exponents we use the Euler formula e i φ = cos φ + i sin φ . As a result, we get
I C ϵ = i δ 1 ρ π δ 2 ρ π ϕ ρ , μ ( ( 1 + ϵ ) e i φ , z ) ( 1 + ϵ ) e i φ 1 ( 1 + ϵ ) e i φ d φ = i δ 1 ρ π δ 2 ρ π ϕ ρ , μ ( ( 1 + ϵ ) e i φ , z ) ( 1 + ϵ ) e i φ ( 1 + ϵ ) e i φ 1 ( 1 + ϵ ) e i φ 1 ( 1 + ϵ ) e i φ 1 d φ = i δ 1 ρ π δ 2 ρ π ϕ ρ , μ ( ( 1 + ϵ ) e i φ , z ) ( 1 + ϵ ) ( 1 + ϵ ) e i φ ( 1 + ϵ ) 2 2 ( 1 + ϵ ) cos φ + 1 d φ = ρ 2 π δ 1 ρ π δ 2 ρ π exp z ( 1 + ϵ ) e i φ ρ z ( 1 + ϵ ) e i φ ρ ( 1 μ ) ( 1 + ϵ ) ( 1 + ϵ ) e i φ ( 1 + ϵ ) 2 2 ( 1 + ϵ ) cos φ + 1 d φ = ρ 2 π δ 1 ρ π δ 2 ρ π exp ( ( 1 + ϵ ) z ) ρ ( cos ( ρ φ ) + i sin ( ρ φ ) ) ( 1 + ϵ ) 2 2 ( 1 + ϵ ) cos φ + 1 ( z ( 1 + ϵ ) ) ρ ( 1 μ ) e i ρ φ ( 1 μ ) ( 1 + ϵ ) ( 1 + ϵ ) e i φ d φ = ρ 2 π δ 1 ρ π δ 2 ρ π exp ( ( 1 + ϵ ) z ) ρ cos ( ρ φ ) ( 1 + ϵ ) 2 2 ( 1 + ϵ ) cos φ + 1 ( z ( 1 + ϵ ) ) ρ ( 1 μ ) e i χ ( ( 1 + ϵ ) , φ , z ) ( 1 + ϵ ) ( 1 + ϵ ) e i φ d φ = δ 1 ρ π δ 2 ρ π P ρ , μ ( 1 + ϵ , φ , z ) d φ ,
where
P ρ , μ ( r , φ , z ) = ρ 2 π exp ( z r ) ρ cos ( ρ φ ) ( z r ) ρ ( 1 μ ) e i χ ( r , φ , z ) r r e i φ r 2 2 r cos φ + 1 , χ ( r , φ , z ) = ( z r ) ρ sin ( ρ φ ) + ρ ( 1 μ ) φ .
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 ρ , μ , δ 1 ρ , δ 2 ρ 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 ρ > 1 / 2 , any real δ 1 ρ , δ 2 ρ 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 S 1 and S 2 of the contour γ ζ and the definite integral corresponds to the integral along the arc of a circle C ϵ . It should be noted that this integral is taken along the arc of a circle of radius 1 + ϵ , where ϵ > 0 . The representation (7) is valid for arbitrary values ρ > 1 / 2 , any μ and any δ 1 ρ and δ 2 ρ that satisfy the condition (5).
However, in general case, for arbitrary values of parameters δ 1 ρ and δ 2 ρ , satisfying the condition (5), the kernel function K ρ , μ ( r , δ 1 ρ , δ 2 ρ , z ) turns to be lengthy. The representation (7) takes the more compact form in case when the half-lines S 1 and S 2 of the contour of integration γ ζ run symmetrically relative to the real axis i.e., when δ 1 ρ = δ 2 ρ = δ ρ . We will formulate the obtained result in the form of a corollary.
Corollary 1.
For any real ϵ > 0 , any complex μ the following integral representations of the Mittag-Leffler function are true:
1. 
at any real ρ > 1 / 2 , any real δ ρ satisfying the condition π 2 ρ < δ ρ min ( π , π ρ ) and any complex z satisfying the condition π 2 ρ δ ρ + π < arg z < π 2 ρ + δ ρ + π
E ρ , μ ( z ) = 1 + ϵ K ρ , μ ( r , δ ρ , z ) d r + δ ρ π δ ρ π P ρ , μ ( 1 + ϵ , φ , z ) d φ ,
where
K ρ , μ ( r , φ , z ) = ρ π ( z r e i π ) ρ ( 1 μ ) exp ( z r e i π ) ρ cos ( ρ φ ) r sin ( η ( r , φ , z ) ) + sin ( η ( r , φ , z ) + φ ) r 2 + 2 r cos φ + 1 ,
η ( r , φ , z ) is defined by (10) and P ρ , μ ( r , φ , z ) has the form (11).
2. 
at any real ρ 1 at δ ρ = π / ρ and any complex z satisfying the condition π 2 ρ + π < arg z < π 2 ρ + π
E ρ , μ ( z ) = 1 + ϵ K ρ , μ ( r , z ) d r + π ρ π π ρ π P ρ , μ ( 1 + ϵ , φ , z ) d φ ,
where
K ρ , μ ( r , z ) = ρ π ( z r e i π ) ρ ( 1 μ ) exp ( z r e i π ) ρ r sin ( π ( 1 μ ) ) + sin ( π ( 1 μ ) + π / ρ ) r 2 + 2 r cos ( π / ρ ) + 1
and P ρ , μ ( r , φ , z ) is defined by (11).
Proof. 
(1) According to Theorem 2, the Mittag-Leffler function can be represented in the form (7). This representation is true for arbitrary δ 1 ρ and δ 2 ρ satisfying the conditions (5). In case if δ 1 ρ = δ 2 ρ = δ ρ the conditions (5) take the form
π 2 ρ < δ ρ min π , π ρ
and the condition (6) can be written in the form
π 2 ρ δ ρ + π < arg z < π 2 ρ + δ ρ + π .
As a result, the Mittag-Leffler function is written in the form
E ρ , μ ( z ) = 1 + ϵ K ρ , μ ( r , δ ρ , δ ρ , z ) d r + δ ρ π δ ρ π P ρ , μ ( 1 + ϵ , φ , z ) d φ .
We consider the integrand of the first integral and denote
K ρ , μ ( r , δ ρ , z ) K ρ , μ ( r , δ ρ , δ ρ , z ) .
Using the definition K ρ , μ ( r , φ 1 , φ 2 , z ) (see (8)) we get
K ρ , μ ( r , δ ρ , z ) = K ρ , μ ( r , δ ρ , δ ρ , z ) = ρ 2 π i ( z r ) ρ ( 1 μ ) e i ρ π ( 1 μ ) A ρ , μ ( r , δ ρ , δ ρ , z ) A ρ , μ ( r , δ ρ , δ ρ , z ) ( r 2 + 2 r cos δ ρ + 1 ) 2 .
Using the definition A ρ , μ ( r , ω 1 , ω 2 , z ) (see (9)) and the fact that the function η ( r , φ , z ) defined by (10) is an odd function according to the variable φ
η ( r , φ , z ) = η ( r , φ , z ) ,
we have
A ρ , μ ( r , δ ρ , δ ρ , z ) A ρ , μ ( r , δ ρ , δ ρ , z ) = exp ( z r ) ρ e i π ρ cos ( ρ δ ρ ) ( r 2 + 2 r cos ( δ ρ ) + 1 ) e i η ( r , δ ρ , z ) ( r + e i δ ρ ) exp ( z r ) ρ e i π ρ cos ( ρ δ ρ ) ( r 2 + 2 r cos δ ρ + 1 ) e i η ( r , δ ρ , z ) ( r + e i δ ρ ) = exp ( z r ) ρ e i π ρ cos ( ρ δ ρ ) ( r 2 + 2 r cos δ ρ + 1 ) r e i η ( r , δ ρ , z ) + e i ( η ( r , δ ρ , z ) + δ ρ ) r e i η ( r , δ ρ , z ) e i ( η ( r , δ ρ , z ) + δ ρ ) = 2 i exp ( z r ) ρ e i π ρ cos ( ρ δ ρ ) ( r 2 + 2 r cos δ ρ + 1 ) r sin ( η ( r , δ ρ , z ) ) + sin ( η ( r , δ ρ , z ) + δ ρ ) .
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 arg z 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 δ 1 ρ = δ 2 ρ = π / ρ . 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 ρ 1 . For the range of values arg z from (23) we get
π 2 ρ + π < arg z < π 2 ρ + π .
Now we consider the representation (18). In the case under consideration it will be written in the form
E ρ , μ ( z ) = 1 + ϵ K ρ , μ r , π / ρ , z d r + π ρ π π ρ π P ρ , μ ( 1 + ϵ , φ , z ) d φ .
We denote
K ρ , μ ( r , z ) K ρ , μ r , π ρ , z .
From (10) it follows that η r , π ρ , z = ( 1 μ ) π . Now using this result in (19) we get (21). ☐
From the proved corollary it follows that if the parameter values δ 1 ρ and δ 2 ρ coincide, then in this the kernel function K ρ , μ ( r , φ 1 , φ 2 , z ) is significantly simplified. Recall that the parameters δ 1 ρ and δ 2 ρ are inclination angles of half-lines S 1 and S 2 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 δ 1 ρ = δ 2 ρ means that half-lines S 1 and S 2 run symmetrically in relation to the real axis. The kernel function K ρ , μ ( r , φ 1 , φ 2 , z ) takes the simplest form in the case δ 1 ρ = δ 2 ρ = π / ρ .
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.
Lemma 1.
For any real ρ > 1 / 2 and any complex values of the parameter μ = μ R + i μ I the integrand of the representation (3)
Φ ρ , μ ( ζ , z ) = ρ 2 π i exp { ( ζ z ) ρ } ( ζ z ) ρ ( 1 μ ) ζ 1 .
relative to the variable ζ has two singular points ζ = 1 and ζ = 0 . The point ζ = 1 is a pole of the first order. The point ζ = 0 is:
1. 
the regular point of the function Φ ρ , μ ( ζ , z ) , with the values of parameters ρ = n , where n = 1 , 2 , 3 , (the positive integer), μ I = 0 and μ R = 1 m 1 / ρ , where m 1 = 0 , 1 , 2 , 3 , (the non-negative integer);
2. 
a pole of the order m 2 , if ρ = n , where n = 1 , 2 , 3 , (the positive integer), μ I = 0 and μ R = 1 + m 2 / ρ , where m 2 = 1 , 2 , 3 , (the positive integer);
3. 
the branch point, for any other values of parameters ρ , μ I , μ R .
The proof of this lemma can be found in the work [27].
We will make the following remark to Corollary 1.
Remark 1.
In Corollary 1 the special case under consideration δ 1 ρ = δ 2 ρ = π / ρ . We will assume that ρ = 1 and study the behavior of the Formula (20) in this case. As a result, we obtain
E 1 , μ ( z ) = 1 + ϵ K 1 , μ ( r , z ) d r + 2 π 0 P 1 , μ ( 1 + ϵ , φ , z ) d φ .
Using (21) for K 1 , μ ( z ) we obtain
K 1 , μ ( r , z ) = 1 π ( z r e i π ) 1 μ exp z r e i π r sin ( π ( 1 μ ) ) + sin ( π ( 1 μ ) + π ) r 2 2 r + 1 = 1 π ( z ) 1 μ e z r sin ( π ( 1 μ ) ) ( r 1 ) ( r 1 ) 2 = 1 π ( z ) 1 μ e z r sin ( π ( 1 μ ) ) r 1 .
From here it is clear, if μ = n , where n = 0 , ± 1 , ± 2 , ± 3 , , then sin ( π ( 1 μ ) ) = 0 . Consequently,
K 1 , ρ ( z ) = 0 , μ = n , 1 π ( z ) 1 μ e z r sin ( π ( 1 μ ) ) r 1 , μ n .
Thus, with integer values of μ the first summand in (28) becomes zero and to calculate the value of E 1 , n ( z ) 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 ρ = 1 . Using the notation (27), the representation (3) takes the form
E 1 , μ ( z ) = γ ζ Φ 1 , μ ( ζ , z ) d ζ ,
where the contour of integration γ ζ , defined by (4), is written in the form
γ ζ = S 1 = { ζ : arg ζ = 2 π , | ζ | 1 + ϵ } , C ϵ = { ζ : 2 π arg ζ 0 , | ζ | = 1 + ϵ } , S 2 = { ζ : arg ζ = 0 , | ζ | 1 + ϵ } .
We represent the complex parameter μ in the form μ = μ R + i μ I and make use of Lemma 1. According to this lemma, the function Φ 1 , μ ( ζ , z ) at values μ I = 0 and μ R = 1 m 1 , where m 1 = 0 , 1 , 2 , 3 , has one singular point ζ = 1 which is a pole of the first order. The point ζ = 0 , in this case, is the regular point. In case, if μ I = 0 and μ R = 1 + m 2 , where m 2 = 1 , 2 , 3 , the function Φ 1 , μ ( ζ , z ) has two singular points: the point ζ = 1 is a pole of the first order and the point ζ = 0 is a pole of the order m 2 . As we can see, in both cases the point ζ = 0 is not a branch point. As a result, in these two cases, the function Φ 1 , μ ( ζ , z ) is the entire function of a complex variable ζ. From here it follows that when μ I = 0 , and μ R = n , where n = 0 , ± 1 , ± 2 , ± 3 , the arc of the circle C ϵ that enters the contour (31) is the closed circle of radius 1 + ϵ . The half-lines S 1 and S 2 pass along the positive part of a real axis in mutually opposite directions. With all other values of the parameter μ (when μ I 0 or μ R n ), according to Lemma 1, the point ζ = 0 is a branch point of the function Φ 1 , μ ( ζ , z ) . In this case, the circle C ϵ of the contour (31) will not close up and half-lines S 1 and S 2 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 C ϵ of the contour (31) is a closed circle and the half-lines S 1 and S 2 run along the positive part of a real axis in mutually opposite directions. Consequently, the sum of the integrals along the half-lines S 1 and S 2 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 S 1 and S 2 . 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.
Corollary 2.
For the values of the parameters ρ = 1 , δ 1 ρ = δ 2 ρ = π , any complex z, satisfying the condition π / 2 < arg z < 3 π / 2 and for integer real values of the parameter μ = n , n = 0 , ± 1 , ± 2 , ± 3 , the Mittag-Leffler function has the form:
1. 
if n 1 (i.e., n = 1 , 0 , 1 , 2 , 3 , ), then
E 1 , n ( z ) = e z z 1 n ,
2. 
if n 2 (i.e., n = 2 , 3 , 4 , ), then
E 1 , n ( z ) = z 1 n e z k = 0 n 2 z k k ! .
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 S 1 and S 2 of the contour γ ζ in the representation (3), a definite integral corresponds to the integral along the arc of the circle C ϵ . As a result, in analytical studies of the function E ρ , μ ( z ) , 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 C ϵ and write the integral representation for the function E ρ , μ ( z ) 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”.
Theorem 3
(Representation “B”). For any real ρ > 1 / 2 and any complex μ satisfying the condition μ < 1 + 1 / ρ for the function E ρ , μ ( z ) , the following integral representations are valid:
1. 
for any real δ 1 ρ , δ 2 ρ satisfying the conditions
π 2 ρ < δ 1 ρ π ρ , π 2 ρ < δ 2 ρ π ρ , i f ρ > 1 , π 2 ρ < δ 1 ρ < π , π 2 ρ < δ 2 ρ < π , i f 1 / 2 < ρ 1 ,
and any complex z satisfying the condition π 2 ρ δ 2 ρ + π < arg z < π 2 ρ + δ 1 ρ + π the Mittag-Leffler function can be represented in the form
E ρ , μ ( z ) = 0 K ρ , μ ( r , δ 1 ρ , δ 2 ρ , z ) d r ,
where K ρ , μ ( r , φ 1 , φ 2 , z ) has the form (8);
2. 
if 1 / 2 < ρ 1 and δ 1 ρ = π , π 2 ρ < δ 2 ρ < π , then for any complex z satisfying the condition π 2 ρ δ 2 ρ + π < arg z < π 2 ρ + 2 π , the Mittag-Leffler function can be represented in the form
E ρ , μ ( z ) = 0 K ρ , μ ( r , δ 2 ρ , z ) d r 0 1 ε 1 K ρ , μ ( r , π , z ) d r 1 + ε 1 K ρ , μ ( r , π , z ) d r + 2 π π P ρ , μ ( ε 1 , ψ , 2 , z ) d ψ ,
where ε 1 is an arbitrary real number satisfying the condition 0 < ε 1 < 1 ,
K ρ , μ ( r , φ , z ) = ρ 2 π i exp ( z r ) ρ e i π ρ cos ( ρ φ ) r 2 + 2 r cos φ + 1 ( z r ) ρ ( 1 μ ) e i [ η ( r , φ , z ) π ρ ( 1 μ ) ] r + e i φ ,
where η ( r , φ , z ) has the form (10) and
P ρ , μ ( τ , ψ , k , z ) = ρ τ 2 π exp ( z r ( τ , ψ ) ) ρ cos ( ρ φ ( τ , ψ , k ) ) ( r ( τ , ψ ) ) 2 2 r ( τ , ψ ) cos ( φ ( τ , ψ , k ) ) + 1 × × ( z r ( τ , ψ ) ) ρ ( 1 μ ) e i [ χ ( τ , ψ , k , z ) + ψ ] r ( τ , ψ ) e i φ ( τ , ψ , k ) 1 ,
where
χ ( τ , ψ , k , z ) = ( z r ( τ , ψ ) ) ρ sin ( ρ φ ( τ , ψ , k ) ) + ρ ( 1 μ ) φ ( τ , ψ , k ) , r ( τ , ψ ) = τ 2 + 2 τ cos ψ + 1 , φ ( τ , ψ , k ) = arctan τ sin ψ τ cos ψ + 1 + k π ;
3. 
if 1 / 2 < ρ 1 and π 2 ρ < δ 1 ρ < π , δ 2 ρ = π , then for any complex z satisfying the condition π 2 ρ < arg z < π 2 ρ + δ 1 ρ + π , the Mittag-Leffler function can be represented in the form
E ρ , μ ( z ) = 0 1 ε 1 K ρ , μ ( r , π , z ) d r 1 + ε 1 K ρ , μ ( r , π , z ) d r + + π 0 P ρ , μ ( ε 1 , ψ , 0 , z ) d ψ 0 K ρ , μ ( r , δ 1 ρ , z ) d r ,
where ε 1 is an arbitrary real number satisfying the condition 0 < ε 1 < 1 ;
4. 
if 1 / 2 < ρ 1 and δ 1 ρ = δ 2 ρ = π , then for any complex z satisfying the condition π 2 ρ < arg z < π 2 ρ + 2 π , the Mittag-Leffler function can be represented in the form
E ρ , μ ( z ) = 0 1 ε 1 K ρ , μ ( r , π , z ) d r 1 + ε 1 K ρ , μ ( r , π , z ) d r + + π 0 P ρ , μ ( ε 1 , ψ , 0 , z ) d ψ + 2 π π P ρ , μ ( ε 1 , ψ , 2 , z ) d ψ ,
where ε 1 is an arbitrary real number satisfying the condition 0 < ε 1 < 1 and K ρ , μ ( r , δ , z ) 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
E ρ , μ ( z ) = γ ζ ϕ ρ , μ ( ζ , z ) ζ 1 d ζ .
The problem consists in calculating this contour integral.
We consider an auxiliary integral
I = Γ ϕ ρ , μ ( ζ , z ) ζ 1 d ζ ,
where the contour Γ (see Figure 2) consists of the arc of the circle C ϵ of radius 1 + ϵ with the center at the origin of coordinates, the segment Γ 2 , the arc of the circle C ε the radius ε with the center at the origin of coordinates and the segment Γ 1 , which are defined in the following way:
Γ = C ϵ = { ζ : δ 1 ρ π arg ζ δ 2 ρ π , | ζ | = 1 + ϵ } , Γ 2 = { ζ : arg ζ = δ 2 ρ π , 1 + ϵ | ζ | ε } , C ε = { ζ : δ 2 ρ π arg ζ δ 1 ρ π , | ζ | = ε } , Γ 1 = { ζ : arg ζ = δ 1 ρ π , ε | ζ | 1 + ϵ } .
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 μ 2 , then there are two poles in the points ζ = 0 and ζ = 1 . If μ 1 , then there is one pole in the point ζ = 1 . In addition, at non-integer values ρ and μ the point ζ = 0 is the branch point. In the case when 1 2 < ρ 1 and δ 1 ρ = π the segment Γ 1 will pass through a singular point ζ = 1 (Figure 2). The similar situation will be at 1 2 < ρ 1 and δ 2 ρ = π . In this case, the segment Γ 2 will pass through the point ζ = 1 . These two cases will be considered separately. The case 1 2 < ρ 1 , δ 1 ρ = δ 2 ρ = π also requires a separate consideration. In this case each of the segments Γ 1 and Γ 2 passes through a singular point ζ = 1 . In view of the foregoing, it is necessary to consider four cases:
  • π 2 ρ < δ 1 ρ π ρ , π 2 ρ < δ 2 ρ π ρ , if ρ > 1 , π 2 ρ < δ 1 ρ < π , π 2 ρ < δ 2 ρ < π , if 1 2 < ρ 1 ,
  • 1 2 < ρ 1 , δ 1 ρ = π , π 2 ρ < δ 2 ρ < π ,
  • 1 2 < ρ 1 , δ 2 ρ = π , π 2 ρ < δ 1 ρ < π ,
  • 1 2 < ρ 1 , δ 1 ρ = π , δ 2 ρ = π .
Case 1. We consider the first case at the beginning. Directly calculating the integral I and substituting the variable of integration ζ = r e i φ we obtain
I = C ϵ ϕ ρ , μ ( ζ , z ) ζ 1 d ζ + Γ 2 ϕ ρ , μ ( ζ , z ) ζ 1 d ζ + C ε ϕ ρ , μ ( ζ , z ) ζ 1 d ζ + Γ 1 ϕ ρ , μ ( ζ , z ) ζ 1 d ζ = = i δ 1 ρ π δ 2 ρ π ϕ ρ , μ ( r e i φ , z ) r e i φ 1 r e i φ d φ r = 1 + ϵ + 1 + ϵ ε ϕ ρ , μ ( r e i φ , z ) r e i φ 1 e i φ d r φ = δ 2 ρ π + i δ 2 ρ π δ 1 ρ π ϕ ρ , μ ( r e i φ , z ) r e i φ 1 r e i φ d φ r = ε + ε 1 + ϵ ϕ ρ , μ ( r e i φ , z ) r e i φ 1 e i φ d r φ = δ 1 ρ π = i δ 1 ρ π δ 2 ρ π ϕ ρ , μ ( ( 1 + ϵ ) e i φ , z ) ( 1 + ϵ ) e i φ 1 ( 1 + ϵ ) e i φ d φ + 1 + ϵ ε ϕ ρ , μ ( r e i ( δ 2 ρ π ) , z ) r e i ( δ 2 ρ π ) 1 e i ( δ 2 ρ π ) d r + i δ 2 ρ π δ 1 ρ π ϕ ρ , μ ( ε e i φ , z ) ε e i φ 1 ε e i φ d φ + ε 1 + ϵ ϕ ρ , μ ( r e i ( δ 1 ρ π ) , z ) r e i ( δ 1 ρ π ) 1 e i ( δ 1 ρ π ) d r = I C ϵ + I Γ 2 + I C ε + I Γ 1 .
Now we let ε 0 in this expression and we will study the behavior of I C ε , I Γ 1 and I Γ 2 . We consider the integral I C ε at the beginning. For this integral the relation is true
lim ε 0 | I C ε | lim ε 0 δ 2 ρ π δ 1 ρ π ϕ ρ , μ ( ε e i φ , z ) ε e i φ 1 ε e i φ | d φ | .
We consider the integrand. To get rid of the complexity in the denominator we multiply and divide the integrand by ε e i φ 1 . We also represent z = | z | e i arg z and μ = μ R + i μ I , and for a power function we will use the representation ξ a = exp { a ln ξ } , ξ > 0 . As a result, we obtain
lim ε 0 ϕ ρ , μ ( ε e i φ , z ) ε e i φ 1 ε e i φ = ρ 2 π lim ε 0 exp z ε e i φ ρ z ε e i φ ρ ( 1 μ ) ε e i φ ε e i φ 1 ε e i φ 1 ε e i φ 1 = ρ 2 π lim ε 0 exp ( | z | ε ) ρ e i ρ ( φ + arg z ) ( | z | ε ) ρ ( 1 μ R i μ I ) e i ρ ( 1 μ R i μ I ) ( φ + arg z ) ε ε e i φ ε 2 2 ε cos φ + 1 = ρ 2 π lim ε 0 1 ε 2 2 ε cos φ + 1 exp ( | z | ε ) ρ e i ρ ( φ + arg z ) + ρ ( 1 μ R i μ I ) ln ( | z | ε ) + + i ρ ( 1 μ R i μ I ) ( φ + arg z ) + ln ε ( ε e i φ ) = ρ 2 π lim ε 0 1 ε 2 2 ε cos φ + 1 exp ( | z | ε ) ρ cos ( ρ ( φ + arg z ) ) + ρ ( 1 μ R ) ln ( | z | ε ) + ρ ( φ + arg z ) μ I + ln ε + i ( | z | ε ) ρ sin ( ρ ( φ + arg z ) ) + ρ ( 1 μ R ) ( φ + arg z ) ρ ln ( | z | ε ) μ I ε e i φ = ρ 2 π lim ε 0 e A + i B ( ε e i φ ) ε 2 2 ε cos φ + 1 ,
where
A = ( | z | ε ) ρ cos ( ρ ( φ + arg z ) ) + ρ ( 1 μ R ) ln ( | z | ε ) + ρ ( φ + arg z ) μ I + ln ε , B = ( | z | ε ) ρ sin ( ρ ( φ + arg z ) ) + ρ ( 1 μ R ) ( φ + arg z ) ρ μ I ln ( | z | ε ) .
For the numerator we have the estimate
e A + i B ( ε e i φ ) e A + i B ε + e A + i B + i φ = e A ( ε + 1 ) .
Thus,
ρ 2 π lim ε 0 e A + i B ( ε e i φ ) ε 2 2 ε cos φ + 1 ρ 2 π lim ε 0 e A ( ε + 1 ) ε 2 2 ε cos φ + 1 = ρ 2 π lim ε 0 1 ε 2 2 ε cos φ + 1 exp ( | z | ε ) ρ cos ( ρ ( φ + arg z ) ) + ρ ( 1 μ R ) ln ( | z | ε ) + ρ ( φ + arg z ) μ I + ln ε + ln ( ε + 1 ) = ρ 2 π lim ε 0 1 ε 2 2 ε cos φ + 1 exp ( | z | ε ) ρ cos ( ρ ( φ + arg z ) ) + + ( ρ ( 1 μ R ) + 1 ) ln ε + ln ( ε + 1 ) + ρ ( 1 μ R ) ln ( | z | ) + ρ μ I ( φ + arg z ) = 0 , μ R < 1 + 1 ρ , ρ 2 π | z | e ρ ( φ + arg z ) μ I , μ R = 1 + 1 ρ .
From this it follows that
lim ε 0 I C ε = 0 , if μ R < 1 + 1 / ρ .
Now we consider the behavior of the integral I Γ 1 at ε 0 . We have the estimate
lim ε 0 ε 1 + ϵ ϕ ρ , μ ( r e i ( δ 1 ρ π ) , z ) r e i ( δ 1 ρ π ) 1 e i ( δ 1 ρ π ) d r lim ε 0 ε 1 + ϵ ϕ ρ , μ ( r e i ( δ 1 ρ π ) , z ) r e i ( δ 1 ρ π ) 1 e i ( δ 1 ρ π ) | d r | = 0 1 + ϵ ϕ ρ , μ ( r e i ( δ 1 ρ π ) , z ) r e i ( δ 1 ρ π ) 1 e i ( δ 1 ρ π ) | d r | .
Consequently, it is necessary to study the behavior of the integrand at r 0 . For convenience, we introduce the notation φ 1 = δ 1 ρ π . 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 ( r e i φ 1 1 ) and represent z = | z | e i arg z , μ = μ R + μ I . For a power function we will make use of the representation ξ a = exp { a ln ξ } . Using (13) we get
lim r 0 ϕ ρ , μ ( r e i ( δ 1 ρ π ) , z ) r e i ( δ 1 ρ π ) 1 e i ( δ 1 ρ π ) = ρ 2 π lim r 0 exp z r e i φ 1 ρ z r e i φ 1 ρ ( 1 μ ) e i φ 1 ( r e i φ 1 1 ) ( r e i φ 1 1 ) ( r e i φ 1 1 ) = ρ 2 π lim r 0 1 r 2 2 r cos φ 1 + 1 exp ( | z | r ) ρ e i ρ ( φ 1 + arg z ) ( | z | r ) ρ ( 1 μ R + i μ I ) e i ρ ( 1 μ R i μ I ) ( φ 1 + arg z ) ( r e i φ 1 ) = ρ 2 π lim r 0 1 r 2 2 r cos φ 1 + 1 exp ( | z | r ) ρ cos ( ρ ( φ 1 + arg z ) ) + ρ μ I ( φ 1 + arg z ) + ρ ( 1 μ R ) ln ( | z | r ) + i ( | z | r ) ρ sin ( ρ ( φ 1 + arg z ) ) + ρ ( 1 μ R ) ( φ 1 + arg z ) ρ μ I ln ( | z | r ) ( r e i φ 1 ) = ρ 2 π lim r 0 e A 1 + i B 1 ( r e i φ 1 ) r 2 2 r cos φ 1 + 1 ,
where
A 1 = ( | z | r ) ρ cos ( ρ ( φ 1 + arg z ) ) + ρ μ I ( φ 1 + arg z ) + ρ ( 1 μ R ) ln ( | z | r ) , B 1 = ( | z | r ) ρ sin ( ρ ( φ 1 + arg z ) ) + ρ ( 1 μ R ) ( φ 1 + arg z ) ρ μ I ln ( | z | r ) .
For the numerator we have the estimate
e A 1 + i B 1 ( r e i φ 1 ) r e A 1 + i B 1 + e A 1 + i B 1 + i φ 1 = e A 1 ( r + 1 ) .
Consequently, we obtain
ρ 2 π lim r 0 e A 1 + i B 1 ( r e i φ 1 ) r 2 2 r cos φ 1 + 1 ρ 2 π lim r 0 e A 1 ( r + 1 ) r 2 2 r cos φ 1 + 1 = ρ 2 π lim r 0 1 r 2 2 r cos φ 1 + 1 exp ( | z | r ) ρ cos ( ρ ( φ 1 + arg z ) ) + ρ μ I ( φ 1 + arg z ) + ρ ( 1 μ R ) ln ( | z | r ) + ln ( r + 1 ) = 0 , μ R < 1 , ρ 2 π e ρ ( φ 1 + arg z ) μ I , μ R = 1 .
We get the similar result for the integral I Γ 2
lim ε 0 ϕ ρ , μ ( r e i ( δ 2 ρ π ) , z ) r e i ( δ 2 ρ π ) 1 e i ( δ 2 ρ π ) = 0 , μ R < 1 , ρ 2 π e ρ ( δ 2 ρ π + arg z ) μ , μ R = 1 .
Thus, at μ 1 the limits lim ε 0 I Γ 1 and lim ε 0 I Γ 2 will converge to the corresponding definite integrals
lim ε 0 I Γ 1 = 0 1 + ϵ ϕ ρ , μ ( r e i ( δ 1 ρ π ) , z ) r e i ( δ 1 ρ π ) 1 e i ( δ 1 ρ π ) d r , lim ε 0 I Γ 2 = 0 1 + ϵ ϕ ρ , μ ( r e i ( δ 2 ρ π ) , z ) r e i ( δ 2 ρ π ) 1 e i ( δ 2 ρ π ) d r .
We will pay attention to the fact that in the case under consideration when the parameters ρ , δ 1 ρ , δ 2 ρ 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
I = I C ϵ + I Γ 2 + I C ε + I Γ 1 = 0
Letting ε 0 in this expression and taking into consideration (42) we obtain
I C ϵ = lim ε 0 ( I Γ 1 + I Γ 2 ) .
Substituting here the expression for I C ϵ and using (43) we obtain
i δ 1 ρ π δ 2 ρ π ϕ ρ , μ ( ( 1 + ϵ ) e i φ , z ) ( 1 + ϵ ) e i φ 1 ( 1 + ϵ ) e i φ d φ = 0 1 + ϵ ϕ ρ , μ ( r e i ( δ 1 ρ π ) , z ) r e i ( δ 1 ρ π ) 1 e i ( δ 1 ρ π ) d r + 0 1 + ϵ ϕ ρ , μ ( r e i ( δ 2 ρ π ) , z ) r e i ( δ 2 ρ π ) 1 e i ( δ 2 ρ π ) d r , μ R < 1 + 1 ρ .
Here it should be pointed out that at values 1 < μ R < 1 + 1 / ρ the integrals (43) have a singularity in the point r = 0 and at values μ R 1 the singularity disappears.
Now let us return to the Mittag-Leffler function. Directly calculating (37) and replacing the variable of integration ζ = r e i φ we get
E ρ , μ ( z ) = S 1 ϕ ρ , μ ( ζ , z ) ζ 1 d ζ + C ϵ ϕ ρ , μ ( ζ , z ) ζ 1 d ζ + S 2 ϕ ρ , μ ( ζ , z ) ζ 1 d ζ = 1 + ϵ ϕ ρ , μ r e i φ , z r e i φ 1 e i φ d r φ = δ 1 ρ π + i δ 1 ρ π δ 2 ρ π ϕ ρ , μ r e i φ , z r e i φ 1 r e i φ φ r = 1 + ϵ + 1 + ϵ ϕ ρ , μ r e i φ , z r e i φ 1 e i φ d r φ = δ 2 ρ π = 1 + ϵ ϕ ρ , μ r e i ( δ 1 ρ π ) , z r e i ( δ 1 ρ π ) 1 e i ( δ 1 ρ π ) d r + i δ 1 ρ π δ 2 ρ π ϕ ρ , μ ( 1 + ϵ ) e i φ , z ( 1 + ϵ ) e i φ 1 ( 1 + ϵ ) e i φ φ + 1 + ϵ ϕ ρ , μ r e i ( δ 2 ρ π ) , z r e i ( δ 2 ρ π ) 1 e i ( δ 2 ρ π ) d r .
Using here (45) we obtain
E ρ , μ ( z ) = 0 ϕ ρ , μ r e i ( δ 2 ρ π ) , z r e i ( δ 2 ρ π ) 1 e i ( δ 2 ρ π ) ϕ ρ , μ r e i ( δ 1 ρ π ) , z r e i ( δ 1 ρ π ) 1 e i ( δ 1 ρ π ) d r , μ R < 1 + 1 ρ .
We consider the first summand in this expression. Getting rid of the complexity in the denominator and using (13) we get
ϕ ρ , μ r e i ( δ 2 ρ π ) , z r e i ( δ 2 ρ π ) 1 e i ( δ 2 ρ π ) = ϕ ρ , μ r e i ( δ 2 ρ π ) , z e i ( δ 2 ρ π ) r e i ( δ 2 ρ π ) 1 r e i ( δ 2 ρ π ) 1 r e i ( δ 2 ρ π ) 1 = ϕ ρ , μ r e i ( δ 2 ρ π ) , z r e i ( δ 2 ρ π ) r 2 + 2 r cos δ 2 ρ + 1 = ρ 2 π i exp z r e i ( δ 2 ρ π ) ρ z r e i ( δ 2 ρ π ) ρ ( 1 μ ) r + e i δ 2 ρ r 2 + 2 r cos δ 2 ρ + 1 = ρ 2 π i exp ( z r ) ρ e i ρ π ( cos ( ρ δ 2 ρ ) + i sin ( ρ δ 2 ρ ) ) ( z r ) ρ ( 1 μ ) e i ρ ( 1 μ ) ( δ 2 ρ π ) r + e i δ 2 ρ r 2 + 2 r cos δ 2 ρ + 1 = ρ 2 π i exp ( z r ) ρ e i ρ π cos ( ρ δ 2 ρ ) ( z r ) ρ ( 1 μ ) e i ρ ( 1 μ ) π e i η ( r , δ 2 ρ , z ) r + e i δ 2 ρ r 2 + 2 r cos δ 2 ρ + 1 ,
where η ( r , φ , z ) has the form (10). Similarly, for the second summand in (46) we have
ϕ ρ , μ r e i ( δ 2 ρ π ) , z r e i ( δ 2 ρ π ) 1 e i ( δ 2 ρ π ) = = ρ 2 π i exp ( z r ) ρ e i ρ π cos ( ρ δ 1 ρ ) ( z r ) ρ ( 1 μ ) e i ρ ( 1 μ ) π e i η ( r , δ 1 ρ , z ) r + e i δ 1 ρ r 2 + 2 r cos δ 1 ρ + 1
Now substituting (47) and (48) in (46) and making simple transformations we finally obtain
E ρ , μ ( z ) = 0 K ρ , μ ( r , δ 1 ρ , δ 2 ρ , z ) d r ,
where K ρ , μ ( r , φ 1 , φ 2 , z ) is defined by the expression (8). The values δ 1 ρ and δ 2 ρ satisfy the conditions (40). The condition for the value arg z 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 1 2 < ρ 1 , δ 1 ρ = π and π / ( 2 ρ ) < δ 2 ρ < π . As it was mentioned earlier, in this case the segment Γ 1 of the auxiliary contour Γ (see Figure 2) will go through the singular point ζ = 1 . 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 C ϵ with the center at the origin of coordinates and radius 1 + ϵ , ϵ > 0 (see Figure 3), the segment A B , the circle C ε with the center at the origin of coordinates and radius ε satisfying the condition 0 < ε < 1 , the segment C D , the arc of the circle C 1 ε 1 with the center in the point ζ = 1 ( arg ζ = 2 π ) and radius ε 1 (here ε 1 < 1 ε and ε 1 < ϵ ) and the segment E F .
For further study we need to parametrize this contour. As a result, in the plane ζ the arc of the circle C ϵ can be written in the form C ϵ = { ζ : 2 π arg ζ δ 2 ρ π , | ζ | = 1 + ϵ } . The segment A B is written in the form A B = { ζ : arg ζ = δ 2 ρ π , ε | ζ | 1 + ϵ } . The arc of the circle C ε has the form C ε = { ζ : δ 2 ρ π arg ζ 2 π , | ζ | = ε } . The segments C D and E F are written in the form C D = { ζ : arg ζ = 2 π , ε | ζ | 1 ε 1 } and E F = { ζ : arg ζ = 2 π , 1 + ε 1 | ζ | 1 + ϵ } .
To parametrize the arc of the circle C 1 ε 1 , we consider the mapping u = ζ 1 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 ζ = 1 in the plane ζ is mapped onto the point u = 0 of the plane u. Thus, the circle C 1 ε 1 with the center in the point ζ = 1 and the radius ε 1 < 1 of the plane ζ will be mapped on the circle C 0 ε 1 u with the center in the point u = 0 and the radius ε 1 on the plane u. Thus, to parametrize the arc of the circle C 1 ε 1 in the plane ζ it is enough to parametrize the arc of the circle C 0 ε 1 u 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 ζ = u + 1 .
We will represent the complex number u in the form u = τ e i ψ . Then, an arc of a circle C 0 ε 1 u in the complex plane u has the form C 0 ε 1 u = { u : 2 π ψ π , τ = ε 1 } , where ε 1 1 . Now we map the complex plane u onto the plane ζ . We have ζ = τ e i ψ + 1 = τ cos ψ + i τ sin ψ + 1 . From this we obtain
| ζ | = τ 2 + 2 τ cos ψ + 1 , arg ζ = arctan τ sin ψ τ cos ψ + 1 .
Here it is necessary to note that arctan ( x ) is a multivalued function. The principal branch of this function takes values in the interval [ π / 2 , π / 2 ] . However, as one can see from the definition of the auxiliary contour Γ , the center of the arc of the circle C 1 ε 1 lies in the point ζ = 1 in which arg ζ = 2 π . Therefore, it is necessary to choose the required branch in arctan ( x ) in such a way that the mapping (49) could map an arc of a circle C 0 ε 1 u of the plane u on to an arc of a circle C 1 ε 1 of the plane ζ with the center in the point ζ = 1 and arg ζ = 2 π . As a result, the mapping (49) will take the form
| ζ | = τ 2 + 2 τ cos ψ + 1 , arg ζ = arctan τ sin ψ τ cos ψ + 1 + k π .
For the arc of the circle C 1 ε 1 we obtain k = 2 . Here it should be pointed out that these formulas produce the mapping of the circle C 0 ε 1 u on to the circle C 1 ε 1 only in the case τ 1 . If τ > 1 , then these formulas will not be true.
Thus, in view of the remarks made the auxiliary contour Γ can be represented in the form
Γ = C ϵ = { ζ : 2 π arg ζ δ 2 ρ π , | ζ | = 1 + ϵ } , A B = { ζ : arg ζ = δ 2 ρ π , ε | ζ | 1 + ϵ } , C ε = { ζ : δ 2 ρ π arg ζ 2 π , | ζ | = ε } , C D = { ζ : arg ζ = 2 π , ε | ζ | 1 ε 1 } , C 1 ε 1 = ζ : arg ζ = arctan τ sin ψ τ cos ψ + 1 + k π , | ζ | = τ 2 + 2 τ cos ψ + 1 , π ψ 2 π , τ = ε 1 , k = 2 , E F = { ζ : arg ζ = 2 π , 1 + ε 1 | ζ | 1 + ϵ } .
The contour is traversed in a positive direction.
Next, we consider an auxiliary integral
I = Γ ϕ ρ , μ ( ζ , z ) ζ 1 d ζ .
By calculating this integral we have
I = Γ ϕ ρ , μ ( ζ , z ) ζ 1 d ζ = C ϵ ϕ ρ , μ ( ζ , z ) ζ 1 d ζ + A B ϕ ρ , μ ( ζ , z ) ζ 1 d ζ + C ε ϕ ρ , μ ( ζ , z ) ζ 1 d ζ + C D ϕ ρ , μ ( ζ , z ) ζ 1 d ζ + C 1 ε 1 ϕ ρ , μ ( ζ , z ) ζ 1 d ζ + E F ϕ ρ , μ ( ζ , z ) ζ 1 d ζ = I C ϵ + I A B + I C ε + I C D + I C 1 ε 1 + I E F .
One should pay attention to the fact that inside the region limited by the contour Γ the integrand ϕ ρ , μ ( ζ , z ) / ( ζ 1 ) of the integral I is an analytical function of the complex variable ζ . Consequently, according to the Cauchy integral theorem
I C ϵ + I A B + I C ε + I C D + I C 1 ε 1 + I E F = 0 .
From here we get that
I C ϵ = I A B I C ε I C D I C 1 ε 1 I E F .
Thus, it is possible to substitute integration on the arc of the circle C ϵ by integration on the remaining parts of the contour Γ .
We consider each integral in the right part of this expression. The integral I C ε has already been considered by us earlier in this lemma. Using (42) we obtain that in the case considered
lim ε 0 I C ε = 0 , μ R < 1 + 1 / ρ .
Consider the integral I A B . Representing the complex number ζ in the form ζ = r e i φ and using (51) we obtain
I A B = A B ϕ ρ , μ ( ζ , z ) ζ 1 d ζ = 1 + ϵ ε ϕ ρ , μ r e i φ , z r e i φ 1 e i φ d r φ = δ 2 ρ π = 1 + ϵ ε ϕ ρ , μ r e i ( δ 2 ρ π ) , z r e i ( δ 2 ρ π ) 1 e i ( δ 2 ρ π ) d r .
We introduce the notation
ϕ ρ , μ r e i ( δ π ) , z r e i ( δ π ) 1 e i ( δ π ) = K ρ , μ ( r , δ , z ) .
Now we transform the function K ρ , μ ( r , δ , z ) . Using the definition for ϕ ρ , μ ( ζ , z ) (see (13)) and getting rid of the complexity in the denominator we obtain
K ρ , μ ( r , δ , z ) = ρ 2 π i exp z r e i ( δ π ) ρ z r e i ( δ π ) ρ ( 1 μ ) e i ( δ π ) r e i ( δ π ) 1 r e i ( δ π ) 1 r e i ( δ π ) 1 = ρ 2 π i exp ( z r ) ρ e i ρ π [ cos ( ρ δ ) + i sin ( ρ δ ) ] z r e i π ρ ( 1 μ ) e i ρ ( 1 μ ) δ r e i ( δ π ) r 2 2 r cos ( δ π ) + 1 = ρ 2 π i exp ( z r ) ρ e i ρ π cos ( ρ δ ) z r e i π ρ ( 1 μ ) r 2 + 2 r cos δ + 1 exp i ( z r ) ρ e i π ρ sin ( ρ δ ) + ρ ( 1 μ ) δ r + e i δ
By introducing the notation η ( r , φ , z ) = ( z r ) ρ e i π ρ sin ( ρ δ ) + ρ ( 1 μ ) δ , that coincides with (10) we obtain that K ρ , μ ( r , δ , z ) can be represented in the form
K ρ , μ ( r , δ , z ) = ρ 2 π i exp ( z r ) ρ e i ρ π cos ( ρ δ ) z r e i π ρ ( 1 μ ) e i η ( r , φ , z ) r + e i δ r 2 + 2 r cos δ + 1 .
Thus, for the integral I A B we get
I A B = 1 + ϵ ε K ρ , μ ( r , δ 2 ρ , z ) d r .
Now consider the integral I C D . Representing the complex number ζ in the form ζ = r e i φ and using (51) we get
I C D = C D ϕ ρ , μ ( ζ , z ) ζ 1 d ζ = ε 1 ε 1 ϕ ρ , μ r e i φ , z r e i φ 1 e i φ d r φ = π π = ε 1 ε 1 ϕ ρ , μ r e i ( π π ) , z r e i ( π π ) 1 e i ( π π ) d r .
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 I C D we obtain
I C D = ε 1 ε 1 K ρ , μ ( r , π , z ) d r .
Similarly, using (51) for the integral I E F we obtain
I E F = E F ϕ ρ , μ ( ζ , z ) ζ 1 d ζ = 1 + ε 1 1 + ϵ ϕ ρ , μ r e i φ , z r e i φ 1 e i φ d r φ = π π = 1 + ε 1 1 + ϵ ϕ ρ , μ r e i ( π π ) , z r e i ( π π ) 1 e i ( π π ) d r .
Using now (56) and (57) the integral I E F takes the form
I E F = 1 + ε 1 1 + ϵ K ρ , μ ( r , π , z ) d r .
In the sum (52) it remains to consider the integral I C 1 ε 1 . This integral is taken along the arc of the circle C 1 ε 1 with the center in the point ζ = 1 (at that arg ζ = 2 π ) and radius ε 1 . 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 C 1 ε 1 in the plane ζ can be given in the form (see (51))
C 1 ε 1 = ζ : arg ζ = arctan τ sin ψ τ cos ψ + 1 + k π , | ζ | = τ 2 + 2 τ cos ψ + 1 , π ψ 2 π , τ = ε 1 , k = 2
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 2 π . Thus, in the case under consideration the contour C 1 ε 1 on the complex plane ζ can be represented in the form
ζ = r ( ε 1 , ψ ) e i φ ( ε 1 , ψ , 2 ) ,
where π ψ 2 π , 0 < ε 1 < 1 and
r ( τ , ψ ) = τ 2 + 2 τ cos ψ + 1 , φ ( τ , ψ , k ) = arctan τ sin ψ τ cos ψ + 1 + k π
Representing now in the integral I C 1 ε 1 the complex number ζ in the form (62) we find
I C 1 ε 1 = C 1 ε 1 ϕ ρ , μ ( ζ , z ) ζ 1 d ζ = i ε 1 2 π π ϕ ρ , μ r ( ε 1 , ψ ) e i φ ( ε 1 , ψ , 2 ) , z r ( ε 1 , ψ ) e i φ ( ε 1 , ψ , 2 ) 1 e i ψ d ψ .
One should pay attention that traversing along the arc of the circle C 1 ε 1 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 ψ = 2 π 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
i τ ϕ ρ , μ r ( τ , ψ ) e i φ ( τ , ψ , k ) , z r ( τ , ψ ) e i φ ( τ , ψ , k ) 1 e i ψ = P ρ , μ ( τ , ψ , k , z ) .
We transform this expression using (13) and getting rid of the complexity in the denominator we get
P ρ , μ ( τ , ψ , k , z ) = ρ τ 2 π exp z r ( τ , ψ ) e i φ ( τ , ψ , k ) ρ z r ( τ , ψ ) e i φ ( τ , ψ , k ) ρ ( 1 μ ) e i ψ r ( τ , ψ ) e i φ ( τ , ψ , k ) 1 = ρ τ 2 π exp z r ( τ , ψ ) ρ ( cos ( ρ φ ( τ , ψ , k ) ) + i sin ( ρ φ ( τ , ψ , k ) ) ) r ( τ , ψ ) e i φ ( τ , ψ , k ) 1 r ( τ , ψ ) e i φ ( τ , ψ , k ) 1 × × z r ( τ , ψ ) ρ ( 1 μ ) e i ρ ( 1 μ ) φ ( τ , ψ , k ) r ( τ , ψ ) e i φ ( τ , ψ , k ) 1 e i ψ = ρ τ 2 π exp ( z r ( τ , ψ ) ) ρ cos ( ρ φ ( τ , ψ , k ) ) ( z r ( τ , φ ) ) ρ ( 1 μ ) ( r ( τ , ψ ) ) 2 2 r ( τ , ψ ) cos φ ( τ , ψ , k ) + 1 × × exp i [ ( z r ( τ , ψ ) ) ρ sin ( ρ φ ( τ , ψ , k ) ) + ρ ( 1 μ ) φ ( τ , ψ , k ) ] r ( τ , ψ ) e i φ ( τ , ψ , k ) 1 e i ψ = ρ τ 2 π exp ( z r ( τ , ψ ) ) ρ cos ( ρ φ ( τ , ψ , k ) ) ( z r ( τ , φ ) ) ρ ( 1 μ ) ( r ( τ , ψ ) ) 2 2 r ( τ , ψ ) cos φ ( τ , ψ , k ) + 1 e i [ χ ( τ , ψ , k , z ) + ψ ] r ( τ , ψ ) e i φ ( τ , ψ , k ) 1 ,
where
χ ( τ , ψ , k , z ) = ( z r ( τ , ψ ) ) ρ sin ( ρ φ ( τ , ψ , k ) ) + ρ ( 1 μ ) φ ( τ , ψ , k ) .
As a result the integral (64) takes the form
I C 1 ε 1 = 2 π π P ρ , μ ( ε 1 , ψ , 2 , z ) d ψ .
Now we get back to the expression (53) and let ε 0 in this expression. It is necessary to point out that the integrals I C ϵ , I C 1 ε 1 and I E F do not depend on ε and, consequently, they will not change with such a passage to the limit. In view of this, we have
I C ϵ = lim ε 0 I A B lim ε 0 I C ε lim ε 0 I C D I C 1 ε 1 I E F .
Taking into account (54) and using the expressions (58), (59), (60) and (67) we get
I C ϵ = lim ε 0 1 + ϵ ε K ρ , μ ( r , δ 2 ρ , z ) d r lim ε 0 ε 1 ε 1 K ρ , μ ( r , π , z ) d r + 2 π π P ρ , μ ( ε 1 , ψ , 2 , z ) d ψ 1 + ε 1 1 + ϵ K ρ , μ ( r , π , z ) d r , μ < 1 + 1 / ρ .
From this we obtain
I C ϵ = 1 + ϵ 0 K ρ , μ ( r , δ 2 ρ , z ) d r 0 1 ε 1 K ρ , μ ( r , π , z ) d r + 2 π π P ρ , μ ( ε 1 , ψ , 2 , z ) d ψ 1 + ε 1 1 + ϵ K ρ , μ ( r , π , z ) d r , μ < 1 + 1 / ρ .
Now we get back to the Mittag-Leffler function. Assuming that δ 1 ρ = π the contour of integration γ ζ in (37) takes the form
γ ζ = S 1 = { ζ : arg ζ = 2 π , | ζ | 1 + ϵ } C ϵ = { ζ : 2 π arg ζ δ 2 ρ π , | ζ | = 1 + ϵ } S 2 = { ζ : arg ζ = δ 2 ρ π , | ζ | 1 + ϵ }
By calculating (37) directly and representing a complex number ζ in the form ζ = r e i φ we obtain
E ρ , μ ( z ) = S 1 ϕ ρ , μ ( ζ , z ) ζ 1 d ζ + C ϵ ϕ ρ , μ ( ζ , z ) ζ 1 d ζ + S 1 ϕ ρ , μ ( ζ , z ) ζ 1 d ζ = 1 + ϵ ϕ ρ , μ r e i φ , z r e i φ 1 e i φ d r φ = π π + I C ϵ + 1 + ϵ ϕ ρ , μ r e i φ , z r e i φ 1 e i φ d r φ = δ 2 ρ π = 1 + ϵ ϕ ρ , μ r e i ( π π ) , z r e i ( π π ) 1 e i ( π π ) d r + I C ϵ + 1 + ϵ ϕ ρ , μ r e i ( δ 2 ρ π ) , z r e i ( δ 2 ρ π ) 1 e i ( δ 2 ρ π ) d r
using the notation (56), this expression can be written in the form
E ρ , μ ( z ) = 1 + ϵ K ρ , μ ( r , π , z ) d r + I C ϵ + 1 + ϵ K ρ , μ ( r , δ 2 ρ , z ) d r ,
where the form of the functions K ρ , μ ( r , δ , z ) is defined by (57). Now we will make use of the representation here (68) for the integral I C ϵ . As a result, we obtain
E ρ , μ ( z ) = 1 + ϵ K ρ , μ ( r , π , z ) d r 1 + ϵ 0 K ρ , μ ( r , δ 2 ρ , z ) d r 0 1 ε 1 K ρ , μ ( r , π , z ) d r + 2 π π P ρ , μ ( ε 1 , ψ , 2 , z ) d ψ 1 + ε 1 1 + ϵ K ρ , μ ( r , π , z ) d r + 1 + ϵ K ρ , μ ( r , δ 2 ρ , z ) d r = 0 K ρ , μ ( r , δ 2 ρ , z ) d r 0 1 ε 1 K ρ , μ ( r , π , z ) d r + 2 π π P ρ , μ ( ε 1 , ψ , 2 , z ) d ψ 1 + ε 1 K ρ , μ ( r , π , z ) d r ,
where μ < 1 + 1 / ρ . It remains to consider how to change the condition (2) in this case. Since, in the case considered 1 2 < ρ 1 , δ 1 ρ = π , π / ( 2 ρ ) < δ 2 ρ π , then the condition (2) takes the form π 2 ρ δ 2 ρ + π < arg z < π 2 ρ + 2 π . Thus, we have obtained the statement of the theorem for the second case.
Case 3. Consider now the case 1 / 2 < ρ 1 , π / ( 2 ρ ) < δ 1 ρ π , δ 2 ρ = π . In this case the segment Γ 2 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 ζ = 1 . It should be pointed out that in this case an argument of points of the segment is equal to arg ζ = 0 . Consequently, in this case arg ζ = 0 . 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 C ϵ = { ζ : δ 1 ρ π arg ζ 0 , | ζ | = 1 + ϵ } with the center in the point ζ = 0 and radius 1 + ϵ , ϵ > 0 , the segment A B = { ζ : arg ζ = 0 , 1 + ϵ | ζ | 1 + ε 1 } , the arc of the circle C 1 ε 1 with the center in the point ζ = 1 , ( arg ζ = 0 ) and radius ε 1 , ( 0 < ε 1 < 1 , ε 1 < ϵ ) , the segment C D = { ζ : arg ζ = 0 , 1 ε 1 | ζ | ε } , the arc of the circle C ε = { ζ : 0 arg ζ δ 1 ρ π , | ζ | = ε } with the center at the origin of coordinates and radius ε > 0 and ε < 1 ε 1 and the segment E F = { ζ : arg ζ = δ 1 ρ π , ε | ζ | 1 + ϵ } .
As one can see, only the arc of the circle C 1 ε 1 remains non-parametric. For its parametrization we will fulfill the procedure as in the previous case. Considering the conformal mapping u = ζ 1 and representing the complex number u in the form u = τ e i ψ we get that the arc of the circle C 1 ε 1 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 ψ = 0 to the point B. It should be pointed out that the center of a circle lies in the point ζ = 1 with arg ζ = 0 . Consequently, in the formulas (50) one should choose the principal branch arctan ( x ) , i.e., to take k = 0 . As a result, an arc of a circle C 1 ε 1 can be represented in the form
C 1 ε 1 = ζ : arg ζ = arctan τ sin ψ τ cos ψ + 1 + k π , | ζ | = τ 2 + 2 τ cos ψ + 1 , 0 ψ π , τ = ε 1 , k = 0 .
In view of the foregoing, the contour Γ is written in the form
Γ = C ϵ = { ζ : δ 1 ρ π arg ζ 0 , | ζ | = 1 + ϵ } , A B = { ζ : arg ζ = 0 , 1 + ϵ | ζ | 1 + ε 1 } , C 1 ε 1 = ζ : arg ζ = arctan τ sin ψ τ cos ψ 1 + k π , | ζ | = τ 2 + 2 τ cos ψ + 1 , 0 arg ζ π , τ = ε 1 , k = 0 , C D = { ζ : arg ζ = 0 , 1 ε 1 | ζ | ε } , C ε = { ζ : 0 arg ζ δ 1 ρ π , | ζ | = ε } , E F = { ζ : arg ζ = δ 1 ρ π , ε | ζ | 1 + ϵ } .
The contour is traversed in a positive direction.
We consider now an auxiliary integral
I = Γ ϕ ρ , μ ( ζ , z ) ζ 1 d ζ .
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
I = C ϵ ϕ ρ , μ ( ζ , z ) ζ 1 d ζ + A B ϕ ρ , μ ( ζ , z ) ζ 1 d ζ + C 1 ε 1 ϕ 1 , μ ( ζ , z ) ζ 1 d ζ + C D ϕ ρ , μ ( ζ , z ) ζ 1 d ζ + C ε ϕ ρ , μ ( ζ , z ) ζ 1 d ζ + E F ϕ ρ , μ ( ζ , z ) ζ 1 d ζ = I C ϵ + I A B + I C 1 ε 1 + I C D + I C ε + I E F .
It is clear from this expression that inside the region limited by the contour Γ the integrand of the integral I is an analytical and continuous function of a variable ζ . Consequently, according to the Cauchy theorem
I C ϵ + I A B + I C 1 ε 1 + I C D + I C ε + I E F = 0 .
From this we obtain that
I C ϵ = I A B I C 1 ε 1 I C D I C ε I E F .
We consider now each of the integrals in the right part separately.
The integral I C ε was considered by us earlier. Using (42), we obtain
lim ε 0 I C ε = 0 , μ < 1 + 1 / ρ .
Now consider the integral I A B . Representing ζ in the form ζ = r e i φ we find
I A B = 1 + ϵ 1 + ε 1 ϕ ρ , μ ( r e i φ , z ) r e i φ 1 e i φ d r π π = 1 + ϵ 1 + ε 1 ϕ ρ , μ ( r e i ( π π ) , z ) r e i ( π π ) 1 e i ( π π ) d r .
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
I A B = 1 + ϵ 1 + ε 1 K ρ , μ ( r , π , z ) d r .
In the same way, using (70), (56) and (57) for the integrals I C D and I E F we find
I C D = 1 ε 1 ε K ρ , μ ( r , π , z ) d r , I E F = ε 1 + ϵ K ρ , μ ( r , π , z ) d r .
It remains to consider the integral I C 1 ε 1 . This integral is taken along the arc of the C 1 ε 1 with the center in the point ζ = 1 , where arg ζ = 0 . This means that the contour is traversed from the point B of the complex plane ζ in which arg ζ B = 0 to the point C in which arg ζ C = 0 . 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 C 1 ε 1 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 C 1 ε 1 can be represented in the form
ζ = r ( ε 1 , ψ ) e i φ ( ε 1 , ψ , 0 ) ,
where r ( τ , ψ ) and φ ( τ , ψ , k ) have the form (63).
Using now the representation (75), in the integral I C 1 ε 1 we obtain
I C 1 ε 1 = 0 π ϕ ρ , μ r ( τ , ψ ) e i φ ( τ , ψ , 0 ) , z i τ e i ψ r ( τ , ψ ) e i φ ( τ , ψ , 0 ) 1 d ψ τ = ε 1 = i ε 1 π 0 ϕ ρ , μ r ( ε 1 , ψ ) e i φ ( ε 1 , ψ , 0 ) , z e i ψ r ( τ , ψ ) e i φ ( τ , ψ , 0 ) 1 d ψ .
Comparing this expression with (64) it is clear that these two integrals are similar. Therefore, using (65) and (66) for the integral I C 1 ε 1 we get
I C 1 ε 1 = π 0 P ρ , μ ( ε 1 , ψ , 0 , z ) d ψ .
We return to the expression (71) and let ε 0 in this expression. Note that the integrals I C ϵ , I A B , I C 1 ε 1 do not depend on ε and, consequently, with such a passage to the limit they will not change. As a result we obtain
I C ϵ = I A B I C 1 ε 1 lim ε 0 I C D lim ε 0 I C ε lim ε 0 I E F .
Using here the expressions (72)–(74), (76) we obtain
I C ϵ = 1 + ϵ 1 + ε 1 K ρ , μ ( r , π , z ) d r + π 0 P ρ , μ ( ε 1 , ψ , 0 , z ) d ψ lim ε 0 1 ε 1 ε K ρ , μ ( r , π , z ) d r lim ε 0 ε 1 + ϵ K ρ , μ ( r , π , z ) d r , μ < 1 + 1 / ρ .
As a result, performing the passage to the limit we have
I C ϵ = 1 + ϵ 1 + ε 1 K ρ , μ ( r , π , z ) d r + π 0 P ρ , μ ( ε 1 , ψ , 0 , z ) d ψ 1 ε 1 0 K ρ , μ ( r , π , z ) d r 0 1 + ϵ K ρ , μ ( r , π , z ) d r , μ < 1 + 1 / ρ .
We return now to the Mittag-Leffler function. Assuming that δ 2 ρ = π the contour of integration γ ζ in (37) takes the form
γ ζ = S 1 = { ζ : arg ζ = δ 1 ρ π , | ζ | 1 + ϵ } C ϵ = { ζ : δ 1 ρ π arg ζ 0 , | ζ | = 1 + ϵ } S 2 = { ζ : arg ζ = 0 , | ζ | 1 + ϵ } ,
and the condition (2) will be written in the form π 2 ρ < arg z < π 2 ρ + δ 1 ρ + π .
By direct calculating (37) and representing the complex number ζ in the form ζ = r e i φ we get
E ρ , μ ( z ) = S 1 ϕ ρ , μ ( ζ , z ) ζ 1 d ζ + C ϵ ϕ ρ , μ ( ζ , z ) ζ 1 d ζ + S 1 ϕ ρ , μ ( ζ , z ) ζ 1 d ζ = 1 + ϵ ϕ ρ , μ r e i φ , z r e i φ 1 e i φ d r φ = δ 1 ρ π + I C ϵ + 1 + ϵ ϕ ρ , μ r e i φ , z r e i φ 1 e i φ d r φ = π π = 1 + ϵ ϕ ρ , μ r e i ( δ 1 ρ π ) , z r e i ( δ 1 ρ π ) 1 e i ( δ 1 ρ π ) d r + I C ϵ + 1 + ϵ ϕ ρ , μ r e i ( π π ) , z r e i ( π π ) 1 e i ( π π ) d r
using the notation (56) this expression can be written in the form
E ρ , μ ( z ) = 1 + ϵ K ρ , μ ( r , δ 1 ρ , z ) d r + I C ϵ + 1 + ϵ K ρ , μ ( r , π , z ) d r ,
where the form of functions K ρ , μ ( r , δ , z ) is defined as (57). We substitute the integral I C ϵ in this expression by the expression (77). As a result, we obtain
E ρ , μ ( z ) = 0 1 ε 1 K ρ , μ ( r , π , z ) d r + π 0 P ρ , μ ( ε 1 , ψ , 0 , z ) d ψ + 1 + ε 1 K ρ , μ ( r , π , z ) d r 0 K ρ , μ ( r , δ 1 ρ , z ) d r , μ < 1 + 1 / ρ .
Thus, we have obtained the statement of the theorem for the third case.
Case 4. Consider now the case 1 / 2 < ρ 1 , δ 1 ρ = δ 2 ρ = π . In this case the segments Γ 1 = { ζ : arg ζ = 2 π , ε | ζ | 1 + ϵ } and Γ 2 = { ζ : arg ζ = 0 , ε | ζ | 1 + ϵ } 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 2 π . Thus, the segment Γ 1 passes through a singular point ζ = 1 ( | ζ | = 1 , arg ζ = 2 π ) , and the segment Γ 2 passes through a singular point ζ = 1 ( | ζ | = 1 , arg ζ = 0 ). 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 C ϵ with the center at the origin of coordinates and radius 1 + ϵ , ϵ > 0 , the segment A B , the arc of the circle C 1 ε 1 with the center in the point ζ = 1 ( | ζ | = 1 , arg ζ = 0 ) and radius 0 < ε 1 < min ( 1 , ϵ ) defined by the expression (69), the segment C D , the arc of the circle C ε with the center at the origin of coordinates and radius 0 < ε < 1 ε 1 , the segment E F , the arc of the circle C 1 ε 1 with the center in the point ζ = 1 ( | ζ | = 1 , arg ζ = 2 π ) and radius 0 < ε 1 < min ( 1 , ϵ ) having the form (61) and the segment G H . The contour is traversed in a positive direction. As a result, the contour Γ is written in the form
Γ = C ϵ = { ζ : 2 π arg ζ 0 , | ζ | = 1 + ϵ } , A B = { ζ : arg ζ = 0 , 1 + ϵ | ζ | 1 + ε 1 } , C 1 ε 1 = ζ : arg ζ = arctan τ sin ψ τ cos ψ + 1 + k π , | ζ | = τ 2 + 2 τ cos ψ + 1 , 0 ψ π , τ = ε 1 , k = 0 , C D = { ζ : arg ζ = 0 , 1 ε 1 | ζ | ε } , C ε = { ζ : 0 arg ζ 2 π , | ζ | = ε } , E F = { ζ : arg ζ = 2 π , ε | ζ | 1 ε 1 } , C 1 ε 1 = ζ : arg ζ = arctan τ sin ψ τ cos ψ + 1 + k π , | ζ | = τ 2 + 2 τ cos ψ + 1 , π ψ 2 π , τ = ε 1 , k = 2 , , G H = { ζ : arg ζ = 2 π , 1 + ε 1 | ζ | 1 + ϵ } .
Next, we consider an auxiliary integral
I = Γ ϕ ρ , μ ( ζ , z ) ζ 1 d ζ ,
where ϕ ρ , μ ( ζ , z ) is defined by (13). By calculating this integral directly we obtain
I = Γ ϕ ρ , μ ( ζ , z ) ζ 1 d ζ = C ϵ ϕ ρ , μ ( ζ , z ) ζ 1 d ζ + A B ϕ ρ , μ ( ζ , z ) ζ 1 d ζ + C 1 ε 1 ϕ ρ , μ ( ζ , z ) ζ 1 d ζ + C D ϕ ρ , μ ( ζ , z ) ζ 1 d ζ + C ε ϕ ρ , μ ( ζ , z ) ζ 1 d ζ + E F ϕ ρ , μ ( ζ , z ) ζ 1 d ζ + C 1 ε 1 ϕ ρ , μ ( ζ , z ) ζ 1 d ζ + G H ϕ ρ , μ ( ζ , z ) ζ 1 d ζ = I C ϵ + I A B + I C 1 ε 1 + I C D + I C ε + I E F + I C 1 ε 1 + I G H .
We immediately note that inside the region bounded by the contour Γ , the function ϕ ρ , μ ( ζ , z ) / ( ζ 1 ) is an analytical function of the complex variable ζ . Consequently, according to the Cauchy integral theorem,
I C ϵ + I A B + I C 1 ε 1 + I C D + I C ε + I E F + I C 1 ε 1 + I G H = 0 .
From this it directly follows that
I C ϵ = I A B I C 1 ε 1 I C D I C ε I E F I C 1 ε 1 I G H .
Next, we calculate each integral on the right-hand side of this expression. However, these integrals were already calculated by us earlier. The integral I A B corresponds to the integral I A B . Therefore, using (73) we obtain
I A B = 1 + ϵ 1 + ε 1 K ρ , μ ( r , π , z ) d r .
The integral I C D corresponds to the integral I C D . Using (74), we obtain
I C D = 1 ε 1 ε K ρ , μ ( r , π , z ) d r .
The integral I C 1 ε 1 corresponds to the integral I C 1 ε 1 . Using (76)
I C 1 ε 1 = π 0 P ρ , μ ( ε 1 , ψ , 0 , z ) d ψ .
The integral I C ε was considered by us when dealing with case 1. Using (42) for the case considered, we obtain
lim ε 0 I C ε = 0 , μ < 1 + 1 / ρ .
The integral I E F corresponds to the integral I C D . Using (59), we obtain
I E F = ε 1 ε 1 K ρ , μ ( r , π , z ) d r .
The integral I G H corresponds to the integral I E F . Therefore, using (60) we obtain
I G H = 1 + ε 1 1 + ϵ K ρ , μ ( r , π , z ) d r .
The integral I C 1 ε 1 corresponds to the integral I C 1 ε 1 . Using (67), we obtain
I C 1 ε 1 = 2 π π P ρ , μ ( ε 1 , ψ , 2 , z ) d ψ .
We get back to the expression (78) and let ε 0 in it. Note that the integrals I A B , I C 1 ε 1 , I C 1 ε 1 , I G H 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
I C ε = 1 + ϵ 1 + ε 1 K ρ , μ ( r , π , z ) d r + π 0 P ρ , μ ( ε 1 , ψ , 0 , z ) d ψ lim ε 0 1 ε 1 ε K ρ , μ ( r , π , z ) d r lim ε 0 ε 1 ε 1 K ρ , μ ( r , π , z ) d r + 2 π π P ρ , μ ( ε 1 , ψ , 2 , z ) d ψ 1 + ε 1 1 + ϵ K ρ , μ ( r , π , z ) d r , μ < 1 + 1 / ρ .
Performing the passage to the limit in this expression ε 0 we obtain
I C ε = 1 + ϵ 1 + ε 1 K ρ , μ ( r , π , z ) d r + π 0 P ρ , μ ( ε 1 , ψ , 0 , z ) d ψ 1 ε 1 0 K ρ , μ ( r , π , z ) d r 0 1 ε 1 K ρ , μ ( r , π , z ) d r + 2 π π P ρ , μ ( ε 1 , ψ , 2 , z ) d ψ 1 + ε 1 1 + ϵ K ρ , μ ( r , π , z ) d r , μ < 1 + 1 / ρ .
Now we consider the Mittag-Leffler (37). Taking into account that in the studied case δ 1 ρ = δ 2 ρ = π we get
E ρ , μ ( z ) = γ ζ ϕ ρ , μ ( ζ , z ) ζ 1 d ζ ,
where the contour γ ζ takes the form
γ ζ = S 1 = { ζ : arg ζ = 2 π , | ζ | 1 + ϵ } , C ε = { ζ : 2 π arg ζ 0 , | ζ | = 1 + ϵ , S 2 = { ζ : arg ζ = 0 , | ζ | 1 + ϵ } ,
and the condition (2) is written in the form π 2 ρ < arg z < π 2 ρ + 2 π .
Directly calculating this integral and representing the complex number ζ in the form ζ = r e i φ we have
E ρ , μ ( z ) = S 1 ϕ ρ , μ ( ζ , z ) ζ 1 d ζ + C ϵ ϕ ρ , μ ( ζ , z ) ζ 1 d ζ + S 2 ϕ ρ , μ ( ζ , z ) ζ 1 d ζ = = 1 + ϵ ϕ ρ , μ r e i φ , z r e i φ 1 e i φ d r φ = π π + I C ϵ + 1 + ϵ ϕ ρ , μ r e i φ , z r e i φ 1 e i φ d r φ = π π = 1 + ϵ ϕ ρ , μ r e i ( π π ) , z r e i ( π π ) 1 e i ( π π ) d r + I C ϵ + 1 + ϵ ϕ ρ , μ r e i ( π π ) , z r e i ( π π ) 1 e i ( π π ) d r .
Now using here (56) we obtain that the expression may be represented in the form
E ρ , μ ( z ) = 1 + ϵ K ρ , μ ( r , π , z ) d r + I C ϵ + 1 + ϵ K ρ , μ ( r , π , z ) d r .
Now here we make use of the representation (79) for the integral I C ϵ . As a result, we obtain
E ρ , μ ( z ) = 1 + ϵ K ρ , μ ( r , π , z ) d r 1 + ϵ 1 + ε 1 K ρ , μ ( r , π , z ) d r + π 0 P ρ , μ ( ε 1 , ψ , 0 , z ) d ψ 1 ε 1 0 K ρ , μ ( r , π , z ) d r 0 1 ε 1 K ρ , μ ( r , π , z ) d r + 2 π π P ρ , μ ( ε 1 , ψ , 2 , z ) d ψ 1 + ε 1 1 + ϵ K ρ , μ ( r , π , z ) d r + 1 + ϵ K ρ , μ ( r , π , z ) d r = 0 1 ε 1 ( K ρ , μ ( r , π , z ) K ρ , μ ( r , π , z ) ) d r + 1 + ε 1 ( K ρ , μ ( r , π , z ) K ρ , μ ( r , π , z ) ) d r + π 0 P ρ , μ ( ε 1 , ψ , 0 , z ) d ψ + 2 π π P ρ , μ ( ε 1 , ψ , 2 , z ) d ψ .
In the expression derived there is the difference of K ρ , μ ( r , π , z ) K ρ , μ ( r , π , z ) . We transform this difference. For this we use the representation (57) for the function K ρ , μ ( r , δ , z ) . As a result, we obtain
K ρ , μ ( r , δ , z ) K ρ , μ ( r , δ , z ) = ρ 2 π i exp ( z r ) ρ e i π ρ cos ( ρ δ ) ( z r e i π ) ρ ( 1 μ ) r 2 + 2 r cos δ + 1 e i η ( r , δ , z ) r + e i δ e i η ( r , δ , z ) r + e i δ .
From the definition for the function η ( r , δ , z ) (10) it follows that it is an odd function relative to the parameter δ
η ( r , δ , z ) = η ( r , δ , z ) ,
In the total, we have
K ρ , μ ( r , δ , z ) K ρ , μ ( r , δ , z ) = ρ 2 π i exp ( z r ) ρ e i π ρ cos ( ρ δ ) ( z r e i π ) ρ ( 1 μ ) r 2 + 2 r cos δ + 1 e i η ( r , δ , z ) r + e i δ e i η ( r , δ , z ) r + e i δ ρ 2 π i exp ( z r ) ρ e i π ρ cos ( ρ δ ) ( z r e i π ) ρ ( 1 μ ) r 2 + 2 r cos δ + 1 r e i η ( r , δ , z ) e i η ( r , δ , z ) + e i ( η ( r , δ , z ) + δ ) e i ( η ( r , δ , z ) + δ ) = ρ 2 π i exp ( z r ) ρ e i π ρ cos ( ρ δ ) ( z r e i π ) ρ ( 1 μ ) r 2 + 2 r cos δ + 1 ( r sin ( η ( r , δ , z ) ) + sin ( η ( r , δ , z ) + δ ) ) = K ρ , μ ( r , δ , z ) .
Here the function K ρ , μ ( r , δ , z ) was obtained by us earlier and has the form (19). As a result, the expression (81) takes the form
E ρ , μ ( z ) = 0 1 ε 1 K ρ , μ ( r , π , z ) d r + 1 + ε 1 K ρ , μ ( r , π , z ) d r + π 0 P ρ , μ ( ε 1 , ψ , 0 , z ) d ψ + 2 π π P ρ , μ ( ε 1 , ψ , 2 , z ) d ψ .
The derived expression proves the theorem completely. ☐
The proved theorem shows that in the representation (3) at the values of parameters ρ , δ 1 ρ , δ 2 ρ satisfying the conditions (32) the integral over the arc of the circle C ϵ of the contour γ ζ can be replaced with integration with respect to the segments Γ 1 and Γ 2 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 S 1 and S 2 and the second summand is the integral along the arc of a circle C ϵ . 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 E ρ , μ ( z ) . As it was shown in Theorem 3, the representation “B” is true only at values μ < 1 + 1 / ρ . If μ 1 + 1 / ρ , 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 1 / 2 < ρ 1 and values δ 1 ρ = π or δ 2 ρ = π the segments Γ 1 or Γ 2 of the auxiliary contour Γ (see Figure 2) pass through the pole ζ = 1 . 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 C 1 ε 1 and C 1 ε 1 . In addition, one has to split the integrals along the half-lines Γ 1 + S 1 and Γ 2 + S 2 into parts to exclude the sections corresponding to the integrals along the arcs of the circle C 1 ε 1 and C 1 ε 1 . All this leads to the appearance of additional terms in the corresponding integral representations. It should also be noted that at parameter values 1 / 2 < ρ 1 , δ 1 ρ = δ 2 ρ = π in the representation (36) the integral 0 1 ε 1 K ρ , μ ( r , π , z ) d r corresponds to the sum of integrals along the segments C D and E F the contour of integration Γ (see Figure 5) the integral 1 + ε 1 K ρ , μ ( r , π , z ) d r corresponds to the sum of integrals along A B + S 2 and G H + S 1 . 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 1 / 2 < ρ < 1 or ρ = 1 and μ 0 or ρ = 1 , μ = 0 and μ is not an integer, the point ζ = 0 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, 0 1 ε 1 K ρ , μ ( r , π , z ) d r 0 and 1 + ε 1 K ρ , μ ( r , π , z ) d r 0 . In case, if ρ = 1 , μ = 0 and μ is an integer, Lemma 1 shows that in this case the point ζ = 0 is a regular point. Consequently, the segments C D and E F , as well as the half-lines A B + S 2 and G H + S 1 will go along one straight line of the complex plane, but in different directions, and arcs of the circles C 1 ε 1 and C 1 ε 1 will close. As a result, we obtain 0 1 ε 1 K 1 , μ ( r , π , z ) d r = 0 and 1 + ε 1 K 1 , μ ( r , π , z ) d r = 0 , and the sum of integrals π 0 P 1 , μ ( ε 1 , ψ , 0 , z ) d ψ + 2 π π P 1 , μ ( ε 1 , ψ , 2 , z ) d ψ corresponds to the integral over the closed contour and will be equal to the residue in the point ζ = 1 . 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 δ 1 ρ and δ 2 ρ coincide, i.e., δ 1 ρ = δ 2 ρ = δ ρ . We will formulate the result obtained in the form of a corollary.
Corollary 3.
For any real ρ > 1 / 2 , any complex μ and z satisfying the conditions μ < 1 + 1 ρ and
π 2 ρ δ ρ + π < arg z < π 2 ρ + δ ρ + π ,
the Mittag-Leffler function can be represented in the form:
1. 
at any real δ ρ satisfying the conditions π 2 ρ < δ ρ π ρ if ρ > 1 and π 2 ρ < δ ρ < π if 1 / 2 < ρ 1
E ρ , μ ( z ) = 0 K ρ , μ ( r , δ ρ , z ) d r ,
where K ρ , μ ( r , δ ρ , z ) has the form (19).
2. 
at δ ρ = π ρ and ρ > 1
E ρ , μ ( z ) = 0 K ρ , μ ( r , z ) d r ,
where K ρ , μ ( r , z ) has the form (21).
Proof. 
(1) According to Theorem 3 the representation (33) is true at values δ 1 ρ and δ 2 ρ satisfying the conditions (32). In this corollary it is assumed that δ 1 ρ = δ 2 ρ = δ ρ . In view of this assumption, the conditions (32) take the form π 2 ρ < δ ρ π ρ if ρ > 1 and π 2 ρ < δ ρ < π if 1 / 2 < ρ 1 . The representation (33) will be written in the form
E ρ , μ ( z ) = 0 K ρ , μ ( r , δ ρ , δ ρ , z ) d r .
From this it is clear that it is necessary to consider how the kernel function K ρ , μ ( r , δ ρ , δ ρ , z ) 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
K ρ , μ ( r , δ ρ , δ ρ , z ) = K ρ , μ ( r , δ ρ , z ) = = ρ π ( z r e i π ) ρ ( 1 μ ) exp ( z r e i π ) ρ cos ( ρ δ ρ ) r sin ( η ( r , δ ρ , z ) ) + sin ( η ( r , δ ρ , z ) + δ ρ ) r 2 + 2 r cos δ ρ + 1 .
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 ρ 1 . However, the value ρ = 1 should be excluded from consideration. In fact, at ρ = 1 the segments Γ 1 and Γ 2 of the auxiliary contour Γ (see Figure 2) pass through the pole ζ = 1 . 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
E ρ , μ ( z ) = 0 K ρ , μ ( r , π / ρ , z ) d r .
Using the definition of the function η ( r , ψ , z ) (see (10)) we obtain η ( r , π / ρ , z ) = ( 1 μ ) π . Now using this result when calculating the kernel function K ρ , μ ( r , π / ρ , z ) and introducing the notation K ρ , μ ( r , π / ρ , z ) K ρ , μ ( r , z ) we obtain
E ρ , μ ( z ) = 0 K ρ , μ ( r , z ) d r ,
where K ρ , μ ( r , z ) has the form (21). ☐

4. Conclusions

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 1 / 2 < ρ 1 there is no need to bypass a pole in the point ζ = 1 . 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 μ < 1 + 1 / ρ . Particular cases also arise at values 1 / 2 < ρ 1 in which one has to bypass a singular point ζ = 1 . 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 ρ , μ , δ 1 ρ , δ 2 ρ . The parameters ρ and μ are referred to the Mittag-Leffler function directly, and parameters δ 1 ρ and δ 2 ρ 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 E ρ , μ ( z ) have a more general form but they are awkward enough. The representations “A” and “B” take a simpler form if δ 1 ρ = δ 2 ρ = δ ρ . 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 δ 1 ρ = δ 2 ρ = π / ρ . 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 δ 1 ρ and δ 2 ρ one can give the geometric interpretation of three introduced parametrizations of the representations “A” and “B”. In fact, the parameters δ 1 ρ and δ 2 ρ describe an inclination angle of half-lines S 1 and S 2 in the contour γ ζ (see Figure 1). Thus, in case if the half-lines S 1 and S 2 independently lie in the range of angles π 2 ρ < δ 1 ρ min π , π ρ , π 2 ρ < δ 2 ρ min π , π ρ , 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 δ 1 ρ = δ 2 ρ = π / ρ .
In conclusion it should be pointed out that the representations for the function E ρ , μ ( z ) , formulated in Theorems 2 and 3 are valid for the vaues arg z satisfying the condition π 2 ρ δ 2 ρ + π < arg z < π 2 ρ + δ 1 ρ + π . 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 E ρ , μ ( z ) . Nevertheless, it is possible to get rid of this constraint and expand the range of admissible values arg z 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

  1. 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]
  2. 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]
  3. 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]
  4. Wiman, A. Über den Fundamentalsatz in der Teorie der Funktionen Ea(x). Acta Math. 1905, 29, 191–201. [Google Scholar] [CrossRef]
  5. Wiman, A. Über die Nullstellen der Funktionen Ea(x). Acta Math. 1905, 29, 217–234. [Google Scholar] [CrossRef]
  6. 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]
  7. Humbert, P. Quelques résultats rélatifs à la fonction de Mittag-Leffler. C. R. Acad. Sci. 1953, 236, 1467–1468. [Google Scholar]
  8. Agarwal, R.P. A propos d’une note de M. Pierre Humbert. C. R. Acad. Sci. 1953, 236, 2031–2032. [Google Scholar]
  9. 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]
  10. 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]
  11. 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]
  12. Popov, A.Y.; Sedletskii, A.M. Distribution of roots of Mittag-Leffler functions. J. Math. Sci. 2013, 190, 209–409. [Google Scholar] [CrossRef]
  13. Rogosin, S. The role of the Mittag-Leffler function in fractional modeling. Mathematics 2015, 3, 368–381. [Google Scholar] [CrossRef]
  14. 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]
  15. Bateman, H. Higher Transcendental Functions; McGraw-Hill Book Company: New York, NY, USA, 1955; Volume 3, p. 292. [Google Scholar]
  16. Podlubny, I. Fractional Differential Equations; Volume 198, Mathematics in Science and Engineering; Academic Press: Amsterdam, The Netherlands, 1999; p. 340. [Google Scholar]
  17. 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]
  18. Seybold, H.J.; Hilfer, R. Numerical results for the generalized Mittag-Leffler function. Fract. Calc. Appl. Anal. 2005, 8, 129–139. [Google Scholar]
  19. 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]
  20. Seybold, H.; Hilfer, R. Numerical Algorithm for Calculating the Generalized Mittag-Leffler Function. SIAM J. Numer. Anal. 2009, 47, 69–88. [Google Scholar] [CrossRef]
  21. 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]
  22. Saenko, V.V. An integral representation of the Mittag-Leffler function. arXiv 2020, arXiv:2001.09606. [Google Scholar]
  23. 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]
  24. Saenko, V.V. Integral Representation of the Fractional Stable Density. J. Math. Sci. 2020, 248, 51–66. [Google Scholar] [CrossRef]
  25. Uchaikin, V.V.; Cahoy, D.O.; Sibatov, R.T. Fractional processes: From Poisson to branching one. Int. J. Bifurc. Chaos 2008, 18, 2717–2725. [Google Scholar] [CrossRef] [Green Version]
  26. Uchaikin, V.; Sibatov, R.; Uchaikin, D. Memory regeneration phenomenon in dielectrics: The fractional derivative approach. Phys. Scr. 2009, T136, 014002. [Google Scholar] [CrossRef]
  27. Saenko, V.V. Singular points of the integral representation of the Mittag-Leffler function. arXiv 2020, 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.
Mathematics 08 01101 g001
Figure 2. The auxiliary contour of integration Γ .
Figure 2. The auxiliary contour of integration Γ .
Mathematics 08 01101 g002
Figure 3. The auxiliary contour of integration Γ .
Figure 3. The auxiliary contour of integration Γ .
Mathematics 08 01101 g003
Figure 4. The auxiliary contour of integration Γ .
Figure 4. The auxiliary contour of integration Γ .
Mathematics 08 01101 g004
Figure 5. The auxiliary contour of integration Γ .
Figure 5. The auxiliary contour of integration Γ .
Mathematics 08 01101 g005

Share and Cite

MDPI and ACS Style

Saenko, V. Two Forms of the Integral Representations of the Mittag-Leffler Function. Mathematics 2020, 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

Back to TopTop