Next Article in Journal
Design of Fractional-Order Lead Compensator for a Car Suspension System Based on Curve-Fitting Approximation
Next Article in Special Issue
Solving a System of Fractional-Order Volterra-Fredholm Integro-Differential Equations with Weakly Singular Kernels via the Second Chebyshev Wavelets Method
Previous Article in Journal
Elastic Rough Surface Contact and the Root Mean Square Slope of Measured Surfaces over Multiple Scales
Previous Article in Special Issue
A Fractional SAIDR Model in the Frame of Atangana–Baleanu Derivative
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On a Five-Parameter Mittag-Leffler Function and the Corresponding Bivariate Fractional Operators

by
Mehmet Ali Özarslan
* and
Arran Fernandez
Department of Mathematics, Faculty of Arts and Sciences, Eastern Mediterranean University, via Mersin-10, Famagusta 99628, Northern Cyprus, Turkey
*
Author to whom correspondence should be addressed.
Fractal Fract. 2021, 5(2), 45; https://doi.org/10.3390/fractalfract5020045
Submission received: 29 March 2021 / Revised: 7 May 2021 / Accepted: 12 May 2021 / Published: 14 May 2021
(This article belongs to the Special Issue Fractional Calculus and Special Functions with Applications)

Abstract

:
Several extensions of the classical Mittag-Leffler function, including multi-parameter and multivariate versions, have been used to define fractional integral and derivative operators. In this paper, we consider a function of one variable with five parameters, a special case of the Fox–Wright function. It turns out that the most natural way to define a fractional integral based on this function requires considering it as a function of two variables. This gives rise to a model of bivariate fractional calculus, which is useful in understanding fractional differential equations involving mixed partial derivatives.

1. Introduction

The original Mittag-Leffler function E α ( z ) , applied to one variable according to one parameter, was first defined and studied by Gösta Mittag-Leffler in the 1900s. In the hundred years since then, many variants and extensions of this function have been defined, including functions of more than one variable and functions with arbitrarily many parameters [1,2].
One of the major motivations for studying Mittag-Leffler functions is their relationship with fractional calculus [2,3,4,5], a field of mathematics which has become very popular due to its many applications in various areas of science [6,7,8]. Mittag-Leffler functions emerge naturally as the solution to some elementary fractional differential equations, and their eigenfunction properties have led them sometimes to be called “fractional exponential functions” [9].
Among the many extensions of the Mittag-Leffler function to more variables and parameters, we mention just a few which are of particular importance or interest in motivating the present work.
  • A Mittag-Leffler function of one variable with three parameters was defined by Prabhakar [10] to solve a certain singular integral equation. Its use as an integral kernel gave rise to a model of fractional calculus which has a semigroup property and which is already broad enough to include many other named fractional-calculus operators [11,12,13], although it is itself a special case of the general Fox–Wright function.
  • A Mittag-Leffler function of n variables with n + 1 parameters was used by Luchko et al. [14,15] to solve multi-term linear fractional differential equations involving multiple independent fractional orders. Note that this is independent from the Mittag-Leffler function of n variables with 2 n + 1 parameters which was defined by Saxena et al. [16] and which gives rise to a model of fractional calculus with a semigroup property [17]: neither of these general functions is a special case of the other.
  • A Mittag-Leffler function of two variables with four parameters and a Mittag-Leffler function of three variables with five parameters were recently defined [18,19] and used to solve multi-order systems of fractional differential equations [20,21]. These are closely related to the general function considered by Luchko et al., but they are not special cases of it. Used as kernels, they also give rise to new models of fractional calculus with a semigroup property [18,21]. For other types of bivariate Mittag-Leffler function that have been defined in the literature, we refer to the papers [22,23,24].
All of these Mittag-Leffler type functions have been connected in some way to fractional calculus and fractional differential equations. Indeed, the strong connection between special functions and fractional calculus has been known for decades and continues to be written about today [25,26,27], while multi-parameter and multi-variable generalisations are being studied both for special functions and for the operators of fractional calculus [17].
Fractional partial differential equations have been an interesting and challenging topic of study [3], with various methods able to be extended (with modifications) from classical partial differential equations in order to solve them, such as the unified transform method [28], distribution theory [29,30], and weak solutions [31,32]. Much attention has been paid to partial differential equations involving the fractional Laplacian operator [33,34], but less attention has gone to differential equations involving mixed partial fractional derivatives. They are mentioned in the classical textbooks [3] (§6.1.1) and [6] (§24.2), and in a few papers (e.g., [35,36,37]), but in general they have attracted little notice in the research literature on fractional partial differential equations.
In the work below, we study a specific type of Mittag-Leffler function, initially a function of one variable with five parameters, and define a double fractional integral operator by converting this function to a bivariate version. In this way, it is possible to involve double integrals and derivatives, with respect to two variables, while preserving the relatively simple structure of a single power series defining the Mittag-Leffler function. These operators are therefore useful in the study of fractional partial differential equations, including those involving mixed partial fractional derivatives.
As an initial motivation, let us consider the following bivariate Abel equation of the second kind:
u ( t , s ) + λ Γ ( α 1 ) Γ ( α 2 ) 0 t 0 s u ( ξ , η ) t τ 1 α 1 s η 1 α 2 d τ d η = f ( t , s ) .
One can reformulate this equation in terms of the standard Riemann–Liouville fractional integral operators as:
1 + λ I t α 1 I s α 2 u ( t , s ) = f ( t , s ) .
Formally, by means of symbolic operational calculus without regard for rigour, this can be solved as follows:
u ( t , s ) = 1 + λ I t α 1 I s α 2 1 f ( t , s ) = n = 0 λ n I t α 1 n I s α 2 n f ( t , s ) .
This last formal result can be written in terms of the following Mittag-Leffler type function:
E α 1 , α 2 ( x ) = n = 0 x n Γ ( α 1 n + 1 ) Γ ( α 2 n + 1 ) ,
namely, in the following way, making use of the double (two-dimensional) Laplace convolution operator denoted here as “∗”:
u ( t , s ) = E α 1 , α 2 λ t α 1 s α 2 f ( t , s ) = 0 t 0 s n = 0 λ n t τ α 1 n s η α 2 n Γ ( α 1 n + 1 ) Γ ( α 2 n + 1 ) f ( τ , η ) d τ d η = n = 0 λ n I t α 1 n I s α 2 n f ( t , s ) .
Therefore, it makes sense to define an integral operator E α 1 , α 2 acting on functions of two variables as follows:
E α 1 , α 2 f ( t , s ) = 0 t 0 s n = 0 λ n t τ α 1 n s η α 2 n Γ ( α 1 n + 1 ) Γ ( α 2 n + 1 ) f ( τ , η ) d τ d η .
This operator emerges naturally from consideration of bivariate Abel equations, but it has the drawback of lacking a semigroup property, or index law, in any of the parameters: if we take a composition of this operator with itself, we do not find another operator in the same form for different values of the parameters.
It should be noted at this point that some useful and important operators, which have been criticised for lacking a semigroup property, can be embedded into a larger class of fractional-calculus operators which has an index law and which therefore contains both the original operators and their compositions. For example, Prabhakar fractional calculus forms a class of operators which contains various useful operators that lack semigroup properties in themselves, their compositions being different elements of the Prabhakar class [11]. We can apply the same way of thinking here, extending the basic two-parameter Mittag-Leffler function (2) and the associated integral operator to more general versions where a semigroup property can be found.
In this case, it is sufficient to add three extra parameters in the definition in order to obtain a generalised operator which has a semigroup property. We consider, throughout this paper, the following five-parameter Mittag-Leffler function of one variable:
E α 1 , α 2 ; β 1 , β 2 γ ( z ) = k = 0 ( γ ) k Γ α 1 k + β 1 Γ α 2 k + β 2 · z k k ! ,
where α 1 , α 2 , β 1 , β 2 , γ are complex constants (with some constraints to be determined later) and ( γ ) k is the Pochhammer symbol defined by
( γ ) 0 = 1 , ( γ ) k = γ ( γ + 1 ) ( γ + k 1 ) , k = 1 , 2 , .
This paper is devoted to a detailed study of the five-parameter Mittag-Leffler function (4), the associated bivariate fractional integral operators, related concepts such as the corresponding fractional derivative operators, and special cases of particular interest.
Specifically, the structure of the paper is as follows. We discuss the five-parameter Mittag-Leffler function and its properties in Section 2, then pass to bivariate fractional calculus in Section 3, firstly integral operators and then derivative operators. In Section 4, we investigate the case where the operators are non-singular and expressible as finite sums of Riemann–Liouville integrals. Section 5 is devoted to discussion of applications and potential future work.

2. The Five-Parameter Mittag-Leffler Function

Definition 1.
Let α 1 , α 2 , β 1 , β 2 , γ C be five parameters satisfying Re ( α 1 + α 2 ) > 0 . The five-parameter Mittag-Leffler function applied to a single variable z is defined by the following power series:
E α 1 , α 2 ; β 1 , β 2 γ ( z ) = k = 0 ( γ ) k Γ α 1 k + β 1 Γ α 2 k + β 2 · z k k ! , z C .
This can be seen as a special case of the Fox–Wright function, which is defined [3] (§1.11) by
Ψ p q ( z ) = Ψ p q ( ( a i , α i ) 1 , p ; ( b i , β i ) 1 , q ; z ) = k = 0 i = 1 p Γ ( a i + α i k ) j = 1 q Γ ( b j + β j k ) · z k k ! , z C ,
where, in [3], the parameters are given as a i , b j C , α i , β j R satisfying α i β j for i = 1 , , p ; j = 1 , , q , and it is proved that this power series is absolutely convergent for all z C in the case Δ : = j = 1 q β j i = 1 p α i > 1 .
The five-parameter Mittag-Leffler function can be written in terms of the Fox–Wright function as follows:
E α 1 , α 2 ; β 1 , β 2 γ ( z ) = 1 Γ ( γ ) Ψ 1 2 ( ( γ , 1 ) ; ( β 1 , α 1 ) , ( β 2 , α 2 ) ; z ) .
As an immediate consequence of the convergence result shown in [3] (Theorem 1.5) and [38], we conclude that, if α 1 , α 2 R and α 1 + α 2 > 0 , then the power series defining E α 1 , α 2 ; β 1 , β 2 γ ( z ) is locally uniformly convergent and therefore it is an entire function.
In fact, there is no need to assume that any of the parameters are real. This assumption is imposed in [3,38] to simplify the calculations, but the same convergence result can be proved for complex values of the parameters by using Stirling’s formula and the ratio test, similarly to the second author’s work in [39]. The condition then to be imposed on the parameters is Re ( α 1 + α 2 ) > 0 , as stated in Definition 1.
It is worth noting that a four-parameter special case of the general Fox–Wright function was recently given particular consideration by Luchko [40], namely the following function:
W ( α 1 , β 1 ) , ( α 2 , β 2 ) ( z ) = k = 0 z k Γ α 1 k + β 1 Γ α 2 k + β 2 ,
which is the case γ = 1 of our five-parameter Mittag-Leffler function. It is also a special case of the vector-index Mittag-Leffler function studied by Al-Bassam and Luchko [41], which is defined as follows:
E ( α 1 , β 1 ) , , ( α n , β n ) ( z ) = k = 0 z k Γ ( α 1 k + β 1 ) Γ ( α n k + β n ) ,
for any n N . The case n = 1 gives the usual two-parameter Mittag-Leffler function, while the case n = 2 gives the four-parameter Wright function (5) considered by Luchko. The function (6) was further extended [42] to include a numerator ( γ ) κ n in terms of two further parameters.
Our function is also a special case of the general Fox–Wright function, but it is different from the special cases (5) and (6) considered previously, because of the Pochhammer symbol appearing in the numerator. This extra parameter and Pochhammer symbol is important because, as we show below, it gives rise to a semigroup property for the resulting fractional integral operators – a property which is lacking for functions such as (5) and (6) which were considered before, and even for functions such as that in [42] which contain a generalised Pochhammer symbol. The semigroup property will arise from the idea of turning this univariate function into a bivariate integral operator, a notion which we justify below from consideration of the gamma functions on the denominator.
Using the Mellin–Barnes integral representation of the Fox–Wright function p Ψ q ( z ) , proved in [3], we have for α 1 > 0 , α 2 > 0 , and γ , β 1 , β 2 C , z 0 that
E α 1 , α 2 ; β 1 , β 2 γ ( z ) = 1 2 π i Γ ( γ ) L Γ ( s ) Γ ( γ s ) Γ ( β 1 α 1 s ) Γ ( β 2 α 2 s ) ( z ) s d s ,
where the integration path L is a Bromwich contour starting at a point C i and terminating at a point C + i and separating the poles m ( m = 0 , 1 , 2 , ) of the gamma function Γ ( s ) to the left and the poles γ + l ( l = 0 , 1 , 2 , ) to the right with the assumption m γ + l for l , m = 0 , 1 , 2 , .
Another complex integral representation is given by the following Theorem.
Theorem 1.
Let α 1 , α 2 , β 1 , β 2 C with Re ( α 1 ) > 0 , Re ( α 2 ) > 0 , and let ε > 0 , π 2 < ϕ π . Let H ( ε ; ϕ ) ( ε > 0 , 0 < ϕ π ) be the contour which is the union of the following three parts oriented according to non-decreasing arg τ :
1. 
the ray arg τ = ϕ , | τ | ε ;
2. 
the arc ϕ arg τ ϕ , | τ | = ε ; and
3. 
the ray arg τ = ϕ , | τ | ε ,
Then, the five-parameter Mittag-Leffler function possesses the following complex integral representations:
E α 1 , α 2 ; β 1 , β 2 γ ( z ) = 1 2 π i H ( ε ; ϕ ) τ β 2 e τ E α 1 , β 1 γ τ α 2 z d τ
= 1 4 π 2 H ( ε ; ϕ ) H ( ε ; ϕ ) ξ β 1 τ β 2 e ξ + τ 1 z ξ α 1 τ α 2 γ d ξ d τ .
Proof. 
Using the known representation [43]
1 Γ ( z ) = 1 2 π i H ( ε ; ϕ ) τ z e τ d τ , z C ,
with z replaced by α 2 k + β 2 , and using the fact that the series represents an entire function, which allows the interchange of the series and the integral, gives
E α 1 , α 2 ; β 1 , β 2 γ ( z ) = k = 0 ( γ ) k Γ α 1 k + β 1 1 2 π i H ( ε ; ϕ ) τ α 2 k β 2 e τ d τ z k k ! = 1 2 π i H ( ε ; ϕ ) τ β 2 e τ k = 0 ( γ ) k Γ α 1 k + β 1 τ α 2 z k k ! d τ ,
which directly yields (7), the first of the desired formulae, in terms of the three-parameter Mittag-Leffler function of Prabhakar.
On the other hand, using (9) for both of the terms 1 Γ α 1 k + β 1 and 1 Γ α 2 k + β 2 , we have
E α 1 , α 2 ; β 1 , β 2 γ ( z ) = 1 4 π 2 k = 0 ( γ ) k H ( ε ; ϕ ) ξ α 1 k β 1 e ξ d ξ · H ( ε ; ϕ ) τ α 2 k β 2 e τ d τ z k k ! = 1 4 π 2 H ( ε ; ϕ ) H ( ε ; ϕ ) ξ β 1 τ β 2 e ξ + τ k = 0 ( γ ) k z ξ α 1 τ α 2 k k ! d ξ d τ ,
which directly yields (8), the second of the desired formulae. This derivation assumes that z ξ α 1 τ α 2 < 1 , which is uniformly true on the given contours provided that Re ( α 1 ) > 0 , Re ( α 2 ) > 0 , and | z | is sufficiently small. The bound on | z | can be removed, by analytic continuation in z, to yield the same result valid for all z. □
Thus far, we have considered the function (4) as a five-parameter function of a single variable z. As a function of z, this function has several properties, such as complex integral representations, which we prove above. However, when we start to construct connections with the topic of fractional calculus, it is more natural to consider a related bivariate function instead.
To see why, consider that fractional differintegral relations for special functions, such as many found in [44,45], usually make use of gamma functions appearing in power series, in the following way:
I a x λ ( x a ) ν Γ ( ν + 1 ) = ( x a ) ν + λ Γ ( ν + λ + 1 ) , D a x λ ( x a ) ν Γ ( ν + 1 ) = ( x a ) ν λ Γ ( ν λ + 1 ) ,
with appropriate choices of ν in order to use such identities for every term of a power series expansion. For many univariate functions defined by power series, the above relations give rise to interesting identities between special functions. In our case, however, the power series (4) has two different gamma functions in the denominator, each of them involving k times a different parameter, α 1 or α 2 . Therefore, to take full advantage of the function’s symmetry in α 1 and α 2 , we should use two different variables: one raised to the power of α 1 , the other raised to the power of α 2 , and then both of them further raised to the power of k.
Motivated by this discussion, we now begin to study, instead of the univariate function E α 1 , α 2 ; β 1 , β 2 γ ( z ) , the closely related bivariate function E α 1 , α 2 ; β 1 , β 2 γ ( x α 1 y α 2 ) . A similar idea, substituting products of fractional powers of new variables instead of old variables, was used in a 2017 paper of the first author [23], but in that case it was used to convert a bivariate function to another bivariate function of different variables. Here, we use it to convert a univariate function to a bivariate function. This seems at first sight as an unnecessary complication, but we see that it makes many things more natural and smooth in the studies related to the five-parameter Mittag-Leffler function.
Definition 2
([3]). The double fractional integrals of a bivariate function f ( x , t ) are defined in the natural way by combining fractional integrals with respect to x and t, namely as follows for λ , μ C with Re ( λ ) > 0 , Re ( μ ) > 0 and for x > a , t > b .
I b t λ I a x μ f ( x , t ) = 1 Γ ( μ ) Γ ( λ ) b t a x ( t τ ) λ 1 ( x ξ ) μ 1 f ( ξ , τ ) d ξ d τ .
Similarly, the partial fractional derivatives of a bivariate function f ( x , t ) are defined as follows, for λ , μ C with Re ( λ ) > 0 , Re ( μ ) > 0 and for n = Re ( μ ) + 1 , m = Re ( λ ) + 1 and x > a , t > b :
D b t β D a x α f ( x , t ) = m + n t m x n I b t n β I a x m α f ( x , t ) .
Lemma 1.
Let α 1 , α 2 , β 1 , β 2 , γ , ω C with Re ( α 1 ) , Re ( α 2 ) , Re ( β 1 ) , Re ( β 2 ) > 0 . Then, for any λ , μ C with Re ( λ ) > 0 , Re ( μ ) > 0 , we have
I a x λ I c y μ ( x a ) β 1 1 ( y c ) β 2 1 E α 1 , α 2 ; β 1 , β 2 γ ω ( x a ) α 1 ( y c ) α 2      = ( x a ) λ + β 1 1 ( y c ) μ + β 2 1 E α 1 , α 2 ; λ + β 1 , μ + β 2 γ ω ( x a ) α 1 ( y c ) α 2 .
Proof. 
We know from above that the infinite series defining the five-parameter Mittag-Leffler function is locally uniformly convergent, so we have the right to interchange the order of this series with fractional integral operators. Since all the exponents of ( x a ) and ( y c ) are greater than 1 , we have:
I a x λ I c y μ ( x a ) β 1 1 ( y c ) β 2 1 E α 1 , α 2 ; β 1 , β 2 γ ω ( x a ) α 1 ( y c ) α 2 = k = 0 ( γ ) k Γ α 1 k + β 1 Γ α 2 k + β 2 · ω k k ! I a x λ ( x a ) β 1 + α 1 k 1 I c y μ ( y c ) β 2 + α 2 k 1 = k = 0 ( γ ) k Γ α 1 k + β 1 Γ α 2 k + β 2 · ω k k ! × Γ ( β 1 + α 1 k ) Γ ( λ + β 1 + α 1 k ) ( x a ) λ + β 1 + α 1 k 1 Γ ( β 2 + α 2 k ) Γ ( μ + β 2 + α 2 k ) ( y c ) μ + β 2 + α 2 k 1 = ( x a ) λ + β 1 1 ( y c ) μ + β 2 1 E α 1 , α 2 ; λ + β 1 , μ + β 2 γ ω ( x a ) α 1 ( y c ) α 2 .
Note how the two gamma functions in the denominator of the series mesh together precisely with the gamma-function quotients arising from the two fractional integrals, in order to achieve the desired result. □
Lemma 2.
Let α 1 , α 2 , β 1 , β 2 , γ , ω C with Re ( α 1 ) , Re ( α 2 ) , Re ( β 1 ) , Re ( β 2 ) > 0 . Then, for any λ , μ C with Re ( λ ) 0 , Re ( μ ) 0 , we have
D a x λ D c y μ ( x a ) β 1 1 ( y c ) β 2 1 E α 1 , α 2 ; β 1 , β 2 γ ω ( x a ) α 1 ( y c ) α 2        = ( x a ) β 1 λ 1 ( y c ) β 2 μ 1 E α 1 , α 2 ; β 1 λ , β 2 μ γ ω ( x a ) α 1 ( y c ) α 2 .
Proof. 
This can be deduced from Lemma 1 by analytic continuation in the variables λ and μ , using the analyticity properties of fractional differintegrals. Alternatively, it can be proved from the series expansion of the function, following similar lines as the proof of Lemma 1. □
Lemma 3.
Setting all parameters to 1 in the bivariate form of the five-parameter Mittag-Leffler function, we can recover a case of the modified Bessel function with parameter 0:
E 1 , 1 ; 1 , 1 1 ( x y ) = I 0 2 x y = k = 0 ( x y ) k ( k ! ) 2 .
Proof. 
This follows immediately from the series definition of the function. □

3. Bivariate Operators with Five-Parameter Mittag-Leffler Kernels

In this section, we move on from functions to operators. Having established the five-parameter Mittag-Leffler function and some of its properties, we now wish to define a fractional integral operator using this function as a kernel, following in the footsteps of other papers [10,18,46] which defined new models of fractional calculus by using various types of Mittag-Leffler functions as kernels.
The five-parameter Mittag-Leffler function (4) is defined by a single power series in terms of a single variable z. When we use it to define a fractional integral operator, however, we must transform it to a bivariate function, again using a single summation, but this time in terms of two independent variables x , y . Correspondingly, we use a double integral and create an operator to be applied to bivariate functions. This is the only natural way to define a model of fractional calculus using the five-parameter Mittag-Leffler function, because the two gamma functions in its denominator will give rise to a double fractional integral in each summand when we wish to write a series representation for the new operator: we need two separate powers in the integrand, one corresponding to each of the gamma functions from the denominator.
Definition 3.
The fractional integral operator based on the five-parameter Mittag-Leffler function as a kernel is given by
I a , b α 1 , α 2 ; β 1 , β 2 γ ; λ ( f ) ( x , y ) : = a x b y f ( t , s ) ( x t ) β 1 1 ( y s ) β 2 1 E α 1 , α 2 ; β 1 , β 2 γ ( λ ( x t ) α 1 ( y s ) α 2 ) d s d t ,
where a , b R are fixed with x > a , y > b and α 1 , α 2 , β 1 , β 2 , γ , λ C are parameters such that Re ( α 1 ) , Re ( α 2 ) , Re ( β 1 ) , Re ( β 2 ) > 0 . Note that these restrictions on the parameters are necessary in order to have a properly convergent singular integral for all reasonably well-behaved functions f (more details later on the function space for f).
It is clear that the two-parameter bivariate integral operator defined in (3) above is a special case of this new fractional integral operator: E α 1 , α 2 f ( x , y ) = I 0 , 0 α 1 , α 2 ; 1 , 1 1 ; λ ( f ) ( x , y ) .
Theorem 2.
Let α 1 , α 2 , β 1 , β 2 , γ , λ C with Re ( α 1 ) , Re ( α 2 ) , Re ( β 1 ) , Re ( β 2 ) > 0 . Then, a , b I α 1 , α 2 ; β 1 , β 2 γ ; λ is a bounded operator from the space L 1 ( [ a , c ] × [ b , d ] ) to itself.
Proof. 
Using Fubini’s theorem, we have
a c b d a , b I α 1 , α 2 ; β 1 , β 2 γ ; λ f ( x , y ) d y d x a c a x b d b y f ( t , s ) ( x t ) β 1 1 ( y s ) β 2 1 E α 1 , α 2 ; β 1 , β 2 γ ( λ ( x t ) α 1 ( y s ) α 2 ) d s d y d t d x = a c b d f ( t , s ) 0 c t 0 d s u β 1 1 w β 2 1 E α 1 , α 2 ; β 1 , β 2 γ ( λ u α 1 w α 2 ) d w d u d s d t .
Since the five-parameter Mittag-Leffler function is an entire function when Re ( α 1 + α 2 ) > 0 , we have a bound of the form E α 1 , α 2 ; β 1 , β 2 γ ( λ u α 1 w β 2 1 ) C on the finite domain [ a , c ] × [ b , d ] . Therefore,
a , b I α 1 , α 2 ; β 1 , β 2 γ ; λ f 1 : = a c b d a , b I α 1 , α 2 ; β 1 , β 2 γ ; λ f ( x , y ) d y d x C ( c a ) β 1 1 d b β 2 1 f 1 .
Since a , b , c , d , β 1 , β 2 are fixed, the result is proved. □
Theorem 3.
Let α 1 , α 2 , β 1 , β 2 , γ , λ C with Re ( α 1 ) , Re ( α 2 ) , Re ( β 1 ) , Re ( β 2 ) > 0 , and let f L 1 ( [ a , c ] × [ b , d ] ) . Then, the operator I a , b α 1 , α 2 ; β 1 , β 2 γ ; λ can be represented as an infinite series of double fractional integrals of Riemann–Liouville type:
I a , b α 1 , α 2 ; β 1 , β 2 γ ; λ f ( x , y ) = k = 0 ( γ ) k λ k k ! I b y β 2 + α 2 k I a x β 1 + α 1 k ( f ) ( x , y ) ,
where the right-hand side is locally uniformly convergent for x , y [ a , c ] × [ b , d ] .
Proof. 
Since E α 1 , α 2 ; β 1 , β 2 γ ( u ) is an entire function defined by a locally uniformly convergent power series, and f L 1 ( [ a , c ] × [ b , d ] ) , we have the right to interchange the order of summation and integration, to yield
I a , b α 1 , α 2 ; β 1 , β 2 γ ; λ ( f ) ( x , y ) = a x b y f ( t , s ) ( x t ) β 1 1 ( y s ) β 2 1 E α 1 , α 2 ; β 1 , β 2 γ ( λ ( x t ) α 1 ( y s ) α 2 ) d t d s = a x b y ( x t ) β 1 1 ( y s ) β 2 1 k = 0 ( γ ) k λ k ( x t ) α 1 k ( y s ) α 2 k Γ α 1 k + β 1 Γ α 2 k + β 2 k ! f ( t , s ) d t d s = k = 0 ( γ ) k λ k k ! · 1 Γ α 1 k + β 1 Γ α 2 k + β 2 a x b y ( x t ) β 1 + α 1 k 1 ( y s ) β 2 + α 2 k 1 f ( t , s ) d t d s = k = 0 ( γ ) k λ k k ! I b y β 2 + α 2 k I a x β 1 + α 1 k ( f ) ( x , y ) .
Whence the result. □
Before the next results, we need to recall the bivariate Laplace transform, or double Laplace transform, which is applied to bivariate functions f ( x , y ) in the following way:
L 2 f ( p , q ) = 0 0 e ( p x + q y ) f ( x , y ) d x d y = L x L y f ( p , q ) , Re ( p ) > 0 , Re ( q ) > 0 ,
provided that this integral is convergent (for example, if f ( x , y ) is exponentially bounded in both variables).
Theorem 4.
Let α 1 , α 2 , β 1 , β 2 , γ , λ C with Re ( α 1 ) , Re ( α 2 ) , Re ( β 1 ) , Re ( β 2 ) > 0 . If f is a bivariate function of exponential order and integrable over [ 0 , ) × [ 0 , ) , then we have
L 2 I a , b α 1 , α 2 ; β 1 , β 2 γ ; λ ( f ) ( p , q ) = 1 p β 1 q β 2 1 λ p α 1 q α 2 γ L 2 [ f ] ( p , q )
for p , q C such that these Laplace transforms exist.
Proof. 
This follows from the series representation of the operator, using the fact that the series is locally uniformly convergent to interchange the Laplace integration and summation:
L 2 I a , b α 1 , α 2 ; β 1 , β 2 γ ; λ ( f ) ( p , q ) = k = 0 ( γ ) k λ k k ! L 2 I b y β 2 + α 2 k I a x β 1 + α 1 k ( f ) ( p , q ) = L 2 [ f ] ( p , q ) k = 0 ( γ ) k λ k k ! q β 2 α 2 k p β 1 α 1 k = 1 p β 1 q β 2 L 2 [ f ] ( p , q ) k = 0 ( γ ) k λ p α 1 q α 2 k k ! = 1 p β 1 q β 2 1 λ p α 1 q α 2 γ L 2 [ f ] ( p , q ) ,
where we have assumed λ p α 1 q α 2 < 1 for convergence of the binomial series, although this condition can be removed by analytic continuation in the variables p and q. □
As an application of the above result, we can use Laplace transforms to quickly learn the result of applying the new fractional operator to the five-parameter Mittag-Leffler function itself, as follows.
Example 1.
Let α 1 , α 2 , β 1 , β 2 , ε 1 , ε 2 , γ , λ , σ C with Re ( α i ) , Re ( β i ) , Re ( ε i ) > 0 for i = 1 , 2 . Then, we have
L 2 I a , b α 1 , α 2 ; β 1 , β 2 γ ; λ x ε 1 1 y ε 2 1 E α 1 , α 2 ; ε 1 , ε 2 σ ( λ x α 1 y α 2 ) ( p , q ) = 1 p β 1 q β 2 1 λ p α 1 q α 2 γ L 2 x ε 1 1 y ε 2 1 E α 1 , α 2 ; ε 1 , ε 2 σ ( λ x α 1 y α 2 ) ( p , q ) = 1 p β 1 q β 2 1 λ p α 1 q α 2 γ 1 p ε 1 q ε 2 1 λ p α 1 q α 2 σ = 1 p β 1 + ε 1 q β 2 + ε 2 1 λ p α 1 q α 2 γ σ .
Taking inverse Laplace transforms on both sides of this equation, we get the following formula showing how the bivariate fractional integral operator with five-parameter Mittag-Leffler function kernel can be applied to this particular type of function:
I a , b α 1 , α 2 ; β 1 , β 2 γ ; λ x ε 1 1 y ε 2 1 E α 1 , α 2 ; ε 1 , ε 2 σ ( λ x α 1 y α 2 ) = x β 1 + ε 1 1 y β 2 + ε 2 1 E α 1 , α 2 ; β 1 + ε 1 , β 2 + ε 2 γ + σ ( λ x α 1 y α 2 ) .
A very important property of the bivariate fractional calculus defined in this paper is that it has a semigroup property in the variables β 1 , β 2 , γ , as expressed by the following theorem. There are several different ways to prove this result, as in [18], and we mention here two of them.
Theorem 5.
Let α 1 , α 2 , β 1 , β 2 , ε 1 , ε 2 , γ , σ , λ C with Re ( α i ) , Re ( β i ) , Re ( ε i ) > 0 for i = 1 , 2 . Then, for any f L 1 ( [ a , c ] × [ b , d ] ) , we have
I a , b α 1 , α 2 ; β 1 , β 2 γ ; λ I a , b α 1 , α 2 ; ε 1 , ε 2 σ ; λ ( f ) = I a , b α 1 , α 2 ; β 1 + ε 1 , β 2 + ε 2 γ + σ ; λ ( f ) .
Proof. 
In the case that a = b = 0 and f is a function whose Laplace transform exists, we know from Theorem 4 that applying the operator I a , b α 1 , α 2 ; β 1 , β 2 γ ; λ corresponds, in the Laplace domain, to multiplication by p β 1 q β 2 1 λ p α 1 q α 2 γ . The latter operation clearly has a semigroup property in the parameters β 1 , β 2 , γ , since these appear only as exponents. Thus, the desired result is clear in this case.
For the general case, we must proceed by manipulation of infinite series and gamma functions, using Theorem 3:
I a , b α 1 , α 2 ; β 1 , β 2 γ ; λ I a , b α 1 , α 2 ; ε 1 , ε 2 σ ; λ ( f ) = k = 0 ( γ ) k λ k k ! I b y β 2 + α 2 k I a x β 1 + α 1 k m = 0 ( σ ) m λ m m ! I b y ε 2 + α 2 m I a x ε 1 + α 1 m ( f ) = k = 0 m = 0 ( γ ) k ( σ ) m λ k + m k ! m ! I b y β 2 + ε 2 + α 2 ( k + m ) I a x β 1 + ε 1 + α 1 ( k + m ) ( f ) = n = 0 k + m = n ( γ ) k ( σ ) m k ! m ! λ n I b y β 2 + ε 2 + α 2 n I a x β 1 + ε 1 + α 1 n ( f ) = n = 0 ( γ + σ ) n n ! λ n I b y β 2 + ε 2 + α 2 n I a x β 1 + ε 1 + α 1 n ( f ) = I a , b α 1 , α 2 ; β 1 + ε 1 , β 2 + ε 2 γ + σ ; λ ( f ) ,
where for the part in square brackets we use a finite-sum identity on gamma functions (see [12], Theorem 2.9)). □
The semigroup property helps us to obtain the left inverse of the bivariate integral operator constructed above, which will motivate the definition of a fractional derivative operator based on the five-parameter Mittag-Leffler function.
Theorem 6.
Let α 1 , α 2 , β 1 , β 2 , γ , λ C with Re ( α 1 ) , Re ( α 2 ) , Re ( β 1 ) , Re ( β 2 ) > 0 . For any ε 1 , ε 2 C with Re ( ε 1 ) , Re ( ε 2 ) > 0 , the following operator is a left inverse to the bivariate fractional integral operator I a , b α 1 , α 2 ; β 1 , β 2 γ ; λ considered above:
D a x β 1 + ε 1 D b y β 2 + ε 2 I a , b α 1 , α 2 ; ε 1 , ε 2 γ ; λ .
Proof. 
Using the semigroup property, we have
I a , b α 1 , α 2 ; ε 1 , ε 2 γ ; λ I a , b α 1 , α 2 ; β 1 , β 2 γ ; λ = I a , b α 1 , α 2 ; β 1 + ε 1 , β 2 + ε 2 0 ; λ = I b y β 2 + ε 2 I a x β 1 + ε 1 ,
using the fact that the Pochhammer symbol ( 0 ) k equals 1 if k = 0 and equals 0 for all integer k 1 . Therefore, applying the double Riemann–Liouville fractional differential operator D a x β 1 + ε 1 D b y β 2 + ε 2 from the left on both sides of the equation, the right-hand side gives the identity operator:
D a x β 1 + ε 1 D b y β 2 + ε 2 I a , b α 1 , α 2 ; ε 1 , ε 2 γ ; λ I a , b α 1 , α 2 ; β 1 , β 2 γ ; λ = D a x β 1 + ε 1 D b y β 2 + ε 2 I b y β 2 + ε 2 I a x β 1 + ε 1 = I .
Thus, we have found the left inverse operator of I a , b α 1 , α 2 ; β 1 , β 2 γ ; λ as required. □
Remark 1.
It should be remarked that the left inverse operator is independent of the parameters ε 1 and ε 2 . The easiest way to see this is by using the series formula:
D a x β 1 + ε 1 D b y β 2 + ε 2 I a , b α 1 , α 2 ; ε 1 , ε 2 γ ; λ ( f ) ( x , y ) = D a x β 1 + ε 1 D b y β 2 + ε 2 k = 0 ( γ ) k λ k k ! I b y ε 2 + α 2 k I a x ε 1 + α 1 k ( f ) ( x , y ) = k = 0 ( γ ) k λ k k ! D a x β 1 + ε 1 D b y β 2 + ε 2 I b y ε 2 + α 2 k I a x ε 1 + α 1 k ( f ) ( x , y ) = k = 0 ( γ ) k λ k k ! I b y β 2 + α 2 k I a x β 1 + α 1 k ( f ) ( x , y ) ,
which is independent of ε 1 and ε 2 . Here, we make use of the semigroup property for Riemann–Liouville fractional calculus, in the case where the inner operator is a fractional integral, and we also use the convention (valid by analytic continuation in the order of integration) that a Riemann–Liouville fractional integral to negative order means a Riemann–Liouville fractional derivative, I c t ν = D c t ν .
Since we can choose any values of ε 1 and ε 2 with positive real part and get the same left inverse operator, we opt for the values which give ordinary (non-fractional) derivatives of order β 1 + ε 1 and β 2 + ε 2 . This means choosing, just like in Riemann–Liouville fractional calculus, ε 1 = N 1 β 1 and ε 2 = N 2 β 2 , to obtain the following definition.
Definition 4.
Let α 1 , α 2 , β 1 , β 2 , γ , λ C with Re ( α 1 ) > 0 , Re ( α 2 ) > 0 , Re ( β 1 ) 0 , Re ( β 2 ) 0 . Define N 1 , N 2 N to be the numbers such that N 1 1 Re ( β 1 ) < N 1 and N 2 1 Re ( β 2 ) < N 2 . Then, the double Riemann–Liouville-type fractional derivative with five-parameter Mittag-Leffler kernel is defined by
D a , b α 1 , α 2 ; β 1 , β 2 γ ; λ ( f ) ( x , y ) = N 1 + N 2 x N 1 y N 2 I a , b α 1 , α 2 ; N 1 β 1 , N 2 β 2 γ ; λ ( f ) ( x , y ) ,
while the double Caputo-type fractional derivative with five-parameter Mittag-Leffler kernel is defined by
D a , b C α 1 , α 2 ; β 1 , β 2 γ ; λ ( f ) ( x , y ) = I a , b α 1 , α 2 ; N 1 β 1 , N 2 β 2 γ ; λ N 1 + N 2 x N 1 y N 2 f ( x , y ) .
Theorem 7.
Let α 1 , α 2 , β 1 , β 2 , γ , λ C with Re ( α 1 ) > 0 , Re ( α 2 ) > 0 , Re ( β 1 ) > 0 , Re ( β 2 ) > 0 . The fractional integrals and derivatives with five-parameter Mittag-Leffler kernel have the following inversion properties:
D a , b α 1 , α 2 ; β 1 , β 2 γ ; λ I a , b α 1 , α 2 ; β 1 , β 2 γ ; λ ( f ) ( x , y ) = f ( x , y ) , I a , b α 1 , α 2 ; β 1 , β 2 γ ; λ D a , b C α 1 , α 2 ; β 1 , β 2 γ ; λ ( f ) ( x , y ) = f ( x , y ) n = 0 N 1 1 ( x a ) n n ! · n x n f ( a , y ) m = 0 N 2 1 ( y b ) m m ! · m y m f ( x , b ) + n = 0 N 1 1 m = 0 N 2 1 ( x a ) n ( y b ) m n ! m ! · n + m x n y m f ( a , b ) .
Proof. 
The first identity is proved during the construction of the left inverse. For the second one, we again use the semigroup property given by Theorem 5:
I a , b α 1 , α 2 ; β 1 , β 2 γ ; λ D a , b C α 1 , α 2 ; β 1 , β 2 γ ; λ ( f ) ( x , y ) = I a , b α 1 , α 2 ; β 1 , β 2 γ ; λ I a , b α 1 , α 2 ; N 1 β 1 , N 2 β 2 γ ; λ N 1 + N 2 x N 1 y N 2 f ( x , y ) = I a , b α 1 , α 2 ; N 1 , N 2 0 ; λ N 1 + N 2 x N 1 y N 2 f ( x , y ) = I a x N 1 I b y N 2 N 1 x N 1 N 2 y N 2 f ( x , y ) = f ( x , y ) n = 0 N 1 1 ( x a ) n n ! · n x n f ( a , y ) m = 0 N 2 1 ( y b ) m m ! · m y m f ( x , b ) + n = 0 N 1 1 m = 0 N 2 1 ( x a ) n ( y b ) m n ! m ! · n + m x n y m f ( a , b ) ,
by using twice the formula for the nth integral of the nth derivative. □
Proposition 1.
Let α 1 , α 2 , β 1 , β 2 , γ , λ C with Re ( α 1 ) > 0 , Re ( α 2 ) > 0 , Re ( β 1 ) 0 , Re ( β 2 ) 0 . Then, the fractional derivatives with five-parameter Mittag-Leffler kernel, of both Riemann–Liouville type and Caputo type, can be expressed by series formulae as follows:
D a , b α 1 , α 2 ; β 1 , β 2 γ ; λ ( f ) ( x , y ) = k = 0 ( γ ) k λ k k ! I b y β 2 + α 2 k I a x β 1 + α 1 k ( f ) ( x , y ) ;
D a , b C α 1 , α 2 ; β 1 , β 2 γ ; λ ( f ) ( x , y ) = k = 0 ( γ ) k λ k k ! I b y N 2 β 2 + α 2 k I a x N 1 β 1 + α 1 k N 1 + N 2 x N 1 y N 2 f ( x , y ) .
Proof. 
This follows immediately from Theorem 3 and the series formula (11). □
Remark 2.
From comparing the series formula (10) for fractional integrals with the series formula (13) for Riemann–Liouville-type fractional derivatives, it is now clear that the latter is an analytic continuation of the former, with β 1 , β 2 , γ replaced by β 1 , β 2 , γ . To see this, we make use of the fact, alluded to in Remark 1 above, that the Riemann–Liouville fractional derivative D c t ν f ( t ) = I c t ν f ( t ) , Re ( ν ) 0 , is the analytic continuation in ν of the fractional integral I c t ν f ( t ) , Re ( ν ) > 0 . Therefore, each term I b y β 2 + α 2 k I a x β 1 + α 1 k ( f ) ( x , y ) appearing in (13) is exactly the same, under the analytic continuation of Riemann–Liouville differintegrals, as the term I b y β 2 + α 2 k I a x β 1 + α 1 k ( f ) ( x , y ) appearing in (10), after replacing β 1 , β 2 , γ by β 1 , β 2 , γ as stated.
Thus, we can adopt the notational convention that
I a , b α 1 , α 2 ; β 1 , β 2 γ ; λ ( f ) ( x , y ) = D a , b α 1 , α 2 ; β 1 , β 2 γ ; λ ( f ) ( x , y ) ,
and use this to extend the meaning of both I a , b α 1 , α 2 ; β 1 , β 2 γ ; λ ( f ) ( x , y ) and D a , b α 1 , α 2 ; β 1 , β 2 γ ; λ ( f ) ( x , y ) to the entire complex plane for the parameters β 1 , β 2 , without any need to impose conditions on the signs of their real parts. This identity is achieved by analytic continuation in the complex variables β 1 , β 2 from one half-plane to the other.
Proposition 2.
Let α 1 , α 2 , β 1 , β 2 , γ , λ C with Re ( α 1 ) > 0 , Re ( α 2 ) > 0 , Re ( β 1 ) 0 , Re ( β 2 ) 0 . Then, the fractional derivatives of Riemann–Liouville type and Caputo type, with five-parameter Mittag-Leffler kernel, are related to each other as follows:
D a , b C α 1 , α 2 ; β 1 , β 2 γ ; λ ( f ) ( x , y ) = D a , b α 1 , α 2 ; β 1 , β 2 γ ; λ ( f ) ( x , y ) k = 0 ( γ ) k λ k k ! n = 0 N 1 1 ( x a ) n β 1 + α 1 k Γ ( n β 1 + α 1 k + 1 ) · I b y β 2 + α 2 k n x n f ( a , y ) k = 0 ( γ ) k λ k k ! m = 0 N 2 1 ( y b ) m β 2 + α 2 k Γ ( m β 2 + α 2 k + 1 ) · I a x β 1 + α 1 k m y m f ( x , b ) + k = 0 ( γ ) k λ k k ! n = 0 N 1 1 m = 0 N 2 1 ( x a ) n β 1 + α 1 k ( y b ) m β 2 + α 2 k Γ ( n β 1 + α 1 k + 1 ) Γ ( m β 2 + α 2 k + 1 ) · n + m x n y m f ( a , b ) .
Proof. 
Starting from the series formulae (13) and (14), we have:
D a , b C α 1 , α 2 ; β 1 , β 2 γ ; λ ( f ) ( x , y ) = k = 0 ( γ ) k λ k k ! I b y β 2 + α 2 k I a x β 1 + α 1 k I b y N 2 I a x N 1 N 1 + N 2 x N 1 y N 2 f ( x , y ) = k = 0 ( γ ) k λ k k ! I b y β 2 + α 2 k I a x β 1 + α 1 k ( f ( x , y ) n = 0 N 1 1 ( x a ) n n ! · n x n f ( a , y ) m = 0 N 2 1 ( y b ) m m ! · m y m f ( x , b ) + n = 0 N 1 1 m = 0 N 2 1 ( x a ) n ( y b ) m n ! m ! · n + m x n y m f ( a , b ) ) = D a , b α 1 , α 2 ; β 1 , β 2 γ ; λ ( f ) ( x , y ) k = 0 ( γ ) k λ k k ! n = 0 N 1 1 ( x a ) n β 1 + α 1 k Γ ( n β 1 + α 1 k + 1 ) · I b y β 2 + α 2 k n x n f ( a , y ) k = 0 ( γ ) k λ k k ! m = 0 N 2 1 ( y b ) m β 2 + α 2 k Γ ( m β 2 + α 2 k + 1 ) · I a x β 1 + α 1 k m y m f ( x , b ) + k = 0 ( γ ) k λ k k ! n = 0 N 1 1 m = 0 N 2 1 ( x a ) n β 1 + α 1 k ( y b ) m β 2 + α 2 k Γ ( n β 1 + α 1 k + 1 ) Γ ( m β 2 + α 2 k + 1 ) · n + m x n y m f ( a , b ) ,
where we use both the formula (12) for the double nth integral of the double nth derivative and also the well-known formulae for Riemann–Liouville fractional differintegrals of power functions. □

4. The Non-Singular Cases of the Operators

Recently, the second author [11] made a detailed study of Prabhakar fractional calculus, separating this class of operators into several subclasses according to their properties. Of particular importance is the consideration of whether an operator is singular or non-singular, and whether the series formula expressing it in terms of Riemann–Liouville integrals is a finite or infinite sum. The same considerations can be applied to other types of fractional calculus, such as the one we are studying here.
In the case of Prabhakar, it was proved [11] (Theorem 4.5) that the most special case is when the operators are non-singular and the sum is finite, and this subclass of Prabhakar fractional calculus consists precisely of the integer powers of the following (inverse) operators:
M f ( t ) = f ( t ) ω I a t α f ( t ) , N f ( t ) = d d t a t E α ω ( t τ ) α f ( τ ) d τ .
These operators are obtained [47] by setting β = 0 (for non-singular operators) and ρ = ± 1 (for finite sums) in the operator I a P t α , β , ρ , ω of Prabhakar integration. We may try to do something similar for the bivariate operator of integration with five-parameter Mittag-Leffler kernel.
Let us firstly compare the Prabhakar fractional integral with the bivariate fractional integral considered in this paper. We have
I a P t α , β , ρ , ω f ( t ) = a t f ( τ ) ( t τ ) β 1 E α , β ρ ω ( t τ ) α d τ = k = 0 ( ρ ) k ω k k ! I a t β + α k f ( t ) ,
versus
I a , b α 1 , α 2 ; β 1 , β 2 γ ; λ ( f ) ( x , y ) = a x b y f ( t , s ) ( x t ) β 1 1 ( y s ) β 2 1 E α 1 , α 2 ; β 1 , β 2 γ ( λ ( x t ) α 1 ( y s ) α 2 ) d s d t = k = 0 ( γ ) k λ k k ! I b y β 2 + α 2 k I a x β 1 + α 1 k ( f ) ( x , y ) .
It is clear that our parameters α 1 , α 2 correspond to the α of Prabhakar, our β 1 , β 2 correspond to the β of Prabhakar, our γ corresponds to the ρ of Prabhakar, and our λ corresponds to the ω of Prabhakar, using all notation as above.
Therefore, for our operator, the process of obtaining a non-singular finite-sum version should involve setting β 1 = β 2 = 0 and γ = ± 1 . Note that, since Definition 3 for the fractional integral operator requires Re ( β 1 ) > 0 and Re ( β 2 ) > 0 , we must use Definition 4 for the case when β 1 = β + 2 = 0 . We find:
I a , b α 1 , α 2 ; 0 , 0 1 ; λ ( f ) ( x , y ) = D a , b α 1 , α 2 ; 0 , 0 1 ; λ ( f ) ( x , y ) = 2 x y I a , b α 1 , α 2 ; 1 , 1 1 ; λ ( f ) ( x , y ) = 2 x y a x b y f ( t , s ) E α 1 , α 2 ; 1 , 1 1 ( λ ( x t ) α 1 ( y s ) α 2 ) d s d t = 2 x y a x b y f ( t , s ) 1 λ ( x t ) α 1 ( y s ) α 2 Γ ( α 1 + 1 ) Γ ( α 2 + 1 ) d s d t = f ( x , y ) λ 2 x y I b y α 2 + 1 I a x α 1 + 1 ( f ) ( x , y ) = f ( x , y ) λ I b y α 2 I a x α 1 ( f ) ( x , y ) .
The same formula can also be obtained more directly from the series formula (13):
I a , b α 1 , α 2 ; 0 , 0 1 ; λ ( f ) ( x , y ) = k = 0 ( 1 ) k λ k k ! I b y α 2 k I a x α 1 k ( f ) ( x , y ) = f ( x , y ) λ I b y α 2 I a x α 1 ( f ) ( x , y ) ,
where again we use the fact that ( 1 ) k equals 1 if k = 0 , 1 if k = 1 , and 0 for all k 2 .
The inverse operator is given by setting γ = 1 instead of γ = 1 , but this one cannot be written in such an elementary way:
I a , b α 1 , α 2 ; 0 , 0 1 ; λ ( f ) ( x , y ) = D a , b α 1 , α 2 ; 0 , 0 1 ; λ ( f ) ( x , y ) = 2 x y I a , b α 1 , α 2 ; 1 , 1 1 ; λ ( f ) ( x , y ) = 2 x y a x b y f ( t , s ) E α 1 , α 2 ; 1 , 1 1 ( λ ( x t ) α 1 ( y s ) α 2 ) d s d t = 2 x y a x b y f ( t , s ) k = 0 λ k ( x t ) α 1 k ( y s ) α 2 k Γ α 1 k + 1 Γ α 2 k + 1 d s d t ,
or equivalently, using the series formula (13),
I a , b α 1 , α 2 ; 0 , 0 1 ; λ ( f ) ( x , y ) = k = 0 λ k I b y α 2 k I a x α 1 k ( f ) ( x , y ) .
The pair of integro-differential operators given by (16) and (17) is of course reminiscent of the so-called Atangana–Baleanu (AB) integral and derivative [48,49]. The latter are defined similarly: the integral by a linear combination of a function with its Riemann–Liouville fractional integral and the derivative by a derivative of an integral transform involving a Mittag-Leffler function kernel. The new development here is that the operators are now bivariate: we can think of the operators (16) and (17) as forming a two-dimensional Atangana–Baleanu calculus.
We also note that there is now a direct connection with the bivariate Abel equation (1) which we considered at first to motivate this paper. The aforementioned Abel equation was rewritten in a form involving the operator (16), and its solution was constructed in a form involving the inverse operator (17).
It is also important to realise that both operators (16) and (17) involve mixed partial integro-differential operators: the AB-type integral operator (16) is defined using fractional integrals with respect to x and y together, and the AB-type derivative operator (17) is defined using mixed partial derivatives with respect to x and y. Thus, these operators may be useful in the study and understanding of fractional PDEs involving mixed partial derivatives, which we note above are under-appreciated in fractional calculus.
To fully recover a bivariate analogue of the AB operators, we need to choose a value for the parameter λ so that appropriate boundary conditions are realised when the parameters α 1 , α 2 go to 0 or 1. This is the reason for the choice of multipliers in the definition of the AB integral and derivative to order α : so that the ordinary first-order integral and derivative are recovered when α 1 and the original function itself when α 0 . Thus, the following definition is motivated.
Definition 5.
The mixed bivariate AB integral is defined to be
I a , b A B x , y α 1 , α 2 f ( x , y ) = ( 1 α 1 ) ( 1 α 2 ) B ( α 1 , α 2 ) f ( x , y ) + α 1 α 2 B ( α 1 , α 2 ) I b y α 2 I a x α 1 ( f ) ( x , y ) ,
where B ( α 1 , α 2 ) is a normalisation function which is assumed to satisfy B ( 0 , 0 ) = 1 and B ( 1 , 1 ) = 1 .
The mixed bivariate AB derivatives, of Riemann–Liouville and Caputo type, respectively, are defined to be
D a , b A B R x , y α 1 , α 2 f ( x , y ) = B ( α 1 , α 2 ) ( 1 α 1 ) ( 1 α 2 ) · 2 x y a x b y f ( t , s ) E α 1 , α 2 ; 1 , 1 1 α 1 α 2 ( 1 α 1 ) ( 1 α 2 ) ( x t ) α 1 ( y s ) α 2 d s d t
and
D a , b A B C x , y α 1 , α 2 f ( x , y ) = B ( α 1 , α 2 ) ( 1 α 1 ) ( 1 α 2 ) a x b y 2 f ( t , s ) t s · E α 1 , α 2 ; 1 , 1 1 α 1 α 2 ( 1 α 1 ) ( 1 α 2 ) ( x t ) α 1 ( y s ) α 2 d s d t ,
where B ( α 1 , α 2 ) is the same normalisation function as above.
These definitions are valid for any α 1 , α 2 C with positive real parts, although the main case of interest is when α 1 , α 2 ( 0 , 1 ) .
Remark 3.
The multiplier function B ( α 1 , α 2 ) is introduced by analogy with the original definition of the AB derivative and AB integral [48], and the restrictions on this function are imposed in order to ensure the following limiting behaviour as the fractional orders α 1 , α 2 approach 0 or 1:
lim α 1 , α 2 0 I a , b A B x , y α 1 , α 2 f ( x , y ) = f ( x , y ) , lim α 1 , α 2 1 I a , b A B x , y α 1 , α 2 f ( x , y ) = I b y 1 I a x 1 ( f ) ( x , y ) , lim α 1 , α 2 0 D a , b A B R x , y α 1 , α 2 f ( x , y ) = f ( x , y ) , lim α 1 , α 2 0 D a , b A B C x , y α 1 , α 2 f ( x , y ) = f ( x , y ) f ( a , y ) f ( b , x ) + f ( a , b ) .
Note that if the pair ( α 1 , α 2 ) takes the values ( 0 , 1 ) or ( 1 , 0 ) , then the mixed bivariate AB integral becomes zero. This is why we call it specifically amixedbivariate operator: in the square set [ 0 , 1 ] × [ 0 , 1 ] for values of the pair ( α 1 , α 2 ) , the values of the mixed bivariate AB integral are weighted towards the diagonal (where both the fractional orders α 1 and α 2 have more similar values) rather than towards the asymmetric corners.
The following result is the bivariate analogue of [47] (Theorem 3.1) or [11] (Proposition 2.4). It follows directly from our consideration above of the non-singular finite-sum special case of the five-parameter Mittag-Leffler kernel operators which culminated in Equations (16) and (17).
Proposition 3.
The mixed bivariate AB integral and derivatives can be written in terms of the five-parameter Mittag-Leffler kernel operators as follows:
I a , b A B x , y α 1 , α 2 f ( x , y ) = ( 1 α 1 ) ( 1 α 2 ) B ( α 1 , α 2 ) · I a , b α 1 , α 2 ; 0 , 0 1 ; α 1 α 2 ( 1 α 1 ) ( 1 α 2 ) ( f ) ( x , y ) = ( 1 α 1 ) ( 1 α 2 ) B ( α 1 , α 2 ) · D a , b α 1 , α 2 ; 0 , 0 1 ; α 1 α 2 ( 1 α 1 ) ( 1 α 2 ) ( f ) ( x , y ) ; D a , b A B R x , y α 1 , α 2 f ( x , y ) = B ( α 1 , α 2 ) ( 1 α 1 ) ( 1 α 2 ) · D a , b α 1 , α 2 ; 0 , 0 1 ; α 1 α 2 ( 1 α 1 ) ( 1 α 2 ) ( f ) ( x , y ) = B ( α 1 , α 2 ) ( 1 α 1 ) ( 1 α 2 ) · I a , b α 1 , α 2 ; 0 , 0 1 ; α 1 α 2 ( 1 α 1 ) ( 1 α 2 ) ( f ) ( x , y ) ; D a , b A B C x , y α 1 , α 2 f ( x , y ) = B ( α 1 , α 2 ) ( 1 α 1 ) ( 1 α 2 ) · D a , b C α 1 , α 2 ; 0 , 0 1 ; α 1 α 2 ( 1 α 1 ) ( 1 α 2 ) ( f ) ( x , y ) .
However, due to the non-singularity properties of the AB-type operators, they can be defined on a larger function space than the general fractional derivatives with five-parameter Mittag-Leffler kernel: no differentiability assumptions are required. The following result is the bivariate analogue in [49] (Lemma 2.1) which was the first result to establish appropriate function spaces for the AB derivatives.
Theorem 8.
The mixed bivariate AB integral I a , b A B x , y α 1 , α 2 f ( x , y ) is defined for any function f L 1 ( [ a , c ] × [ b , d ] ) . The mixed bivariate ABR derivative D a , b A B R x , y α 1 , α 2 f ( x , y ) can also be defined for any function f L 1 ( [ a , c ] × [ b , d ] ) , while its ABC counterpart D a , b A B C x , y α 1 , α 2 f ( x , y ) is defined for any twice differentiable function f on [ a , c ] × [ b , d ] such that 2 x y f L 1 ( [ a , c ] × [ b , d ] ) .
Proof. 
For the mixed bivariate AB integral, the result is clear, since this is just a linear combination of f ( x , y ) with its double Riemann–Liouville integral.
For the mixed bivariate ABC derivative, the result follows from Theorem 2, since this operator is defined by applying a special case of the five-parameter Mittag-Leffler kernel operator to the function 2 x y f ( x , y ) .
It remains to consider the mixed bivariate ABR derivative, which we can simplify in the following way since the kernel is non-singular, as a bivariate version of the arguments used in [11,49]:
( 1 α 1 ) ( 1 α 2 ) B ( α 1 , α 2 ) · D a , b A B R x , y α 1 , α 2 f ( x , y ) = 2 x y a x b y f ( t , s ) E α 1 , α 2 ; 1 , 1 1 α 1 α 2 ( 1 α 1 ) ( 1 α 2 ) ( x t ) α 1 ( y s ) α 2 d s d t = f ( x , y ) + a x f ( t , s ) x E α 1 , α 2 ; 1 , 1 1 α 1 α 2 ( 1 α 1 ) ( 1 α 2 ) ( x t ) α 1 ( y s ) α 2 s = y d t + b y f ( t , s ) y E α 1 , α 2 ; 1 , 1 1 α 1 α 2 ( 1 α 1 ) ( 1 α 2 ) ( x t ) α 1 ( y s ) α 2 t = x d s + a x b y f ( t , s ) 2 x y E α 1 , α 2 ; 1 , 1 1 α 1 α 2 ( 1 α 1 ) ( 1 α 2 ) ( x t ) α 1 ( y s ) α 2 d s d t = f ( x , y ) + a x b y f ( t , s ) 2 x y E α 1 , α 2 ; 1 , 1 1 α 1 α 2 ( 1 α 1 ) ( 1 α 2 ) ( x t ) α 1 ( y s ) α 2 d s d t ,
where the new kernel function is
2 x y E α 1 , α 2 ; 1 , 1 1 α 1 α 2 ( 1 α 1 ) ( 1 α 2 ) ( x t ) α 1 ( y s ) α 2 = k = 1 ( x t ) α 1 k 1 ( y s ) α 2 k 1 Γ ( α 1 k ) Γ ( α 2 k ) ( x t ) α 1 1 ( y s ) α 2 1 Γ ( α 1 ) Γ ( α 2 ) ,
and therefore integrable, as x t 0 and y s 0 . Thus, the ABR derivative of f ( x , y ) equals the function f ( x , y ) plus an integral term which behaves as a double Riemann–Liouville integral of f ( x , y ) in the singular limit x t 0 , y s 0 . This means the operator is well-defined for any f L 1 ( [ a , c ] × [ b , d ] ) , as required. □
The following result is the bivariate analogue of the original series formulae for AB derivatives, [49] (Theorems 2.1 and 2.2). It follows directly from the work we did above to obtain Equation (18).
Proposition 4.
The mixed bivariate AB derivatives can be given by the following infinite series of double Riemann–Liouville fractional integrals:
D a , b A B R x , y α 1 , α 2 f ( x , y ) = B ( α 1 , α 2 ) ( 1 α 1 ) ( 1 α 2 ) k = 0 α 1 α 2 ( 1 α 1 ) ( 1 α 2 ) k I b y α 2 k I a x α 1 k ( f ) ( x , y ) , D a , b A B C x , y α 1 , α 2 f ( x , y ) = B ( α 1 , α 2 ) ( 1 α 1 ) ( 1 α 2 ) k = 0 α 1 α 2 ( 1 α 1 ) ( 1 α 2 ) k I b y α 2 k + 1 I a x α 1 k + 1 2 f x y ( x , y ) .
where this series is locally uniformly convergent for any f L 1 ( [ a , c ] × [ b , d ] ) .
Corollary 1.
From the series formulae above, it is clear that the mixed bivariate ABC derivative can be written in terms of its ABR counterpart as follows:
D a , b A B C x , y α 1 , α 2 f ( x , y ) = D a , b A B R x , y α 1 , α 2 f ( x , y ) f ( a , y ) f ( b , x ) + f ( a , b ) .
Therefore, the domain of this operator can be extended from the one mentioned in Theorem 8, and the mixed bivariate ABC derivative can also be defined on the whole space L 1 ( [ a , c ] × [ b , d ] ) , without any differentiability conditions.
Theorem 9.
The mixed bivariate AB integral is both a left and right inverse of the mixed bivariate ABR derivative on the space L 1 ( [ a , c ] × [ b , d ] ) , and it shares the following relationship with the mixed bivariate ABC derivative:
I a , b A B x , y α 1 , α 2 D a , b A B C x , y α 1 , α 2 f ( x , y ) = f ( x , y ) f ( a , y ) f ( b , x ) + f ( a , b ) .
Proof. 
This is immediate from Proposition 4 and Corollary 1. Note that the series for the ABR derivative consists only of Riemann–Liouville fractional integrals, which have a semigroup property unlike their fractional derivative counterparts. □
The following result is the bivariate analogue of the Laplace transform for AB derivatives [48,49].
Proposition 5.
The double Laplace transforms of the mixed bivariate AB integral and derivatives, in the case a = b = 0 , can be expressed as follows:
L 2 I 0 , 0 A B x , y α 1 , α 2 f ( x , y ) ( p , q ) = ( 1 α 1 ) ( 1 α 2 ) B ( α 1 , α 2 ) 1 + α 1 α 2 ( 1 α 1 ) ( 1 α 2 ) p α 1 q α 2 L 2 [ f ] ( p , q ) ; L 2 D 0 , 0 A B R x , y α 1 , α 2 f ( x , y ) ( p , q ) = B ( α 1 , α 2 ) ( 1 α 1 ) ( 1 α 2 ) 1 + α 1 α 2 ( 1 α 1 ) ( 1 α 2 ) p α 1 q α 2 1 L 2 [ f ] ( p , q ) ; L 2 D 0 , 0 A B C x , y α 1 , α 2 f ( x , y ) ( p , q ) = B ( α 1 , α 2 ) ( 1 α 1 ) ( 1 α 2 ) 1 + α 1 α 2 ( 1 α 1 ) ( 1 α 2 ) p α 1 q α 2 1 × L 2 [ f ] ( p , q ) 1 p L y q [ f ] ( 0 , q ) 1 q L x p [ f ] ( p , 0 ) + 1 p q f ( 0 , 0 ) .
Proof. 
For the mixed bivariate AB integral, the result follows directly from the definition and the known facts on Laplace transforms of Riemann–Liouville fractional integrals.
For the mixed bivariate AB derivative of Riemann–Liouville type, the result follows from the series formula of Proposition 4 and a simple application of the binomial theorem. An important fact to notice here is that, via the series formula, this Riemann–Liouville type operator can be written solely in terms of Riemann–Liouville integrals, so there is no need to involve initial conditions here.
For the mixed bivariate AB derivative of Caputo type, derivatives and therefore initial conditions become involved. We use the correct form for the double Laplace transform of a mixed partial derivative, given in [46,50], to achieve the result:
L 2 D a , b A B C x , y α 1 , α 2 f ( x , y ) ( p , q ) = B ( α 1 , α 2 ) ( 1 α 1 ) ( 1 α 2 ) k = 0 α 1 α 2 ( 1 α 1 ) ( 1 α 2 ) k q α 2 k 1 p α 1 k 1 L 2 2 f x y ( x , y ) ( p , q ) = B ( α 1 , α 2 ) ( 1 α 1 ) ( 1 α 2 ) 1 + α 1 α 2 ( 1 α 1 ) ( 1 α 2 ) p α 1 q α 2 1 p 1 q 1 ( p q L 2 [ f ] ( p , q ) q L y q [ f ] ( 0 , q ) p L x p [ f ] ( p , 0 ) + f ( 0 , 0 ) ) ,
which trivially rearranges to the stated result. □

5. Discussion and Conclusions

In this paper, we construct a type of bivariate fractional calculus based originally on a function of a single variable defined by a single power series. We prove some of the fundamental properties of this type of fractional calculus: boundedness of operators, series formulae, Laplace transforms, semigroup and inversion properties, etc. As a special case, we also consider the non-singular versions of our bivariate operators, in order to construct a two-dimensional version of the well-known Atangana–Baleanu calculus.
Some of the ideas used in the above work are similar to those for previously existing models of fractional calculus with various Mittag-Leffler type kernels; however, this is the first time that a univariate function has been used in this way to construct bivariate fractional-calculus operators. Previous contributions in this direction have included using univariate single series to construct univariate operators [10], using bivariate double series to construct univariate operators [18], and using bivariate double series to construct bivariate operators [23,24].
At first glance, it may seem that our construction is unnatural: Given a function of one variable defined by a single series, why would one make a simple thing more complicated by replacing z with x α 1 y α 2 and introducing a bivariate integral operator? The answer is that the underlying mathematical structure is still simple, because everything is defined by a single series whose convergence is easy to describe, but the range of problems that can be tackled is now richer and more diverse, because more variables and parameters allow for more flexibility in adapting the operators to particular scenarios. Our “trick” of turning a univariate function into a bivariate operator gives a shortcut, a possibility of modelling complicated problems using simpler tools.
The bivariate operators that we define lend themselves more to partial differential equations than to ordinary differential equations. Even our initial motivation for proposing them, in the first section of this paper, arose from a bivariate Abel-type equation for a function of two variables. In the literature thus far, Mittag-Leffler kernel operators have mostly been applied to modelling problems with ordinary differential equations. While the existing operators could of course be combined to trivially create bivariate derivatives and integrals with Mittag-Leffler behaviour, we believe that a richer structure will emerge from considering operators such as those here, where the behaviours with respect to x and y are intertwined so that the operator cannot simply be broken down into one with respect to x and another with respect to y. We believe that our work here may have useful ramifications in the study of fractional partial differential equations in 2 or 2 + 1 dimensions.
For higher dimensions, it will be possible to construct similarly a model of trivariate or multivariate fractional calculus based on a univariate function defined by a single series with 3 or n different gamma functions on the denominator. We give here, without justification, the trivariate version.
Starting from the following function, defined by a series convergent for parameters α i , β i , γ C ( i = 1 , 2 , 3 ) with Re ( α 1 + α 2 + α 3 ) > 0 :
E α 1 , α 2 , α 3 ; β 1 , β 2 , β 3 γ ( z ) = k = 0 ( γ ) k Γ α 1 k + β 1 Γ α 2 k + β 2 Γ α 3 k + β k · z k k ! ,
it is possible to define the following trivariate fractional integral operator, convergent for any function f L 1 ( [ a 1 , b 1 ] × [ a 2 , b 2 ] × [ a 3 , b 3 ] ) and parameters α i , β i , γ , λ C with Re ( α i ) > 0 and Re ( β i ) > 0 for i = 1 , 2 , 3 :
I a 1 , a 2 , a 3 α 1 , α 2 , α 3 ; β 1 , β 2 , β 3 γ ; λ ( f ) ( x , y )       : = a 1 x a 2 y a 3 z f ( s , t , u ) ( x s ) β 1 1 ( y t ) β 2 1 ( z u ) β 3 1            × E α 1 , α 2 , α 3 ; β 1 , β 2 , β 3 γ λ ( x s ) α 1 ( y t ) α 2 ( z u ) α 3 d u d t d s ,
and then the corresponding trivariate fractional derivative operator are defined as follows for parameters α i , β i , γ , λ C with Re ( α i ) > 0 and Re ( β i ) 0 for i = 1 , 2 , 3 :
D a 1 , a 2 , a 3 α 1 , α 2 , α 3 ; β 1 , β 2 , β 3 γ ; λ ( f ) ( x , y ) = N 1 + N 2 + N 3 x N 1 y N 2 z N 3 D a 1 , a 2 , a 3 α 1 , α 2 , α 3 ; N 1 β 1 , N 2 β 2 , N 3 β 3 γ ; λ f ( x , y ) .
The detailed analysis and investigation of these functions and operators, including verification of their essential properties and discussion of special cases and applications, is left for a future research project.
Other related directions of research may include a deeper study of the function spaces on which the integral and derivative operators considered here can be defined. For example, Theorem 2 can be extended easily to show that the bivariate fractional integral operator is bounded on an L p space. Then, a measure-theoretical approach may yield extensions to Morrey spaces [51], or a distributional approach may yield extensions to still larger spaces related to generalised integral operators [52]. Describing more deeply the various relevant function spaces may be useful in the qualitative theory of partial differential equations related to these operators, for example regularity theory [29].

Author Contributions

All authors contributed equally to this article. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gorenflo, R.; Kilbas, A.A.; Mainardi, F.; Rogosin, S.V. Mittag-Leffler Functions, Related Topics and Applications, 2nd ed.; Springer: Berlin, Germany, 2020. [Google Scholar]
  2. Paneva-Konovska, J. From Bessel to Multi-Index Mittag-Leffler Functions: Enumerable Families, Series in Them and Convergence; World Scientific: Singapore, 2016. [Google Scholar]
  3. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar]
  4. Shishkina, E.; Sitnik, S. Transmutations, Singular and Fractional Differential Equations with Applications to Mathematical Physics; Academic Press: Cambridge, MA, USA, 2020. [Google Scholar]
  5. Mathai, A.M.; Haubold, H.J. Mittag-Leffler Functions and Fractional Calculus. In Special Functions for Applied Scientists; Mathai, A.M., Haubold, H.J., Eds.; Springer: Berlin, Germany, 2008; pp. 79–134. [Google Scholar]
  6. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives: Theory and Applications; Gordon & Breach Science Publishers: Yverdon, Switzerland, 1993; Originally in Russian: Nauka i Tekhnika, Minsk, 1987. [Google Scholar]
  7. Baleanu, D.; Diethelm, K.; Scalas, E.; Trujillo, J.J. Fractional Calculus: Models and Numerical Methods, 2nd ed.; World Scientific: New York, NY, USA, 2017. [Google Scholar]
  8. Hilfer, R. (Ed.) Applications of Fractional Calculus in Physics; World Scientific: Singapore, 2000. [Google Scholar]
  9. Rabotnov, Y.N. Elements of Hereditary Solid Mechanics; Mir Publishers: Moscow, Russia, 1980. [Google Scholar]
  10. Prabhakar, T.R. A singular integral equation with a generalized Mittag Leffler function in the kernel. Yokohama Math. J. 1971, 19, 7–15. [Google Scholar]
  11. Fernandez, A.; Baleanu, D. Classes of Operators in Fractional Calculus: A Case Study. Math. Methods Appl. Sci. 2020, 1–20. [Google Scholar] [CrossRef]
  12. Fernandez, A.; Baleanu, D.; Srivastava, H.M. Series representations for models of fractional calculus involving generalised Mittag-Leffler functions. Commun. Nonlinear Sci. Numer. Simul. 2019, 67, 517–527. [Google Scholar] [CrossRef] [Green Version]
  13. Kilbas, A.A.; Saigo, M.; Saxena, R.K. Generalized Mittag-Leffler function and generalized fractional calculus operators. Integral Transform. Spec. Funct. 2004, 15, 31–49. [Google Scholar] [CrossRef]
  14. Hadid, S.B.; Luchko, Y.F. An operational method for solving fractional differential equations of an arbitrary real order. Panam. Math. J. 1996, 6, 57–73. [Google Scholar]
  15. Luchko, Y.F.; Gorenflo, R. An operational method for solving fractional differential equations. Acta Math. Vietnam. 1999, 24, 207–234. [Google Scholar]
  16. Saxena, R.K.; Kalla, S.L.; Saxena, R. Multivariate analogue of generalised Mittag-Leffler function. Integral Transform. Spec. Funct. 2011, 22, 533–548. [Google Scholar] [CrossRef]
  17. Özarslan, M.A.; Fernandez, A. On the fractional calculus of multivariate Mittag-Leffler functions. Int. J. Comput. Math. 2021, 1–27. [Google Scholar] [CrossRef]
  18. Fernandez, A.; Kürt, C.; Özarslan, M.A. A naturally emerging bivariate Mittag-Leffler function and associated fractional-calculus operators. Comput. Appl. Math. 2020, 39, 200. [Google Scholar] [CrossRef]
  19. Huseynov, I.T.; Ahmadova, A.; Ojo, G.O.; Mahmudov, N.I. A natural extension of Mittag-Leffler function associated with a triple infinite series. arXiv 2020, arXiv:2011.03999. [Google Scholar]
  20. Huseynov, I.T.; Ahmadova, A.; Fernandez, A.; Mahmudov, N.I. Explicit analytic solutions of incommensurate fractional differential equation systems. Appl. Math. Comput. 2021, 390, 125590. [Google Scholar]
  21. Ahmadova, A.; Huseynov, I.T.; Fernandez, A.; Mahmudov, N.I. Trivariate Mittag-Leffler functions used to solve multi-order systems of fractional differential equations. Commun. Nonlinear Sci. Numer. Simul. 2021, 97, 105735. [Google Scholar] [CrossRef]
  22. Garg, M.; Manohar, P.; Kalla, S.L. A Mittag-Leffler-type function of two variables. Integral Transform. Spec. Funct. 2013, 24, 934–944. [Google Scholar] [CrossRef]
  23. Özarslan, M.A.; Kürt, C. On a double integral equation including a set of two variables polynomials suggested by Laguerre polynomials. J. Comput. Anal. Appl. 2017, 22, 1198–1207. [Google Scholar]
  24. Özarslan, M.A.; Kürt, C. Bivariate Mittag-Leffler functions arising in the solutions of convolution integral equation with 2D-Laguerre-Konhauser polynomials in the kernel. Appl. Math. Comput. 2019, 347, 631–644. [Google Scholar] [CrossRef]
  25. Kiryakova, V. The multi-index Mittag-Leffler functions as an important class of special functions of fractional calculus. Comput. Math. Appl. 2010, 59, 1885–1895. [Google Scholar] [CrossRef] [Green Version]
  26. Kiryakova, V. Unified approach to fractional calculus images of special functions—A survey. Mathematics 2020, 8, 2260. [Google Scholar] [CrossRef]
  27. Kiryakova, V. A guide to special functions in fractional calculus. Mathematics 2021, 9, 106. [Google Scholar] [CrossRef]
  28. Fernandez, A.; Baleanu, D.; Fokas, A.S. Solving PDEs of fractional order using the unified transform method. Appl. Math. Comput. 2018, 339, 738–749. [Google Scholar] [CrossRef] [Green Version]
  29. Fernandez, A. An elliptic regularity theorem for fractional partial differential operators. Comput. Appl. Math. 2018, 37, 5542–5553. [Google Scholar] [CrossRef] [Green Version]
  30. Williams, P. Fractional Calculus of Schwartz Distributions. Bachelor’s Thesis, University of Melbourne, Melbourne, Australia, 2007. [Google Scholar]
  31. Djida, J.-D.; Area, I.; Nieto, J.J. Nonlocal time porous medium equation with fractional time derivative. Rev. Mat. Complut. 2019, 32, 273–304. [Google Scholar] [CrossRef] [Green Version]
  32. Djida, J.-D.; Area, I.; Nieto, J.J. Nonlocal time-porous medium equation: Weak solutions and finite speed of propagation. Discret. Contin. Dyn. Syst. Ser. B 2019, 24, 4031–4053. [Google Scholar] [CrossRef] [Green Version]
  33. Lischke, A.; Pang, G.; Gulian, M.; Song, F.; Glusa, C.; Zheng, X.; Mao, Z.; Cai, W.; Meerschaert, M.M.; Ainsworth, M.; et al. What is the fractional Laplacian? A comparative review with new results. J. Comput. Phys. 2020, 404, 109009. [Google Scholar] [CrossRef]
  34. Stinga, P.R. User’s guide to the fractional Laplacian and the method of semigroups. arXiv 2018, arXiv:1808.05159. [Google Scholar]
  35. Conlan, J. Hyperbolic Differential Equations of Generalized Order. Appl. Anal. 1983, 14, 167–177. [Google Scholar] [CrossRef]
  36. Jiang, J.; Feng, Y.; Li, S. Exact Solutions to the Fractional Differential Equations with Mixed Partial Derivatives. Axioms 2018, 7, 10. [Google Scholar] [CrossRef] [Green Version]
  37. Kilbas, A.A.; Repin, O.A. An Analog of the Tricomi Problem for a Mixed Type Equation with a Partial Fractional Derivative. Fract. Calc. Appl. Anal. 2010, 13, 69–84. [Google Scholar]
  38. Kilbas, A.A.; Saigo, M.; Trujillo, J.J. On the generalized Wright function. Fract. Calc. Appl. Anal. 2002, 5, 437–460. [Google Scholar]
  39. Fernandez, A.; Husain, I. Modified Mittag-Leffler functions with applications in complex formulae for fractional calculus. Fractal Fract. 2020, 4, 45. [Google Scholar] [CrossRef]
  40. Luchko, Y. The Four-Parameters Wright Function of the Second kind and its Applications in FC. Mathematics 2020, 8, 970. [Google Scholar] [CrossRef]
  41. Al-Bassam, M.A.; Luchko, Y.F. On generalized fractional calculus and its application to the solution of integro-differential equations. J. Fract. Calc. 1995, 7, 69–88. [Google Scholar]
  42. Choi, J.; Agarwal, P. A note on fractional integral operator associated with multiindex Mittag-Leffler functions. Filomat 2016, 30, 1931–1939. [Google Scholar] [CrossRef]
  43. Marichev, O.I. Handbook of Integral Transforms of Higher Transcendental Functions, Theory and Algorithmic Tables; Ellis Horwood: Chichester, UK, 1983. [Google Scholar]
  44. Miller, K.S.; Ross, B. An Introduction to the Fractional Calculus and Fractional Differential Equations; Wiley: New York, NY, USA, 1993. [Google Scholar]
  45. Oldham, K.B.; Spanier, J. The Fractional Calculus; Academic Press: San Diego, CA, USA, 1974. [Google Scholar]
  46. Kürt, C.; Özarslan, M.A.; Fernandez, A. On a certain bivariate Mittag-Leffler function analysed from a fractional-calculus point of view. Math. Methods Appl. Sci. 2021, 44, 2600–2620. [Google Scholar] [CrossRef]
  47. Fernandez, A.; Abdeljawad, T.; Baleanu, D. Relations between fractional models with three-parameter Mittag-Leffler kernels. Adv. Differ. Equations 2020, 2020, 186. [Google Scholar] [CrossRef]
  48. Atangana, A.; Baleanu, D. New fractional derivative with non-local and non-singular kernel. Therm. Sci. 2016, 20, 757–763. [Google Scholar] [CrossRef] [Green Version]
  49. Baleanu, D.; Fernandez, A. On some new properties of fractional derivatives with Mittag-Leffler kernel. Commun. Nonlinear Sci. Numer. Simul. 2018, 59, 444–462. [Google Scholar] [CrossRef] [Green Version]
  50. Anwar, A.M.O.; Jarad, F.; Baleanu, D.; Ayaz, F. Fractional Caputo heat equation within the double Laplace transform. Rom. J. Phys. 2013, 58, 15–22. [Google Scholar]
  51. Eridani, A.; Kokilashvili, V.; Meskhi, A. Morrey spaces and fractional integral operators. Expo. Math. 2009, 27, 227–239. [Google Scholar] [CrossRef] [Green Version]
  52. Morales, M.G.; Došla, Z.; Mendoza, F.J. Riemann–Liouville derivative over the space of integrable distributions. Electron. Res. Arch. 2020, 28, 567–587. [Google Scholar] [CrossRef] [Green Version]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Özarslan, M.A.; Fernandez, A. On a Five-Parameter Mittag-Leffler Function and the Corresponding Bivariate Fractional Operators. Fractal Fract. 2021, 5, 45. https://doi.org/10.3390/fractalfract5020045

AMA Style

Özarslan MA, Fernandez A. On a Five-Parameter Mittag-Leffler Function and the Corresponding Bivariate Fractional Operators. Fractal and Fractional. 2021; 5(2):45. https://doi.org/10.3390/fractalfract5020045

Chicago/Turabian Style

Özarslan, Mehmet Ali, and Arran Fernandez. 2021. "On a Five-Parameter Mittag-Leffler Function and the Corresponding Bivariate Fractional Operators" Fractal and Fractional 5, no. 2: 45. https://doi.org/10.3390/fractalfract5020045

APA Style

Özarslan, M. A., & Fernandez, A. (2021). On a Five-Parameter Mittag-Leffler Function and the Corresponding Bivariate Fractional Operators. Fractal and Fractional, 5(2), 45. https://doi.org/10.3390/fractalfract5020045

Article Metrics

Back to TopTop