A Survey of Some Recent Developments on Higher Transcendental Functions of Analytic Number Theory and Applied Mathematics

Abstract

Often referred to as special functions or mathematical functions, the origin of many members of the remarkably vast family of higher transcendental functions can be traced back to such widespread areas as (for example) mathematical physics, analytic number theory and applied mathematical sciences. Here, in this survey-cum-expository review article, we aim at presenting a brief introductory overview and survey of some of the recent developments in the theory of several extensively studied higher transcendental functions and their potential applications. For further reading and researching by those who are interested in pursuing this subject, we have chosen to provide references to various useful monographs and textbooks on the theory and applications of higher transcendental functions. Some operators of fractional calculus, which are associated with higher transcendental functions, together with their applications, have also been considered. Many of the higher transcendental functions, especially those of the hypergeometric type, which we have investigated in this survey-cum-expository review article, are known to display a kind of symmetry in the sense that they remain invariant when the order of the numerator parameters or when the order of the denominator parameters is arbitrarily changed.

1. Introduction and Motivation

Throughout this survey-cum-expository review article, we use the following standard notations:

$\mathbb{N}:=\{1,2,3,\cdots \},{\mathbb{N}}_0:=\{0,1,2,3,\cdots \}=\mathbb{N}\cup \left\{0\right\}$

and

${\mathbb{Z}}^{-}:=\{-1,-2,-3,\cdots \}={\mathbb{Z}}_0^{-}\setminus \left\{0\right\}.$

Additionally, as usual, $\mathbb{Z}$ denotes the set of integers, $\mathbb{R}$ denotes the set of real numbers and $\mathbb{C}$ denotes the set of complex numbers.

By far the most useful and the most fundamental special function of mathematical analysis and applied mathematics happens to be the familiar (Euler’s) Gamma function $\Gamma \left(z\right)$. Essentially, it stemmed from an attempt by Leonhard Euler to render a meaning to $x!$ when x is any positive real number. In the year 1729, Euler undertook the problem of interpolating $n!$ between the positive integer values of n. It led to what is widely known as the (Euler’s) Gamma function $\Gamma \left(z\right)$ which is defined, for $z\in \mathbb{C}\setminus {\mathbb{Z}}_0^{-}$, by

(1)$\Gamma \left(z\right)=\left\{\begin{array}{cc}{\displaystyle {\int }_0^{\infty }{e}^{-t}{t}^{z-1}dt}\hfill & \left(\mathfrak{R}\left(z\right)>0\right)\hfill \\ \\ {\displaystyle \frac{\Gamma (z+n)}{{\displaystyle \prod _{j=0}^{n-1}}(z+j)}}\hfill & \left(z\in \mathbb{C}\setminus {\mathbb{Z}}_0^{-};n\in \mathbb{N}\right).\hfill \end{array}\right.$

Historically, of course, the origin of the Gamma function $\Gamma \left(z\right)$ defined by (1) can be traced back to two letters from Leonhard Euler (1707–1783) to Christian Goldbach (1690–1764), elaborating upon a simple desire to extend factorials to values between the integers. It is regrettable to see that, in several recent amateurish-type publications, some authors have trivially changed the variable t of integration in the integral definition in (1) and have thereby claimed to have produced a “generalization” of this classical Gamma function $\Gamma \left(z\right)$ (see, for details, [1], Section 3, pp. 1505–1506).

Among the numerous special functions, which are related rather closely to the Gamma function, we choose to mention the incomplete Gamma functions, the Beta and the incomplete Beta functions, the Error functions and the Probability integral, the Digamma and Polygamma functions, and so on.

Our main objective, in this survey-cum-expository review article, is to present a brief introductory overview and survey of some of the recent developments in the theory of a large variety of extensively-studied higher transcendental functions and their potential applications. The detailed plan of this review is as follows. In the next section (Section 2), we present a discussion of the hypergeometric type functions and their multivariate extensions and generalizations. In Section 3, we describe the Polygamma and related functions of analytic number theory. Section 4 is concerned with the Mittag-Leffler type functions and their applications to the modeling and analysis of applied problems by means of fractional calculus are presented in Section 5. Finally, in the concluding section (Section 6), we present a number of further remarks and observations.

We remark in passing that we have made this review article as much complete, comprehensive and self-contained as possible. Of course, wherever applicable, we have chosen to clearly and fully give due credits to the authors of the earlier publications and also included them in the list of references. This has, naturally, resulted in some necessary overlaps with earlier, but clearly and explicitly cited, publications.

2. The Hypergeometric Functions and Their Extensions and Multivariate Generalizations

We begin this section by recalling the following second-order homogeneous linear ordinary differential equation (popularly known as the Gauss hypergeometric equation):

(2)$z(1-z)\frac{d^{2}w}{d{z}^2}+[c-(a+b+1)z]\frac{dw}{dz}-abw=0,$

with three singularities at $z=0,1,\infty$, each of which is a regular singular point of (2). In fact, it happens to be the most celebrated differential equation of the Fuchsian class consisting of differential equations whose only singularities (including the point at ∞) are regular singular points. The importance of the differential Equation (2) arises from the following known observation in the theory of differential equations.

Every homogeneous linear differential equation of the second order, whose singularities (including the point at infinity) are regular and at most three in number, can be transformed into the hypergeometric Equation (2).

The general solution of the hypergeometric Equation (2), which is valid in a neighborhood of the origin ($z=0$), can be expressed in terms of the Gauss hypergeometric function ${}_2{F}_1,$ which is named after the famous German mathematician, Carl Friedrich Gauss (1777–1855), who introduced this function in the year 1812 and gave the F-notation for it. In terms of the general Pochhammer symbol or the shifted factorial ${\left(\lambda \right)}_{\nu }$, since

${\left(1\right)}_n=n!(n\in {\mathbb{N}}_0),$

which is defined (for $\lambda ,\nu \in \mathbb{C})$, in terms of the Gamma function in (1), by

(3)$\begin{array}{c}\hfill {\left(\lambda \right)}_{\nu }:=\frac{\Gamma (\lambda +\nu )}{\Gamma \left(\lambda \right)}=\left\{\begin{array}{cc}1\hfill & (\nu =0;\lambda \in \mathbb{C}\setminus \{0\left\}\right)\hfill \\ \\ \lambda (\lambda +1)\cdots (\lambda +n-1)\hfill & (\nu =n\in \mathbb{N};\lambda \in \mathbb{C}),\hfill \end{array}\right.\end{array}$

it being understood conventionally that ${\left(0\right)}_0:=1$ and assumed tacitly that the $\Gamma$-quotient exists, the Gauss hypergeometric function ${}_2{F}_1$ is defined by

(4)$\begin{array}{cc}\hfill {}_2{F}_1(a,b;c;z)& :=1+\frac{ab}{1·c}z+\frac{a(a+1)b(b+1)}{1·2·c(c+1)}{z}^2+\cdots \hfill \\ \hfill & =\sum _{n=0}^{\infty }\frac{{\left(a\right)}_n{\left(b\right)}_n}{{\left(c\right)}_n}\frac{z^n}{n!}(a,b\in \mathbb{C};c\in \mathbb{C}\setminus {\mathbb{Z}}_0^{-})\hfill \end{array}$

or by its analytic continuation for all $z\in \mathbb{C}\setminus [1,\infty )$. The function ${}_2{F}_1$ does indeed provide the fundamental solution of the hypergeometric differential Equation (2) in the neighborhood of the singular point $z=0$.

In the general case, which is usually credited to Bishop Ernest William Barnes (1874–1953) of the Church of England in Birmingham, who used the notation ${}_p{F}_q(p,q\in {\mathbb{N}}_0)$ analogously in the year 1907 for a generalized hypergeometric function, with p numerator parameters ${\alpha }_j\in \mathbb{C}(j=1,\cdots ,p)$ and q denominator parameters ${\beta }_j\in \mathbb{C}\setminus {\mathbb{Z}}_0^{-}$$(j=1,\cdots ,q)$, given by

(5)$\begin{array}{cc}\hfill {}_p{F}_q\left[\begin{array}{c}\hfill {\alpha }_1,\cdots ,{\alpha }_p;\\ \\ \hfill {\beta }_1,\cdots ,{\beta }_q;\end{array}z\right]& :=\sum _{n=0}^{\infty }\frac{{\left({\alpha }_1\right)}_n\cdots {\left({\alpha }_p\right)}_n}{{\left({\beta }_1\right)}_n\cdots {\left({\beta }_q\right)}_n}\frac{z^n}{n!}\hfill \\ \hfill & =:{}_p{F}_q\left({\alpha }_1,\cdots ,{\alpha }_p;{\beta }_1,\cdots ,{\beta }_q;z\right),\hfill \end{array}$

in which the infinite series

• (i). converges absolutely for $\left|z\right|<\infty$ if $p\leqq q,$

• (ii). converges absolutely for $\left|z\right|<1$ if $p=q+1,$ and

• (iii). diverges for all z$(z\ne 0)$ if $p>q+1.$

Furthermore, if we set

$\omega :=\sum _{j=1}^q{\beta }_j-\sum _{j=1}^p{\alpha }_j,$

then it is known that the generalized hypergeometric ${}_p{F}_q$ series in (5) (with $p=q+1$) is

• I.. absolutely convergent for $\left|z\right|=1$ if $\Re \left(\omega \right)>0,$

• II.. conditionally convergent for $\left|z\right|=1$ $(z\ne 1)$ if $-1<\Re \left(\omega \right)\leqq 0,$ and

• III.. divergent for $\left|z\right|=1$ if $\Re \left(\omega \right)\leqq -1.$

Remarkably, from the definition (5) of the generalized hypergeometric function ${}_p{F}_q,$ it is clearly seen that the sets of the p numerator parameters $a_1,\cdots ,a_p$ and the q denominator parameters $b_1,\cdots ,b_q,$ are symmetric individually. Such potentially useful facts apply equally well to some of the further generalizations and extensions of the generalized hypergeometric function ${}_p{F}_q,$ which is defined above by (5).

In the particular case when $p-1=q=2$, the function ${}_3{F}_2$ is known as the Clausen hypergeometric function in honor of the Danish mathematician and astronomer, Thomas Clausen (1801–1885) who, in the year 1828, established the following hypergeometric identity:

(6)${\left({}_2{F}_1\left[\begin{array}{c}\hfill a,b;\\ \\ \hfill a+b+\frac{1}{2};\end{array}z\right]\right)}^2={}_3{F}_2\left[\begin{array}{c}\hfill 2a,a+b,2b;\\ \\ \hfill a+b+\frac{1}{2},2(a+b);\end{array}z\right].$

For $a,b,c,z\in \mathbb{R}$, the positivity of the function ${}_3{F}_2$ on the right-hand side of Clausen’s identity (6) was instrumental in de Branges’ proof of the 68-year-old Bieberbach conjecture that

$\left|a_n\right|\leqq n(n\in \mathbb{N}\setminus \left\{1\right\})$

for functions f in the normalized analytic and univalent function class $\mathcal{S}$ given by

$\mathcal{S}:=\left\{f:f\in \mathcal{A},f\left(z\right)=z+\sum _{n=0}^{\infty }{a}_n{z}^n(z\in \mathbb{U})\mathrm{and}f\mathrm{is}\mathrm{univalent}\mathrm{in}\mathbb{U}\right\},$

where $\mathcal{A}$ denotes the class of normalized functions which are univalent in the open unit disk $\mathbb{U}$ (see, for details, [2]).

Remarkably, all of such widely used special functions of mathematical physics and applied mathematics as (for example) the Bessel, Legendre and Lommel functions, Whittaker functions, Elliptic integrals, and so on, are expressible in terms of the generalized hypergeometric function ${}_p{F}_q$. Moreover, the Jacobi, Laguerre, Hermite and other associated orthogonal and non-orthogonal polynomials (such as, for example, the Bessel polynomials) are essentially hypergeometric functions ${}_p{F}_q(p,q\in {\mathbb{N}}_0)$ in which one or more of the numerator parameters ${\alpha }_j(j=1,\cdots ,p)$ is a negative integer. In particular, the special orthogonal polynomials ${\left\{{\mathcal{S}}_n\left(x\right)\right\}}_{n\in {\mathbb{N}}_0}$, which are related rather closely to the Bessel polynomials as well as the Laguerre polynomials, occurred in the investigation of energy spectral functions for a certain family of isotropic turbulence fields (see, for details, [3] and the references cited therein). Quite recently in [4], a modified form of these special orthogonal polynomials ${\left\{{\mathcal{S}}_n\left(x\right)\right\}}_{n\in {\mathbb{N}}_0}$ has provided a novel set of (orthogonal) basis functions along with some suitable collocation points in a certain matrix technique for computationally treating a class of multi-order fractional pantograph differential equations by using matrix techniques. Thus, clearly, the potential for the usefulness of the generalized hypergeometric function ${}_p{F}_q(p,q\in {\mathbb{N}}_0)$ and its considerably wide variety of special cases cannot be overemphasized. Some of the books, monographs and tables along these lines include those by Abramowitz and Stegun [5], Andrews [6], Andrews et al. [7], Carlson [8], (Erdélyi et al. [9], Volumes I and II), Rainville [10], Lebedev [11], Luke (see [12,13]), Magnus et al. [14], Miller [15], Srivastava et al. (see [16,17]), Temme [18], Szegö [19], and others.

Various extensions and multivariate generalizations of hypergeometric functions have stemmed from, and are somewhat motivated by, the potential for their applications in solutions of systems of partial differential equations. Some of the multivariate generalizations of hypergeometric functions include the two-variable Appell and Kampé de Fériet functions and the Lauricella functions in several variables, the Srivastava-Daoust hypergeometric functions in two and more variables, and so on. For the theory and applications of many of these univariate and multivariate hypergeometric functions, the interested reader is referred to the monographs by (for example) Appell and Kampé de Fériet [20], Bailey [21], (Erdélyi et al. [9], Vol. I), Slater [22], Seaborn [23], Srivastava et al. (see, for details, [16]), and other authors. In particular, for the theory and applications of the substantially more general G and H functions in one, two and more variables, the reader is referred to the monographs by Mathai et al. [24] and Srivastava et al. (see [17,25]), and others. It should be remarked in passing that a nontrivial and non-vacuous generalization of the H-function was encountered in a systematic study of a class of Feynman integrals, the $\overline{H}$-function was introduced and investigated rather widely and extensively (see, for details, [26,27,28,29,30]).

We choose to conclude this section by recalling the following notable remark which was once made by Gian-Carlo Rota (1932–1999):

The families of univariate and multivariate hypergeometric functions are capable of encompassing just about everything in sight.

3. The Polygamma and Related Functions of Analytic Function Theory

The Polygamma functions ${\psi }^{\left(n\right)}\left(z\right)(n\in {\mathbb{N}}_0)$ are defined by

(7)${\psi }^{\left(n\right)}\left(z\right):=\frac{d^{n+1}}{d{z}^{n+1}}\left\{log\Gamma \left(z\right)\right\}=\frac{d^n}{d{z}^n}\left\{\psi \left(z\right)\right\}(n\in {\mathbb{N}}_0;z\in \mathbb{C}\setminus {\mathbb{Z}}_0^{-}),$

of which ${\psi }^{\left(0\right)}\left(z\right)=\psi \left(z\right)$ is called the Digamma function. In terms of the Hurwitz (or generalized) Zeta function $\zeta (s,a)$, which is defined by

(8)$\zeta (s,a)=\sum _{k=0}^{\infty }\frac{1}{{(k+a)}^s}\left(\Re \left(s\right)>1;a\in \mathbb{C}\setminus {\mathbb{Z}}_0^{-}\right)$

and by its meromorphic continuation to the whole complex s-plane except for a simple pole at $s=1$ with the residue 1, it is easily seen that

(9)${\psi }^{\left(n\right)}\left(z\right)={(-1)}^{n+1}n!\sum _{k=0}^{\infty }\frac{1}{{(k+z)}^{n+1}}={(-1)}^{n+1}n!\zeta (n+1,z)$

$(n\in \mathbb{N};z\in \mathbb{C}\setminus {\mathbb{Z}}_0^{-}),$

which may be used to deduce the properties of the Polygamma functions ${\psi }^{\left(n\right)}\left(z\right)(n\in \mathbb{N})$ from those of the Hurwitz (or generalized) Zeta function $\zeta (s,z)$ $(s=n+1;n\in \mathbb{N}).$ The special case of the Hurwitz (or generalized) Zeta function $\zeta (s,z)$ when $z=1$ is the familiar Riemann Zeta function $\zeta \left(s\right)$, which is involved in the following Taylor-Maclaurin series expansion of $log\Gamma (1+z)$ about the origin $z=0$:

(10)$log\Gamma (1+z)=-\gamma z+\sum _{n=2}^{\infty }{(-1)}^n\zeta \left(n\right)\frac{z^n}{n}\left(\right|z|<1),$

where $\gamma$ denotes the Euler-Mascheroni constant. The Riemann Zeta function $\zeta \left(s\right)$ is defined by

(11)$\zeta \left(s\right):=\left\{\begin{array}{cc}{\displaystyle \sum _{n=1}^{\infty }}{\displaystyle \frac{1}{n^s}}={\displaystyle \frac{1}{1-2^{-s}}}{\displaystyle \sum _{n=1}^{\infty }}{\displaystyle \frac{1}{{\left(2n-1\right)}^s}}\hfill & \left(\Re \left(s\right)>1\right)\hfill \\ \\ {\displaystyle \frac{1}{1-2^{1-s}}}{\displaystyle \sum _{n=1}^{\infty }}{\displaystyle \frac{{\left(-1\right)}^{n-1}}{n^s}}\hfill & \left(\Re \left(s\right)>0;s\ne 1\right)\hfill \end{array}\right.$

and by its meromorphic continuation to the whole complex s-plane except for a simple pole at $s=1$ with the residue 1.

Just as in its obvious special cases $\zeta \left(s\right)$ and $\zeta (s,a)$, the general Hurwitz-Lerch Zeta function $\mathsf{\Phi }(z,s,a)$ defined by

(12)$\mathsf{\Phi }(z,s,a):=\sum _{n=0}^{\infty }\frac{z^n}{{(n+a)}^s}$

$\left(a\in \mathbb{C}\setminus {\mathbb{Z}}_0^{-};s\in \mathbb{C}\mathrm{when}\left|z\right|<1;\Re \left(s\right)>1\mathrm{when}\left|z\right|=1\right),$

can be continued meromorphically to the whole complex s-plane except for a simple pole at $s=1$ with the residue 1. In fact, the general Hurwitz-Lerch Zeta function $\mathsf{\Phi }(z,s,a)$ defined by (12) contains, as its special cases, not only the Riemann Zeta function $\zeta \left(s\right)$ and the Hurwitz (or generalized) Zeta function $\zeta (s,a)$ and the Lerch Zeta function ${\ell }_s\left(\xi \right)$ defined by

(13)${\ell }_s\left(\xi \right):=\sum _{n=1}^{\infty }\frac{e^{2n\pi \mathrm{i}\xi }}{n^s}=e^{2\pi \mathrm{i}\xi }\mathsf{\Phi }\left(e^{2\pi \mathrm{i}\xi },s,1\right)\left(\xi \in \mathbb{R};\Re \left(s\right)>1\right),$

but also such other important functions of Analytic Number Theory as the Polylogarithmic function ${\mathrm{Li}}_s\left(z\right)$:

(14)${\mathrm{Li}}_s\left(z\right):=\sum _{n=1}^{\infty }\frac{z^n}{n^s}=z\mathsf{\Phi }(z,s,1)$

$\left(s\in \mathbb{C}\mathrm{when}\left|z\right|<1;\Re \left(s\right)>1\mathrm{when}\left|z\right|=1\right)$

and the Lipschitz-Lerch Zeta function $\varphi (\xi ,a,s)$:

(15)$\varphi (\xi ,a,s):=\sum _{n=0}^{\infty }\frac{e^{2n\pi \mathrm{i}\xi }}{{(n+a)}^s}=\mathsf{\Phi }\left(e^{2\pi \mathrm{i}\xi },s,a\right)=:L\left(\xi ,s,a\right)$

$\left(a\in \mathbb{C}\setminus {\mathbb{Z}}_0^{-};\Re \left(s\right)>0\mathrm{when}\xi \in \mathbb{R}\setminus \mathbb{Z};\Re \left(s\right)>1\mathrm{when}\xi \in \mathbb{Z}\right),$

which was first studied by Rudolf Lipschitz (1832–1903) and Matyáš Lerch (1860–1922) in connection with Dirichlet’s famous theorem on primes in arithmetic progressions (see, for details, [31]).

For various further extensions of the Hurwitz-Lerch Zeta function $\mathsf{\Phi }(s,z,a)$, including its multi-parameter extensions, the reader is referred to the recent monograph by Srivastava and Choi [32] (see also the other more recent developments reported in [33,34]). In the monograph [32], one can also find a detailed systematic presentation of the double, triple and multiple Gamma functions which were studied by Barnes and others (see also (Whittaker and Watson [35], p. 264) and (Gradshteyn and Ryzhik [36], p. 661, Entry 6.441(4); p. 937, Entry 8.333)).

An interesting and potentially useful family of the $\lambda$-generalized Hurwitz-Lerch Zeta functions, which further extend the multi-parameter Hurwitz-Lerch Zeta function

${\mathsf{\Phi }}_{{\lambda }_1,\cdots ,{\lambda }_p;{\mu }_1,\cdots ,{\mu }_q}^{({\rho }_1,\cdots ,{\rho }_p;{\sigma }_1,\cdots ,{\sigma }_q)}(z,s,\kappa )$

defined by ([37], p. 503, Equation (6.2)) (see also [32])

(16)$\begin{array}{cc}\hfill & {\mathsf{\Phi }}_{{\lambda }_1,\cdots ,{\lambda }_p;{\mu }_1,\cdots ,{\mu }_q}^{({\rho }_1,\cdots ,{\rho }_p,{\sigma }_1,\cdots ,{\sigma }_q)}(z,s,\kappa )\hfill \\ \hfill & :=\sum _{n=0}^{\infty }\frac{{\displaystyle \prod _{j=1}^p}{\left({\lambda }_j\right)}_{n{\rho }_j}}{n!·{\displaystyle \prod _{j=1}^q}{\left({\mu }_j\right)}_{n{\sigma }_j}}\frac{z^n}{{(n+\kappa )}^s}\hfill \end{array}$

$(p,q\in {\mathbb{N}}_0;{\lambda }_j\in \mathbb{C}(j=1,\cdots ,p);\kappa ,{\mu }_j\in \mathbb{C}\setminus Z_0^{-}(j=1,\cdots ,q);$

${\rho }_j,{\sigma }_k\in {\mathbb{R}}^{+}(j=1,\cdots ,p;k=1,\cdots ,q);\Delta >-1\mathrm{when}s,z\in \mathbb{C};$

$\Delta =-1\mathrm{and}s\in \mathbb{C}\mathrm{when}\left|z\right|<{\nabla }^{*};$

$\Delta =-1\mathrm{and}\Re \left(\mathsf{\Xi }\right)>\frac{1}{2}\mathrm{when}\left|z\right|={\nabla }^{*}),$

where, for convenience,

(17)$\Delta :=\sum _{j=1}^q{\sigma }_j-\sum _{j=1}^p{\rho }_j\mathrm{and}\mathsf{\Xi }:=s+\sum _{j=1}^q{\mu }_j-\sum _{j=1}^p{\lambda }_j+\frac{p-q}{2}$

and

(18)${\nabla }^{*}:=\left(\prod _{j=1}^p{\rho }_j^{-{\rho }_j}\right)·\left(\prod _{j=1}^q{\sigma }_j^{{\sigma }_j}\right),$

It was introduced and investigated systematically in a recent paper by Srivastava [38], who also discussed their potential application in Number Theory by appropriately constructing a presumably new continuous analogue of Lippert’s Hurwitz measure and, in addition, considered some other statistical applications of these families of the $\lambda$-generalized Hurwitz-Lerch Zeta functions in probability distribution theory (see also the references to several related earlier works cited by Srivastava [38]). For the convenience of the interested reader in pursuing some of the related open problems, we choose to reproduce here the definition of the $\lambda$-generalized Hurwitz-Lerch Zeta function whose investigation was initiated by Srivastava [38]:

(19)$\begin{array}{cc}\hfill & {\mathsf{\Phi }}_{{\lambda }_1,\cdots ,{\lambda }_p;{\mu }_1,\cdots ,{\mu }_q}^{({\rho }_1,\cdots ,{\rho }_p,{\sigma }_1,\cdots ,{\sigma }_q)}(z,s,a;b,\lambda ):=\frac{1}{\Gamma \left(s\right)}{\int }_0^{\infty }{t}^{s-1}exp\left(-at-\frac{b}{t^{\lambda }}\right)\hfill \\ \hfill & ·{}_p{\mathsf{\Psi }}_q^{*}\left[\begin{array}{c}\hfill ({\lambda }_1,{\rho }_1),\cdots ,({\lambda }_p,{\rho }_p);\\ \\ \hfill ({\mu }_1,{\sigma }_1),\cdots ,({\mu }_q,{\sigma }_q);\end{array}z{e}^{-t}\right]\mathrm{d}t\hfill \end{array}$

$\left(min\{\Re (a),\Re (s\left)\right\}>0;\Re \left(b\right)\geqq 0;\lambda \geqq 0\right),$

so that, obviously, we have the following relationship with the multi-parameter Hurwitz-Lerch Zeta function defined by (16):

(20)$\begin{array}{cc}\hfill {\mathsf{\Phi }}_{{\lambda }_1,\cdots ,{\lambda }_p;{\mu }_1,\cdots ,{\mu }_q}^{({\rho }_1,\cdots ,{\rho }_p,{\sigma }_1,\cdots ,{\sigma }_q)}(z,s,a;0,\lambda )& ={\mathsf{\Phi }}_{{\lambda }_1,\cdots ,{\lambda }_p;{\mu }_1,\cdots ,{\mu }_q}^{({\rho }_1,\cdots ,{\rho }_p,{\sigma }_1,\cdots ,{\sigma }_q)}(z,s,a)\hfill \\ \hfill & =e^b{\mathsf{\Phi }}_{{\lambda }_1,\cdots ,{\lambda }_p;{\mu }_1,\cdots ,{\mu }_q}^{({\rho }_1,\cdots ,{\rho }_p,{\sigma }_1,\cdots ,{\sigma }_q)}(z,s,a;b,0).\hfill \end{array}$

4. The Mittag-Leffler Type Functions: Extensions and Generalizations

We begin this section by recalling the familiar Mittag-Leffler function $E_{\alpha }\left(z\right)$ and its two-parameter version $E_{\alpha ,\beta }\left(z\right)$, which are defined, respectively, by (see [39,40,41])

(21)$E_{\alpha }\left(z\right):=\sum _{k=0}^{\infty }\frac{z^k}{\Gamma (\alpha k+1)}\mathrm{and}{E}_{\alpha ,\beta }\left(z\right):=\sum _{k=0}^{\infty }\frac{z^k}{\Gamma (\alpha k+\beta )}$

$\left(z,\alpha ,\beta \in \mathbb{C};\Re \left(\alpha \right)>0\right).$

The one-parameter function $E_{\alpha }\left(z\right)$ was first considered by Magnus Gustaf (Gösta) Mittag-Leffler (1846–1927) in 1903 and its two-parameter version $E_{\alpha ,\beta }\left(z\right)$ was introduced by Anders Wiman (1865–1959) in 1905 (see also [42]).

The Mittag-Leffler functions $E_{\alpha }\left(z\right)$ and $E_{\alpha ,\beta }\left(z\right)$ are natural extensions of the exponential, hyperbolic and trigonometric functions. In fact, it is easily observed that

$E_1\left(z\right)=e^z,E_2\left(z^2\right)=coshz,E_2\left(-z^2\right)=cosz,$

$E_{1,2}\left(z\right)=\frac{e^z-1}{z}\mathrm{and}{E}_{2,2}\left(z^2\right)=\frac{sinhz}{z}.$

Some of the potentially useful properties, extensions and applications of the Mittag-Leffler functions $E_{\alpha }\left(z\right)$ and $E_{\alpha ,\beta }\left(z\right)$ were presented by Gorenflo et al. [43], Haubold et al. [44] and Kilbas et al. ([45,46] and ([47], Chapter 1)). The Mittag-Leffler function $E_{\alpha }\left(z\right)$ in (21) and its several generalized forms were computed numerically in the recent works in [48,49].

It should be remarked in passing that the families of Mittag-Leffler type functions have the demonstrated potential for application in a wide range of problems in applied as well as engineering sciences. Moreover, their several extended and generalized forms, including those involving multiple indices, are known to be needed in solving fractional differential and integro-differential equations (see, for example, [50]; see also [51,52,53], as well as the references to the related earlier works which are cited therein). For a detailed presentation of some of these recent developments, one is referred to a recent survey-cum-expository review article [54].

Most (if not all) of the above-mentioned multi-parameter generalizations and extensions of the Mittag-Leffler function $E_{\alpha }\left(z\right)$ in (21) are contained, as obvious special cases, of the general Fox-Wright function ${}_p{\mathsf{\Psi }}_q(p,q\in {\mathbb{N}}_0)$ or ${}_p{\mathsf{\Psi }}_q^{*}(p,q\in {\mathbb{N}}_0)$, with p numerator parameters $a_1,\cdots ,a_p$ and q denominator parameters $b_1,\cdots ,b_q$ such that

$a_j\in \mathbb{C}(j=1,\cdots ,p)\mathrm{and}{b}_j\in \mathbb{C}\setminus {\mathbb{Z}}_0^{-}(j=1,\cdots ,q).$

These general Fox-Wright functions ${}_p{\mathsf{\Psi }}_q$ and ${}_p{\mathsf{\Psi }}_q^{*}$ are defined by (see, for details, ([9], Volume I, p. 183) and ([16], p. 21); see also ([47], p. 56), ([55], p. 65) and ([25], p. 19))

(22)$\begin{array}{cc}\hfill & {}_p{\mathsf{\Psi }}_q^{*}\left[\begin{array}{c}\hfill \left(a_1,A_1\right),\cdots ,\left(a_p,A_p\right);\\ \\ \hfill \left(b_1,B_1\right),\cdots ,\left(b_q,B_q\right);\end{array}z\right]\hfill \\ \hfill & :=\sum _{n=0}^{\infty }\frac{{\left(a_1\right)}_{A_{1}n}\cdots {\left(a_p\right)}_{A_{p}n}}{{\left(b_1\right)}_{B_{1}n}\cdots {\left(b_q\right)}_{B_{q}n}}\frac{z^n}{n!}\hfill \\ \hfill & =:\frac{\Gamma \left(b_1\right)\cdots \Gamma \left(b_q\right)}{\Gamma \left(a_1\right)\cdots \Gamma \left(a_p\right)}{}_p{\mathsf{\Psi }}_q\left[\begin{array}{c}\hfill \left(a_1,A_1\right),\cdots ,\left(a_p,A_p\right);\\ \\ \hfill \left(b_1,B_1\right),\cdots ,\left(b_q,B_q\right);\end{array}z\right]\hfill \end{array}$

$\left(\Re \left(A_j\right)>0\left(j=1,\cdots ,p\right);\Re \left(B_j\right)>0\left(j=1,\cdots ,q\right);\Re \left(\sum _{j=1}^q{B}_j-\sum _{j=1}^p{A}_j\right)\geqq -1\right),$

where, and in what follows, ${\left(\lambda \right)}_{\nu }$ denotes the general Pochhammer symbol or the shifted factorial defined already by (3), and the equality in the convergence condition holds true only for suitably bounded values of $\left|z\right|$ given by

$\left|z\right|<\nabla :=\left(\prod _{j=1}^p{A}_j^{-A_j}\right)·\left(\prod _{j=1}^q{B}_j^{B_j}\right).$

Clearly, the generalized hypergeometric function ${}_p{F}_q(p,q\in {\mathbb{N}}_0)$ in (23), with p numerator parameters $a_1,\cdots ,a_p$ and q denominator parameters $b_1,\cdots ,b_q$, is a widely and extensively investigated and potentially useful special case of the general Fox-Wright function ${}_p{\mathsf{\Psi }}_q(p,q\in {\mathbb{N}}_0)$ when

$A_j=1(j=1,\cdots ,p)\mathrm{and}{B}_j=1(j=1,\cdots ,q),$

given by

(23)$\begin{array}{cc}\hfill {}_p{F}_q\left[\begin{array}{c}\hfill a_1,\cdots ,a_p;\\ \\ \hfill b_1,\cdots ,b_q;\end{array}z\right]& :=\sum _{n=0}^{\infty }\frac{{\left(a_1\right)}_n\cdots {\left(a_p\right)}_n}{{\left(b_1\right)}_n\cdots {\left(b_q\right)}_n}\frac{z^n}{n!}\hfill \\ \hfill & ={}_p{\mathsf{\Psi }}_q^{*}\left[\begin{array}{c}\hfill \left(a_1,1\right),\cdots ,\left(a_p,1\right);\\ \\ \hfill \left(b_1,1\right),\cdots ,\left(b_q,1\right);\end{array}z\right]\hfill \\ \hfill & =\frac{\Gamma \left(b_1\right)\cdots \Gamma \left(b_q\right)}{\Gamma \left(a_1\right)\cdots \Gamma \left(a_p\right)}{}_p{\mathsf{\Psi }}_q\left[\begin{array}{c}\hfill \left(a_1,1\right),\cdots ,\left(a_p,1\right);\\ \\ \hfill \left(b_1,1\right),\cdots ,\left(b_q,1\right);\end{array}z\right].\hfill \end{array}$

Just as we have already remarked above, various families of special functions of the Mittag-Leffler and Fox-Wright types are known to play important roles in the theory of fractional calculus and operational calculus as well as in their applications in the basic processes of evolution, relaxation, diffusion, oscillation and wave propagation. Furthermore, the Mittag-Leffler type functions have only recently been calculated numerically in the whole complex plane (see, for example, [48,49]; see also [5,56]). Some other general families of Mittag-Leffler type functions have been investigated and applied recently (see, for details, [54] and the references to earlier works cited therein).

It was my proud privilege to have met many times and discussed mathematical researches, especially on various families of higher transcendental functions and related topics, with my Canadian colleague, Charles Fox (1897–1977) of birth and education in England, both at McGill University and Sir George Williams University (now Concordia University) in Montréal, mainly during the 1970s (see, for details, [57]). Another remarkable mathematical scientist of modern times happens to be Sir Edward Maitland Wright (1906–2005), with whom I had the privilege to meet and discuss research emerging from his publications on hypergeometric and related higher transcendental functions during my visit to the University of Aberdeen in Scotland in the year 1976. We recall here a series of monumental works by Wright (see, for example, [58,59,60]), in which he introduced and systematically studied the asymptotic expansion of the following Taylor-Maclaurin series (see [58], p. 424):

(24)${\mathfrak{E}}_{\alpha ,\beta }(\varphi ;z):=\sum _{n=0}^{\infty }\frac{\varphi \left(n\right)}{\Gamma (\alpha n+\beta )}{z}^n\left(\alpha ,\beta \in \mathbb{C};\Re \left(\alpha \right)>0\right),$

where $\varphi \left(t\right)$ is a function satisfying suitable conditions. Wright’s work was motivated substantially by the earlier developments involving simpler cases, which were reported by Mittag-Leffler in 1905, Wiman in 1905, Ernest William Barnes (1874–1953) in 1906, Godfrey Harold Hardy (1877–1947) in 1905, George Neville Watson (1886–1965) in 1913, Fox in 1928, and many other authors. In particular, the aforementioned work [61] by Bishop Ernest William Barnes (1874–1953) of the Church of England in Birmingham considered the asymptotic expansions of functions in the class defined below:

(25)$E_{\alpha ,\beta }^{\left(\kappa \right)}(s;z):=\sum _{n=0}^{\infty }\frac{z^n}{{(n+\kappa )}^s\Gamma (\alpha n+\beta )}\left(\alpha ,\beta \in \mathbb{C};\Re \left(\alpha \right)>0\right)$

for suitably restricted parameters $\kappa$ and s. It is easy to deduce, from the definition (25), the following relationships with the Mittag-Leffler type function $E_{\alpha ,\beta }^{\left(\kappa \right)}(s;z)$ of Barnes [61]:

(26)$E_{\alpha }\left(z\right)=\underset{s\to 0}{lim}\left\{E_{\alpha ,1}^{\left(\kappa \right)}(s;z)\right\}\mathrm{and}{E}_{\alpha ,\beta }\left(z\right)=\underset{s\to 0}{lim}\left\{E_{\alpha ,\beta }^{\left(\kappa \right)}(s;z)\right\}.$

More interestingly, we also have the following relationship:

$\underset{\alpha \to 0}{lim}\left\{E_{\alpha ,\beta }^{\left(\kappa \right)}(s;z)\right\}=\frac{1}{\Gamma \left(\beta \right)}\mathsf{\Phi }(z,s,\kappa )$

with the classical Lerch transcendent (or the Hurwitz-Lerch Zeta function) $\mathsf{\Phi }(z,s,\kappa )$ defined by (12).

As we have indicated above, Barnes [61] systematically presented asymptotic expansions of many functions such as those in the class of the Mittag-Leffler type function $E_{\alpha ,\beta }^{\left(\kappa \right)}(s;z)$ defined by (25), and the classical Mittag-Leffler functions $E_{\alpha }\left(z\right)$ and $E_{\alpha ,\beta }\left(z\right)$ defined by (21). On the other hand, the multi-parameter Hurwitz-Lerch Zeta function

${\mathsf{\Phi }}_{{\lambda }_1,\cdots ,{\lambda }_p;{\mu }_1,\cdots ,{\mu }_q}^{({\rho }_1,\cdots ,{\rho }_p;{\sigma }_1,\cdots ,{\sigma }_q)}(z,s,\kappa ),$

defined by (16), obviously provides a natural unification and generalization of the Fox-Wright function ${}_p{\mathsf{\Psi }}_q^{*}$ defined by (22) as well as the Hurwitz-Lerch Zeta function $\mathsf{\Phi }(z,s,\kappa )$ defined by (12).

We conclude this section by recalling the following interesting unification of the definitions in (24) and (16) for suitably restricted function $\phi \left(\tau \right)$ (see, for details, [1,54]):

(27)${\mathcal{E}}_{\alpha ,\beta }(\phi ;z,s,\kappa ):=\sum _{n=0}^{\infty }\frac{\phi \left(n\right)}{{(n+\kappa )}^s\Gamma (\alpha n+\beta )}{z}^n\left(\alpha ,\beta \in \mathbb{C};\Re \left(\alpha \right)>0\right),$

where the parameters $\alpha ,$$\beta ,$s and $\kappa$ are appropriately constrained as above. If, upon replacing the sequence ${\left\{\phi \left(n\right)\right\}}_{n=0}^{\infty }$ in the definition (27) by the sequence ${\left\{\varphi \left(n\right)\right\}}_{n=0}^{\infty },$ we readily observe that

(28)${\mathfrak{E}}_{\alpha ,\beta }(\varphi ;z)=\underset{s\to 0}{lim}\left\{{\mathcal{E}}_{\alpha ,\beta }(\phi ;z,s,\kappa )\right\}{|}_{\phi \equiv \varphi }.$

Moreover, if we put $\alpha =\beta =1$ and

(29)$\phi \left(n\right)={\displaystyle \frac{{\displaystyle \prod _{j=1}^p}{\left({\lambda }_j\right)}_{n{\rho }_j}}{{\displaystyle \prod _{j=1}^q}{\left({\mu }_j\right)}_{n{\sigma }_j}}}(n\in {\mathbb{N}}_0)$

in the definition (27), then the definition (27) will immediately yield the definition (16) of the extended Hurwitz-Lerch Zeta function

${\mathsf{\Phi }}_{{\lambda }_1,\cdots ,{\lambda }_p;{\mu }_1,\cdots ,{\mu }_q}^{({\rho }_1,\cdots ,{\rho }_p;{\sigma }_1,\cdots ,{\sigma }_q)}(z,s,\kappa ).$

Alternatively, in the special case of (27) when $\alpha \to 0,$ $\beta =1$ and

(30)$\phi \left(n\right)={\displaystyle \frac{{\displaystyle \prod _{j=1}^p}{\left({\lambda }_j\right)}_{n{\rho }_j}}{n!·{\displaystyle \prod _{j=1}^q}{\left({\mu }_j\right)}_{n{\sigma }_j}}}(n\in {\mathbb{N}}_0)$

or, more simply, by setting

(31)$\phi \left(n\right)={\displaystyle \frac{\Gamma (\alpha n+\beta ){\displaystyle \prod _{j=1}^p}{\left({\lambda }_j\right)}_{n{\rho }_j}}{n!·{\displaystyle \prod _{j=1}^q}{\left({\mu }_j\right)}_{n{\sigma }_j}}}(n\in {\mathbb{N}}_0),$

we are led to the extended Hurwitz-Lerch Zeta function

${\mathsf{\Phi }}_{{\lambda }_1,\cdots ,{\lambda }_p;{\mu }_1,\cdots ,{\mu }_q}^{({\rho }_1,\cdots ,{\rho }_p;{\sigma }_1,\cdots ,{\sigma }_q)}(z,s,\kappa )$

defined by (16).

5. Operators of Fractional Calculus and Their Applications

The current literature is loaded by fractional-order modeling and analysis of real-world and other problems which arise in the mathematical, physical, biological, statistical and engineering sciences. It is truly amazing to see how fast the subject of fractional calculus has grown in recent years, the concept of which is rooted essentially in a question raised in the year 1695 by Marquis de l’Hôpital (1661–1704) to Gottfried Wilhelm Leibniz (1646–1716), which sought the meaning of Leibniz’s (currently popular) notation

$\frac{d^{n}y}{d{x}^n}$

for the derivative of order $n\in {\mathbb{N}}_0$ when $n=\frac{1}{2}$, and indeed also in Leibniz’s reply dated 30 September 1695 to l’Hôpital as follows:

“⋯ This is an apparent paradox from which, one day, useful consequences will be drawn. ⋯”

In widespread applications of fractional calculus, use is made of fractional-order derivatives of different (and, occasionally, ad hoc) kinds (see, for example, [62,63,64,65,66,67,68,69,70,71,72,73,74]). It is fairly traditional to define the fractional-order integrals and fractional-order derivatives by means of the following right-sided Riemann-Liouville fractional integral operator ${}^{\mathrm{RL}}{I}_{a+}^{\mu }$ and the left-sided Riemann-Liouville fractional integral operator ${}^{\mathrm{RL}}{I}_{a-}^{\mu }$, and the corresponding Riemann-Liouville fractional derivative operators ${}^{\mathrm{RL}}{D}_{a+}^{\mu }$ and ${}^{\mathrm{RL}}{D}_{a-}^{\mu }$, as follows (see, for example, ([9], Volume II, Chapter 13), ([47], pp. 69–70) and [75]):

(32)$\left({}^{\mathrm{RL}}{I}_{a+}^{\mu }f\right)\left(x\right)=\frac{1}{\Gamma \left(\mu \right)}{\int }_a^x{\left(x-t\right)}^{\mu -1}f\left(t\right)dt\left(x>a;\Re \left(\mu \right)>0\right),$

(33)$\left({}^{\mathrm{RL}}{I}_{a-}^{\mu }f\right)\left(x\right)=\frac{1}{\Gamma \left(\mu \right)}{\int }_x^a{\left(t-x\right)}^{\mu -1}f\left(t\right)dt\left(x<a;\Re \left(\mu \right)>0\right)$

and

(34)$\left({}^{\mathrm{RL}}{D}_{a\pm }^{\mu }f\right)\left(x\right)={\left(\pm \frac{d}{dx}\right)}^n\left(I_{a\pm }^{n-\mu }f\right)\left(x\right)\left(\Re \left(\mu \right)\geqq 0;n=\left[\Re \left(\mu \right)\right]+1\right),$

where the function f is assumed to be at least locally integrable and $\left[\Re \left(\mu \right)\right]$ indicates the greatest integer in $\Re \left(\mu \right)$. It is worthwhile to remark in passing that, in the current literature, there is an unfortunate trend to trivially and inconsequentially translate known theory and known results based upon the classical Riemann-Liouville and other types of fractional integrals and fractional derivatives in terms of the corresponding theory and the corresponding results by forcing in some obviously redundant (or superfluous) parameters and variables in the widely and extensively investigated definitions and results (see, for details, [1], Section 2, pp. 1504–1505).

In a series of papers, Hilfer et al. (see [62,63,64]; see also [66,67,70]) introduced and studied an interesting family of generalized Riemann-Liouville fractional derivatives of order $\mu (0<\mu <1)$ and type $\nu (0\leqq \nu \leqq 1)$. The right-sided Hilfer fractional derivative operator ${}^{\mathrm{H}}{D}_{a+}^{\mu ,\nu }$ and the left-sided Hilfer fractional derivative operator ${}^{\mathrm{H}}{D}_{a-}^{\alpha ,\beta }$ of order $\mu (0<\mu <1)$ and type $\nu (0\leqq \beta \leqq 1)$ with respect to x are defined by

(35)$\left({}^{\mathrm{H}}{D}_{a\pm }^{\mu ,\nu }f\right)\left(x\right)=\left(\pm {}^{\mathrm{H}}{I}_{a\pm }^{\nu \left(1-\mu \right)}\frac{d}{dx}\left({}^{\mathrm{H}}{I}_{a\pm }^{\left(1-\nu \right)\left(1-\mu \right)}f\right)\right)\left(x\right),$

in which the second member is tacitly assumed to exist.

For $\nu =0$, (35) reduces to the familiar Riemann-Liouville fractional derivative operator. When $\nu =1$, (35) yields the fractional derivative operator that was introduced by (Liouville [76], p. 10) which, quite frequently, is attributed obviously incorrectly to Caputo [77] now-a-days, but which should more appropriately be called the Liouville-Caputo fractional derivative, giving due credit to the originator, Joseph Liouville (1809–1882), who considered such fractional derivatives many decades earlier in 1832 (see [76]). The general operators in (35) are referred to as the Hilfer fractional derivative operators (see, for example, [68]). The Hilfer fractional derivative operator $D_{a\pm }^{\alpha ,\beta }$ was applied in [64] (see also [50,78]).

The Fox-Wright hypergeometric function ${}_p{\mathsf{\Psi }}_q\left(z\right)$ defined by (22), and also such more general functions as Meijer’s G-function and Fox’s H-function, were used as kernels of many different classes of operators of fractional calculus (see, for details, [25,72,79], as well as the references cited therein). It should be noted here that Srivastava et al. [72] used fractional integrals of the Riemann-Liouville type with kernels involving the Fox H-function and the Fox-Wright hypergeometric function ${}_p{\mathsf{\Psi }}_q\left(z\right)$. Moreover, they also considered applications of their results to the general $\overline{H}$-function (see, for example, [29]).

The Wright function ${\mathfrak{E}}_{\alpha ,\beta }(\phi ;z)$ in (24), which was introduced and studied in [58] as long ago as 1940, has appeared recently in [80] in connection with fractional calculus, but without giving due credit to Wright [58]. Closely following the recent works [1,54], the general right-sided fractional integral operator ${\mathcal{I}}_{a+}^{\mu }(\phi ;s,\kappa )$ and the general left-sided fractional integral operator ${\mathcal{I}}_{a-}^{\mu }(\phi ;z,s,\kappa ,\nu )$, and the corresponding fractional derivative operators ${\mathcal{D}}_{a+}^{\mu }(\phi ;z,s,\kappa ,\nu )$ and ${\mathcal{D}}_{a-}^{\mu }(\phi ;z,s,\kappa ,\nu ),$ each of the Riemann-Liouville type, are defined by

(36)$\left({\mathcal{I}}_{a+}^{\mu }(\phi ;z,s,\kappa ,\nu )f\right)\left(x\right)=\frac{1}{\Gamma \left(\mu \right)}{\int }_a^x{\left(x-t\right)}^{\mu -1}{\mathcal{E}}_{\alpha ,\beta }\left(\phi ;z{(x-t)}^{\nu },s,\kappa \right)f\left(t\right)dt$

$\left(x>a;\Re \left(\mu \right)>0\right),$

(37)$\left({\mathcal{I}}_{a-}^{\mu }(\phi ;z,s,\kappa ,\nu )f\right)\left(x\right)=\frac{1}{\Gamma \left(\mu \right)}{\int }_x^a{\left(t-x\right)}^{\mu -1}{\mathcal{E}}_{\alpha ,\beta }\left(\phi ;z{(t-x)}^{\nu },s,\kappa \right)f\left(t\right)dt$

$\left(x<a;\Re \left(\mu \right)>0\right)$

and

(38)$\left({\mathcal{D}}_{a\pm }^{\mu }(\phi ;z,s,\kappa ,\nu )f\right)\left(x\right)={\left(\pm \frac{d}{dx}\right)}^n\left({\mathcal{I}}_{a\pm }^{n-\mu }(\phi ;z,s,\kappa ,\nu )f\right)\left(x\right)$

$\left(\Re \left(\mu \right)\geqq 0;n=\left[\Re \left(\mu \right)\right]+1\right),$

where f is in the space $L(\mathfrak{a},\mathfrak{b})$ of Lebesgue integrable functions on an interval $[\mathfrak{a},\mathfrak{b}]$ $(\mathfrak{b}>\mathfrak{a})$ of $\mathbb{R}$, which is given explicitly by

(39)$L(\mathfrak{a},\mathfrak{b})=\left\{f:{∥f∥}_1={\int }_{\mathfrak{a}}^{\mathfrak{b}}\left|f\left(x\right)\right|dx<\infty \right\}.$

For potential applications based upon the general fractional-calculus operators defined by the Equations (36) to (38), we list below several useful properties of the kernel ${\mathcal{E}}_{\alpha ,\beta }\left(\phi ;z{x}^{\nu },s,\kappa \right)$ involved therein.

(40)$\begin{array}{cc}\hfill & \frac{d^n}{d{x}^n}\left\{x^{\mu -1}{\mathcal{E}}_{\alpha ,\beta }\left(\phi ;z{x}^{\nu },s,\kappa \right)\right\}\hfill \\ \hfill & =x^{\mu -n-1}\sum _{j=0}^{\infty }\frac{\phi \left(j\right)}{{(j+\kappa )}^s\Gamma (\alpha j+\beta )}\frac{\Gamma (\nu j+\mu )}{\Gamma \left(\nu j+\mu -n\right)}{\left(z{x}^{\nu }\right)}^j\hfill \end{array}$

$\left(n\in {\mathbb{N}}_0;\Re \left(\mu \right)>0;\Re \left(\nu \right)>0;\Re \left(\alpha \right)>0\right),$

which, in the special case when $\mu =\beta$ and $\nu =\alpha ,$ yields

(41)$\frac{d^n}{d{x}^n}\left\{x^{\beta -1}{\mathcal{E}}_{\alpha ,\beta }\left(\phi ;z{x}^{\alpha },s,\kappa \right)\right\}=x^{\beta -n-1}{\mathcal{E}}_{\alpha ,\beta -n}\left(\phi ;z{x}^{\alpha },s,\kappa \right)$

$\left(n\in {\mathbb{N}}_0;\Re \left(\alpha \right)>0;\Re \left(\beta \right)>0\right),$

provided that each member of the Equations (40) and (41) exists.

(42)$\begin{array}{cc}\hfill & \mathfrak{I}\left\{t^{\mu -1}{\mathcal{E}}_{\alpha ,\beta }\left(\phi ;z{t}^{\nu },s,\kappa \right)\right\}\left(x\right)\hfill \\ \hfill & :={\int }_0^x{t}^{\mu -1}{\mathcal{E}}_{\alpha ,\beta }\left(\phi ;z{t}^{\nu },s,\kappa \right)dt\hfill \\ \hfill & =x^{\mu }\sum _{j=0}^{\infty }\frac{\phi \left(j\right)}{{(j+\kappa )}^s\Gamma (\alpha j+\beta )}\frac{\Gamma (\nu j+\mu )}{\Gamma \left(\nu j+\mu +1\right)}{\left(z{x}^{\nu }\right)}^j\hfill \end{array}$

$\left(\Re \left(\mu \right)>0;\Re \left(\nu \right)>0;\Re \left(\alpha \right)>0\right),$

provided that the integral exists. More generally, we have

(43)$\begin{array}{cc}\hfill & {\mathfrak{I}}^n\left\{t^{\mu -1}{\mathcal{E}}_{\alpha ,\beta }\left(\phi ;z{t}^{\nu },s,\kappa \right)\right\}\left(x\right)\hfill \\ \hfill & =x^{\mu +n-1}\sum _{j=0}^{\infty }\frac{\phi \left(j\right)}{{(j+\kappa )}^s\Gamma (\alpha j+\beta )}\frac{\Gamma (\nu j+\mu +n-1)}{\Gamma \left(\nu j+\mu +n\right)}{\left(z{x}^{\nu }\right)}^j\hfill \end{array}$

$\left(n\in \mathbb{N};\Re \left(\mu \right)>0;\Re \left(\nu \right)>0;\Re \left(\alpha \right)>0\right),$

which, in the special case when $\mu =\beta$ and $\nu =\alpha ,$ yields

(44)${\mathfrak{I}}^n\left\{t^{\beta -1}{\mathcal{E}}_{\alpha ,\beta }\left(\phi ;z{t}^{\alpha },s,\kappa \right)\right\}\left(x\right)=x^{\beta +n-1}{\mathcal{E}}_{\alpha ,\beta +n}\left(\phi ;z{x}^{\alpha },s,\kappa \right)$

$\left(n\in \mathbb{N};\Re \left(\alpha \right)>0;\Re \left(\beta \right)>0\right),$

provided that each member of the Equations (43) and (44) exists.

For the operator $\mathcal{L}$ of the Laplace transform given by

(45)$\mathcal{L}\left\{f\left(\tau \right):\mathfrak{s}\right\}:={\int }_0^{\infty }{\mathrm{e}}^{-\mathfrak{s}\tau }f\left(\tau \right)d\tau =:F\left(\mathfrak{s}\right)\left(\Re \left(\mathfrak{s}\right)>0\right),$

where the function $f\left(\tau \right)$ is so constrained that the integral exists, we get

(46)$\begin{array}{c}\hfill \mathcal{L}\left\{{\tau }^{\mu -1}{\mathcal{E}}_{\alpha ,\beta }\left(\phi ;z{\tau }^{\nu },s,\kappa \right):\mathfrak{s}\right\}=\frac{1}{{\mathfrak{s}}^{\mu }}\sum _{j=0}^{\infty }\frac{\phi \left(j\right)\Gamma (\nu j+\mu )}{{(j+\kappa )}^s\Gamma (\alpha j+\beta )}{\left(\frac{z}{{\mathfrak{s}}^{\nu }}\right)}^j\end{array}$

$\left(\Re \left(\mathfrak{s}\right)>0;\Re \left(\mu \right)>0;\Re \left(\nu \right)>0;\Re \left(\alpha \right)>0\right),$

provided that each member of (46) exists. For $\mu =\beta$ and $\nu =\alpha ,$ the Laplace transform Formula (46) assumes the following simpler form:

(47)$\begin{array}{c}\hfill \mathcal{L}\left\{{\tau }^{\beta -1}{\mathcal{E}}_{\alpha ,\beta }\left(\phi ;z{\tau }^{\alpha },s,\kappa \right):\mathfrak{s}\right\}=\frac{1}{{\mathfrak{s}}^{\mu }}\sum _{j=0}^{\infty }\frac{\phi \left(j\right)}{{(j+\kappa )}^s}{\left(\frac{z}{{\mathfrak{s}}^{\alpha }}\right)}^j\end{array}$

$\left(\Re \left(\mathfrak{s}\right)>0;\Re \left(\alpha \right)>0;\Re \left(\beta \right)>0\right).$

Unfortunately, by trivially changing, in the definition of the classical Laplace transform (45), the index $\mathfrak{s}$ or the integration variable t or both the index $\mathfrak{s}$ and the integration variable t, many mainly amateurish-type researchers have made and continue to make the obviously false claim to have “generalized” the classical Laplace transform (45) itself. Some of the examples in this connection include, but are not limited to, the so-called Sumudu transform, the so-called natural transform, the so-called Shehu transform, the so-called ${\mathcal{P}}_{\delta }$-transform, the so-called k-Laplace transform, and so on (see, for details, [1], Section 4, pp. 1508–1510).

Various special cases and consequences of the above key results, involving simpler functions of the types which we have considered in this and the preceding sections, can be deduced fairly easily. We turn instead to the following examples of applications involving fractional-order derivatives of one kind or the other. In this connection, we note the following relationship between the Riemann-Liouville fractional derivative ${}^{\mathrm{RL}}{D}_{0+}^{\mu }$ and the Liouville-Caputo fractional derivative ${}^{\mathrm{LC}}{D}_{0+}^{\mu }$ of order $\mu$ (see, for example, ([47], p. 91, Equation (2.4.1)) with $a=0$):

(48)$\left({}^{\mathrm{LC}}{D}_{0+}^{\mu }f\right)\left(x\right)=\left({}^{\mathrm{RL}}{D}_{0+}^{\mu }\left\{f\left(t\right)-\sum _{k=0}^{n-1}\frac{f^{\left(k\right)}\left(0\right)}{k!}{t}^k\right\}\right)\left(x\right),$

where n is given by

(49)$n=\left\{\begin{array}{cc}[\Re (\mu \left)\right]+1\hfill & (\mu \ne {\mathbb{N}}_0)\hfill \\ \\ \mu \hfill & (\mu \in {\mathbb{N}}_0).\hfill \end{array}\right.$

Equivalently, since

(50)$\left({}^{\mathrm{RL}}{D}_{0+}^{\mu }\left\{t^{\lambda -1}\right\}\right)\left(x\right)=\frac{\Gamma \left(\lambda \right)}{\Gamma (\lambda -\mu )}{x}^{\lambda -\mu -1}\left(\Re \left(\lambda \right)>0;\Re \left(\mu \right)\geqq 0\right),$

the relationship (48) can be written as follows:

(51)$\left({}^{\mathrm{LC}}{D}_{0+}^{\mu }f\right)\left(x\right)=\left({}^{\mathrm{RL}}{D}_{0+}^{\mu }f\right)\left(x\right)-\sum _{k=0}^{n-1}\frac{f^{\left(k\right)}\left(0\right)}{\Gamma (k-\mu +1)}{x}^{k-\mu },$

where n is given, as before in (48), by (49).

I. The basic processes of relaxation, diffusion, oscillations and wave propagation were revisited by Gorenflo et al. [43] who introduced the Liouville-Caputo type fractional-order derivatives in the governing (ordinary or partial) differential equations and considered each of the following fractional differential equations:

(52)$\frac{d^{\alpha }u}{d{t}^{\alpha }}+c^{\alpha }u\left(t;\alpha \right)=0\left(c>0;0<\alpha \leqq 2\right)$

and

(53)$\frac{{\partial }^{2\beta }u}{\partial t^{2\beta }}=k\frac{{\partial }^{2}u}{\partial x^2}\left(k>0;-\infty <x<\infty ;0<\beta \leqq 1\right),$

where the aforementioned Liuoville-Caputo fractional derivative of order $\mu >0$ of a causal function $f\left(t\right)$, that is,

$f\left(t\right)=0(t<0),$

is given by

(54)$\begin{array}{cc}\hfill \frac{d^{\mu }}{d{x}^{\mu }}\left\{f\left(x\right)\right\}& =\left({}^{\mathrm{LC}}{D}_{0+}^{\mu }f\right)\left(x\right)\hfill \\ \hfill & :=\left\{\begin{array}{cc}{f}^{\left(n\right)}\left(x\right)\hfill & (\mu =n\in {\mathbb{N}}_0)\hfill \\ \\ {\displaystyle {\displaystyle \frac{1}{\Gamma (n-\mu )}}{\int }_0^x{\displaystyle \frac{f^{\left(n\right)}\left(t\right)}{{\left(x-t\right)}^{\mu -n+1}}}dt}\hfill & \left(n-1<\Re \left(\mu \right)<n;n\in \mathbb{N}\right).\hfill \end{array}\right.\hfill \end{array}$

Here, as usual, n is given by (49) and $f^{\left(n\right)}\left(t\right)$ denotes the ordinary derivative of $f\left(t\right)$ of order n.

The Equation (52) represents fractional relaxation when $0<\alpha \leqq 1$ under the initial condition $u\left(0+;\alpha \right)=u_0$. Moreover, it can be viewed as fractional oscillation when $1<\alpha \leqq 2$ under the following initial conditions:

$u\left(0+;\alpha \right)=u_0\mathrm{and}\dot{u}\left(0+;\alpha \right)=v_0,$

where $v_0\equiv 0$ for continuous dependence of the solution on the parameter $\alpha$ also in the transition from $\alpha =1-$ to $\alpha =1+$, and $\dot{u}$ is the time-derivative of u.

In the Equation (53), $u=u\left(x,t;\beta \right)$ is assumed to be a causal function of time $\left(t>0\right)$ such that

$u\left(\mp \infty ,t;\beta \right)=0.$

Clearly, the Equation (53) represents fractional diffusion when $0<\beta \leqq \frac{1}{2}$ under the initial condition $u\left(x,0+;\beta \right)=f\left(x\right)$. It can also be viewed as fractional wave equation when $\frac{1}{2}<\beta \leqq 1$ under the following initial conditions:

$u\left(x,0+;\beta \right)=f\left(x\right)\mathrm{and}\dot{u}\left(x,0+;\beta \right)=g\left(x\right),$

where $g\left(x\right)\equiv 0$ for continuous dependence of the solution on the parameter $\beta$ also in the transition from $\beta =\frac{1}{2}-$ to $\beta =\frac{1}{2}+$.

In terms of the Mittag-Leffler function $E_{\alpha }\left(z\right)$ defined by (21), the explicit solution of the initial-value problem involving the fractional differential Equation (52) is given by

(55)$\begin{array}{c}\hfill u\left(t;\alpha \right)=u_0{E}_{\alpha }\left(-{\left(ct\right)}^{\alpha }\right).\end{array}$

On the other hand, the explicit solution of the initial-value problem involving the fractional differential Equation (53) can be expressed as follows:

(56)$u\left(x,t;\beta \right)={\int }_{-\infty }^{\infty }{\mathcal{G}}_c\left(\xi ,t;\beta \right)f\left(x-\xi \right)d\xi ,$

where ${\mathcal{G}}_c\left(x,t;\beta \right)$ denotes the Green function given by

(57)$\left|x\right|{\mathcal{G}}_c\left(x,t;\beta \right)=\frac{z}{2}\sum _{n=0}^{\infty }\frac{{\left(-z\right)}^n}{n!\Gamma \left(1-\beta -\beta n\right)}\left(z=\frac{\left|x\right|}{\sqrt{k}{t}^{\beta }};0<\beta <1\right),$

which, in turn, can easily be expressed in terms of Wright’s generalized Bessel function or the Bessel-Wright function $J_{\nu }^{\mu }\left(z\right)$ defined by (see, for example, [17], p. 42, Equation II.5 (22))

(58) $\begin{array}{c}\hfill J_{\nu }^{\mu }\left(z\right):=\sum _{n=0}^{\infty }\frac{{\left(-z\right)}^n}{n!\Gamma \left(\mu n+\nu +1\right)}={}_0\mathsf{\Psi }{}_2\left[\begin{array}{c}\hfill \overline{};\\ \\ \hfill (1,1),(\nu +1,\mu );\end{array}-z\right].\end{array}$

II. In recent years, various different forms of kinetic equations of fractional order have been widely used, especially in the modeling and analysis of a number of important problems of physics and astrophysics (see, for details, [81]). Particularly, in the past decade or so, kinetic equations of fractional order have apparently gained popularity, mainly because of the discovery of their relation with the theory of CTRW (Continuous-Time Random Walks) (see [66]). These equations have been and are being investigated with the aim to first determine and then interpret certain physical phenomena which are known to govern such processes as (for example) diffusion in porous media, reaction and relaxation in complex systems, anomalous diffusion, and so on (see also [63,82]).

Theorems 1 to 3 below, each of which was established in [54], are sufficiently general key results, which are capable of being appropriately and suitably specialized with a view to including solutions of the corresponding (known or new) fractional-order kinetic equations associated with a large variety of simpler functions than those involved herein. With a view to making our presentation accessible readily to the interested reader, we have chosen to essentially reproduce each of these fundamental results (Theorems 1 to 3 below) from our earlier work [54].

Theorem 1.

Let $c,\mu ,\nu ,\rho ,\sigma \in {\mathbb{R}}^{+}$. Suppose also that the general function ${\mathcal{E}}_{\alpha ,\beta }(\phi ;z,s,\kappa ),$ defined by $\left(27\right),$ exists. Then the solution of the following generalized fractional kinetic equation:

(59) $N\left(t\right)-N_0{t}^{\mu -1}{\mathcal{E}}_{\alpha ,\beta }\left(\phi ;z{t}^{\nu },s,\kappa \right)=-c^{\rho }\left({}^{\mathrm{RL}}{I}_{0+}^{\sigma }N\right)\left(t\right)$

is given by

(60) $\begin{array}{cc}\hfill N\left(t\right)& =N_0{t}^{\mu -1}\sum _{r=0}^{\infty }{\left(-c^{\rho }{t}^{\sigma }\right)}^r\hfill \\ \hfill & ·\sum _{j=0}^{\infty }\frac{\phi \left(j\right)\Gamma (\nu j+\mu )}{{(j+\kappa )}^s\Gamma (\alpha j+\beta )\Gamma (\nu j+\sigma r+\mu )}{\left(z{t}^{\nu }\right)}^j(t>0),\hfill \end{array}$

provided that the right-hand side of the solution asserted by $\left(60\right)$ exists.

Theorem 2.

Let $c,\mu ,\nu ,\rho ,\sigma \in {\mathbb{R}}^{+}$. Suppose also that the general function ${\mathfrak{E}}_{\alpha ,\beta }(\varphi ;z),$ defined by $\left(24\right),$ exists. Then the solution of the following generalized fractional kinetic equation:

(61) $N\left(t\right)-N_0{t}^{\mu -1}{\mathfrak{E}}_{\alpha ,\beta }\left(\varphi ;z{t}^{\nu }\right)=-c^{\rho }\left({}^{\mathrm{RL}}{I}_{0+}^{\sigma }N\right)\left(t\right)$

is given by

(62) $\begin{array}{cc}\hfill N\left(t\right)& =N_0{t}^{\mu -1}\sum _{r=0}^{\infty }{\left(-c^{\rho }{t}^{\sigma }\right)}^r\hfill \\ \hfill & ·\sum _{j=0}^{\infty }\frac{\varphi \left(j\right)\Gamma (\nu j+\mu )}{\Gamma (\alpha j+\beta )\Gamma (\nu j+\sigma r+\mu )}{\left(z{t}^{\nu }\right)}^j(t>0),\hfill \end{array}$

provided that the right-hand side of the solution asserted by $\left(62\right)$ exists.

Theorem 3.

For $c,\mu ,\nu ,\rho ,\sigma \in {\mathbb{R}}^{+},$ let the extended Hurwitz-Lerch Zeta function:

${\mathsf{\Phi }}_{{\lambda }_1,\cdots ,{\lambda }_p;{\mu }_1,\cdots ,{\mu }_q}^{({\rho }_1,\cdots ,{\rho }_p;{\sigma }_1,\cdots ,{\sigma }_q)}\left(z,s,\kappa \right),$

defined by $\left(16\right),$ exist. Then the solution of the following generalized fractional kinetic equation:

(63) $N\left(t\right)-N_0{t}^{\mu -1}{\mathsf{\Phi }}_{{\lambda }_1,\cdots ,{\lambda }_p;{\mu }_1,\cdots ,{\mu }_q}^{({\rho }_1,\cdots ,{\rho }_p;{\sigma }_1,\cdots ,{\sigma }_q)}\left(z{t}^{\nu },s,\kappa \right)=-c^{\rho }\left({}^{\mathrm{RL}}{I}_{0+}^{\sigma }N\right)\left(t\right)$

is given by

(64) $\begin{array}{cc}\hfill N\left(t\right)& =N_0{t}^{\mu -1}\sum _{r=0}^{\infty }{\left(-c^{\rho }{t}^{\sigma }\right)}^r\frac{\Gamma \left(\mu \right)}{\Gamma (\sigma r+\mu )}\hfill \\ \hfill & ·{\mathsf{\Phi }}_{\mu ,{\lambda }_1,\cdots ,{\lambda }_p;\sigma r+\mu ,{\mu }_1,\cdots ,{\mu }_q}^{(\nu ,{\rho }_1,\cdots ,{\rho }_p;\nu ,{\sigma }_1,\cdots ,{\sigma }_q)}\left(z{t}^{\nu },s,\kappa \right)(t>0),\hfill \end{array}$

provided that the right-hand side of the solution asserted by $\left(64\right)$ exists.

Remarkably, the general function ${\mathcal{E}}_{\alpha ,\beta }(\phi ;z,s,\kappa ),$ defined by (27), occurring in the non-homogeneous term of the kinetic Equation (59) of fractional order is distinctly advantageous because of its generality. Naturally, therefore, solutions of other kinetic equations involving other simpler non-homogeneous terms can be derived as appropriate special cases of the solution (60) given by Theorem 1. Similar remarks would apply equally strongly to the results (62) and (64), which are provided by Theorems 2 and 3, respectively.

III. The unified fractional derivative operator ${}^{\mathrm{H}}{D}_{0+}^{\alpha ,\beta }$ of order $\alpha (0<\alpha <1)$ and type $\beta (0\leqq \beta \leqq 1)$, defined by the Equation (35), was considered by Hilfer (see [63]) in order to derive the solution of the following general fractional differential equation:

(65)$\left({}^{\mathrm{H}}{D}_{0+}^{\alpha ,\beta }f\right)\left(x\right)=\lambda f\left(x\right)(x>0)$

under the initial condition given, in terms of the corresponding two-parameter fractional integral operator ${}^{\mathrm{H}}{I}_{0+}^{\alpha ,\beta }$, by

(66)$\left({}^{\mathrm{H}}{I}_{0+}^{\left(1-\beta \right)\left(1-\alpha \right)}f\right)(0+)=c_0,$

where it is tacitly assumed that

$\left({}^{\mathrm{H}}{I}_{0+}^{\left(1-\beta \right)\left(1-\alpha \right)}f\right)(0+):=\underset{x\to 0+}{lim}\left\{\left({}^{\mathrm{H}}{I}_{0+}^{\left(1-\beta \right)\left(1-\alpha \right)}f\right)\left(x\right)\right\},$

$c_0$ is a given constant and the parameter $\lambda$ is the eigenvalue. Hilfer’s solution of the above initial-value problem is given by (see [63], p. 115, Equation (124)):

(67)$f\left(x\right)=c_0{x}^{\left(1-\beta \right)\left(\alpha -1\right)}{E}_{\alpha ,\alpha +\beta \left(1-\alpha \right)}\left(\lambda x^{\alpha }\right),$

where $E_{\alpha ,\beta }\left(z\right)$ denotes the two-parameter Mittag-Leffler function defined by (21).

If we put $\beta =0$ and $c_0=1$ in Hilfer’s solution (66), we can deduce the corrected version of the claimed solution (see [74], p. 802, Equation (3.1); see also [78]) of the following initial-value problem:

(68)$\left({}^{\mathrm{RL}}{D}_{0+}^{\alpha }f\right)\left(x\right)=\lambda f\left(x\right)(x>0),$

under the initial condition given by

(69)$\left({}^{\mathrm{RL}}{I}_{0+}^{1-\alpha }f\right)(0+)=1,$

where, as also in the equation (66),

$\left({}^{\mathrm{RL}}{I}_{0+}^{1-\alpha }f\right)(0+):=\underset{x\to 0+}{lim}\left\{\left({}^{\mathrm{RL}}{I}_{0+}^{1-\alpha }f\right)\left(x\right)\right\}$

in the form given by

(70)$f\left(x\right)=x^{\alpha -1}{E}_{\alpha ,\alpha }\left(\lambda x^{\alpha }\right),$

in terms of the two-parameter Mittag-Leffler function defined by (21).

In concluding this section, we remark that, in recent years, various real-world problems and issues in many areas of mathematical, physical and engineering sciences have been modeled and analyzed by appealing to several powerful tools, one of which involves applications of the operators of fractional calculus. Notably, a number of important and potentially useful definitions have been introduced and used for fractional-order derivatives. These include, among others, the fractional-order derivative operators which stem from the Riemann-Liouville, the Grünwald-Letnikov, the Liouville-Caputo, the Caputo-Fabrizio and the Atangana-Baleanu fractional-order derivatives (see, for example, [47,83,84,85]).

The Riemann-Liouville fractional derivative is known to involve the convolution of a given function and a power-law kernel (see, for details, [47,85]). On the other hand, the Liouville-Caputo (LC) fractional derivative involves the convolution of the local derivative of a given function with a power-law function [86]. The fractional-order derivatives proposed by Caputo and Fabrizio [84] and Atangana and Baleanu [83] are based upon the exponential decay law which is a generalized power-law function (see [87,88,89,90,91,92]). The Caputo-Fabrizio (CFC) fractional-order derivative as well as the Atangana-Baleanu (ABC) fractional-order derivative allow us to describe complex physical problems that follow, at the same time, the power law and the exponential decay law (see, for details, [87,88,89,90,91,92]). Such items of information as those presented in many of these works are believed to have the potential to generate further developments on fractional-order modeling and analysis of interesting applied problems. Many experiments and theories have shown that a fairly large number of abnormal phenomena, which occur in the engineering and applied sciences, can be described well by using discrete fractional calculus (see, for example, [93]). In particular, fractional difference equations have been found to provide powerful tools in the modeling and analysis of various phenomena in many different fields of science and engineering including those in mathematical physics, fluid mechanics, heat conduction, and so on.

6. Concluding Remarks and Observations

Commonly known as special functions or mathematical functions, the origin of many members of the remarkably vast family of higher transcendental functions can be traced back to such widespread areas as (for example) mathematical physics, analytic number theory and applied mathematical sciences. Here, in this survey-cum-expository review article, our objective has been to briefly present an introductory overview and survey of some important recent developments in the theory of several extensively studied families of higher transcendental functions (or, more popularly, special functions) and their potential applications in (for example) mathematical physics, analytic number theory and applied mathematical sciences. For further reading and researching by those who are interested in pursuing this subject, we have chosen to provide references to various useful monographs and textbooks on the theory and applications of higher transcendental functions. We have also considered several operators of fractional calculus, which are associated with higher transcendental functions, briefly indicating their applications as well.

The bibliography in this review article includes a number of recently published journal articles which have dealt with the extensively investigated subject of fractional calculus and its widespread applications. In fact, having an overview of the on-going contributions to the theory and applications of fractional calculus, which are continually appearing in some of the leading journals devoted to mathematical and physical sciences, biological sciences, statistical sciences, engineering sciences, and so on, the subject-matter, which we have dealt with in this review article, is remarkably important and potentially useful. Moreover, the interested future researchers will surely benefit from the listing of references to some of the other applications of various fractional-calculus operators in the mathematical and other sciences, which we have not considered in the preceding sections (see, for example, [94,95,96,97,98,99,100,101,102,103,104,105,106,107,108,109,110,111]).

There is considerable literature investigating and applying the quantum or basic (or $\mathfrak{q}$-) calculus not only in the area of higher transcendental functions, which we have presented in Section 2, and Geometric Function Theory of Complex Analysis (see, for a detailed historical and introductory overview, the recently published survey-cum-expository review article [102]), but also in the modeling and analysis of applied problems as well as in extending the well-established theory and applications of various rather classical mathematical functions and mathematical inequalities. It is regretful, however, to see that a large number of mostly amateurish-type researchers on these and other related topics continue to produce and publish obvious and inconsequential variations and straightforward translations of the known $\mathfrak{q}$-results in terms of the so-called $(\mathfrak{p},\mathfrak{q})$-calculus by unnecessarily forcing-in an obviously superfluous (or redundant) parameter $\mathfrak{p}$ into the classical $\mathfrak{q}$-calculus and thereby falsely claiming “generalization” (see [102], p. 340 and [1], Section 5, pp. 1511–1512). Such tendencies to produce and flood the literature with trivialities should be discouraged by all means.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

It is my pleasure to express my appreciation and sincere thanks to the colleagues in the Symmetry Editorial Office (especially to the Managing Editor, Celina Ci) for kindly inviting me to contribute this survey-cum-expository review article to *Symmetry*.

Conflicts of Interest

The author declares no conflict of interest.

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

References

1. Srivastava, H.M. Some parametric and argument variations of the operators of fractional calculus and related special functions and integral transformations. J. Nonlinear Convex Anal.; 2021; 22, pp. 1501-1520.

2. De Branges, L. A proof of the Bieberbach conjecture. Acta Math.; 1985; 154, pp. 137-152. [DOI: https://dx.doi.org/10.1007/BF02392821]

3. Srivastava, H.M. Some orthogonal polynomials representing the energy spectral functions for a family of isotropic turbulence fields. Z. Angew. Math. Mech.; 1984; 64, pp. 255-257. [DOI: https://dx.doi.org/10.1002/zamm.19840640612]

4. Izadi, M.; Srivastava, H.M. A novel matrix technique for multi-order pantograph differential equations of fractional order. Proc. R. Soc. Lond. Ser. A Math. Phys. Engergy Sci.; 2021; 477, 2021031. [DOI: https://dx.doi.org/10.1098/rspa.2021.0321]

5. Abramowitz, M.; Stegun, I.A. Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables; Applied Mathematics Series; National Bureau of Standards: Washington, DC, USA, 1964; Reprinted by Dover Publications: New York, NY, USA, 1965

6. Andrews, L.C. Special Functions for Engineers and Applied Mathematicians; Macmillan Company: New York, NY, USA, 1984.

7. Andrews, G.E.; Askey, R.; Roy, R. Special Functions; Encyclopedia of Mathematics and Its Applications Cambridge University Press: Cambridge, UK, London, UK, New York, NY, USA, 1999; Volume 71.

8. Carlson, B.C. Special Functions of Applied Mathematics; Academic Press: New York, NY, USA, San Francisco, CA, USA, London, UK, 1977.

9. Erdélyi, A.; Magnus, W.; Oberhettinger, F.; Tricomi, F.G. Higher Transcendental Functions; McGraw-Hill Book Company: New York, NY, USA, Toronto, ON, Canada, 1953; London, UK, 1955; Volumes I, II and III

10. Rainville, E.D. Special Functions; Macmillan Company: New York, NY, USA, 1960; Reprinted by Chelsea Publishing Company: Bronx, NY, USA, 1971

11. Lebedev, N.N. Special Functions and Their Applications; Translated from the Russian by R. A. Silverman Prentice-Hall: Englewood Cliffs, NJ, USA, 1965.

12. Luke, Y.L. The Special Functions and Their Approximations; Mathematics in Science and Engineering A Series of Monographs and Textbooks; Academic Press: New York, NY, USA, London, UK, 1969; Volumes I (53-I) and II (53-II)

13. Luke, Y.L. Mathematical Functions and Their Approximations; Academic Press: New York, NY, USA, San Francisco, CA, USA, London, UK, 1975.

14. Magnus, W.; Oberhettinger, F.; Soni, R.P. Formulas and Theorems for the Special Functions of Mathematical Physics; 3rd ed. Die Grundlehren der Mathematischen Wissenschaften in Einzeldarstellungen mit besonderer Berücksichtingung der Anwendungsgebiete, Bd Springer: Berlin/Heidelberg, Germany, New York, NY, USA, 1966; Volume 52.

15. Miller, W., Jr. Lie Theory and Special Functions; Mathematics in Science and Engineering A Series of Monographs and Textbooks; Academic Press: London, UK, New York, NY, USA, 1968; Volume 43.

16. Srivastava, H.M.; Karlsson, P.W. Multiple Gaussian Hypergeometric Series; Halsted Press: Chichester, UK, John Wiley and Sons: New York, NY, USA, Chichester, UK, Brisbane, Australia, Toronto, ON, Canada, 1985.

17. Srivastava, H.M.; Kashyap, B.R.K. Special Functions in Queuing Theory and Related Stochastic Processes; Academic Press: New York, NY, USA, London, UK, 1982.

18. Temme, N.M. Special Functions: An Introduction to Classical Functions of Mathematical Physics; A Wiley-Interscience Publication, John Wiley and Sons: New York, NY, USA, Chichester, UK, Brisbane, Australia, Toronto, ON, Canada, 1996.

19. Szegö, G. Orthogonal Polynomials; 4th ed. American Mathematical Society: Providence, RI, USA, 1975; Volume 23.

20. Appell, P.; Kampé de Fériet, J. Fonctions Hypergéométriques et Hypersphériques; Polynômes d’Hermite Gauthier-Villars: Paris, France, 1926.

21. Bailey, W.N. Generalized Hypergeometric Series; Cambridge Tracts in Mathematics and Mathematical Physics Cambridge University Press: Cambridge, UK, London, UK, New York, NY, USA, 1935; Volume 32, Reprinted by Stechert-Hafner Service Agency: New York, NY, USA; London, UK, 1964

22. Slater, L.J. Generalized Hypergeometric Functions; Cambridge University Press: Cambridge, UK, London, UK, New York, NY, USA, 1966.

23. Seaborn, J.B. Hypergeometric Functions and Their Applications; Springer: Berlin/Heidelberg, Germany, New York, NY, USA, 1991.

24. Mathai, A.M.; Saxena, R.K.; Haubold, H.J. The H-Function: Theory and Applications; Springer: Berlin/Heidelberg, Germany, New York, NY, USA, 2010.

25. Srivastava, H.M.; Gupta, K.C.; Goyal, S.P. The H-Functions of One and Two Variables with Applications; South Asian Publishers: New Delhi, India, Madras, India, 1982.

26. Inayat-Hussain, A.A. New properties of hypergeometric series derivable from Feynman integrals. I: Transformation and reduction formulae. J. Phys. A Math. Gen.; 1987; 20, pp. 4109-4118. [DOI: https://dx.doi.org/10.1088/0305-4470/20/13/019]

27. Inayat-Hussain, A.A. New properties of hypergeometric series derivable from Feynman integrals. II: A generalisation of the H-function. J. Phys. A Math. Gen.; 1987; 20, pp. 4119-4128. [DOI: https://dx.doi.org/10.1088/0305-4470/20/13/020]

28. Buschman, R.G.; Srivastava, H.M. The H function associated with a certain class of Feynman integrals. J. Phys. A Math. Gen.; 1990; 23, pp. 4707-4710. [DOI: https://dx.doi.org/10.1088/0305-4470/23/20/030]

29. Srivastava, H.M.; Lin, S.-D.; Wang, P.-Y. Some fractional-calculus results for the H-function associated with a class of Feynman integrals. Russ. J. Math. Phys.; 2006; 13, pp. 94-100. [DOI: https://dx.doi.org/10.1134/S1061920806010092]

30. Srivastava, H.M.; Saxena, R.K.; Parmar, R.K. Some families of the incomplete H-functions and the incomplete H-functions and associated integral transforms and operators of fractional calculus with applications. Russ. J. Math. Phys.; 2018; 25, pp. 116-138. [DOI: https://dx.doi.org/10.1134/S1061920818010119]

31. Titchmarsh, E.C. The Theory of the Riemann Zeta-Function; Clarendon (Oxford University) Press: Oxford, UK, London, UK, New York, NY, USA, 1951; 2nd ed.; (Revised by D. R. Heath-Brown), 1986

32. Srivastava, H.M.; Choi, J. Zeta and q-Zeta Functions and Associated Series and Integrals; Elsevier Science Publishers: Amsterdam, The Netherlands, London, UK, New York, NY, USA, 2012.

33. Srivastava, H.M. The Zeta and related functions: Recent developments. J. Adv. Eng. Comput.; 2019; 3, pp. 329-354. [DOI: https://dx.doi.org/10.25073/jaec.201931.229]

34. Srivastava, H.M. Some general families of the Hurwitz-Lerch Zeta functions and their applications: Recent developments and directions for further researches. Proc. Inst. Math. Mech. Nat. Acad. Sci. Azerbaijan; 2019; 45, pp. 234-269. [DOI: https://dx.doi.org/10.29228/proc.7]

35. Whittaker, E.T.; Watson, G.N. A Course of Modern Analysis: An Introduction to the General Theory of Infinite Processes and of Analytic Functions; with an Account of the Principal Transcendental Functions; 4th ed. Cambridge University Press: Cambridge, UK, London, UK, New York, NY, USA, 1927.

36. Gradshteyn, I.S.; Ryzhik, I.M. Tables of Integrals, Series, and Products; Corrected and Enlarged Edition Jeffrey, A. Academic Press: New York, NY, USA, London, UK, 1980; 6th ed.; 2000

37. Srivastava, H.M.; Saxena, R.K.; Pogány, T.K.; Saxena, R. Integral and computational representations of the extended Hurwitz-Lerch Zeta function. Integr. Transf. Spec. Funct.; 2011; 22, pp. 487-506. [DOI: https://dx.doi.org/10.1080/10652469.2010.530128]

38. Srivastava, H.M. A new family of the λ-generalized Hurwitz-Lerch zeta functions with applications. Appl. Math. Inform. Sci.; 2014; 8, pp. 1485-1500. [DOI: https://dx.doi.org/10.12785/amis/080402]

39. Mittag-Leffler, G.M. Sur la nouvelle fonction E_α(x). C. R. Acad. Sci. Paris; 1903; 137, pp. 554-558.

40. Wiman, A. Über den Fundamentalsatz in der Theorie der Funcktionen E_α(x). Acta Math.; 1905; 29, pp. 191-201. [DOI: https://dx.doi.org/10.1007/BF02403202]

41. Wiman, A. Über die Nullstellen der Funktionen E_α(x). Acta Math.; 1905; 29, pp. 217-234. [DOI: https://dx.doi.org/10.1007/BF02403204]

42. Mainardi, F. Why the Mittag-Leffler function can be considered the queen function of the fractional calculus?. Entropy; 2020; 22, 1359. [DOI: https://dx.doi.org/10.3390/e22121359] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33266284]

43. Gorenflo, R.; Mainardi, F.; Srivastava, H.M. Special functions in fractional relaxation-oscillation and fractional diffusion-wave phenomena. Proceedings of the Eighth International Colloquium on Differential Equations; Plovdiv, Bulgaria, 18–23 August 1997; Bainov, D. VSP Publishers: Utrecht, The Netherlands, Tokyo, Japan, 1998; pp. 195-202.

44. Haubold, H.J.; Mathai, A.M.; Saxena, R.K. Mittag-Leffler functions and their applications. J. Appl. Math.; 2011; 2011, 298628. [DOI: https://dx.doi.org/10.1155/2011/298628]

45. Kilbas, A.A.; Saigo, M.; Saxena, R.K. Generalized Mittag-Leffler function and generalized fractional calculus operators. Integr. Transf. Spec. Funct.; 2004; 15, pp. 31-49. [DOI: https://dx.doi.org/10.1080/10652460310001600717]

46. Kilbas, A.A.; Saigo, M.; Saxena, R.K. Solution of Volterra integro-differential equations with generalized Mittag-Leffler functions in the kernels. J. Integral Equ. Appl.; 2002; 14, pp. 377-396. [DOI: https://dx.doi.org/10.1216/jiea/1181074929]

47. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; North-Holland Mathematical Studies Elsevier (North-Holland) Science Publishers: Amsterdam, The Netherlands, London, UK, New York, NY, USA, 2006; Volume 204.

48. Hilfer, R.; Seybold, H.J. Computation of the generalized Mittag-Leffler function and its inverse in the complex plane. Integr. Transf. Spec. Funct.; 2006; 17, pp. 637-652. [DOI: https://dx.doi.org/10.1080/10652460600725341]

49. Seybold, H.J.; Hilfer, R. Numerical results for the generalized Mittag-Leffler function. Fract. Calc. Appl. Anal.; 2005; 8, pp. 127-139.

50. Srivastava, H.M. Some families of Mittag-Leffler type functions and associated operators of fractional calculus. TWMS J. Pure Appl. Math.; 2016; 7, pp. 123-145.

51. Srivastava, H.M.; Fernandez, A.; Baleanu, D. Some new fractional-calculus connections between Mittag-Leffler functions. Mathematics; 2019; 7, 485. [DOI: https://dx.doi.org/10.3390/math7060485]

52. Srivastava, H.M.; Bansal, M.K.; Harjule, P. A study of fractional integral operators involving a certain generalized multi-index Mittag-Leffler function. Math. Methods Appl. Sci.; 2018; 41, pp. 6108-6121. [DOI: https://dx.doi.org/10.1002/mma.5122]

53. Srivastava, H.M.; Bansal, M.K.; Harjule, P. A class of fractional integral operators involving a certain general multi-index Mittag-Leffler function. Ukr. Math. J.; 2021; in press

54. Srivastava, H.M. An introductory overview of fractional-calculus operators based upon the Fox-Wright and related higher transcendental functions. J. Adv. Eng. Comput.; 2021; 5, pp. 135-166.

55. Kilbas, A.A.; Saigo, M. H-Transforms: Theory and Applications; Analytical Methods and Special Functions: An International Series of Monographs in Mathematics; Chapman and Hall (A CRC Press Company): Boca Raton, MA, USA, London, UK, New York, NY, USA, 2004; Volume 9.

56. Olver, F.W.J.; Lozier, D.W.; Boisvert, R.F.; Clark, C.W. NIST Handbook of Mathematical Functions; With 1 CD-ROM (Windows, Macintosh and UNIX) U. S. Department of Commerce, National Institute of Standards and Technology: Washington, DC, USA, Cambridge University Press: Cambridge, UK, London, UK, New York, NY, USA, 2010.

57. Srivastava, H.M. Charles Fox. Bull. Lond. Math. Soc.; 1980; 12, pp. 67-70. [DOI: https://dx.doi.org/10.1112/blms/12.1.67]

58. Wright, E.M. The asymptotic expansion of integral functions defined by Taylor series. I. Philos. Trans. R. Soc. Lond. Ser. A Math. Phys. Sci.; 1940; 238, pp. 423-451.

59. Wright, E.M. The asymptotic expansion of integral functions defined by Taylor series. II. Philos. Trans. R. Soc. Lond. Ser. A Math. Phys. Sci.; 1941; 239, pp. 217-232.

60. Wright, E.M. The asymptotic expansion of integral functions and of the coefficients in their Taylor series. Trans. Am. Math. Soc.; 1948; 64, pp. 409-438. [DOI: https://dx.doi.org/10.1090/S0002-9947-1948-0028405-5]

61. Barnes, E.W. The asymptotic expansion of integral functions defined by Taylor’s series. Philos. Trans. R. Soc. Lond. Ser. A Math. Phys. Sci.; 1906; 206, pp. 249-297.

62. Hilfer, R. Applications of Fractional Calculus in Physics; World Scientific Publishing Company: Singapore, 2000.

63. Hilfer, R. Fractional time evolution. Applications of Fractional Calculus in Physics; World Scientific Publishing Company: Singapore, New Jersey, NJ, USA, London, UK, Hong Kong, China, 2000; pp. 87-130.

64. Hilfer, R. Experimental evidence for fractional time evolution in glass forming materials. J. Chem. Phys.; 2002; 284, pp. 399-408. [DOI: https://dx.doi.org/10.1016/S0301-0104(02)00670-5]

65. Hilfer, R. Threefold introduction to fractional derivatives. Anomalous Transport: Foundations and Applications; Klages, R.; Radons, G.; Sokolov, I.M. Wiley-VCH Verlag: Weinheim, Germany, 2008; pp. 17-73.

66. Hilfer, R.; Anton, L. Fractional master equations and fractal time random walks. Phys. Rev. E; 1995; 51, pp. R848-R851. [DOI: https://dx.doi.org/10.1103/PhysRevE.51.R848]

67. Hilfer, R.; Luchko, Y.; Tomovski, Ž. Operational method for solution of the fractional differential equations with the generalized Riemann-Liouville fractional derivatives. Fract. Calc. Appl. Anal.; 2009; 12, pp. 299-318.

68. Mainardi, F.; Gorenflo, R. Time-fractional derivatives in relaxation processes: A tutorial survey. Fract. Calc. Appl. Anal.; 2007; 10, pp. 269-308.

69. Mokhtary, P.; Ghoreishi, F.; Srivastava, H.M. The Müntz-Legendre Tau method for fractional differential equations. Appl. Math. Model.; 2016; 40, pp. 671-684. [DOI: https://dx.doi.org/10.1016/j.apm.2015.06.014]

70. Sandev, T.; Tomovski, Ž. General time fractional wave equation for a vibrating string. J. Phys. A Math. Theoret.; 2010; 43, 055204. [DOI: https://dx.doi.org/10.1088/1751-8113/43/5/055204]

71. Saxena, R.K.; Kalla, S.L. Solution of Volterra-type integro-differential equations with a generalized Lauricella confluent hypergeometric function in the kernels. Int. J. Math. Math. Sci.; 2005; 2005, pp. 1155-1170. [DOI: https://dx.doi.org/10.1155/IJMMS.2005.1155]

72. Srivastava, H.M.; Harjule, P.; Jain, R. A general fractional differential equation associated with an integral operator with the H-function in the kernel. Russ. J. Math. Phys.; 2015; 22, pp. 112-126. [DOI: https://dx.doi.org/10.1134/S1061920815010124]

73. Srivastava, H.M.; Saxena, R.K. Some Volterra-type fractional integro-differential equations with a multivariable confluent hypergeometric function as their kernel. J. Integral Equ. Appl.; 2005; 17, pp. 199-217. [DOI: https://dx.doi.org/10.1216/jiea/1181075324]

74. Tomovski, Ž.; Hilfer, R.; Srivastava, H.M. Fractional and operational calculus with generalized fractional derivative operators and Mittag-Leffler type functions. Integr. Transf. Spec. Funct.; 2010; 21, pp. 797-814. [DOI: https://dx.doi.org/10.1080/10652461003675737]

75. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives: Theory and Applications; Gordon and Breach Science Publishers: Yverdon, Switzerland, 1993.

76. Liouville, J. Mémoire sur quelques quéstions de géometrie et de mécanique, et sur un nouveau genre de calcul pour résoudre ces quéstions. J. École Polytech.; 1832; 13, pp. 1-69.

77. Caputo, M. Elasticità e Dissipazionne; Zanichelli: Bologna, Italy, 1969.

78. Srivastava, H.M. Remarks on some fractional-order differential equations. Integr. Transf. Spec. Funct.; 2017; 28, pp. 560-564. [DOI: https://dx.doi.org/10.1080/10652469.2017.1317248]

79. Srivastava, H.M.; Saxena, R.K. Operators of fractional integration and applications. Appl. Math. Comput.; 2001; 118, pp. 1-52. [DOI: https://dx.doi.org/10.1016/S0096-3003(99)00208-8]

80. Raina, R.K. On generalized Wright’s hypergeometric functions and fractional calculus operators. East Asian Math. J.; 2005; 21, pp. 191-203.

81. Haubold, H.J.; Mathai, A.M. The fractional kinetic equation and thermonuclear functions. Astrophys. Space Sci.; 2000; 273, pp. 53-63. [DOI: https://dx.doi.org/10.1023/A:1002695807970]

82. Srivastava, H.M. Fractional-order derivatives and integrals: Introductory overview and recent developments. Kyungpook Math. J.; 2020; 60, pp. 73-116.

83. Atangana, A.; Baleanu, D. New fractional derivatives with nonlocal and non-singular kernel: Theory and application to heat transfer model. Thermal Sci.; 2016; 20, pp. 763-769. [DOI: https://dx.doi.org/10.2298/TSCI160111018A]

84. Caputo, M.; Fabrizio, M. A new definition of fractional derivative without singular kernel. Progr. Fract. Differ. Appl.; 2015; 1, pp. 73-85.

85. Podlubny, I. Fractional Differential Equations: An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications; Mathematics in Science and Engineering Academic Press: New York, NY, USA, London, UK, Sydney, Australia, Tokyo, Japan, Toronto, ON, Canada, 1999; Volume 198.

86. Caputo, M.; Mainardi, F. A new dissipation model based on memory mechanism. Pure Appl. Geophys.; 1971; 91, pp. 134-147. [DOI: https://dx.doi.org/10.1007/BF00879562]

87. Alsaedi, A.; Baleanu, D.; Etemad, S.; Rezapour, S. On coupled systems of time-fractional differential problems by using a new fractional derivative. J. Funct. Spaces; 2016; 2016, 4626940. [DOI: https://dx.doi.org/10.1155/2016/4626940]

88. Atangana, A. On the new fractional derivative and application to nonlinear Fisher’s reaction diffusion equation. Appl. Math. Comput.; 2016; 273, pp. 948-956. [DOI: https://dx.doi.org/10.1016/j.amc.2015.10.021]

89. Atangana, A.; Alkahtani, B.S.T. Extension of the RLC electrical circuit to fractional derivative without singular kernel. Adv. Mech. Eng.; 2015; 7, 1687814015591937.

90. Atangana, A.; Alkahtani, B.S.T. Analysis of the Keller-Segel model with a fractional derivative without singular kernel. Entropy; 2015; 17, pp. 4439-4453. [DOI: https://dx.doi.org/10.3390/e17064439]

91. Atangana, A.; Alkahtani, B.S.T. New model of groundwater flowing within a confine aquifer: Application of Caputo-Fabrizio derivative. Arab. J. Geosci.; 2016; 9, 8. [DOI: https://dx.doi.org/10.1007/s12517-015-2060-8]

92. Atangana, A.; Nieto, J.J. Numerical solution for the model of RLC circuit via the fractional derivative without singular kernel. Adv. Mech. Eng.; 2015; 7, 1687814015613758. [DOI: https://dx.doi.org/10.1177/1687814015613758]

93. Goodrich, C.; Peterson, A.C. Discrete Fractional Calculus; Springer: Berlin/Heidelberg, Germany, New York, NY, USA, 2015.

94. Ahmad, B.; Henderson, J.; Luca, R. Boundary Value Problems for Fractional Differential Equations and Systems; Trends in Abstract and Applied Analysis World Scientific Publishing Company: Singapore, Hoboken, NJ, USA, London, UK, Hong Kong, China, 2021; Volume 9.

95. 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. [DOI: https://dx.doi.org/10.1016/j.cnsns.2021.105735]

96. Anastassiou, G.A. Generalized Fractional Calculus: New Advancements and Applications; Studies in Systems, Decision and Control Springer: Cham, Switzerland, 2021; Volume 305.

97. Badík, A.; Fečkan, M. Applying fractional calculus to analyze final consumption and gross investment influence on GDP. J. Appl. Math. Stat. Inform.; 2021; 17, pp. 65-72. [DOI: https://dx.doi.org/10.2478/jamsi-2021-0004]

98. Chen, S.-B.; Jahanshahi, H.; Abba, O.A.; Solís-Pérez, J.E.; Bekiros, S.; Gómez-Aguilar, J.F.; Yousefpour, A.; Chu, Y.-M. The effect of market confidence on a financial system from the perspective of fractional calculus: Numerical investigation and circuit realization. Chaos Solitons Fractals; 2020; 140, 110223. [DOI: https://dx.doi.org/10.1016/j.chaos.2020.110223]

99. Garrappa, R.; Giusti, A.; Mainardi, F. Variable-order fractional calculus: A change of perspective. Commun. Nonlinear Sci. Numer. Simul.; 2021; 102, 105904. [DOI: https://dx.doi.org/10.1016/j.cnsns.2021.105904]

100. Gómez-Aguilar, J.F.; Atangana, A. New chaotic attractors: Application of fractal-fractional differentiation and integration. Math. Methods Appl. Sci.; 2020; 130, pp. 3036-3065. [DOI: https://dx.doi.org/10.1002/mma.6432]

101. Guariglia, E. Fractional calculus, zeta functions and Shannon entropy. Open Math.; 2021; 19, pp. 87-100. [DOI: https://dx.doi.org/10.1515/math-2021-0010]

102. Srivastava, H.M. Operators of basic (or q-) calculus and fractional q-calculus and their applications in geometric function theory of complex analysis. Iran. J. Sci. Technol. Trans. A Sci.; 2020; 44, pp. 327-344. [DOI: https://dx.doi.org/10.1007/s40995-019-00815-0]

103. Saad, K.M. New fractional derivative with non-singular kernel for deriving Legendre spectral collocation method. Alex. Energy J.; 2020; 59, pp. 1909-1917. [DOI: https://dx.doi.org/10.1016/j.aej.2019.11.017]

104. Saad, K.M.; Al-Shareef, E.H.F.; Alomari, A.K.; Baleanu, D.; Gómez-Aguilar, J.F. On exact solutions for time-fractional Korteweg-de Vries and Korteweg-de Vries-Burgers equations using homotopy analysis transform method. Chin. J. Phys.; 2020; 63, pp. 149-162. [DOI: https://dx.doi.org/10.1016/j.cjph.2019.11.004]

105. Saad, K.M.; Gómez-Aguilar, J.F.; Almadiy, A.A. A fractional numerical study on a chronic hepatitis C virus infection model with immune response. Chaos Solitons Fractals; 2020; 139, 110062. [DOI: https://dx.doi.org/10.1016/j.chaos.2020.110062]

106. Shishkina, E.; Sitnik, S. Transmutations, Singular and Fractional Differential Equations with Applications to Mathematical Physics; Mathematics in Science and Engineering Academic Press (Elsevier Science Publishers): New York, NY, USA, London, UK, Toronto, ON, Canada, 2020.

107. Srivastava, H.M.; Alomari, A.-K.N.; Saad, K.M.; Hamanah, W.M. Some dynamical models involving fractional-order derivatives with the Mittag-Leffler type kernels and their applications based upon the Legendre spectral collocation method. Fractal Fract.; 2021; 5, 131. [DOI: https://dx.doi.org/10.3390/fractalfract5030131]

108. Srivastava, H.M.; Deniz, S.; Saad, K.M. An efficient semi-analytical method for solving the generalized regularized long wave equations with a new fractional derivative operator. J. King Saud Univ. Sci.; 2021; 33, 101345. [DOI: https://dx.doi.org/10.1016/j.jksus.2021.101345]

109. Srivastava, H.M.; Jan, R.; Jan, A.; Deebani, W.; Shutaywi, M. Fractional-calculus analysis of the transmission dynamics of the dengue infection. Chaos; 2021; 31, 53130. [DOI: https://dx.doi.org/10.1063/5.0050452]

110. Wang, W.-T.; Khan, M.A. Analysis and numerical simulation of fractional model of bank data with fractal-fractional Atangana-Baleanu derivative. J. Comput. Appl. Math.; 2020; 369, 112646. [DOI: https://dx.doi.org/10.1016/j.cam.2019.112646]

111. You, Z.-L.; Fečkan, M.; Wang, J.-R. Relative controllability of fractional delay differential equations via delayed perturbation of Mittag-Leffler functions. J. Comput. Appl. Math.; 2020; 378, 112939. [DOI: https://dx.doi.org/10.1016/j.cam.2020.112939]