This work is licensed under the Creative Commons Attribution 4.0 International License.
Introduction and Preliminaries
Fractional calculus is a very rapidly growing subject of mathematics which deals with the study of fractional order derivatives and integrals. Fractional calculus is an efficient tool to study many complex real world systems [1]. It is demonstrated that fractional order representation of complex processes appearing in various fields of science, engineering and finance, provides a more realistic approach with memory effects to study these problems. (see e.g. [2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14] and [15]). Among the research work developing the theory of fractional calculus and presenting some applications, we have to point out some literature (see [16, 17, 18, 19, 20]). Kumar et al. [21] analyzed the fractional model of modified Kawahara equation by using newly introduced Caputo-Fabrizio fractional derivative. One also et al. [22] studied a heat transfer problem and presented a new non-integer model for convective straight fins with temperature-dependent thermal conductivity associated with Caputo-Fabrizio fractional derivative. Recently, one et al. [23] presented a new fractional extension of regularized long wave equation by using Atangana-Baleano fractional operator. In et al. [24] one introduced a new numerical scheme for fractional Fitzhugh-Nagumo equation arising in transmission of new impulses. In et al. [25] one constituted a modified numerical scheme to study fractional model of Lienard’s equations. Hajipour et al. [26] in their work formulated a new scheme for class of fractional chaotic systems. Baleanu et al. [27] proposed a new formulation of the fractional control problems involving Mittag-Leffler non-singular kernel. In another work, Baleanu et al. [28] studied the motion of a Bead sliding on a wire in fractional analysis. Jajarmi et al. [29] analyzed a hyperchaotic financial system and its chaos control and synchronization by using fractional calculus.
For mathematical modeling of many complex problems appearing in various fields of science and engineering such as fluid dynamics, plasma physics, astrophysics, image processing, stochastic dynamical system, controlled thermonuclear fusion, nonlinear control theory, nonlinear biological systems, quantum physics and heat transfer problems, the fractional calculus operators involving various special functions have been used successfully. There is rich literature available revealing the notable development in fractional order derivatives and integrals (see, [1, 10, 11, 18, 19, 20, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39]). Recently, Caputo and Fabrizio [40] introduced a new fractional derivative which is more suitable than the classical Caputo fractional derivative for many engineering and thermodynamical processes. Atangana [41] used a new fractional derivative to study the nature of Fisher’s reaction diffusion equation. Riemann and Caputo fractional derivative operators both have a singular kernel which cannot exactly represent the complete memory effect of the system. To overcome these limitations of the old derivatives, very recently Atangana and Baleanu [42] presented a new non-integer order derivative having a non-local, non-singular and Mittag-Leffler type kernel.
In recent years, many researchers have extensively studied the properties, applications and extensions of various fractional integral and differential operators involving the various special functions. (for detail see McBride [43], Kalla [44, 45], Kalla and Saxena [46, 47], Saigo [48, 49, 50], Saigo and Maeda [51], Kiryakova [32, 52], [53] etc).
For our present study, we recall the following pair of Saigo hypergeometric fractional integral operators.
For x > 0, λ, σ, ϑ ∈ ℂ and ℛ(λ) > 0, we have
\left( {I_{0,x}^{\lambda ,\sigma ,\vartheta } f(t)} \right)(x) = \frac{{x^{ - \lambda - \sigma } }}{{\Gamma (\lambda )}}\int_0^x (x - t)^{\lambda - 1} _2 F_1 \left( {\lambda + \sigma , - \vartheta ;\lambda ;1 - \frac{t}{x}} \right)f(t){\kern 1pt} dtand
\left( {I_{x,\infty }^{\lambda ,\sigma ,\vartheta } f(t)} \right)(x) = \frac{1}{{\Gamma (\lambda )}}\int_x^\infty (t - x)^{\lambda - 1} t^{ - \lambda - \sigma } _2 F_1 \left( {\lambda + \sigma , - \vartheta ;\lambda ;1 - \frac{x}{t}} \right)f(t){\kern 1pt} dtwhere the 2F1(.), a special case of the generalized hypergeomteric function, is the Gauss hypergeometric function.
The operator I_{0,x}^{\lambda ,\sigma ,\vartheta } (.) contains the Riemann-Liouville R_{0,x}^\lambda (.) fractional integral operators by means of the following relationships:
\left( {R_{0,x}^\lambda f(t)} \right)(x) = \left( {I_{0,x}^{\lambda , - \lambda ,\vartheta } f(t)} \right)(x) = \frac{1}{{\Gamma (\lambda )}}\int_0^x (x - t)^{\lambda - 1} f(t){\kern 1pt} dt\left( {W_{x,\infty }^\lambda f(t)} \right)(x) = \left( {J_{x,\infty }^{\lambda , - \lambda ,\vartheta } f(t)} \right)(x) = \frac{1}{{\Gamma (\lambda )}}\int_x^\infty (t - x)^{\lambda - 1} f(t){\kern 1pt} dt
The image formulas for special functions of one or more variables are very useful in the evaluation and solution of differential and integral equations. Motivating by the above discussion, we developed new fractional calculus formulas involving extended generalized Mathieu series.
The following familiar infinite series
S(r) = \sum\limits_{n = 1}^\infty \frac{{2n}}{{{{({n^2} + {r^2})}^2}}},\,\,\,\,\,\,\,\,\,\,\,\,(r \in {\mathbb{R}^ + }),is called a Mathieu series. It was introduced and studied by Émile Leonard Mathieu in his book [55] devoted to the elasticity of solid bodies. Bounds for this series are needed for the solution of boundary value problems for the biharmonic equations in a two dimensional rectangular domain, see [56, Eq. (54), p. 258]. Several interesting problems and solutions dealing with integral representations and bounds for the following generalization of the Mathieu series, which is so-called generalized Mathieu series with a fractional power can be found in [57, 58, 60]:
S_\mu (r) = \sum\limits_{n = 1}^\infty \frac{{2n}}{{(n^2 + r^2 )^{\mu + 1} }},\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,(\mu > 0,r > 0).
In [59], the authors derived the following new Laplace type integral representation series
S_\mu (r) = \frac{{\sqrt \pi }}{{2^{\mu - \frac{1}{2}} \Gamma (\mu + 1)}}\int_0^\infty e^{ - rt} \kappa _\mu (t)dt,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\left( {\mu > \frac{3}{2}} \right)\kappa _\mu (t) = t^{\mu + \frac{1}{2}} \sum\limits_{k = 1}^\infty \frac{{J_{\mu + \frac{1}{2}} (kt)}}{{k^{\mu - \frac{1}{2}} }}
and Jμ(z) is the Bessel function. Motivated essentially by the works of Cerone and Lenard [61], Srivastava and Tomovski in [62] defined a family of generalized Mathieu series
S_\mu ^{(\alpha ,\beta )} (r;a) = S_\mu ^{(\alpha ,\beta )} (r;\{ a_n \} _{n = 1}^\infty ) = \sum\limits_{n = 1}^\infty \frac{{2a_n^\beta }}{{(a_n^\alpha + r^2 )^\mu }},\,\,\,\,\,\,\,\,\,\,\,\,\,(\alpha ,\beta ,\mu ,r > 0),where it is tacitly assumed that the positive sequence
a = \left\{ {a_n } \right\} = \left\{ {a_1 ,a_2 ,...} \right\}
such that
lim_{n \to \infty } a_n = \infty
is so chosen that the infinite series in definition (1.11) converges, that is, that the following auxiliary series
\sum\limits_{n = 1}^\infty \frac{1}{{a_n^{\mu \alpha - \beta } }}
is convergent.
Definition 1
(see [63, Eq. (6.1), p. 256] ) The extended Beta function Bp;q(x; y) is defined by\begin{gathered}
{B_{p,q}}(x,y) = \int_0^1 {t^{x - 1}}{(1 - t)^{y - 1}}{E_{p,q}}dt, \hfill \\
\left( {x,y,p,q, \in \mathbb{C};\min \left\{ {\Re (x),\Re (y)} \right\} > 0,\,\,\min \left\{ {\Re (p),\Re (q)) \geqslant 0} \right\}} \right. \hfill \\
\end{gathered}where Ep,q(t) is defined by
E_{p,q} (t) = \exp \left( { - \frac{p}{t} - \frac{q}{{1 - t}}} \right)
(p, q ∈ ℂ and min{ℛ(p), ℛ(q)} > 0).
In particular, Chaudhry et al. [64, p. 591, Eq. (1.7)], introduced the p–extension of Euler’s Beta function B(x, y) :
B_p (x,y) = \int_0^1 t^{x - 1} (1 - t)^{y - 1} e^{ - \frac{p}{{t(1 - t)}}} dt
(ℛ(p) > 0) whose special case, when p = 0 ( or p = q = 0 in (1.12)), is the familiar Beta integral
B(x,y) = \int_0^1 t^{x - 1} (1 - t)^{y - 1} dt
(ℛ(x), ℛ(y) > 0).
Recently, Mehrez and Tomovski [65] introduces the (p, q)-Mathieu-type power series in terms of the extended Beta function (1.12), which is defined as:
\begin{array}{l} S_{\mu ,\vartheta ,\tau ,\xi }^{\alpha ,\beta } (r;a;p,q;z) = \sum\limits_{n = 1}^\infty \frac{{2a_n^\beta (\vartheta )_n B_{p,q} (\tau + n,\xi - \tau )}}{{B(\tau ,\xi - \tau )(a_n^\alpha + r^2 )^\mu }}\frac{{z^n }}{{n!}} \\ (r,\alpha ,\beta ,\nu > 0;\xi > \tau > 0;p,q \in \mathbb{C};\min \{ \Re (p),\Re (q)\} \ge 0;|z| < 1) \\ \end{array}
In particular case when p = q; we define the p-Mathieu-type power series defined by
\begin{array}{*{20}c} {S_{\mu ,\vartheta ,\tau ,\xi }^{\alpha ,\beta } (r;a;p;z) = \sum\limits_{n = 1}^\infty \frac{{2a_n^\beta (\vartheta )_n B_p (\tau + n,\xi - \tau )}}{{B(\tau ,\xi - \tau )(a_n^\alpha + r^2 )^\mu }}\frac{{z^n }}{{n!}}} \\ {(\tau ,\alpha ,\beta ,\vartheta ,\xi ,\tau > 0,p \in C,|z| \le 1)} \\\end{array}
The function S_{\mu ,\vartheta ,\tau ,\xi }^{\alpha ,\beta } (\tau ;a;p;z) has many other special cases. If we set p = q = 0; we get
S_{\mu ,\vartheta ,\tau ,\xi }^{\alpha ,\beta } (r;a;z) = \sum\limits_{n = 1}^\infty \frac{{2a_n^\beta (\vartheta )_n }}{{(a_n^\alpha + r^2 )^\mu (\xi )_n }}\frac{{z^n }}{{n!}}
(τ, α, β, ϑ, ξ > 0, |z| ≤ 1)
On the other hand, by letting τ = ω in (1.15) we obtain [66, Eq. 5, p. 974]:
S_{\mu ,\vartheta }^{\alpha ,\beta } (r;a;z) = \sum\limits_{n = 1}^\infty \frac{{2a_n^\beta (\vartheta )_n }}{{(a_n^\alpha + r^2 )^\mu }}\frac{{z^n }}{{n!}},\,\,\,\,\,\,\,\,\,\,\,(\tau ,\alpha ,\beta ,\vartheta > 0,|z| \le 1).
The concept of the Hadamard product (or the convolution) of two analytic functions is very useful in our present study. It can help us to decompose a newly emerging function into two known functions. Let
f(z): = \sum\limits_{n = 0}^\infty a_n z^n ,(|z| < R_f )and
g(z): = \sum\limits_{n = 0}^\infty b_n z^n ,(|z| < R_g )
be two power series whose radii of convergence are denoted by Rf and Rg, respectively.
Then their Hadamard product is the power series defined by
\begin{array}{l} (f*g)(z): = \sum\limits_{n = 0}^\infty a_n b_n z^n = (g*f)(z) \\ (|z| < R), \\ \end{array}where
R = \mathop {\lim }\limits_{n \to \infty } \left| {\frac{{a_n b_n }}{{a_{n + 1} b_{n + 1} }}} \right| = \mathop {\lim }\limits_{n \to \infty } \left| {\frac{{a_n }}{{a_{n + 1} }}} \right|.\mathop {\lim }\limits_{n \to \infty } \left| {\frac{{b_n }}{{b_{n + 1} }}} \right| = R_f R_g
Therefore, in general, we have R ≥ Rf.Rg [67, 68].
For various investigations involving the Hadamard product (or the convolution), the interested reader may refer to several recent papers on the subject (see, for example, [69, 70] and the references cited therein).
Fractional integration
In this section, we will establish some fractional integral formulas for the generalized (p, q)-Mathieu-type power series. Then their special cases also introduced here.
In this section, we establish certain theorems involving the results obtained in previous section associated with the integral transforms like, Beta transform, Laplace transform and Whittaker transform.
Beta Transform
The Beta transform of f(z) is defined as [71]:B\{ f(z):a,b\} = \int_0^1 z^{a - 1} (1 - z)^{b - 1} f(z)dz
The proof of this theorem would run parallel as those of Theorem 7.
Fractional Kinetic Equations
The importance of fractional differential equations in the field of applied science has gained more attention not only in mathematics but also in physics, dynamical systems, control systems and engineering, to create the mathematical model of many physical phenomena. Especially, the kinetic equations describe the continuity of motion of substance. The extension and generalization of fractional kinetic equations involving many fractional operators were found in [72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85].
In view of the effectiveness and a great importance of the kinetic equation in certain astrophysical problems the authors develop a further generalized form of the fractional kinetic equation involving generalized k-Mittag-Leffler function.
The fractional differential equation between rate of change of the reaction, the destruction rate and the production rate was established by Haubold and Mathai [78] given as follows:\frac{{dN}}{{dt}} = - d(N_t ) + p(N_t ),where N = N(t) the rate of reaction, d = d(N) the rate of destruction, p = p(N) the rate of production and Nt denotes the function defined by Nt(t*) = N(t − t*), t* > 0.
The special case of (4.1) for spatial fluctuations and inhomogeneities in N(t) the quantities are neglected, that is the equation\frac{{dN}}{{dt}} = - c_i N_i (t),
with the initial condition that Ni(t = 0) = N0 is the number density of the species i at time t = 0 and ci > 0. If we remove the index i and integrate the standard kinetic equation (4.2), we haveN(t) - N_0 = - c_0 D_t^{ - 1} N(t)where _0 D_t^{ - 1} is the special case of the Riemann-Liouville integral operator _0 D_t^{ - \nu } defined as{}_0D_t^{ - \nu } f(t) = \frac{1}{{\Gamma (\nu )}}\int_0^t \left( {t - s} \right)^{\nu - 1} f(s)ds,\,\,\,\,\,\,\,\,\,\,\,\,(t > 0,R(\nu ) > 0)
The fractional generalization of the standard kinetic equation(4.3) is given by Haubold and Mathai [78] as follows:N(t) - N_0 = - c^\nu {}_0D_t^{ - 1} N(t)and obtained the solution of (4.5) as follows:N(t) = N_0 \sum\limits_{k = 0}^\infty \frac{{( - 1)^k }}{{\Gamma \left( {\nu k + 1} \right)}}\left( {ct} \right)^{\nu k}
Further, (Saxena and Kalla [83]) considered the the following fractional kinetic equation:N(t) - N_0 f(t) = - c^\nu {}_0D_t^{ - \nu } N(t),\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,(\Re (v) > 0),where N(t) denotes the number density of a given species at time t, N0 = N(0) is the number density of that species at time t = 0, c is a constant and f ∈ ℒ (0, ∞).
By applying the Laplace transform to (4.7) (see [79]),\begin{array}{l} L\left\{ {N(t);p} \right\} = N_0 \frac{{F(p)}}{{1 + c^\nu p^{ - \nu } }} = N_0 \left( {\sum\limits_{n = 0}^\infty ( - c^\nu )^n p^{ - \nu n} } \right)F(p), \\ \,\,\,\,\,\,\,\,\,\,\,\,\left( {n \in N_0 ,\left| {\frac{c}{p}} \right| < 1} \right) \\ \end{array}where the Laplace transform [86] is given by
F(p) = L\left\{ {N(t);p} \right\} = \int_0^\infty e^{ - pt} f(t)dt,\,\,\,\,\,\,\,\,\,({\cal R}(p) > 0).
Solution of generalized fractional kinetic equations
In this section, we investigated the solutions of the generalized fractional kinetic equations by considering generalized (p, q)-Mathieu Type Series
Theorem 9
If a > 0, d > 0, ν > 0, r, α, β, ϑ > 0; ξ > τ > 0; p, q ∈ ℂ;min{ℛ(p), ℛ(q)} ≥ 0, then the solution of the fractional kinetic equationN(t) - N_0 S_{\mu ,\vartheta ,\tau ,\xi }^{\alpha ,\beta } (r;a;p,q;d^\nu t^\nu ) = - a^\nu {}_0D_t^{ - \nu } N(t)is given by the following formula\begin{array}{l} N(t) = N_0 \sum\limits_{n = 1}^\infty \frac{{2a_n^\beta (\vartheta )_n B_{p,q} (\tau + n,\xi - \tau )}}{{B(\tau ,\xi - \tau )(a_n^\alpha + r^2 )^\mu }}\frac{{(a^\nu t^\nu )^n }}{{n!}} \\ \,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\, \times \Gamma (\nu n + 1)E_{\nu ,\nu n + 1} ( - a^\nu t^\nu ). \\ \end{array}
The equation (5.12) can be written asN(t) = N_0 \sum\limits_{n = 1}^\infty \frac{{2a_n^\beta (\vartheta )_n B_{p,q} (\tau + n,\xi - \tau )}}{{B(\tau ,\xi - \tau )(a_n^\alpha + r^2 )^\mu }}\Gamma (\nu n + 1)\frac{{(a^\nu t^\nu )^n }}{{n!}}E_{\nu ,\nu n + 1} ( - a^\nu t^\nu ).
Theorem 10
If d > 0, ν > 0, r, α, β, ϑ > 0; ξ > τ > 0; p, q ∈ ℂ; min{ℛ(p), ℛ(q)} ≥ 0, then the solution of the fractional kinetic equationN(t) - N_0 S_{\mu ,\vartheta ,\tau ,\xi }^{\alpha ,\beta } (r;a;p,q;d^\nu t^\nu ) = - d^\nu {}_0D_t^{ - \nu } N(t)is given by the following formula\begin{array}{l} N(t) = N_0 \sum\limits_{n = 1}^\infty \frac{{2a_n^\beta (\vartheta )_n B_{p,q} (\tau + n,\xi - \tau )}}{{B(\tau ,\xi - \tau )(a_n^\alpha + r^2 )^\mu }}\frac{{(d^\nu t^\nu )^n }}{{n!}} \\ \,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\, \times \Gamma (\nu n + 1)E_{\nu ,\nu n + 1} ( - d^\nu t^\nu ). \\ \end{array}
Theorem 11
If d > 0, ν > 0, r, α, β, ϑ > 0; ξ > τ > 0; p, q ∈ ℂ; min{ℛ(p), ℛ(q)} ≥ 0, then the solution of the fractional kinetic equationN(t) - N_0 S_{\mu ,\vartheta ,\tau ,\xi }^{\alpha ,\beta } (r;a;p,q;t) = - d^\nu {}_0D_t^{ - \nu } N(t)is given by the following formula\begin{array}{l} N(t) = N_0 \sum\limits_{n = 1}^\infty \frac{{2a_n^\beta (\vartheta )_n B_{p,q} (\tau + n,\xi - \tau )}}{{B(\tau ,\xi - \tau )(a_n^\alpha + r^2 )^\mu }}\frac{{(t^\nu )^n }}{{n!}} \\ \,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\, \times \Gamma (\nu n + 1)E_{\nu ,\nu n + 1} ( - d^\nu t^\nu ). \\ \end{array}
Proof
The proof of the Theorem 10 and Theorem 11 are same as that of Theorem 9, so we would like to skip here.
Special cases
Here we introduce some special cases of our results established in this section.
Case 4.
If p = q, then Theorem 9, Theorem 10 and Theorem 11 reduces to
Corollary 7
If a > 0, d > 0, ν > 0, r, α, β, ϑ > 0; ξ > τ > 0; p ∈ ℂ;ℛ(p) ≥ 0, then the solution of the fractional kinetic equationN(t) - N_0 S_{\mu ,\vartheta ,\tau ,\xi }^{\alpha ,\beta } (r;a;p;d^\nu t^\nu ) = - a^\nu {}_0D_t^{ - \nu } N(t)is given by the following formulaN(t) = N_0 \sum\limits_{n = 1}^\infty \frac{{2a_n^\beta (\vartheta )_n B_p (\tau + n,\xi - \tau )}}{{B(\tau ,\xi - \tau )(a_n^\alpha + r^2 )^\mu }}\frac{{(a^\nu t^\nu )^n }}{{n!}}\Gamma (\nu n + 1)E_{\nu ,\nu n + 1} ( - a^\nu t^\nu ).
Corollary 8
If d > 0, ν > 0, r, α, β, ϑ > 0; ξ > τ > 0; p ∈ ℂ; ℛ(p) ≥ 0, then the solution of the fractional kinetic equationN(t) - N_0 S_{\mu ,\vartheta ,\tau ,\xi }^{\alpha ,\beta } (r;a;p;d^\nu t^\nu ) = - d^\nu {}_0D_t^{ - \nu } N(t)is given by the following formulaN(t) = N_0 \sum\limits_{n = 1}^\infty \frac{{2a_n^\beta (\vartheta )_n B_p (\tau + n,\xi - \tau )}}{{B(\tau ,\xi - \tau )(a_n^\alpha + r^2 )^\mu }}\frac{{(d^\nu t^\nu )^n }}{{n!}}\Gamma (\nu n + 1)E_{\nu ,\nu n + 1} ( - d^\nu t^\nu ).
Corollary 9
If d > 0, ν > 0, r, α, β, ϑ > 0; ξ > τ > 0; p ∈ ℂ;ℛ(p) ≥ 0, then the solution of the fractional kinetic equationN(t) - N_0 S_{\mu ,\vartheta ,\tau ,\xi }^{\alpha ,\beta } (r;a;p;t) = - d^\nu {}_0D_t^{ - \nu } N(t)is given by the following formulaN(t) = N_0 \sum\limits_{n = 1}^\infty \frac{{2a_n^\beta (\vartheta )_n B_p (\tau + n,\xi - \tau )}}{{B(\tau ,\xi - \tau )(a_n^\alpha + r^2 )^\mu }}\frac{{(t^\nu )^n }}{{n!}}\Gamma (\nu n + 1)E_{\nu ,\nu n + 1} ( - d^\nu t^\nu ).
Case 5.
If p = q = 0, then Theorem 9, Theorem 10 and Theorem 11 reduces to
Corollary 10
If a > 0, d > 0, ν > 0, r, α, β, ϑ > 0; ξ > τ > 0, then the solution of the fractional kinetic equationN(t) - N_0 S_{\mu ,\vartheta ,\tau ,\xi }^{\alpha ,\beta } (r;a;d^\nu t^\nu ) = - a^\nu {}_0D_t^{ - \nu } N(t)is given by the following formulaN(t) = N_0 \sum\limits_{n = 1}^\infty \frac{{2a_n^\beta (\vartheta )_n }}{{(a_n^\alpha + r^2 )^\mu (\xi )_n }}\frac{{(a^\nu t^\nu )^n }}{{n!}}\Gamma (\nu n + 1)E_{\nu ,\nu n + 1} ( - a^\nu t^\nu ).
Corollary 11
If d > 0, ν > 0, r, α, β, ϑ > 0; ξ > τ > 0, then the solution of the fractional kinetic equationN(t) - N_0 S_{\mu ,\vartheta ,\tau ,\xi }^{\alpha ,\beta } (r;a;p;d^\nu t^\nu ) = - d^\nu {}_0D_t^{ - \nu } N(t)is given by the following formulaN(t) = N_0 \sum\limits_{n = 1}^\infty \frac{{2a_n^\beta (\vartheta )_n }}{{(a_n^\alpha + r^2 )^\mu (\xi )_n }}\frac{{(d^\nu t^\nu )^n }}{{n!}}\Gamma (\nu n + 1)E_{\nu ,\nu n + 1} ( - d^\nu t^\nu ).
Corollary 12
If d > 0, ν > 0, r, α, β, ϑ > 0; ξ > τ > 0, then the solution of the fractional kinetic equationN(t) - N_0 S_{\mu ,\vartheta ,\tau ,\xi }^{\alpha ,\beta } (r;a;p;t) = - d^\nu {}_0D_t^{ - \nu } N(t)is given by the following formulaN(t) = N_0 \sum\limits_{n = 1}^\infty \frac{{2a_n^\beta (\vartheta )_n }}{{(a_n^\alpha + r^2 )^\mu (\xi )_n }}\frac{{(t^\nu )^n }}{{n!}}\Gamma (\nu n + 1)E_{\nu ,\nu n + 1} ( - d^\nu t^\nu ).
Conclusion
In the present work, fractional integral formulae involving (p, q)-Mathieu Type series has established. The image formulae of our findings by employing integral transform has been also introduced. Further in this work we gave the solution of fractional kinetic equation in terms of Mittag-Leffler function. All the results are general in nature and give numerous results as their special cases.