Abstract

The aim of this research paper is to obtain explicit expressions of (1𝑥)22𝑎𝐹1[𝑎,𝑏;4𝑥/(1𝑥)22𝑏+𝑗]for 𝑗=0,±1,±2. For 𝑗=0, we have the well-known transformation formula due to Gauss. The results are derived with the help of generalized Watson's theorem. Some known results obtained earlier follow special cases of our main findings.

1. Introduction

The generalized hypergeometric function with 𝑝 numerator and 𝑞 denominator parameters is defined by [1, page 73, equation (2)] 𝑝𝐹𝑞𝛼1,,𝛼𝑝𝛽;𝑧1,,𝛽𝑞,=𝑝𝐹𝑞𝛼1,,𝛼𝑝;𝛽1,,𝛽𝑞=;𝑧𝑛=0𝛼1𝑛𝛼𝑝𝑛𝛽1𝑛𝛽𝑞𝑛𝑧𝑛,𝑛!(1.1) where (𝛼)𝑛 denotes the Pochhammer symbol (or the shifted factorial, since (1)𝑛=𝑛!) defined, for any complex number 𝛼, by (𝛼)𝑛=𝛼(𝛼+1)(𝛼+𝑛1),𝑛,1,𝑛=0.(1.2) Using the fundamental property Γ(𝛼+1)=𝛼Γ(𝛼), (𝛼)𝑛 can be written in the form (𝛼)𝑛=Γ(𝛼+𝑛),Γ(𝛼)(1.3) where Γ() is the well known Gamma function.

The special case of (1.1) for 𝑝=2 and 𝑞=1, namely 2𝐹1𝑐𝑎,𝑏;𝑧=1+𝑎𝑏1𝑐𝑧+𝑎(𝑎+1)𝑏(𝑏+1)𝑧12𝑐(𝑐+1)2=+𝑛=0(𝑎)𝑛(𝑏)𝑛(𝑐)𝑛𝑧𝑛,𝑛!(1.4) was systematically studied by Gauss [2] in 1812.

The series (1.1) is of great importance to mathematicians and physicists. All the elements 𝑎, 𝑏, and 𝑐 (similarly for (1.1)) in (1.4) are called the parameters of the series and 𝑧 is called the variable of the series. All four quantities 𝑎, 𝑏, 𝑐, and 𝑧 may be real or complex with one exception that the denominator parameter 𝑐 should not be zero or a negative integer. Also it can easily been seen that if any one of the numerator parameters 𝑎 or 𝑏 or both is a negative integer, the series terminates that is, reduces to a polynomial.

The series (1.4) is known as Gauss series or the ordinary hypergeometric series and may be regarded as a generalization of the elementary geometric series. In fact (1.4) reduces to the elementary geometric series in two cases, when 𝑎=𝑐 and 𝑏=1 and also when 𝑏=𝑐 and 𝑎=1.

For convergence (including absolute convergence) we refer the reader to the standard texts [3] and [1].

It is interesting to mention here that in (1.4), if we replace 𝑧 by 𝑧/𝑏 and let 𝑏, then since ((𝑏)𝑛𝑧𝑛)/𝑏𝑛𝑧𝑛 we arrive at the following series: 1𝐹1𝑎𝑐𝑎;𝑧=1+1𝑐𝑧+𝑎(𝑎+1)𝑧12𝑐(𝑐+1)2=+𝑛=0(𝑎)𝑛(𝑐)𝑛𝑧𝑛,𝑛!(1.5) which is called the Kummer’s series or the confluent hypergeometric series.

Gauss’s hypergeometric function 2𝐹1 and its confluent case 1𝐹1 form the core of the special functions and include, as their special cases, most of the commonly elementary functions.

It should be remarked here that whenever hypergeometric and generalized hypergeometric functions reduce to gamma functions, the results are very important from an application point of view. Only a few summation theorems for the series 2𝐹1 and 3𝐹2 are available.

In this context, it is well known that the classical summation theorems such as of Gauss, Gauss second, Kummer and Bailey for the series 2𝐹1; Watson, Dixon, Whipple, and Saalschütz for the series 3𝐹2 play an important rule in the theory of hypergeometric and generalized series.

Several formulae were given by Gauss [2] and Kummer [4] expressing the product of the hypergeometric series as a hypergeometric series, such as 𝑒1𝑥𝐹1(𝑥) as a series of the type 1𝐹1(𝑥) and (1+𝑥)2𝑝𝐹1[4𝑥/(1+𝑥)2] as a series of the type 2𝐹1(𝑥). In 1927, Whipple [5] has obtained a formula expressing (1𝑥)3𝑝𝐹2[4𝑥/(1𝑥)2] as a series of the type 3𝐹2.

By employing the above mentioned classical summation theorems for the series 2𝐹1 and 3𝐹2, Bailey [6] in his well known, interesting and popular research paper made a systematic study and obtained a large number of such formulas.

Gauss [2] obtained the following quadratic transformation formula, namely (1𝑥)22𝑎𝐹1𝑎,𝑏;4𝑥(1𝑥)2=2𝑏2𝐹11𝑎,𝑎𝑏+2;𝑥21𝑏+2,(1.6) which is also contained in [7, entry (8.1.1.41), page 573].

Berndt [8] pointed out that the result (1.6) is precisely (5) of Erdèlyi treatise [9, page 111], and is the Entry 3 of the Chapter 11 of Ramanujan’s Notebooks [8, page 50] (of course, by replacing 𝑥 by 𝑥).

Bailey [6] established the result (1.6) with the help of the following classical Watson’s summation theorem [3], namely 3𝐹21𝑎,𝑏,𝑐;12=(𝑎+𝑏+1),2𝑐Γ(1/2)Γ(𝑐+1/2)Γ(𝑐(1/2)𝑎(1/2)𝑏+1/2)Γ((1/2)𝑎+(1/2)𝑏+1/2),Γ((1/2)𝑎+1/2)Γ((1/2)𝑏+1/2)Γ(𝑐(1/2)𝑎+1/2)Γ(𝑐(1/2)𝑏+1/2)(1.7) provided that Re(2𝑐𝑎𝑏)>1.

The proof of (1.7) when one of the parameters 𝑎 or 𝑏 is a negative integer was given in Watson [10]. Subsequently, it was established more generally in the nonterminating case by Whipple [5]. The standard proof of the nonterminating case was given in Bailey’s tract [3] by employing the fundamental transformation due to Thomae combined with the classical Dixon’s theorem of the sum of a 3𝐹2. For a very recent proof of (1.7), see [11].

It is not out of place to mention here that in (1.6), if we replace 𝑥 by 𝑥/𝑎 and let 𝑎, then after a little simplification, we get the following well-known Kummer’s second theorem [4, page 140] [12, page 132], namely 𝑒1𝑥/2𝐹1𝑏=;𝑥2𝑏0𝐹1;𝑥2116𝑏+2,(1.8) which also appeared as Entry 7 of the chapter 11 of Ramanujan’s Notebooks [8, page 50] (of course, by replacing 𝑥 by 𝑥/2).

Very recently, Kim et al. [13] have obtained sixty six results closely related to (1.8) out of which four results are given here. These are 𝑒1𝑥/2𝐹1(𝛼;2𝛼+1;𝑥)=0𝐹11;𝛼+2;𝑥2𝑥162(2𝛼+1)0𝐹13;𝛼+2;𝑥2𝑒16,(1.9)1𝑥/2𝐹1(𝛼;2𝛼1;𝑥)=0𝐹11;𝛼2;𝑥2+𝑥162(2𝛼1)0𝐹11;𝛼+2;𝑥2𝑒16,(1.10)1𝑥/2𝐹1(𝛼;2𝛼+2;𝑥)=0𝐹13;𝛼+2;𝑥2𝑥162(𝛼+1)0𝐹13;𝛼+2;𝑥2+𝑥1624(𝛼+1)(2𝛼+3)0𝐹15;𝛼+2;𝑥2,𝑒16(1.11)1𝑥/2𝐹1(𝛼;2𝛼2;𝑥)=0𝐹11;𝛼2;𝑥2+𝑥162(𝛼1)0𝐹13;𝛼2;𝑥2+𝑥1624(𝛼1)(2𝛼1)0𝐹11;𝛼+2;𝑥2.16(1.12)

We remark in passing that the results (1.9) and (1.10) are also recorded in [14].

Recently, a good progress has been made in generalizing the classical Watson’s theorem (1.7) on the sum of a 3𝐹2. In 1992, Lavoie et al. [15] have obtained explicit expressions of 3𝐹21𝑎,𝑏,𝑐;12(𝑎+𝑏+𝑖+1),2𝑐+𝑗for𝑖,𝑗=0,±1,±2.(1.13)

For 𝑖=𝑗=0, we get Watson’s theorem (1.7). In the same paper [15], they have also obtained a large number of very interesting limiting and special cases of their main findings.

In [16], a summation formula for (1.7) with fixed 𝑗 and arbitrary 𝑖(𝑖,𝑗) was given. This result generalizes the classical Watson’s summation theorem with the case 𝑖=𝑗=0.

For the a recent generalization of Watson’s summation theorems and other classical summation theorems for the series 2𝐹1 and 3𝐹2 in the most general case, see [17].

The aim of this research paper is to obtain the explicit expressions of (1𝑥)22𝑎𝐹1𝑎,𝑏;4𝑥(1𝑥)22𝑏+𝑗for𝑗=0,±1,±2.(1.14)

In order to derive our main results, we shall require the following.(1)The following special cases of (1.13) for 𝑖=0, recorded in [15]: 3𝐹212𝑛,2𝑎+2𝑛,𝑐;1𝑎+2,2𝑐+𝑗=𝐷𝑗(1/2)𝑛([])𝑎𝑐+(1/2)𝑗/2𝑛(𝑎+(1/2))𝑛[])(𝑐+(1/2)+𝑗/2𝑛,(1.15)3𝐹212𝑛1,2𝑎+2𝑛+1,𝑐;1𝑎+2,2𝑐+𝑗=𝐸𝑗(3/2)𝑛[])(𝑎𝑐+(3/2)(𝑗+1)/2𝑛(𝑎+(1/2))𝑛[])(𝑐+(1/2)+(𝑗+1)/2𝑛(1.16) each for 𝑗=0, ±1, ±2. Also, as usual, [𝑥] denotes the greatest integer less than or equal to 𝑥, and its modulus is defined by |𝑥|. The coefficients 𝐷𝑗 and 𝐸𝑗 are given in Table 1.(2)The known identities [1, page 22, lemma 5; page 58, equation 1; page 52, equation 2; page 58, equation 3] (𝜆)2𝑛=22𝑛12𝜆𝑛121𝜆+2𝑛((𝑛{0}),(1.17)1𝑥)𝑎=𝑛=0(𝑎)𝑛𝑥𝑛!𝑛,(1.18)(𝜆)𝑛𝑘=(1)𝑘(𝜆)𝑛(1𝜆𝑛)𝑘(0𝑘𝑛;𝑛{0}),(1.19)(𝑛𝑘)!=(1)𝑘𝑛!(𝑛)𝑘(0𝑘𝑛;𝑛{0}).(1.20)

2. Main Transformation Formulae

The generalization of the quadratic transformation (1.6) due to Gauss to be established is (1𝑥)22𝑎𝐹1𝑎,𝑏;4𝑥(1𝑥)2=2𝑏+𝑗𝑛=0𝐷𝑗(𝑎)𝑛[])(𝑎𝑏+(1/2)𝑗/2𝑛[])(𝑏+(1/2)+𝑗/2𝑛𝑥𝑛!2𝑛+2𝑎𝑛=0𝐸𝑗(𝑎+1)𝑛([(])𝑎𝑏+(3/2)𝑗+1)/2𝑛[])(𝑏+(1/2)+(𝑗+1)/2𝑛𝑥𝑛!2𝑛+1for𝑗=0,±1,±2.(2.1)

Also, as usual, [𝑥] represents the greatest integer less than or equal to 𝑥, and its modulus is denoted by |𝑥|. The coefficients 𝐷𝑗 and 𝐸𝑗 are given in Table 1.

2.1. Derivation

In order to derive our main transformation (2.1), we proceed as follows.

Proof. Denoting the left-hand side of (2.1) by 𝑆, we have 𝑆=(1𝑥)22𝑎𝐹1𝑎,𝑏;4𝑥(1𝑥)22𝑏+𝑗.(2.2) Expressing 2𝐹1 as a series and after a little simplification =𝑘=0(𝑎)𝑘(𝑏)𝑘(1)𝑘22𝑘𝑥𝑘𝑘!(2𝑏+𝑗)𝑘(1𝑥)(2𝑎+2𝑘).(2.3) Using Binomial theorem (1.18), we have =𝑘=0(𝑎)𝑘(𝑏)𝑘(1)𝑘22𝑘𝑥𝑘𝑘!(2𝑏+𝑗)𝑘𝑛=0(2𝑎+2𝑘)𝑛𝑥𝑛!𝑛,(2.4) which on simplification gives =𝑘=0𝑛=0(2𝑎+2𝑘)𝑛(𝑎)𝑘(𝑏)𝑘(1)𝑘22𝑘(2𝑏+𝑗)𝑘𝑥𝑛!𝑘!𝑛+𝑘.(2.5) Changing 𝑛 to 𝑛𝑘 and using the result [1, page 57, lemma 11] 𝑛=0𝑘=0𝐴(𝑘,𝑛)=𝑛𝑛=0𝑘=0𝐴(𝑘,𝑛𝑘),(2.6) we have 𝑆=𝑛𝑛=0𝑘=0(2𝑎+2𝑘)𝑛𝑘(𝑎)𝑘(𝑏)𝑘(1)𝑘22𝑘(2𝑏+𝑗)𝑘𝑥(𝑛𝑘)!𝑘!𝑛=𝑛𝑛=0𝑘=0Γ(2𝑎+𝑛+𝑘)Γ(2𝑎+2𝑘)(𝑎)𝑘(𝑏)𝑘(1)𝑘22𝑘(2𝑏+𝑗)𝑘(𝑥𝑛𝑘)!𝑘!𝑛.(2.7) Using (1.20) and after a little algebra =𝑛=0(2𝑎)𝑛𝑥𝑛!𝑛𝑛𝑘=0(𝑛)𝑘(2𝑎+𝑛)𝑘(𝑏)𝑘(𝑎+(1/2))𝑘(2𝑏+𝑗)𝑘.𝑘!(2.8)
Summing up the inner series, we have 𝑆=𝑛=0(2𝑎)𝑛𝑥𝑛!𝑛3𝐹21𝑛,2𝑎+𝑛,𝑏;1𝑎+2,2𝑏+𝑗,(2.9) separating into even and odd powers of 𝑥, we have 𝑆=𝑛=0(2𝑎)2𝑛𝑥(2𝑛)!32𝑛𝐹212𝑛,2𝑎+2𝑛,𝑏;1𝑎+2+,2𝑏+𝑗𝑛=0(2𝑎)2𝑛+1𝑥(2𝑛+1)!32𝑛+1𝐹212𝑛1,2𝑎+2𝑛+1,𝑏;1𝑎+2.,2𝑏+𝑗(2.10)
Finally, using (1.17), (1.15), and (1.16) and after a little algebra, we easily arrive at the right-hand side of (2.1).
This completes the proof of (2.1).

3. Special Cases

In (2.1), if we put 𝑗=0, ±1, ±2, we get, after summing up the series in terms of generalized hypergeometric function, the following interesting results:(i)For 𝑗=0, (1𝑥)22𝑎𝐹1𝑎,𝑏;4𝑥(1𝑥)2=2𝑏2𝐹11𝑎,𝑎𝑏+2;𝑥21𝑏+2.(3.1)(ii)For 𝑗=1, (1𝑥)22𝑎𝐹1𝑎,𝑏;4𝑥(1𝑥)2=2𝑏+12𝐹11𝑎,𝑎𝑏+2;𝑥21𝑏+2+2𝑎𝑥(2𝑏+1)2𝐹11𝑎+1,𝑎𝑏+2;𝑥23𝑏+2.(3.2)(iii)For 𝑗=1, (1𝑥)22𝑎𝐹1𝑎,𝑏;4𝑥(1𝑥)2=2𝑏12𝐹13𝑎,𝑎𝑏+2;𝑥21𝑏22𝑎𝑥(2𝑏1)2𝐹13𝑎+1,𝑎𝑏+2;𝑥21𝑏+2.(3.3)(iv)For 𝑗=2, (1𝑥)22𝑎𝐹1𝑎,𝑏;4𝑥(1𝑥)2=2𝑏+22𝐹11𝑎,𝑎𝑏+2;𝑥23𝑏+2+2𝑎𝑥(𝑏+1)2𝐹11𝑎+1,𝑎𝑏+2;𝑥23𝑏+2+4𝑎(𝑎+1)𝑥(𝑏+1)(2𝑏+3)22𝐹11𝑎+2,𝑎𝑏+2;𝑥25𝑏+2.(3.4)(v)For 𝑗=2, (1𝑥)22𝑎𝐹1𝑎,𝑏;4𝑥(1𝑥)2=2𝑏22𝐹13𝑎,𝑎𝑏+2;𝑥21𝑏22𝑎𝑥(𝑏1)2𝐹15𝑎+1,𝑎𝑏+2;𝑥21𝑏2+4𝑎(𝑎+1)𝑥(𝑏1)(2𝑏1)22𝐹15𝑎+2,𝑎𝑏+2;𝑥21𝑏+2.(3.5)

Clearly, the result (3.1) is the well-known quadratic transformation due to Gauss (1.6) and the results (3.2) to (3.5) are closely related to (3.1).

Remark 3.1. The results (3.2) and (3.3) are also recorded in [18].
In (3.1), (3.2), and (3.4), if we take 𝑏=1/2, we get the following results:(1)For 𝑗=0, (1𝑥)22𝑎𝐹11𝑎,2;4𝑥(1𝑥)21=2𝐹1𝑎,𝑎;𝑥21.(3.6)(2)For 𝑗=1, (1𝑥)22𝑎𝐹11𝑎,2;4𝑥(1𝑥)22=2𝐹1𝑎,𝑎;𝑥21+𝑎𝑥2𝐹1𝑎+1,𝑎;𝑥22.(3.7)(3)For 𝑗=2, (1𝑥)22𝑎𝐹11𝑎,2;4𝑥(1𝑥)23=2𝐹1𝑎,𝑎;𝑥22+4𝑎𝑥32𝐹1𝑎+1,𝑎;𝑥22+2𝑎(𝑎+1)3𝑥22𝐹1𝑎+2,𝑎;𝑥23.(3.8)
We remark in passing that the result (3.6) is the Entry 5 of Chapter 11 in Ramanujan’s Notebooks [8, page 50] (with 𝑥 replaced by 𝑥), and the results (3.7) and (3.8) are closely related to (3.6).

3.1. Limiting Cases

In the special cases (3.1) to (3.5), if we replace 𝑥 by 𝑥/𝑏 and let 𝑏, we get, after a little simplification, the known results (1.8) and (1.9) to (1.12), respectively.

4. Application

In this section, we shall first establish the following result, which is given as Entry 4 in the Ramanujan’s Notebooks [8, page 50]: 2𝐹1121𝑎,21𝑎+2;4𝑥(1+𝑥)21𝑏+2=(1+𝑥)𝑎2𝐹11𝑎,𝑎𝑏+21;𝑥𝑏+2,(4.1) by employing (3.1).

Proof. In order to prove (4.1), we require the following result due to Kummer [4]: 2𝐹1;𝑎,𝑏2𝑥=1+𝑥2𝑏(1+𝑥)𝑎2𝐹1121𝑎,21𝑎+2;𝑥21𝑏+2.(4.2) Equation (4.2) is a well-known quadratic transformation recorded in Erdèlyi et al. [9, equation 4, page 111] and also recorded as an Entry 2 in the Ramanujan’s Notebooks [8, page 50]. In (4.2), if we replace 𝑥 by 2𝑥/(1+𝑥), then we have 2𝐹1;4𝑎,𝑏𝑥1+𝑥2=2𝑏1+𝑥2𝑎(1+𝑥)𝑎2𝐹1121𝑎,21𝑎+2;4𝑥(1+𝑥)21𝑏+2.(4.3) Transposing the above equation, we have 2𝐹1121𝑎,21𝑎+2;4𝑥(1+𝑥)21𝑏+2=(1+𝑥)𝑎1+𝑥22𝑎𝐹1;4𝑎,𝑏𝑥1+𝑥22𝑏.(4.4) Now, in (3.1) first replacing 𝑥 by 𝑥 and then replacing 𝑥 by 𝑥 and using on the right-hand side of (4.4), we get 2𝐹1121𝑎,21𝑎+2;4𝑥(1+𝑥)21𝑏+2=(1+𝑥)𝑎2𝐹11𝑎,𝑎𝑏+21;𝑥𝑏+2.(4.5) This completes the proof of (4.1).

Remark 4.1. (1) The result (4.1) can also be established by employing Gauss’s summation theorem.
(2) For generalization of (4.2), see a recent paper by Kim et al. [13].

In our next application, we would like to mention here that in 1996, there was an open problem posed by Baillon and Bruck [19, equation (5.17)] who needed to verify the following hypergeometric identity: 2𝐹11𝑚,22;4𝑧(1𝑧)=(𝑚+1)(1𝑧)𝑧22𝑚1𝐹1;𝑚,𝑚1𝑧𝑧22+(2𝑧1)𝑧22𝑚1𝐹1;𝑚,𝑚1𝑧𝑧21,(4.6) in order to derive a quantitative form of the Ishikawa-tdelstin-ó Brain asymptotic regularity theorem. Using Zeilberger’s algorithm [20], Baillon and Bruck [19] gave a computer proof of this identity which is the key to the integral representation [19, equation  (2.1)] of their main theorem.

Soon after, Paule [21] gave the proof of (4.6) by using classical hypergeometric machinery by means of contiguous functions relations and Gauss’s quadratic transformation (3.1).

Our objective of this section is to obtain first three results from (3.1), (3.2), and (3.4) and then establish again three new results out of which one will be the natural generalization of the Baillon-Bruck identity (4.6).

For this, in our results (3.1), (3.2), and (3.4), if we replace 𝑥 by (𝑧1)/𝑧, we get after a little simplification the following results: 2𝐹1𝑎,𝑏;4𝑧(1𝑧)2𝑏=𝑧22𝑎𝐹11𝑎,𝑎𝑏+2;1𝑧𝑧21𝑏+2,(4.7)2𝐹1𝑎,𝑏;4𝑧(1𝑧)2𝑏+1=𝑧2𝑎2𝐹11𝑎,𝑎𝑏+2;1𝑧𝑧21𝑏+22𝑎2𝑏+11𝑧𝑧2𝐹11𝑎+1,𝑎𝑏+2;1𝑧𝑧23𝑏+2,(4.8)2𝐹1𝑎,𝑏;4𝑧(1𝑧)2𝑏+2=𝑧2𝑎2𝐹11𝑎,𝑎𝑏+2;1𝑧𝑧21𝑏+22𝑎2𝑏+11𝑧𝑧2𝐹11𝑎+1,𝑎𝑏+2;1𝑧𝑧23𝑏+2+4𝑎(𝑎+1)(𝑏+1)(2𝑏+3)1𝑧𝑧22𝐹11𝑎+2,𝑎𝑏+2;1𝑧𝑧25𝑏+2.(4.9)

Finally, in (4.7), (4.8), and (4.9) if we take 𝑎=𝑚 and 𝑏=1/2, we get the following very interesting results: 2𝐹1121,𝑚;4𝑧(1𝑧)=𝑧22𝑚𝐹1;𝑚,𝑚1𝑧𝑧21,(4.10)2𝐹1121,𝑚;4𝑧(1𝑧)=𝑧2𝑚2𝐹1;𝑚,𝑚1𝑧𝑧21+𝑚1𝑧𝑧2𝐹1;𝑚,𝑚+11𝑧𝑧22,(4.11)2𝐹1122,𝑚;4𝑧(1𝑧)=𝑧2𝑚2𝐹1;𝑚,𝑚1𝑧𝑧21+4𝑚31𝑧𝑧2𝐹1;𝑚,𝑚+11𝑧𝑧222𝑚(1𝑚)31𝑧𝑧22𝐹1;𝑚,𝑚+21𝑧𝑧23.(4.12)

Equation (4.7) is a natural generalization of Baillon-Bruck result (4.6). The result (4.11) is an alternate form of the Baillon-Bruck result (4.6). Its exact form can be obtained from (4.11) by using the contiguous function relation 𝑎𝑏𝑧𝑐(𝑐1)2𝐹1=𝑎+1,𝑏+1;𝑧𝑐+12𝐹1𝑎,𝑏;𝑧𝑐+12𝐹1𝑐𝑎,𝑏;𝑧(4.13) with 𝑎=𝑏=𝑚 and 𝑐=1.

We conclude this section by remarking that the result (4.7) is also recorded in [22] by Rathie and Kim who obtained it by other means and the results (4.8) and (4.9) are believed to be new.

Acknowledgments

The author is highly grateful to the referee for carefully reading the manuscript and providing certain very useful suggestions which led to a better presentation of this research article. He also, would like to express his thanks to Professor A. K. Rathie (Center for Mathematical Sciences, Pala, Kerala-India) for all suggestions and for his encouraging and fruitful discussions during the preparation of this research article. The author was supported by the research Grant (IG/SCI/DOMS/12/05) funded by Sultan Qaboos University-Oman.