Arxiv:2008.09357V2 [Math.CA] 8 Sep 2020

Home , Appell series

arXiv:2008.09357v2 [math.CA] 8 Sep 2020 a 2018 May 17-219 DESY -eiaie fmliaibeqhpremti functio q-hypergeometric multivariable of q-derivatives 2 1 -al [email protected] E-mail: [email protected] E-mail: aeo umto nee ihpstv elcecet.A nexam function an As hypergeometric equ type coeﬃcients. Horn Explicit real non-conﬂuent positive variables. with of indexes several embodie summation and function of one case This of parameters. the series to q-hypergeometric respect with function ric ecnie h -eiaie fteSiatv n autbscm basic Daoust and Srivastava the of q-derivatives the consider We b colo hsc n srnm,SnYtsnUniversity, Yat-sen Sun Astronomy, and Physics of School ldmrV Bytev, V. Vladimir 490Dba(ocwRgo) Russia Region), (Moscow Dubna 141980 a on nttt o ula Research, Nuclear for Institute Joint 102Zua,China, Zhuhai, 519082 parameters Abstract a 1 egigZag, Zhang Pengming H 3 spresented. is nienme fvarious of number entire a s tosaegvnfrgeneral for given are ations l eiaie fq-analog of derivatives ple ihrsett their to respect with n liaibehypergeomet- ultivariable b, 2 SN0418-9833 ISSN 1 Introduction

Basic hypergeometric functions, namely q-analogue generalizations of ordinary hypergeometric series, has a long history [2]. Contrary to the generalized hypergeometric series, where the ratio of successive terms is a rational function of n, the ratio of successive terms in basic hypergeometric series is a rational function of qn. The basic hypergeometric series 2φ1(qα, qβ; qγ; q,x), q-extension of the Gauss hypergeometric function 2F1 was first considered by Eduard Heine in 1846. The last could be restored from q- analogue by applying parameters equal qαi and considering the limit when the base q tends to unity. Renewed interest in hypergeometric-type series containing digamma functions has emerged in connection with derivatives of hypergeometric functions with respect to their parameters. The first derivatives for some special values of the parameters were already known a long time ago [12–14]. Later on, Ancarani et al. have found in a series of papers [3, 4, 7] the derivatives of Gaussian hypergeometric functions and some derivatives of two-variable series, namely the Appell series and four degenerate confluent series [5]. Moreover, it has been shown that the first derivatives of generalized hypergeometric functions are expressible in terms of Kampéde Fériet functions, and, with the same technique, derivatives of the Appell hypergeometric function have been obtained in Ref. [6]. In the paper [19] the general case of hypergeometric function derivative with respect to the parameters was considered. Namely, the case of multi-summation index with arbitrary natural coefficients was calculated. It was shown, that the derivatives with respect to the upper and lower parameters of the (Srivastava-Daoust) generalized Lauricella series are expressed as a final sum of the generalized Lauricella series with shifted parameters. In quantum field theory, one has to calculate higher-order Feynman diagrams for quantum corrections to scattering processes. For the evaluation of the dimensionally regularized Feynman integrals in D = 4 − 2ε dimensional space time, one has to construct the expansion in ε. In the [19] it was concluded that the ε expansions of Feynman integrals at any order are expressible in terms of Horn-type hypergeometric functions: the n-th term of the ε series can be expressed as a Horn-type hypergeometric function in n + m variables, where m is the number of summations in the Horn- type representation of the Feynman integral. The region of convergence of any of these parameter derivatives, i.e., the coefficients in the ε expansion, and the initial Feynman integral are the same. Classes of the Feynman integrals that could be expressed in the terms of multiple polylogaritms (MP) [11] are mainly related to the massless cases. More interesting massive case, for example massive sunrise and kite integrals could be expressed through elliptic generalization of MP (EP), see [8], [9]. At [10] it was shown, that EP can be written as a multiple integral of basic hypergeometric function convolution. The utilization of q-difference equations and q-contiguous relations gives one the possibility to establish a connection with basic hypergeometric functions and find efficient algorithm for analytical continuation and numerical evaluation of EP. By analogy with paper [19] the knowledge of of q-basic function derivatives over parameters can help us to find out the properties of Feynman integral expansion over dimensional parameter and make stable numerical estimations. The generalization of hypergeometric series derivatives, the q-derivatives of basic hypergeometric series with respect to the parameters was considered in the papers [20–24]. The case of q-hypergeometric series rΦs (q-extension of generalized hypergeometric series) was studied at [20], later the derivatives of Appell functions with respect to parameters and q-derivatives of 3-variable

2 q-Lauricella functions [21] as well as k- variable q-Lauricella functions and q-Kampe de Feriet function [22, 23] was calculated. Some special cases of q-derivative of Srivastava’s general triple q-hypergeometric series with respect to its parameters HA,q,HB,q,HC,q are considered in [24]. These multiple basic hypergeometric series considered above are special cases of the q-extension of generalized Lauricella series in n variables (1). At this paper we would like to present the q- derivatives of the basic (Srivastava-Daoust) generalized Lauricella series with respect to the upper and lower parameters and generalize results obtained at [20–24]. Here we consider the case when (1) (1) (1) (1) (n) (n) (n) (n) coefficients θ1 , ψ1 , φ1 , δ1 ,...,θA , ψC , φB(n) , δD(n) of multi summation indexes in (2) are positive reals, contrary to the considered special cases, where these coefficient are implied equal zero or unity. This paper is organized as follows. We begin in Sec. 2 by considering definitions for q-extension of generalized Lauricella hypergeometric function and some basic definitions of the q-special function theory, namely Jackson derivative, q-Pochhammer symbol and q-bracket. Next, Sec. 3 is devoted to the q-derivatives in the case of double summation index parameters with real positive coefficients, namely as an example the q-analog of non-confluent Horn type hypergeometric function H3(a, b, c, x, y) q-derivatives with respect to the upper and lower parameters are considered. Our main result are presented in Sec. 4, where we consider the q-derivatives of q-extension of the (Srivastava-Daoust) generalized Lauricella series in n variables with respect to the upper and lower parameters. The cases of multiple summation indexes with positive real coefficient are considered.

2 Deﬁnitions

Q-extension of the (Srivastava-Daoust) generalized Lauricella series in n variables, given by Srivas- tava [1] is deﬁned:

x1 (1) (n) F A:B ;...;B  .  C:D(1);...;D(n) .  xn    (1) (n) 1 (1) n (n) A:B(1);...;B(n) [(a) : θ ,...,θ ] : [(b ) : φ ]; . . . ; [(b ) : φ ] = F q; x ,...,xn C:D(1);...;D(n) [(c) : ψ(1),...,ψ(n)] : [(d1) : δ(1)]; . . . ; [(dn) : δ(n)] 1 ∞ s1 sn x1 xn = Ω(s1,...,sn) . . . , (1) (q,q)s1 (q,q)sn s1,...,sXn=0 where A B(1) (1) · · · B(n) (n) j=1(aj,q) (1) ··· (n) j=1 (bj ,q) (1) j=1 (bj ,q) (n) s1θj + +snθj s1φj snφj Ω(s1,...,sn)= Q Q Q n , (2) C D(1) (1) · · · D( ) (n) j=1(cj,q) (1) ··· (n) j=1 (dj ,q) (1) j=1 (dj ,q) (n) s1ψj + +snψj s1δj snδj Q Q Q and all parameters

(1) (1) (1) (1) (n) (n) (n) (n) θ1 , ψ1 , φ1 , δ1 ,...,θA , ψC , φB(n) , δD(n) (3)

3 are so constrained that multiple series (1) converges. Here we impose additional condition that parameters (3) are non-negative real numbers. Here (a, q)n is q-shifted factorial and is defined by [2]: n−1 − m (a, q)n = Πm=0(1 aq ), (a, q)0 = 1, (4) (a, q)∞ (a, q)n = (5) (aqn,q)∞ The last definition of q-shifted factorial helps us to extend the definition area of variable n to negative integers and real values: − n n(n 1) 1 1 (−q/a) q 2 (a, q)− = = = (6) n n − −m −n Πm=1(1 aq ) (aq ,q)n (q/a, q)n The q-analog of n, also known as the q-bracket or q-number of n, is defined to be 2 n−1 [n]q =1+ q + · · · + q , (7) The equivalent definition as 1 − qn [n] = , (8) q 1 − q gives us possibility to extend the definition over real values of n. The Jackson [16] q-derivative is a q-analog of the ordinary derivative of a function f(x) defined as: f(qx) − f(x) D f(x)= . (9) x,q (q − 1)x It has product rule analogous to the ordinary derivative product, we could produce from above equation the q-derivative for the ratio of two functions:

f(x) g(x)Dx,qf(x) − f(x)Dx,qg(x) Dx,q = = , (10) g(x) g(qx)g(x) more properties could be find in [17], [2], [18]. For our purposes we shorten notation for (1,12) functions to names without parameters and variables and explicitly written out only those parameters that are different from those in definition (1,12). For example, F (λx2) means that instead of variable x2 we will use variable λx2 in definition (1) of F .

3 Q-diﬀerence of q-analog of non-conﬂuent Horn type hypergeometric function H3

The q-analog of non-conﬂuent Horn-type hypergeometric function H3(a, b, c, z1, z2)

(a)2n1+n2(b)n2 n1 n2 H3(a, b, c, q, z1, z2)= z1 z2 (11) n1!n2!(c)n1+n2 nX1,n2 could be written in the form:

(a, q)2n1+n2 (b, q)n2 n1 n2 Hq,3(a, b, c, q, z1, z2)= z1 z2 . (12) (q,q)n1 (q,q)n2 (c, q)n1+n2 nX1,n2

4 3.1 Q-derivative with respect to the upper parameters Derivative over upper one-summation parameter b could be calculated by the same algorithm as in [15]: − (bq, q)n2 (b, q)n2 (a, q)2n1+n2 n1 n2 Db,qHq,3(a, b, c, q, z1, z2)= − z1 z2 (q 1)b (q,q)n1 (q,q)n2 (c, q)n1+n2 nX1,n2 1 (bq, q)∞ (b, q)∞ (a, q) = − 2n1+n2 zn1 zn2 − n2+1 n2 1 2 (q 1)b (bq ,q)∞ (bq ,q)∞ (q,q)n1 (q,q)n2 (c, q)n1+n2 nX1,n2 n 1 (b, q)∞ 1 − bq 2 (a, q) = − 1 2n1+n2 zn1 zn2 − n2 − 1 2 (q 1)b (bq ,q)∞ 1 b (q,q)n1 (q,q)n2 (c, q)n1+n2 nX1,n2 − n2 − 1 q 1 (a, q)2n1+n2 (b, q)n2 n1 n2 = − − z1 z2 . (13) (1 b) (q 1) (q,q)n1 (q,q)n2 (c, q)n1+n2 nX1,n2

By using the derivative of H3 over variable zi

1 (a, q)2n1+n2 (b, q)n2 n1 n2 Dzi,qHq,3(a, b, c, q, z1, z2)= [ni]q z1 z2 , (14) zi (q,q)n1 (q,q)n2 (c, q)n1+n2 nX1,n2 we could write (13) as follows: z D H (a, b, c, q, z , z )= − 2 D H (a, b, c, q, z , z ). (15) b,q q,3 1 2 1 − b z2,q q,3 1 2 It is not so straightforward calculation of q-derivative with respect to the double-summation parameter a as in the case of one-summation index b, eq. (13): − (aq, q)2n1+n2 (a, q)2n1+n2 (b, q)n2 n1 n2 Da,qHq,3(a, b, c, q, z1, z2)= − z1 z2 (q 1)a (q,q)n1 (q,q)n2 (c, q)n1+n2 nX1,n2 1 (aq, q)∞ (a, q)∞ (b, q) = − n2 zn1 zn2 − 2n1+n2+1 2n1+n2 1 2 (q 1)a (aq ,q)∞ (aq ,q)∞ (q,q)n1 (q,q)n2 (c, q)n1+n2 nX1,n2 n n 1 (a, q)∞ 1 − aq2 1+ 2 (b, q) = − 1 n2 zn1 zn2 − 2n1+n2 − 1 2 (q 1)a (aq ,q)∞ 1 a (q,q)n1 (q,q)n2 (c, q)n1+n2 nX1,n2 − 2n1+n2 − 1 q 1 (a, q)2n1+n2(b, q)n2 n1 n2 = − − z1 z2 . (16) (1 a) (q 1) (q,q)n1 (q,q)n2 (c, q)n1+n2 nX1,n2

Unlike to the eq. (13) we have here q-number [2n1 +n2]q, which contains double summation indexes with the integer coeﬃcient 2. By means of identity that splits double summation indexes to two items with t one summation indexes 1 [2n + n ] = (1 + q2n1 )[n ] +(1+ qn2 )[2n ] , (17) 1 2 q 2 2 q 1 and the expression for the q-derivative of H3 function with squared variable z2,

2 1 (a, q)2n1+n2 (b, q)n2 n1 2n2 Dz2,qHq,3(a, b, c, q, z1, z2 )= [2n2]q z1 z2 , (18) z2 (q,q)n1 (q,q)n2 (c, q)n1+n2 nX1,n2

5 we get the following answer for q-derivative with respect to the double-summation index parameter a: −1 Da,qHq, (a, b, c, q, z , z )= 3 1 2 2(1 − a) × 2 2 2 z2Dz2,q(Hq,3 + Hq,3(q z1)) + z1Dz1,q(Hq,3(z1 )+ Hq,3(z1 , qz2)) .(19) 3.2 Q-derivative with respect to the lower parameters The similar procedure could be applied for derivation with respect to the double-summation lower parameter c: − − 1 (cq,q)n1+n2 (c, q)n1+n2 Dc,qHq,3(a, b, c, z1, z2) = − (c, q)n1+n2 (q 1)c nX1,n2 (a, q) (b, q) × 2n1+n2 n2 n1 n2 z1 z2 (q,q)n1 (q,q)n2 (cq,q)n1+n2 −1 (cq,q)∞ (c, q)∞ = − − n1+n2+1 n1+n2 c(q 1)(c, q)n1+n2 (cq ,q)∞ (cq ,q)∞ nX1,n2 (a, q) (b, q) × 2n1+n2 n2 n1 n2 z1 z2 (q,q)n1 (q,q)n2 (cq,q)n1+n2 −1 1 − cqn1+n2 = − 1 (q − 1)c 1 − c nX1,n2 (a, q) (b, q) × 2n1+n2 n2 n1 n2 z1 z2 (q,q)n1 (q,q)n2 (cq,q)n1+n2 n1+n2 − 1 q 1 (a, q)2n1+n2(b, q)n2 n1 n2 = − − z1 z2 . (20) (1 c) (q 1) (q,q)n1 (q,q)n2 (cq,q)n1+n2 nX1,n2

With the help of identity similar to the eq. (17)

1 n1 n2 [n + n ]q = ((1 + q )[n ]q +(1+ q )[n ]) , (21) 1 2 2 2 1 and eq. (14) we could write the ﬁnal answer for the q-derivative with the respect to the lower parameter c: 1 D H (a, b, c, z , z )= [z D (H + H (qz )) + z D (H + H (qz ))] . (22) c,q q,3 1 2 2(1 − c) 2 z2,q q,3 q,3 1 1 z1,q q,3 q,3 2

6 4 Q-derivative of q-extension of the generalized Lauricella series with respect to the parameters

4.1 Derivative with respect to the upper parameter For calculation of the derivative we need and extension of identity (17) to the case of multiple summation index with real non-negative coeﬃcients θi. We will put it here explicitly:

1 θ2m2 θ2m2+θ3m3 θ2m2+...θkmk [θ m + ...θ m ]q = ((1 + q + q + ... + q )[θ m ]q 1 1 k k k 1 1 θ3m3 θ3m3+θ4m4 θ3m3+...θkmk+θ1m1 +(1 + q + q + ... + q )[θ2m2]q · · ·

θ1m1 θ1m1+θ2m2 θ1m1+...θk−1mk−1 +(1 + q + q + ... + q )[θkmk]q), (23) here all mj are non-negative integers. The derivative over variable zk of hypergeometric function (1) reads: ∞ s1 sn 1 z1 zn Dzk,qF (z1, ..., zn)= [sk]qΩ(s1,...,sn) . . . . (24) zk (q,q)s1 (q,q)sn s1,...,sXn=0

Now we could obtain the explicit formula for q-derivative of hypergeometric function F with respect to the upper parameter aj:

∞ − (ajq,q) (1) ··· (n) (aj,q) (1) ··· (n) s1θj + +snθj s1θj + +snθj Daj ,qF = (q − 1)aj s1,...,sXn=0 s1 sn 1 x1 xn × Ω(s1,...,sn) . . . (aj,q) (1) ··· (n) (q,q)s1 (q,q)sn s1θj + +snθj n ∞ s θ(1)+···+s θ( ) 1 (a ,q)∞ 1 − a q 1 j n j = j j − 1 − (1) ··· (n)  −  (q 1)aj s1θj + +snθj 1 aj s1,...,sXn=0 (ajq ,q)∞   s1 sn 1 x1 xn × Ω(s1,...,sn) . . . (aj,q) (1) ··· (n) (q,q)s1 (q,q)sn s1θj + +snθj ∞ s θ(1) ··· s θ(n) −1 q 1 j + + n j − 1 = (1 − aj) (q − 1) s1,...,sXn=0 s1 sn x1 xn ×Ω(s1,...,sn) . . . , (25) (q,q)s1 (q,q)sn by applying the multivariable versions (23, 24) of (17, 14) we obtain the main result for q-derivation

7 of the generalized Lauricella q-extension series with respect to the upper parameter: −1 D F = z D (F (zθ1 )+ F (zθ1 ,qθ2 z )+ F (zθ1 ,qθ2 z ,qθ3 z ) aj ,q − 1 z1,q 1 1 2 1 2 3 k(1 aj) h · · · θ1 θ2 θk + + F (z1 ,q z2,...,q zk)) θ2 θ2 θ3 θ2 θ3 θ4 +z2Dz2,q(F (z2 )+ F (z2 ,q z3)+ F (z2 ,q z3,q z4) · · · θ1 θ2 θ3 θk + + F (q z1, z2 ,q z3,...,q zk)) · · ·

θk θ1 θk θ1 θ2 θk +zkDzk,q(F (zk )+ F (q z1, zk )+ F (q z1,q z2, zk )

· · · θ1 θk−1 θk + + F (q z1,...,q zk−1, zk )) . (26) i (1) For the one-summation upper parameter index q-derivation bj = bj we could produce:

∞ − (bjq,q) (1) (bj,q) (1) s1 sn s1φ s1φ 1 x x D F = j j Ω(s ,...,s ) 1 . . . n bj ,q − 1 n (q 1)bj (bj,q) (1) (q,q)s1 (q,q)sn s1,...,sXn=0 s1φj ∞ s φ(1) 1 (b ,q)∞ 1 − b q 1 j = j j − 1 − (1)  −  (q 1)bj s1φj 1 bj s1,...,sXn=0 (bjq ,q)∞   s1 sn 1 x1 xn × Ω(s1,...,sn) . . . (bj,q) (1) (q,q)s1 (q,q)sn s1φj ∞ (1) − s1φj − s1 sn 1 q 1 x1 xn = − − Ω(s1,...,sn) . . . , (27) (1 bj) (q 1) (q,q)s1 (q,q)sn s1,...,sXn=0 and with the help of eq. (24) we obtain:

−1 (1) φj Dbj ,qF = z1Dz1,qF (q z1). (28) (1 − bj)

→ (1) → (1) (2) (n) → The eq. (28) also could be obtained by applying aj bj, θj φj , θj ...θj 0. (2) (n) The derivative with respect to other upper parameters bj . . . bj could be obtained from eq. → (i) (1) → (i) → (28) by substitution bj bj , φj φj , z1 zi:

−1 (i) φj Db(i),qF = (i) ziDzi,qF (q zi). (29) j − (1 bj )

8 4.2 Q-derivative with respect to the lower parameter With the same procedure we consider the q-derivative of the Srivastava and Daoust multivariable hypergeometric series with respect to the lower parameter cj:

∞ − (cjq,q) (1) ··· (n) (aj,q) (1) ··· (n) s1ψj + +snψj s1ψj + +snψj Dcj ,qF = (q − 1)cj s1,...,sXn=0 − s1 sn 1 x1 xn × Ω(s1,...,sn) . . . (qcj,q) (1) ··· (n) (q,q)s1 (q,q)sn s1ψj + +snψj ∞ s ψ(1)+···+s ψ(n) −1 (c ,q)∞ 1 − c q 1 j n j = j j − 1 − (1) ··· (n)  −  (q 1)cj s1ψj + +snψj 1 aj s1,...,sXn=0 (cjq ,q)∞   s1 sn 1 x1 xn × Ω(s1,...,sn) . . . (qcj,q) (1) ··· (n) (q,q)s1 (q,q)sn s1ψj + +snψj ∞ (1) ··· (n) s1ψj + +snψj − (cj,q) (1) ··· (n) 1 q 1 s1ψj + +snψj = − − (1 cj) (q 1) (qcj,q) (1) ··· (n) s1,...,sXn=0 s1ψj + +snψj s1 sn x1 xn ×Ω(s1,...,sn) . . . . (30) (q,q)s1 (q,q)sn By applying (23, 24) we obtain: 1 D F = z D (F (qc , zθ1 )+ F (qc , zθ1 ,qθ2 z )+ F (qc , zθ1 ,qθ2 z ,qθ3 z ) cj ,q − 1 z1,q j 1 j 1 2 j 1 2 3 k(1 cj) h · · · θ1 θ2 θk + + F (qcj, z1 ,q z2,...,q zk)) θ2 θ2 θ3 θ2 θ3 θ4 +z2Dz2,q(F (qcj , z2 )+ F (qcj, z2 ,q z3)+ F (qcj, z2 ,q z3,q z4) · · · θ1 θ2 θ3 θk + + F (qcj,q z1, z2 ,q z3,...,q zk)) · · ·

θk θ1 θk θ1 θ2 θk +zkDzk,q(F (qcj , zk )+ F (qcj,q z1, zk )+ F (qcj,q z1,q z2, zk )

· · · θ1 θk−1 θk + + F (qcj,q z1,...,q zk−1, zk )) . (31) i

9 (i) For the one-summation index lower parameter q-derivation dj we could produce:

(i) (i) ∞ (i) − (i) s (dj q,q)s δ (dj ,q)s δ − 1 sn i j i j 1 x1 xn D (i) F = Ω(s1,...,sn) . . . dj ,q (i) (i) − i (q,q)s1 (q,q)sn s1,...,sXn=0 (q 1)dj (dj q,q) ( ) siδj ∞ (i) (i) (i) − siδj 1 (dj ,q)∞ 1 dj q =  − 1 (i) i (i) (i) − ( ) siδj − s ,...,sXn=0 (q 1)dj (d q ,q)∞ 1 dj 1 j   −1 xs1 xsn × Ω(s ,...,s ) 1 . . . n (i) 1 n (q,q)s (q,q)s (dj q,q) (i) 1 n siδj

∞ (i) (i) siδ (d ,q) (i) s1 sn 1 q j − 1 j siδ x x = j Ω(s ,...,s ) 1 . . . n , (32) (i) − (i) 1 n − (q 1) i (q,q)s1 (q,q)sn s1,...,sXn=0 (1 dj ) (dj q,q) ( ) siδj and with the help of eq. (24) we obtain: 1 (i) δj D (i) F = z D F (q z ). (33) d ,q (i) i zi,q i j − (1 dj ) The equations (26, 28, 31, 33) show the ﬁnal expressions for any q-derivative with respect to the parameters of the q-extension of the generalized Lauricella series.

5 Conclusion

We extend the algorithm [19] for derivation of generalized Lauricella hypergeometric series to the case of basic hypergeometric function q-derivatives and consider the general case of q-derivatives with respect to the parameters for q-extension of the generalized Lauricella series. Namely, the derivation formula for upper parameter is presented in (26,28), derivatives over lower parameter could be find at (31, 33). By that equations we cover all possible cases parameter derivation of Srivastava and Daoust basic multivariable hypergeometric function with positive summation index and explicitly show that the parameter derivative could be expressed as a finite sum of derivative over variables with shifted values. As an example we produce q-derivative of non-confluent Horn-type hypergeometric function H3(a, b, c, z1, z2) q-analog, over upper (15) and lower (22) parameters correspondingly.

Acknowledgments

This work was partially supported by the National Key Research and Development Program of China (No. 2016YFE0130800), the National Natural Science Foundation of China (Grants No. 11975320) and by the Heisenberg–Landau Program.

References

[1] H. M. Srivastava, Certain q-polynomial expansions for functions of several variables, IMA J. App. Math. 30 315 (1983).

10 [2] G. Gasper, M. Rahman Basic Hypergeometric Series, Cambridge, Encyclopedia of Mathe- matics and its Applications) (2004).

[3] L. U. Ancarani and G. Gasaneo, Derivatives of any order of the Gaussian hypergeometric function 2F1(a, b, c; z) with respect to the parameters a, b and c, Journal of Physics A: Math- ematical and Theoretical 42, 395208 (2009).

[4] L. U. Ancarani and G. Gasaneo, Derivatives of any order of the hypergeometric function pFq(a1, ..., ap; b1, ..., bq; z) with respect to the parameters ai and bi, Journal of Physics A: Math- ematical and Theoretical 43, 085210 (2010).

[5] L. U. Ancarani, J. A. D. Punta, and G. Gasaneo, Derivatives of Horn hypergeometric functions with respect to their parameters, J. Math. Phys. 58, 073504 (2017).

[6] V. Sahai and A. Verma, Derivatives of Appell functions with respect to parameters, Journal of Inequalities and Special Functions 6, 1 (2015).

[7] B. Xh. Fejzullahu, Parameter derivatives of the generalized hypergeometric function, Integral Transforms and Special Functions 28, 781 (2017).

[8] L. Adams, C. Bogner and S. Weinzierl, J. Math. Phys. 57 (2016) no.3, 032304 doi:10.1063/1.4944722 [arXiv:1512.05630 [hep-ph]].

[9] L. Adams, C. Bogner, A. Schweitzer and S. Weinzierl, J. Math. Phys. 57 (2016) no.12, 122302 doi:10.1063/1.4969060 [arXiv:1607.01571 [hep-ph]].

[10] G. Passarino, Eur. Phys. J. C 77 (2017) no.2, 77 doi:10.1140/epjc/s10052-017-4623-1 [arXiv:1610.06207 [math-ph]].

[11] A. LEVIN, Compositio Mathematica, 106 (1997) no.3, 267 doi:10.1023/A:1000193320513

[12] M. Abramowitz, Handbook of mathematical functions, with formulas, graphs, and mathematical tables, (Dover Publications Inc., New York, N.Y., 1974).

[13] Yu. A. Brychkov, Handbook of special functions: derivatives, integrals, series and other formulas, (Chapman and Hall/CRC Press, New York, 2008).

[14] J. Froehlich, Parameter derivatives of the Jacoby polynomials and the Gaussian hypergeometric function, Integral Transforms and Special Functions 2, 253 (1994).

[15] H. Ghany, Q-derivative of basic hypergeometric series with respect to parameters, Int. J. of Math. Analysis, 3 1617 (2009).

[16] F. N. Jackson, On q-functions and certain diﬀerence operator, Trans. Roy. Soc. Edin, 46 (1908).

[17] H. Exton, Q-hypergeometric functions and applications, Ellis Horwood Limited, Chichester, (1983).

[18] F. N. Jackson, On basic hypergeometric functions, Quart. J. Math. 13 69 (1942).

11 [19] V. V. Bytev and B. A. Kniehl, Derivatives of Horn-type hypergeometric functions with respect to their parameters, arXiv:1712.07579 [math-ph]. Nucl. Phys. B 952 (2020)

[20] H. A. Ghany, Q-derivative of basic hypergoemtric series with respect to parameters, Int. J. Math. Anal 3 33 (2009).

[21] V. B. Kuznetsov and E. K. Sklyanin, Factorisation of Macdonald polynomials, Symmetries and integrability of diﬀerence equations London Math. Soc. Lecture Note Ser., 255 370 (1999).

[22] V. Sahai and A. Verma, nth-order q-derivatives of multivariable q-hypergeometric series with respect to parameters, Asian-Eur. J. Math. 7 29 (2014).

[23] V. Sahai and A. Verma, Derivatives of Appell functions with respect to parameters, J. Inequal. Spec. Funct., 6 1 (2015)

[24] V. Sahai and A. Verma, Derivatives of Appell functions with respect to parameters, KYUNG- POOK Math. J. 56 911 (2016)