Certain Fractional Integral Formulas Involving the Product of Generalized Bessel Functions

We apply generalized operators of fractional integration involving Appell's function F 3(·) 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.


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][19][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][24][25][26][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][31][32][33][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 2 The Scientific World Journal 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 > 0 and , , , , ∈ C by ( , , , , 0,+ These operators (integral transforms) were introduced by Marichev [35] as Mellin type convolution operators with a special function 3 (⋅) 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 3 (⋅) 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 2 1 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 ∈ C\{0} and , , ] ∈ C with R(]) > −1 by the following series (see, e.g., [30, page 10, equation (1.15)]; for a very recent work, see also [31][32][33], [ where C denotes the set of complex numbers and Γ( ) is the familiar Gamma function (see [ where The coefficients ( = 1, . . and Z − 0 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 Scientific World Journal 3 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.

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
The Marichev-Saigo-Maeda fractional integration (1) of product of generalized Bessel functions (5) and Binomial function is given by the following result.
On setting = 0 in Theorem 2, we get the following Saigo hypergeometric fractional image of the product of generalized Bessel function of the first formula kind ] ( ). , , , , , , ] , , , ∈ C and > 0,
Again, on letting = 0 and = 0 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.

Right-Side Fractional Integration of Generalized Bessel Functions
In the sequel, we use the following results.
The Marichev-Saigo-Maeda fractional integration (2) of product of generalized Bessel functions (5) and Binomial function is given by the following result.
Proof. For convenience, let the left-hand side of the (31) be denoted by J. Applying definition (5) and the binomial 6 The Scientific World Journal 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.
Again, on setting = 0 and = 0 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.

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 ( ≥ 1) the function √ /2 ] ( ) is well approximated by the trigonometric function cos( − ]/2− /4). 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 = 1 and = 1, the Generalized Bessel function of the first kind (5) reduces to ] ( ), and when = −1 and = 1 the function ] ( ) becomes ] ( ). Similarly, when = 1 and = 2, the function ] ( ) reduces to 2 ] /√ , while if = −1 and = 2, then ] ( ) becomes 2 ] /√ . In the sequel from (5), we have ] (0) = 0. Therefore, the results presented in this paper are easily converted in terms of the various special Bessel functions after some suitable parametric replacement.
Further, it is interesting to observe that, if we set at = 0, 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 = 0 and = = 1, then our main results (11) and (31) follow the kown results due to [28]. Finally, if we set = 0, = = = 1 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, radiophysics, 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.

8
The Scientific World Journal