- Research
- Open Access
- Published:

# Computable solution of fractional kinetic equations using Mathieu-type series

*Advances in Difference Equations*
**volume 2019**, Article number: 234 (2019)

## Abstract

The Mathieu series appeared in the study of elasticity of solid bodies in the work of Émile Leonard Mathieu. Since then numerous authors have studied various problems arising from the Mathieu series in several diverse ways. In this line, our aim is to study the solution of fractional kinetic equations involving generalized Mathieu-type series. The generality of this series will help us to deduce results related to a fractional kinetic equation involving another form of Mathieu series. To obtain the solution, we use the Laplace transform technique. Besides, a graphical representation is given to observe the behavior of the obtained solutions.

## Introduction and preliminaries

Fractional calculus (FC) is a useful mathematical tool to study fractional order integrals and derivatives. The fractional calculus has been developed and used in different areas of science and engineering. The concept of fractional differential equations and their applications have played a significant role in many diverse fields such as applied science, physics, biology, chemistry, and engineering. The kinetic equations designate a system of differential equations, which describe the rate of change of the chemical composition of a star for each order in terms of the reaction rates for production and destruction. During the last several decades, fractional kinetic equations in various forms have been broadly and usefully employed when describing and solving several important problems of physics and astrophysics (see, e.g., [2, 7,8,9,10,11, 13,14,15,16, 20, 24,25,26,27,28,29, 31] and the references therein). The special functions and their applications appear in the solutions of fractional integral and differential equations and are related to comprehensive problems in the several areas of mathematics and mathematical physics. In view of the effectiveness and great importance of the fractional kinetic equations in certain astrophysical problems, the authors develop a generalized form of the fractional kinetic equations along with Mathieu-type series. The broad generality of Mathieu-type series will allow us to deduce many special cases of the main results.

In [14], Haubold and Mathai explored the kinetic equation describing the rate of change of detraction, production, and reaction, which is given as

where \(d=d(N)\) is the rate of destruction, \(p=p(N)\) the rate of production, \(N=N(t)\) the rate of reaction, and \(N_{t}\) denotes the function given by \(N_{t}(t^{*})=N(t-t^{*})\), \(t>0\). A particular case of (1), when special fluctuations or inhomogeneities in quantity \(N(t)\) are decayed, is given by the equation

with basic condition \(N_{i}(t=0)=N_{0}\) describing the variety of density of species *i* at time \(t=0\) and \(c_{i}>0\). Neglecting index *i* and integrating, (2) becomes

where \({}_{0}D_{t}^{-1}\) is the standard integral operator.

Further, a fractional generalization of (3) was given by Haubold and Mathai [14] in the following form:

where \({}_{0}D_{t}^{-\nu }\) represents the Riemann–Liouville integral operator defined as

Saxena et al. [26] investigated solutions of three generalized forms of (4) in terms of the following generalized Mittag-Leffler function (see, e.g., [40]):

We here recall just one of them (see [26, Theorem 2]). The equation

is solved by

where \(\mu ,\nu ,c,d \in \mathbb{R^{+}}\); \(c\neq d\).

Sexena et al. [28] investigated the following generalized fractional kinetic equation:

and its solution is

where \(E_{\nu ,\mu }^{\gamma }(z)\) is an extension of the generalized Mittag-Leffler function (6) (see [23]):

and \((\eta )_{\varepsilon }\) denotes the Pochhammer symbol which is defined (for \(\eta ,\varepsilon \in \mathbb{C}\)) in terms of the gamma function *Γ* as

Saxena and Kalla [25] investigated the following fractional kinetic equation:

where \(f\in {L(0,\infty )}\). Applying the Laplace transform (LT) to (13), we have

Suppose that \(f(t)\) is a real- (or complex-) valued function of the (time) variable \(t>0\) and *s* is real or complex parameter. The Laplace transform of the function \(f(t)\) is defined by

whenever the limit exits (as a finite number). The convolution of two functions \(f(t)\) and \(g(t)\), which are defined for \(t>0\), plays an important role in a number of different physical applications. The Laplace convolution of functions \(f(t)\) and \(g(t)\) is given by the following integral:

which exists if the functions *f* and *g* are at least piecewise continuous. One of the most significant properties possessed by the convolution in connection with the Laplace transform is that the Laplace transform of the convolution of two functions is the product of their transforms (see, e.g., [20, 30, 32, 33]).

### The Laplace Convolution Theorem

*If*
*f*
*and*
*g*
*are piecewise continuous on*
\([0,\infty )\)
*and of exponential order*
*α*
*when*
\(t \rightarrow \infty \), *then*

*We find*

*by using the following well*-*known identity*:

*in which*
\(L^{-1}\)
*denotes the inverse Laplace transform*.

Recently, many papers investigated the solution of generalized fractional kinetic equations involving a variety of special functions by using Laplace and Sumudu transformations. For instance, for generalized fractional kinetic equations with \(f(t)\) in (13) replaced by certain forms of generalized Mittag-Leffler function (11), Saxena et al. [29] presented their corresponding solutions. Moreover, Chaurasia and Kumar [7] investigated the solution of generalized fractional kinetic equations involving M-series, and Choi and Kumar [9] gave their corresponding solutions involving Aleph function. Furthermore, Agarwal et al. [1, 2] investigated solutions of fractional kinetic equations with \(f(t)\) in (13) replaced by certain forms of *k*-Mittag-Leffler function and *k*-Bessel function. Subsequently, Nisar and Qi [22] presented solution of fractional kinetic equations with \(f(t)\) in (13) replaced by certain forms of generalized *k*-Bessel function.

The following familiar infinite series was introduced by Mathieu [21]:

An integral representation of (21) is given by (see [12])

Numerous interesting problems and solutions handling integral representations and related to the following generalization of the Mathieu series with fractional power

can be found in the works by Cerone and Lenard [6] as well as Tomovski and Trencevski [37]. Extending the work of Cerone and Lenard [6], Srivastava and Tomovski [34] defined a family of generalized Mathieu series

where the positive sequence \(d= \lbrace d_{k} \rbrace _{k=1}^{\infty }= \lbrace d_{1},d_{2},\dots \rbrace \) (\(\lim_{k\rightarrow \infty }d_{k}=\infty \)) is such that the infinite series \(\sum_{k=1}^{\infty }\frac{1}{d^{\mu {\rho }-\sigma } _{k}}\) is convergent.

In the sequel, Tomovski and Mehrez [35] proposed a generalization of definition (24) in the following power series:

Evidently, the case \(d_{k} = k\), \(\rho =2\), \(\sigma =1=\tau \) and *μ* replaced with \(\mu +1\) corresponds to the Mathieu series defined by Tomovski and Pogány [36]:

Due to the great importance of fractional kinetic equations involving special functions, in this paper, we aim to investigate solution of generalized fractional kinetic equations (13) with \(f(t)\) replaced by several generalized Mathieu-type series, by mainly using the Laplace and Sumudu transforms. The results presented here, being general, are also shown to reduce to fractional kinetic equations involving simpler special functions. The manuscript is organized as follows. In Sect. 2, the solution of fractional kinetic equations involving Mathieu-type series is established, and in Sect. 3, a graphical interpretation and nature of the solution are discussed. Section 4 is focused on some examples, and in Sect. 5 concluding remarks are given.

## Solution of fractional kinetic equations

In this section, we obtain the solution of generalized fractional kinetic equation (13) involving generalized Mathieu-type series (25) by applying the Laplace transform technique.

### Theorem 2.1

*For all*
\(\nu , c, h >{0}\); \(c\neq {h}\), *if*
\(\mu , \rho , \sigma , \tau , l \in {\mathbb{R}^{+}}\), \(|ht|\leq {1}\), *then equation*

*has the solution*

### Proof

Applying the Laplace transform to (27) and using (25) and (18) gives

where \(N^{*}(s)=L \lbrace N(t);s \rbrace \).

Under the given assumptions, computing the integral in (29) term by term and using \(L \lbrace t^{\lambda };s \rbrace =\frac{ \varGamma (\lambda +1)}{s^{\lambda +1}}\), we have

Employing the geometric series expansion of \(( 1+ ( \frac{c}{s} )^{\nu } )^{-1}\) for \(c<|s|\), we have

Taking the inverse Laplace transform and applying the relation \(L^{-1}\{s^{-\nu };t\}=\frac{t^{\nu -1}}{\varGamma (\nu )}\), \(\Re ( \nu )>0\), we get

□

In view of relation (26), we state the following consequence of Theorem 2.1.

### Corollary

*For all*
\(\nu , c, h >{0}\); \(c\neq {h}\), *if*
\(\mu , l \in {\mathbb{R} ^{+}}\), \(|ht|\leq {1}\), *then equation*

*has the solution*

### Theorem 2.2

*Assume that*
\(\nu , c, w >{0}\)
*and*
\(\mu , \rho , \sigma , \tau , l \in {\mathbb{R}^{+}}\), \(|wt|\leq {1}\), *then equation*

*has the solution*

### Proof

Applying LT to both sides of (33) and then using (25) and (15), we get

which, upon solving for \(N^{*}(s)\), yields

Employing the binomial formula

which obviously converges for all \(|u|<1\), we have

If we now take the inverse LT of (35) and use (11), after easy simplification, the result (34) readily follows. □

In view of relation (26), we establish the following consequence of Theorem 2.2.

### Corollary

*Assume that*
\(\nu , \alpha , c, w >{0}\), *and let*
\(\mu , l \in {\mathbb{R} ^{+}}\), \(|wt|\leq {1}\), *then equation*

*has the solution*

### Theorem 2.3

*For all*
\(\nu >0\), \(c>{0}\), \(\mu , \rho , \sigma , \tau , l \in {\mathbb{R} ^{+}}\), \(|t|\leq {1}\), *the equation*

*has the solution*

### Proof

This can be proved by the same procedure as in the proof of Theorem 2.1. So we omit all details. □

Using relation (26), we state the following consequence of Theorem 2.3.

### Corollary

*For all*
\(\nu >0\), \(c>{0}\), \(\mu , l \in {\mathbb{R}^{+}}\), \(|t|\leq {1}\), *the equation*

*has the solution*

## Graphical representations

Figures 1–5 represents the graph of solution (28) by considering some fractional values, and observe that, by choosing small values for *ν*, the graph is increasing in nature, and if we choose larger values for *ν*, then it is gradually decreasing. The graphical results confirm that the region of convergence of solutions depends continuously on the fractional parameter *ν*. Hence, by observing the nature of the solutions for various parameters and time interval, it is concluded that \(N(t)\) can be negative as well as positive. A similar observation can be obtained for the solutions (34) and (39).

## Example

A detailed account of Mathieu-type series and their applications can be found in the monographs by various authors [21, 35,36,37]. An integral transform, called the Sumudu transform, was defined and studied by Watugala [38] to facilitate the process of solving differential and integral equations in the time domain, and for the use in various applications of system control engineering and applied mathematics. In [3, 4], some fundamental properties of the Sumudu transform are investigated. It turns out that the Sumudu transform has very special and useful properties, and it is useful in solving problems of science and control engineering governing kinetic equations. The Sumudu transform is defined on the set of functions

by

where *M* is a real finite number and \(\delta _{1}\) and \(\delta _{2}\) can be finite or infinite. For more details, see [38, 39].

By using the convolution theorem for the Sumudu transform [3,4,5], (5) can be written in the following form:

It is easy to see that the Sumudu transform of the function \(f(t)=t^{\varepsilon }\) is given by

The Laplace transform is a very potent mathematical tool which is used to solve diverse engineering and science problems. The Sumudu transform (which is not so well-known as the Laplace transform) was proposed in the early 1990s, and it has some interesting advantages over other integral transforms, especially regarding the “unity” feature which could become convenient when solving differential equations. As a comparison, the Sumudu transform is a simple variant of the Laplace transform. This paper focuses on the effectiveness of both transforms to obtain the solutions of fractional kinetic equations. The Sumudu transform is essentially identical with the Laplace transform. Given an initial \(f(t)\), its Laplace transform \(F(s)\) can be translated into the Sumudu transform \(F_{s}(u)\) of *f* by means of the relation

and its inverse is

For more details about the Sumudu transform and its properties in comparison to the Laplace transform, the interested readers can check [17,18,19].

Due to the importance of the above observation, in this section, we first define the Sumudu transform of Mathieu-type series (25) and further determine the solution of fractional kinetic equations by applying Sumudu transform technique as given in the following Examples 4.1, 4.2, and 4.3.

The following integral gives the Sumudu transform of Mathieu-type series (25):

where \(\Re (u)>|w|^{-1/\Re (\nu )}\), \(\Re (\varepsilon )>0\), \(u\in (-\delta _{1},\delta _{2})\), \(|f(t)|< Me^{|t|/\delta _{j}}\) and \(l,d,\rho , \sigma , \mu \in \mathbb{R^{+}}\).

Using the same procedure of analysis as in Theorems 2.1, 2.2, and 2.3, we can find the solutions of generalized fractional kinetic equations involving Mathieu-type series, which are given in the following three examples.

### Example 4.1

If \(\varepsilon , \nu , c,h >{0}\), \(\Re (s)>0\) with \(|s|< c^{-1}\), \(c\neq {h}\) and \(\mu , \rho , \sigma , \tau , l \in {\mathbb{R}^{+}}\), then the solution of the generalized fractional kinetic equation

is given by

### Example 4.2

If \(\nu , c, w >{0}\) and \(\mu , \rho , \sigma , \tau , l \in {\mathbb{R} ^{+}}\); \(|s|< c^{-1}\), \(c\neq {w}\), then the solution of the generalized fractional kinetic equation

is given by

### Example 4.3

If \(\nu >0\), \(c>{0}\), \(\mu , \rho , \sigma , \tau , l \in {\mathbb{R} ^{+}}\) and \(|s|< {c^{-1}}\), then the solution of the generalized fractional kinetic equation

is given by

## Conclusion

In the current study, we have obtained the solution of generalized fractional kinetic equations involving generalized Mathieu-type series with the help of Laplace and Sumudu transform techniques. The results of Sect. 2 are general in character and are likely to find certain applications in the theory of fractional calculus and special functions. From the graphical analysis, we conclude that \(N(t) > 0\) or \(N(t) < 0\) for distinct values of the parameters. By suitably specializing the values of the parameters of generalized Mathieu-type series (25), our main results can yield several new solutions of generalized fractional kinetic equations corresponding to the Mathieu series given by various authors [6, 21, 34, 37]. Therefore, the investigated results in this paper would at once give many results involving diverse special functions occurring in the problems of astrophysics, mathematical physics, and engineering.

## References

- 1.
Agarwal, P., Chand, M., Baleanu, D., O’Regan, D., Jain, S.: On the solution of certain fractional kinetic equations involving

*k*-Mittag-Leffler function. Adv. Differ. Equ.**2018**, 249 (2018) https://doi.org/10.1186/s13662-018-1694-8 - 2.
Agarwal, P., Chand, M., Singh, G.: Fractional kinetic equations involving generalized

*k*-Bessel function via Sumudu transform. Alex. Eng. J.**55**(4), 3053–3059 (2016) - 3.
Asiru, M.A.: Sumudu transform and the solution of integral equations of convolution type. Int. J. Math. Educ. Sci. Technol.

**32**(6), 906–910 (2001) - 4.
Belgacem, F.B.M., Karaballi, A.A.: Sumudu transform fundamental properties investigations and applications. J. Appl. Math. Stoch. Anal.

**2006**, Article ID 91083 (2006) https://doi.org/10.1155/JAMSA/2006/91083 - 5.
Belgacem, F.B.M., Karaballi, A.A., Kalla, S.L.: Analytical investigations of the Sumudu transform and applications to integral production equations. Math. Probl. Eng.

**3**, 103–118 (2003) - 6.
Cerone, P., Lenard, C.T.: On integral forms of generalized Mathieu series. JIPAM. J. Inequal. Pure Appl. Math.

**4**, Article ID 100 (2003) - 7.
Chaurasia, V.B.L., Kumar, D.: On the solution of generalized fractional kinetic equation. Adv. Stud. Theor. Phys.

**4**, 773–780 (2010) - 8.
Chaurasia, V.B.L., Pandey, S.C.: On the new computable solution of the generalized fractional kinetic equations involving the generalized function for the fractional calculus and related functions. Astrophys. Space Sci.

**317**, 213–219 (2008) - 9.
Choi, J., Kumar, D.: Solutions of generalized fractional kinetic equations involving Aleph functions. Math. Commun.

**20**, 113–123 (2015) - 10.
Dorrego, G., Kumar, D.: A generalization of the kinetic equation using the Prabhakar-type operators. Honam Math. J.

**39**(3), 401–416 (1017) - 11.
Dutta, B.K., Arora, L.K., Borah, J.: On the solution of fractional kinetic equation. Gen. Math. Notes

**6**, 40–48 (2011) - 12.
Emersleben, O.: Über die Reihe. Math. Ann.

**125**, 165–171 (1952) - 13.
Gupta, V.G., Sharma, B., Belgacem, F.B.M.: On the solutions of generalized fractional kinetic equations. Appl. Math. Sci.

**5**(17–20), 899–910 (2011) - 14.
Haubold, H.J., Mathai, A.M.: The fractional kinetic equation and thermonuclear functions. Astrophys. Space Sci.

**273**, 53–63 (2000) - 15.
Kamarujjama, M., Khan, N.U., Khan, O.: The generalized

*p*-*k*-Mittag-Leffler function and solution of fractional kinetic equations. J. Anal. in press. https://doi.org/10.1007/s41478-018-0160-z - 16.
Kamarujjama, M., Khan, O.: Computation of new class of integrals involving generalized Galue type Struve function. J. Comput. Appl. Math.

**351**, 228–236 (2019) - 17.
Kilicman, A., Eltayeb, H.: Some remarks on the Sumudu and Laplace transforms and applications to differential equations. ISRN Appl. Math.

**2012**, Article ID 591517 (2012) https://doi.org/10.5402/2012/591517 - 18.
Kilicman, A., Eltayeb, H., Atan, K.A.M.: A note on the comparison between Laplace and Sumudu transforms. Bull. Iran. Math. Soc.

**37**(1), 131–141 (2011) - 19.
Kilicman, A., Gadain, H.E.: On the applications of Laplace and Sumudu transforms. J. Franklin Inst.

**347**, 848–862 (2010) - 20.
Kumar, D., Choi, J., Srivastava, H.M.: Solution of a general family of fractional kinetic equations associated with the generalized Mittag-Leffler function. Nonlinear Funct. Anal. Appl.

**23**(3), 455–471 (2018) - 21.
Mathieu, E.L.: Traité de Physique Mathematique, VI-VII, Théorie de l’élasticité des Corps Solids. Gauthier-Villars, Paris (1980)

- 22.
Nisar, K.S., Qi, F.: On salutation of fractional kinetic equations involving the generalized

*k*-Bessel function. Note Mat.**37**(2), 11–20 (2017) - 23.
Prabhakar, T.R.: A singular integral equation with a generalized Mittag-Leffler function in the kernel. Yokohama Math. J.

**19**, 7–15 (1971) - 24.
Saichev, A.I., Zaslavsky, G.M.: Fractional kinetic equations: solutions and applications. Chaos

**7**(4), 753–764 (1997) - 25.
Saxena, R.K., Kalla, S.L.: On the solutions of certain fractional kinetic equations. Appl. Math. Comput.

**199**, 504–511 (2008) - 26.
Saxena, R.K., Mathai, A.M., Haubold, H.J.: On fractional kinetic equations. Astrophys. Space Sci.

**282**, 281–287 (2002) - 27.
Saxena, R.K., Mathai, A.M., Haubold, H.J.: Unified fractional kinetic equations and a fractional diffusion equation. Astrophys. Space Sci.

**290**, 299–310 (2004) - 28.
Saxena, R.K., Mathai, A.M., Haubold, H.J.: On generalized fractional kinetic equations. Physica A

**344**, 657–664 (2004) - 29.
Saxena, R.K., Ram, J., Kumar, D.: Alternative derivation of generalized fractional kinetic equations. J. Fract. Calc. Appl.

**4**(2), 322–334 (2013) - 30.
Schiff, J.L.: The Laplace Transform: Theory and Applications. Springer, Berlin (1999)

- 31.
Shaktawat, B.S., Gupta, R.K., Kumar, D.: Generalized fractional kinetic equations and its solutions involving generalized Mittag-Leffler function. J. Rajasthan Acad. Phys. Sci.

**16**, 63–74 (2017) - 32.
Spiegel, M.R.: Theory and Problems of Laplace Transforms. Schaums Outline Series. McGraw-Hill, New York (1965)

- 33.
Srivastava, H.M., Saxena, R.K.: Operators of fractional integration and their applications. Appl. Math. Comput.

**118**, 1–52 (2001) - 34.
Srivastava, H.M., Tomovski, Z.: Some problems and solutions involving Mathieu’ series and its generalizations. J. Inequal. Pure Appl. Math.

**5**(2), Article ID 45 (2004) - 35.
Tomovski, Z., Mehrez, M.: Some families of generalized Mathieu-type power series. Appl. Anal. Discrete Math.

**20**(4), 973–986 (2017) - 36.
Tomovski, Z., Pogány, T.K.: Integral expressions for Mathieu-type power series and for the Butzer–Flocke–Hauss function. Fract. Calc. Appl. Anal.

**14**(4), 623–634 (2011) - 37.
Tomovski, Z., Trencevski, K.: On an open problem of Bai-Ni Guo and Feng Qi. J. Inequal. Pure Appl. Math.

**4**(2), Article ID 29 (2003) - 38.
Watugala, G.K.: Sumudu transform: a new integral transform to solve differential equations and control engineering problems. Int. J. Math. Educ. Sci. Technol.

**24**(1), 35–43 (1998) - 39.
Watugala, G.K.: The Sumudu transform for function of two variables. Math. Eng. Ind.

**8**, 293–302 (2002) - 40.
Wiman, A.: Über den fundamental Satz in der Theorie der Funktionen \(E_{\alpha }(x)\). Acta Math.

**29**, 191–201 (1995)

### Acknowledgements

Not applicable.

### Availability of data and materials

Not applicable.

## Funding

Not applicable.

## Author information

### Affiliations

### Contributions

All authors contributed equally. All authors read and approved the final manuscript.

### Corresponding author

## Ethics declarations

### Competing interests

The authors declare that they have no competing interest.

## Additional information

### Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Rights and permissions

**Open Access** This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

## About this article

### Cite this article

Khan, O., Khan, N., Baleanu, D. *et al.* Computable solution of fractional kinetic equations using Mathieu-type series.
*Adv Differ Equ* **2019, **234 (2019). https://doi.org/10.1186/s13662-019-2167-4

Received:

Accepted:

Published:

### MSC

- 26A33
- 33A10
- 33E20
- 44A10

### Keywords

- Generalized fractional kinetic equation
- Mathieu-type series
- Laplace transform
- Sumudu transform