#### Abstract

In investigation of boundary-value problems for certain partial differential equations arising in applied mathematics, we often need to study the solution of system of partial differential equations satisfied by hypergeometric functions and find explicit linearly independent solutions for the system. Here we choose the Exton function among his 21 functions to show how to find the linearly independent solutions of partial differential equations satisfied by this function . Based upon the classical derivative and integral operators, we introduce a new operational images for hypergeometric function . By means of these operational images, a number of finite series and decomposition formulas are then found.

#### 1. Introduction and Preliminaries

Solutions of many applied problems involving thermal conductivity and dynamics, electromagnetic oscillation and aerodynamics, and quantum mechanics and potential theory are obtainable with the help of hypergeometric (higher and special or transcendent) functions [1–4]. Such kinds of functions are often referred to as special functions of mathematical physics. They mainly appear in the solution of partial differential equations which are dealt with harmonic analysis method (see [5]). In view of various applications, it is interesting in itself and seems to be very important to conduct a continuous research of multiple hypergeometric functions. For instance, in [6], a comprehensive list of hypergeometric functions of three variables as many as 205 is recorded, together with their regions of convergence. It is noted that Riemann functions and the fundamental solutions of the degenerate second-order partial differential equations are expressible by means of hypergeometric functions of several variables (see [7–19]). Therefore, in investigation of boundary-value problems for these partial differential equations, we need to study the solution of the system of hypergeometric functions and find explicit linearly independent solutions (see [12–17, 19]). Exton [20, pages 78-79] introduced 21 complete hypergeometric functions of four variables. In [21] Sharma and Parihar introduced 83 complete hypergeometric functions of four variables. It is remarked that, out of these 83 functions, the following 19 functions had already appeared in the work of Exton [20] in the different notations: Each quadruple hypergeometric function is of the form where is a certain sequence of complex parameters and there are twelve parameters in each function . Here, for example, we choose the Exton function among his twenty-one functions to find the linearly independent solutions of partial differential equations satisfied by this function .

#### 2. The System of Partial Differential Equations for

According to the theory of multiple hypergeometric functions (see [22]), the system of partial differential equations for Exton hypergeometric function is readily seen to be given as follows: where .

Now by making use of some elementary calculations, we find the following system of second-order partial differential equations: It is noted that three equations of the system (5) are linearly dependent, because the hypergeometric function satisfies the system. Now, in order to find the linearly independent solutions of the system (5), we consider as in the form an unknown function, and , , , and are constants which are to be determined. So, substituting into the system (5), we obtain where It is noted that the system (6) is analogical to the system (5). Therefore, it is required that the following conditions should be satisfied: It is not difficult to see that the system (8) satisfies the following solutions:

Finally, substituting all eight solutions (9) into (6), we find the following linearly independent solutions of the system (5): Moreover, it is seen that the global solution of the system (5) is combined to be in the form where are constants.

#### 3. Operational Representations for and Applications

In this section, we apply the concept of the right Riemann-Liouville fractional derivative to obtain operational images for hypergeometric functions . Indeed, in this section we build up the right Riemann-Liouville fractional derivative operator, which plays the role of augmenting parameters in the hypergeometric functions involved (see [23, Chapters 4 and 5]). In particular, we will deal with operational definitions ruled by the operators and where denotes the derivative operator and defines the inverse of the derivative and once acting on unity, it yields It is evident that is essentially an integral operator and the lower limit has been assumed to be zero. The following two formulas are well-known consequences of the derivative operator and the integral operator (see [24]). Consider where . Based on the operational relations (13), we have In view of the definition of the Gaussian hypergeometric function [6], Bin-Saad [25] introduced the following operational representations for Gaussian hypergeometric .

Lemma 1. *Let and ; then
**
For the quadruple hypergeometric function , one introduces the following operational representations.*

Theorem 2. *Let , , , , , , and ; then
*

*Proof. *Denote, for convenience, the left-hand side of assertion (17) by . Then, as a consequence of the binomial theorem, it is easily seen that
Upon using (14) and considering definition (3), we are led finally to right-hand side of assertion (17). The proof of assertion (18) runs parallel to that of assertion (17); then we skip the details.

One of the advantages offered by the use of the operational images (17) and (18) is the possibility of establishing finite sums for the quadruple hypergeometric series . First, let, for convenience, By employing (18) and taking advantage from the identity and the definition of the general triple hypergeometric series (see [6]), we can state that Similarly, in view of (17), the definition of Appell hypergeometric series (see [6, page 23, equation (5)]) and the use of the following identity allow to conclude that Also, by virtue of the identity we find the following decomposition formula: where is Lauricella function of four variables [6, page 33, equation (3)].

Moreover, in view of the composition rule and the operational image (18), we easily find for the function that where

Finally, let us stress that the schema suggested in Sections 2 and 3 can be applied to find linear independent solutions, symbolic operational images, finite sums, and decomposition formulas for other quadruple hypergeometric functions.

#### Conflict of Interests

The authors declare that they have no competing interests.