1. Introduction and Preliminaries
Despite being a field of study for four centuries, special functions has seen tremendous growth in the last fifty years due to its applications in science and engineering. Bessel functions, generalised hypergeometric functions, generalised Mittag-Leffler functions, Wright hypergeometric functions, H-function, Legendre polynomials, Laguerre polynomials, Hermite polynomials, Srivastava polynomials, and many more special functions are included in the list.
In this paper, we shall be working with the well-known Mittag-Leffler function, which was first presented by Mittag-Leffler [
1] in relation to a summation of some divergent series. The Mittag-Leffler function appears in the solution of fractional order differential equations and integral equations, particularly in the study of complex systems, random walks, Levy flights, superdiffusive transport, and the fractional generalisation of the kinetic equation. Numerous works, including Haubold , Mathai and Saxena [
2] and Hilfer [
3] are available for use as references.
The works of Ankit pal, AK Shukla, JC Prajapati, Hari M Srivastava, Arran Fernandez, F Ghanim, Hiba F Al-Janaby, Omar Bazighifan, Anatoly Kilbas, Anna Koroleva, and Sergei Rogosin [
4,
5,
6,
7,
8,
9], and many others contain numerous integral formulas connected to the function of Mittag-Leffler. Inspired by their efforts, we provided unified integrals involving the Mittag-Leffler function as defined in (
1) in this study, and we take into consideration and comment on several intriguing specific instances as a result.
For our purpose, we begin with the Mittag-Leffler function,
, and its two-parameter version,
, which are introduced in Mittag-Leffler and Wiman[
1,
10] in 1903 and 1905 respectively.
where
The Mittag-Leffler functions,
and
, are considered to be a generalization of well-known mathematical functions such as the exponential, hyperbolic and trigonometric functions. This can be seen from the following equalities:
The properties, extensions, and applications of Mittag-Leffler functions
and
have been explored and documented in various works such as Giusti [
11], Gorenflo [
12], Mainardi and Srivastava [
12]. These authors present a range of useful properties with additions to the Mittag-Leffler
function in (
1) and its generalized forms, including numerical calculations.
An increase in the number of parameters is associated with the next phase of development within the Mittag-Leffler function theory. Accordingly, Prabhakar [
13] presented a three-parametric Mittag-Leffler type function,
and continuing with this function, we get
as mentioned by Shukla and Prajapati [
4] in their work, where
and
. It is evident that the aforementioned functions exhibit the following relationships:
For all z, if the function converges, and for , if the function absolutely converges.
The various multi-parameter extensions and generalizations of the Mittag-Leffler function
mentioned in (
1) can all be seen as special cases of the general Fox-Wright function
or its modified version
, which have
p numerator parameters
and
q denominator parameters
such that
Further details can be found in several references, which include work of Erdelyi, Gorenflo, Mainardi, Srivastava, Pal, Khan, Ahmed, Mathai and Haubold [
9,
14,
15,
16].
For Fox-Wright functions,
and
, we referred to Srivastava work [
17] defined as
where,
and the notation
used includes Pochhammer symbols, which are shifted factorials, defined as for
,
The convergence of series (
4) is subject to the the value of
which must be suitably bounded.
The special case of the general Fox-Wright function
is the generalized hypergeometric function
with
p numerator parameters
and
q denominator parameters
,where
, which has been extensively studied and has potential applications. This only happens when
This will provide us with
To set up our results, we will use a result found in Prudnikov, 2018 [
18], which we would like to highlight.
for
and we choose
such that
are not zero where
.
Furthermore, we will look at the Lavoie-Trottier integral formula [
19], which is as follows:
for
.
By using Generalized hypergeometric function representation of
as refered in Shukla and Prajapati [
4]
where
is a
-tuple
is a
-tuple
. Refer to [
4] for more information on the convergence criterion.