#### Abstract

We apply generalized operators of fractional integration involving Appell’s function due to Marichev-Saigo-Maeda, to the product of the generalized Bessel function of the first kind due to Baricz. The results are expressed in terms of the multivariable generalized Lauricella functions. Corresponding assertions in terms of Saigo, Erdélyi-Kober, Riemann-Liouville, and Weyl type of fractional integrals are also presented. Some interesting special cases of our two main results are presented. We also point out that the results presented here, being of general character, are easily reducible to yield many diverse new and known integral formulas involving simpler functions.

#### 1. Introduction, Definitions, and Preliminaries

The fractional calculus is nowadays one of the most rapidly growing subject of mathematical analysis. It is a field of applied mathematics that deals with derivatives and integrals of arbitrary orders. The fractional integral operator involving various special functions has found significant importance and applications in various subfields of applicable mathematical analysis. Many applications of fractional calculus can be found in turbulence and fluid dynamics, stochastic dynamical system, plasma physics and controlled thermonuclear fusion, nonlinear control theory, image processing, nonlinear biological systems, astrophysics, and in quantum mechanics. Since last four decades, a number of workers like Love [1], McBride [2], Kalla [3, 4], Kalla and Saxena [5], Saigo [6, 7], Kilbas [8], Kilbas and Sebastian [9], Kiryakova [10, 11], Baleanu [12], Baleanu and Mustafa [13], Baleanu et al. [14], Baleanu et al. [15], Purohit and Kalla [16], Purohit [17], Agarwal [18–20] and Agarwal and Jain [21], and so forth have studied, in depth, the properties, applications, and different extensions of various operators of fractional calculus. A detailed account of fractional calculus operators along with their properties and applications can be found in the research monographs by Miller and Ross [22], Kiryakova [10], and so forth.

The computation of fractional derivatives (and fractional integrals) of special functions of one and more variables is important from the point of view of the usefulness of these results in the evaluation of generalized integrals and the solution of differential and integral equations. Motivated by these avenues of applications, a remarkably large number of fractional integral formulas involving a variety of special functions have been developed by many authors (see, e.g., [23–27]). Fractional integration formulae for the Bessel function and generalized Bessel functions are given recently by Kilbas and Sebastian [9], Saxena et al. [28], Malik et al. [29] and Purohit et al. [27]. A useful generalization of the Bessel function has been introduced and studied in [30–34]. Here we aim at presenting composition formula of Marichev-Saigo-Maeda fractional integral operators and the product of generalized Bessel function, which are expressed in terms of the multivariable generalized Lauricella functions. Some interesting special cases of our main results are also considered.

For our purpose, we begin by recalling some known functions and earlier works. This paper deals with two integral transforms defined for and by These operators (integral transforms) were introduced by Marichev [35] as Mellin type convolution operators with a special function in the kernel. These operators were rediscovered and studied by Saigo in [36] as generalization of the so-called Saigo fractional integral operators (see [37]). Such further properties as, for example, their relations with the Mellin transform and with the hypergeometric operators (or the Saigo fractional integral operators), together with their decompositional, operational and other properties in the McBride space (see [38]) were studied by Saigo and Maeda [39].

In (1) and (2) the symbol denotes so-called 3rd Appell function (known also as Horn function) (see [40, page 413]): The properties of this function are discussed in [40, pages 412–415]. In particular, its relation to the Gauss hypergeometric function is presented as follows: It is known that the 3rd Appell function cannot be expressed as a product of two functions that satisfy pairs of linear partial differential equations of the second order.

We investigate compositions of integral transforms (1) and (2) with the product of generalized Bessel function of the first kind, , which is defined for and with by the following series (see, e.g., [30, page 10, equation (1.15)]; for a very recent work, see also [31–33], [34, page 182, equation (2.2)], and [29, page 2, equation (8)]): where denotes the set of complex numbers and is the familiar Gamma function (see [41, Section 1.1]).

Then we show that the composition is expressed in terms of the multivariable generalized Lauricella functions due to Srivastava and Daoust (see, [42, page 454]) is a generalization of the Wright function in several variables and defined by (see e.g., [43, page 107]) where

The coefficients (), (), (), and (), for all , are real and positive and means the array of -parameters ; with similar interpretations for , , , and so forth, and where is the Pochhammer symbol defined (for ) by (see [41, page 2 and pages 4–6]) and denotes the set of nonpositive integers.

The paper is organized as follows. The composition formula of Marichev-Saigo-Maeda fractional integral operators (1) and (2) with the product of generalized Bessel function of the first kind (5) is proved in terms of the multivariable generalized Lauricella functions (6) in Sections 2 and 3, respectively. The corresponding results for the corresponding assertions in terms of Saigo, Erdélyi-Kober, Riemann-Liouville and Weyl type of fractional integrals are also presented in Sections 2 and 3. Special cases giving compositions of fractional integrals with the product of trigonometric functions and concluding remarks are considered in Section 4.

#### 2. Left-Side Fractional Integration of Generalized Bessel Functions

Our results in this Section are based on the preliminary assertions giving composition formula of fractional integral (1) with a power function.

Lemma 1 (Saigo and Maeda [39, Lemma 1]). *Let and if , , then
*

The Marichev-Saigo-Maeda fractional integration (1) of product of generalized Bessel functions (5) and Binomial function is given by the following result.

Theorem 2. *Let and satisfying the inequalities
**
then for , , one has
**
where , , , , , , , and are given by the following
*

*Proof. *For sake of convenience, let the left-hand side of the (11) be denoted by . Using definition (5) and the binomial expansion, namely,
we find

Following the convergence condition of Theorem 2, for any ,
and then changing the order of integration and summation, we obtain
Now on applying Lemma 1 and using (9) with replaced by , we obtain

This, in accordance with (6), gives the required result (11). This completed the proof of the Theorem 2.

On setting in Theorem 2, we get the following Saigo hypergeometric fractional image of the product of generalized Bessel function of the first formula kind .

Corollary 3. *If and , **
then there hold the results:
**
where , , , , , and are given by (12), (14), (15), (17), (18), and (19), respectively.*

Again, on letting and in Theorem 2, we get the Riemann-Liouville fractional image of the product of generalized Bessel function of the first formula kind; asserted by the following corollary.

Corollary 4. *Let and , , and
**
Then, one has
**
where , , , and are given by (12), (15), (18), and (19), respectively.*

#### 3. Right-Side Fractional Integration of Generalized Bessel Functions

In the sequel, we use the following results.

Lemma 5 (Saigo and Maeda [39, Lemma 2]). *Let and if , , then
*

The Marichev-Saigo-Maeda fractional integration (2) of product of generalized Bessel functions (5) and Binomial function is given by the following result.

Theorem 6. *Let and , such that
**
Then, one has
**
where , , , , , and are given by the following:
**
and are given by (18) and (19), respectively.*

*Proof. *For convenience, let the left-hand side of the (31) be denoted by . Applying definition (5) and the binomial expansion for and then changing the order of integration and summation, we find
Following the convergence condition of Theorem 6, for any

Now, on making use of Lemma 5 and (9), we obtain

This, in accordance with (6), gives the required result (31). This completed the proof of Theorem 6.

On taking in Theorem 6, we get the right-sided Saigo hypergeometric fractional image of the product of generalized Bessel function of the first formula kind, as follows.

Corollary 7. *Let and , , and
**
Then
**
where , , , , and are given by (32), (33), (34), (35), (18), and (19), respectively.*

Again, on setting and in Theorem 6, we get the following image formula for the product of generalized Bessel function of the first kind involving Riemann-Liouville type fractional integral operator.

Corollary 8. *Let and , , and
**
Then
**
where , , , and are given by (33), (34), (18), and (19), respectively.*

#### 4. Consequence Results and Concluding Remarks

In this section, we briefly consider another variation of the results derived in the preceding sections. Bessel functions are important special functions that appear widely in science and engineering. Bessel functions of the first kind are oscillatory and may be regarded as generalizations of trigonometric functions. Indeed, for large argument the function is well approximated by the trigonometric function . Similarly, modified Bessel functions of the first kind , which are Bessel functions of imaginary argument, may be regarded as generalization of exponentials. Exponential functions have the unique and special property that they are particularly easy to multiply and to raise to powers: and . Further, it can be easily seen that for and , the Generalized Bessel function of the first kind (5) reduces to , and when and the function becomes . Similarly, when and , the function reduces to , while if and , then becomes . In the sequel from (5), we have . Therefore, the results presented in this paper are easily converted in terms of the various special Bessel functions after some suitable parametric replacement.

For example, if we set , then the generalized Bessel function in (5) have following relation with cosine function when is replaced by (see, e.g., [29]):

While if and , in (5) have relation with hyperbolic cosine function as (see, e.g., [29]) In the sequel, if we set , then the generalized Bessel function in (5) have following relation with sine function when is replaced by (see, e.g., [29]):

Further, if we put and in (5), we have the following relation with and hyperbolic sine function as (see, e.g., [29])

Now, by virtue of the relations (45) to (48), one can easily derive the Marichev-Saigo-Maeda type fractional integrals involving product of trigonometric functions. Therefore, we omit the details of these results.

Further, it is interesting to observe that, if we set at , the results given by Malik et al. [29, pages 5 and 8, Theorems 4.1 and 5.1] follow as special cases of our main results. Again, if we set and , then our main results (11) and (31) follow the kown results due to [28]. Finally, if we set , and give some suitable parametric replacement in the main results (11) and (31), we can arrive at the known results due to Purohit et al. [27] in a slightly different notation.

Fractional integral formulas involving products of Bessel functions have been developed and play an important role in several physical problems. In fact, Bessel functions are playing the important role in studying solutions of differential equations, integral equations [44], and they are associated with a wide range of problems in important areas of mathematical physics, like problems of acoustics, radio-physics, hydrodynamics, and atomic and nuclear physics. These considerations have led various workers in the field of special functions for exploring the possible extensions and applications for the Bessel functions. Among many properties of Bessel functions, they also have investigated some possible extensions of the Bessel functions. The generalized Bessel function defined by (5), possess the advantage that a number of Bessel functions, trigonometric functions and hyperbolic functions happen to be the particular cases of this function. Therefore, we conclude this paper with the remark that, the results deduced above are significant and can lead to yield numerous other fractional integrals involving various Bessel functions and trigonometric functions by the suitable specializations of arbitrary parameters in the theorems. More importantly, they are expected to find some applications to the solutions of fractional differential and integral equations. The results thus derived in this paper are general in character and likely to find certain applications in the theory of special functions.

#### Acknowledgments

The authors take this opportunity to express their deepest thanks to the members of the editorial and reviewer board for their valuable comments and essential suggestions to improve this paper as in the present form.