arxiv:math/ v4 [math.ca] 29 Jul 2006

Transcription

1 arxiv:math/0109v4 [math.ca] 9 Jul 006 Contiguous relations of hypergeometri series Raimundas Vidūnas University of Amsterdam Abstrat The 15 Gauss ontiguous relations for F 1 hypergeometri series imply that any three F 1 series whose orresponding parameters differ by integers are linearly related over the field of rational funtions in the parameters. We prove several properties of oeffiients of these general ontiguous relations, and use the results to propose effetive ways to ompute ontiguous relations. We also disuss ontiguous relations of generalized and basi hypergeometri funtions, and several appliations of them. Keywords: Hypergeometri funtion, ontiguous relations. 1 Contiguous relations of F 1 series In this paper, let a, b a, b F := F 1 z = 1 + a b a a + 1bb + 1 z + z denote the Gauss hypergeometri funtion with the argument z. Two hypergeometri funtions with the same argument z are ontiguous if their parameters a, b and differ by integers. For example, F a, b and F a+10, b 7 are ontiguous. As is known [AAR99,.5], for any three ontiguous F 1 funtions there is a ontiguous relation, whih is a linear relation, with oeffiients being rational funtions in the parameters a, b, and the argument z. For example, a + 1, b a, b a z 1F + a az + bzf a 1, b + af = 0, a + 1, b a, b a, b af 1F + a 1F = 0, 3 1 a + 1, b a, b + 1 af bf Supported by NWO, projet number b af +3 a, b =

2 In these relations two hypergeometri series differ just in one parameter from the third hypergeometri series, and the differene is 1. The relations of this kind were found by Gauss; there are 15 of them. A ontiguous relation between any three ontiguous hypergeometri funtions an be found by ombining linearly a sequene of Gauss ontiguous relations. In the next setion we disuss this and other ways to ompute ontiguous relations. The following theorem summarizes some properties of oeffiients of ontiguous relations. These results are useful in omputations and appliations of ontiguous relations. We assume that the parameters a, b, and z are not related, and by S a, S b, S we denote the shift operators a a + 1, b b + 1 and + 1 respetively. Theorem 1.1 For any integers k, l, m there are unique funtions Pk, l, m and Qk, l, m, rational in the parameters a,b, and z, suh that a + k, b + l a, b a + 1, b F = Pk, l, m F + Qk, l, mf. 5 + m These funtions satisfy the same ontiguous relations as F a, b, that is, if 3 a + αn, b + β n f n F = 0, f 1, f, f 3 Ca, b,, z, 6 + γ n n=1 is a ontiguous relation, then 3 n=1 S k a S l b Sm f n Pk+αn, l+β n, m+γ n = 0, and similarly for Qk, l, m. Besides, the following expressions hold: P0, 0, 0 = 1, Q0, 0, 0 = 0, P1, 0, 0 = 0, Q1, 0, 0 = 1, Pk, l, m = a 1 a + 11 z S aqk 1, l, m, 8 Qk + k, l + l, m + m = Sa k SbS l m Qk, l, m Qk+1, l, m + Sa k+1 Sb l a Sm a 1 z Qk 1, l, m Qk, l, m, 9 Pk, l, m Pk, l, m Qk, l, m Qk, l, m = m m z m z 1 m k l a+1 k b l a m k b m l 7 S k a S l b S m Qk k, l l, m m, 10 Q k, l, m = 1m+1 a k 1 b l z m 1 z k+l m 1 m 1 m a k m b l m S k a S l b S m Qk, l, m. 11

3 Here a k = Γa + k/γa is the Pohhammer symbol; for positive k it is equal to a a a + n 1. This theorem is proved in Setion 3. In the following setion, we disuss omputational aspets and appliations of ontiguous relations, inluding ontiguous relations of generalized and basi hypergeometri series. Computational aspets and appliations Computational aspets. To ompute a ontiguous relation 6 between three F 1 series, one an take one of the series as F a, b. If neessary, one an apply a suitable shift operator to the ontiguous relation afterwards. Then the other two hypergeometri series are F a+k, b+l +m and F a+k, b+l +m for some integers k, l, m and k, l, m. Aording to Theorem 1.1, these two funtions an be expressed linearly in F a, b and F a+1, b as in 5. Elimination of F a+1, b gives: a + k Qk, l, m, b + l a + k F + m Qk, l, m, b + l F + m = Pk, l, m Qk, l, m Pk, l, m Qk, l, m a, b F. 1 Due to formula 8 it is enough to be able to ompute Qk, l, m for any integers k, l, m. For this one an take a finite sequene of integer triples k i, l i, m i whih starts with 0, 0, 0, 1, 0, 0, ends with k, l, m, and any two neighboring triples differ just in one omponent preisely by 1. Then all Qk i, l i, m i inluding Qk, l, m an be omputed onsequently using the simplest Gauss ontiguous relations. A possible hoie of the sequene is to take all triples with intermediate integer values in the first omponent between 0, 0, 0 and k, 0, 0, then all intermediate triples between k, 0, 0 and k, l, 0, and finally all triples between k, l, 0 and k, l, m. Formula 9 allows us to ompute Qk, l, m by divide and onquer tehniques, that is, by reduing the shift k, l, m reursively by half more or less at eah step; ompare with [Knu71, 4.6.3]. A straightforward algorithm of this type is the following: ompute Qk, 0, 0 and Qk+1, 0, 0 using 9 with shifts in the first parameter only so there would be Olog k intermediate Q-terms, perhaps linearly growing in size; then use 4 to ompute Qk, 1, 0; and then alternate similar appliation of 9 and Gauss ontiguity relations to ompute Qk, l, 0, Qk, l+1, 0, Qk, l, 1 and Qk, l, m. Notie that most formulas of Theorem 1.1 annot be used when a, b or are speialized, sine the orresponding shift operators would not be defined then. 3

4 Returning to the omputation of a ontiguous relation 6, note that one an take any other of the given three funtions in 1 as F a, b. Proeeding in the same way one would obtain a ontiguous relation involving, say, Q k, l, m and Qk k, l l, m m. The two ontiguous relations must be the same up to a rational fator and orresponding shifts in the parameters. Formulas 10 and 11 are the expliit relations between the oeffiients implied by this fat. Notie that ombination of 1 and 10 gives an expression of a ontiguous relation with oeffiients linear in Q s. Contiguous relations 5 and 1 an be also omputed by ombining Gauss ontiguous relations themselves in similar ways: after hoosing a sequene of ontiguous hypergeometri series onneting F a+k, b+l +m and F a, b so that neighboring series differ just in one parameter by 1, or by using a similar formula to 9. But this is double work ompared to omputing the Q s. Takayama [Tak89] showed how Gröbner basis an be used to ompute ontiguous relalations. Paule [Pau01] has shown that ontiguous relations an be omputed by a generalized version of Zeilberger s algorithm. However, for large shifts k, l, m this method is not effiient, beause the growing degree in the disrete parameters of oeffiients of ontiguous relations imply larger linear problems. On the other hand, reurrene relations for hypergeometri funtions an be omputed as speial ases of ontiguous relations by the methods desribed here, alternatively to Zeilberger s algorithm. Similar remarks hold also for omputation of ontiguous relations of generalized and basi hypergeometri series. On the website [VK00] there is a link to a Maple pakage for omputing ontiguous relations for F 1 series. Generalization. The generalized hypergeometri series is defined as: a1...,a r rf s b 1,..., b s z a 1 k... a r k := b 1 k... b s k k! zk. 13 Two p F q series are ontiguous if their orresponding upper and lower parameters differ by integers. As for F 1 series, there are linear relations between ontiguous generalized hypergeometri series. In general they relate 1 + maxp, q + 1 hypergeometri series, see [Rai60, 48] and [AAR99, 3.7]. In partiular, relations 3 and 4 also hold for hypergeometri series with more upper and lower parameters. These relations allow to transform any differene operator to a differene operator with different shifts in one upper parameter only by lowering all other upper parameters and raising the lower ones. It follows that general ontiguous relations for p F q series are generated by the relations of type 3 4, and a reurrene relation with the shifts in one upper parameter. Like for F 1 series, ontiguous relations for a lass of p F q funtions an be omputed by linearly ombining a sequene of simplest ontiguous relations. k=0 4

5 One an derive a formula analogous to 5 with at most maxp, q + 1 fixed hypergeometri funtions on the right-hand side. The oeffiients of the fixed hypergeometri funtions would satisfy the ontiguity relations of p F q funtions, with orresponding initial onditions like 7. Those oeffiients usually are not all related by a formula like 8, unless the hypergeometri funtions under onsideration satisfy three-term ontiguous relations say, 3 F 1 funtions. Similarly, there are ontiguous relations for basi hypergeometri or q- hypergeometri series; see [AAR99, 10.9] for the definition of r φ s series. Two suh series are ontiguous if their orresponding upper and lower parameters differ by a power of the base q. Moreover, multipliative q-shifts in the argument z of these funtions an also be allowed, sine the q-shift in the argument an be expressed in q-shifts of the parameters, say: a φ 1 a, b q; qz = a 1 φ 1 aq, b b q; z + φ 1 a, b q; z. For any three ontiguous φ 1 series, where also q-shifts in the argument z are allowed, there is a ontiguous relation. Allowing q-shifts in the argument is natural, beause many transformations of basi hypergeometri series mix the parameters a, b, and the argument z, see [AAR99, ]. A q-differential equation for r φ s series an be interpreted as a ontiguous relation in this more general sense, sine it an be seen as a q-differene equation. Evaluation. Contiguous relations an be used to evaluate a hypergeometri funtion whih is ontiguous to a hypergeometri series whih an be satisfatorily evaluated. For example, Kummer s identity [AAR99, Cor. 3.1.] a, b F a b 1 an be generalized to a + n, b F 1 a b 1 = Pn = Γ1 + a b Γ1 + a Γ1 + a Γ1 + a 14 b. Γa bγ a+1 ΓaΓ a+1 b + Qn Γa bγ a ΓaΓ a 15 b. Here n is an integer, the fators Pn and Qn an be expressed for n 0 as 1 n Pn = n+1 3 F, n+1, a b 1, a 1, 16 Qn = n+1 n 1 n+1 3 F, n, a+1 b 1, 17 3, a+1 and similarly for n < 0, see [Vid03]. In fat, formula 15 is a ontiguous relation between F a+n, b 1 a b 1, F a 1, b 1 a b 1 and F a, b 1 1+a b 1, where the 5

6 last two terms are evaluated in terms of Γ-funtions using Kummer s identity 14, and the oeffiients are expressed as terminating 3 F 1 series. Formulas 35 and 36 in [Vid03] present similar evaluations of a, 1/ F 1 a + 3/ + n 1 4 a + n, b, and 3 F a b, a 1 when n is an integer. These series are ontiguous to the F 1 1/4 and 3 F 1 series evaluable by Gosper s or Dixon s respetively identities. Both new formulas are also three-term ontiguous relations, where two hypergeometri terms are evaluated and the oeffiients are written as terminating hypergeometri series. In both ases all series ontiguous to Gosper s F 1 1/4 or well-posed 3 F 1 series an be evaluated by ontiguous relations, but general expressions for the oeffiients in the final three-term expression like 15 are not known. Formulas like 15 may speialize to one-term evaluations. For example, P 5 = 0 if a 4 a b + b 1 a + 17 b + 1 = 0, and under this ondition formula 15 may be brought to F 1 a 5, b a b 1 = a b 1 a b Γa b Γ a 1 Γa 3 Γ a. 18 b By parameterizing the urve given by the relation between a and b one gets an exoti formula, see [Vid03, 33]. Apparently, this kind of formula an be obtained only by using ontiguous relations and a known evaluation. Transformations. General transformations of hypergeometri series an be derived from the symmetries of their ontiguous relations. For example, all terms in the relations between the 4 Kummer s F 1 funtions see [Erd53,.9] satisfy not only the same hypergeometri differential equation, but also the same ontiguous relations with the same shifts in the parameters a, b and. To show this statement, one an hek that the two funtions a, b u 1 = F 1 z ΓΓ a b a, b, u = Γ aγ b F 1 a+b+1 1 z 19 satisfy basi ontiguous relations 4. Then the funtions ΓΓb a a, a + 1 u 3 = Γ aγb z a F 1 1 a + 1 b, 0 z ΓΓa b b + 1, b u 4 = Γ bγa z b F 1 1 b + 1 a, 1 z Γ Γ1 a Γ1 b a+1, b+1 u 5 = Γ Γ bγ a z1 F 1 z, ΓΓa + b a, b u 6 = 1 z a b F 1 Γa Γb +1 a b, 1 z 3 6

7 satisfy the same ontiguous relations as well, sine the oeffiients in the relations u 6 = u 1 u, u 4 = u 1 u 3, e iπa sinπb a sinπb u 3 = u 1 eiπ a sinπa+b u 6, 4 sinπb u 5 = u 1 sinπsinπa+b u 6, 5 sinπa sinπb are onstants with respet to integral shifts in the parameters a, b,. Other 18 hypergeometri funtions are alternative representations of u 1,...,u 6, see [Erd53,.9]. All relations between the 4 Kummer s funtions are generated by 4 5, so the statement follows. Mind that different expressions of u i s as hypergeometri funtions may have different arguments. In priniple, relations 4 5 an be found by showing that the three terms satisfy the same seond order reurrene relation with respet to integral shifts in some parameter, and omparing their asymptotis as the orresponding parameter approahes and. However, this approah may be umbersome. Quadrati or higher order algebrai transformations an also be found in this way. An interesting question is whether the symmetries of ontiguous relations an lassify all identities between hypergeometri series. For F 1 series this boils down to the symmetries of Qk, l, m. In the remainder of this setion, we demonstrate a way to investigate the symmetries of reurrene relations a speial ase of ontiguous relations of hypergeometri series. Symmetries of reurrene relations. To obtain reurrene relations for some hypergeometri series one an introdue the disrete parameter n by replaing one or several ontinuous parameters as a a + n and similarly in different ways. Then a reurrene is a ontiguous relation between suffiiently many hypergeometri funtions with suessive n. We are interested in situations when two hypergeometri series with a ommon disrete parameter n satisfy the same reurrene relation, perhaps after multiplying the funtions by a hypergeometri term that is, a solution of a first order reurrene relation with oeffiients rational in n, see [PWZ96, 5.1]. Therefore we all two hypergeometri funtions with disrete parameter n equivalent if they differ by suh a fator. More speifially, suppose that a hypergeometri series satisfies a seond order reurrene relation for n 0 AnSn BnSn + CnSn 1 = 0, 6 with An, Bn and Cn being rational funtions in n. For onveniene, we assume that Sn does not satisfy a first order reurrene relation, and that 7

8 Bn 0 for all n 0. Then Sn is equivalent to: n A0...An 1 Zn = 1 B0...Bn 1 Sn. The sequene Zn satisfies the reurrene: Zn + 1 Zn + CnAn 1 Zn 1 = 0. 7 BnBn 1 In fat, Zn is the unique sequene equivalent to Sn and satisfying a reurrene relation of form Zn + 1 Zn + HnZn 1 = 0. Note that the seond order reurrene for Zn is unique up to a fator, and that we want to keep the oeffiients of two terms. It follows that an equivalene lass of hypergeometri funtions is determined by the rational funtion BnBn 1 CnAn 1. 8 Hypergeometri series with a disrete parameter n an be lassified by their equivalene lass funtion. One an ompute this funtion for various types of hypergeometri funtions and different ways of introduing the disrete parameter n. Then one an look for the ases when different equivalene lass funtions are equal perhaps under some relations between their ontinuous parameters. For example, the equivalene lass funtion for 0 F 1 ẑ is the polynomial n+ĉ 1n+ĉ /ẑ. This is also the equivalene lass funtion for 0F 1 ẑ. The equivalene lass funtion for 1 F a+n 1 z is ĉ n ĉ+n +n n+ 1n+ 34n +4n 4n+ +az z 4n +4n 1n+... z. n + a 1n + a 1n + n + 4 When = a, this rational funtion is also a polynomial of degree two, namely 4n+a 1n+a 3/z. By omparing these two polynomials and orresponding transformations of the reurrenes to normalized form 7, we onlude that there must be a linear relation between the three funtions a + n 1F 1 a + n 4z, 0F 1 a+ 1 +n z and Γa+ 1 +n Γ 3 a n z1 a n 0F 1 3 a n z, where oeffiients are onstants with respet to the shift n n+1. Note that the first two funtions are bounded as n, whereas the last one is unbounded for fixed generi a and z. Hene the oeffiient of the last funtion is zero. By omparing the limits of the first two funtions one onludes that a+n 1F 1 4z = e z 0F 1 z ; see also [AAR99, 4.1.1]. a+n a+n+ 1 8

9 The author omputed a number of equivalene lass funtions for r F 1 with r = 0, 1, and 3 F 1 funtions. They do not imply new identities, exept plenty of exoti two-term identities suh as 18, and some interesting onsequenes of known identities. Straightforward omputations are too umbersome when the degree of numerators or denominators of equivalene lass funtions exeeds 10. However, this method generalizes easily to q-hypergeometri funtions. A more intelligent onsideration of symmetries of reurrene and ontiguous relations ould be helpful in finding relations between hypergeometri funtions when one annot use symmetries of differential equations for them, say when the variable z is speialized, or when the parameters and z are related. 3 Proof of Theorem 1.1 To show the existene of 5, observe that ontiguous relations 4 express F a 1, b, F a, b 1, F a, b+1 in terms of F a, b and F a+1, b. There are similar expressions for F a, b +1 and F a, b 1. Using shifted versions of these relations one an express F a+k, b+l +m linearly in hypergeometri funtions without shifts in the parameters b and, and then leave only terms F a, b and F a+1, b. / If expression 5 is not unique, then F a+1, b F a, b is a rational funtion in the parameters, just as F a+1, b+1 / F a+1, b by the symmetry of / the upper parameters. Then F a+1, b a b 1 a, b F 1 1+a b 1 must be a rational funtion in a and b. But this funtion has an unbounded set of poles aording to Kummer s identity 14. Hene a ontradition. The uniqueness of 5 implies that Pk, l, m and Qk, l, m satisfy the same ontiguous relations as F 1 series; hek the triple substitution of 5 into 6. The initial onditions 7 are obvious. To prove 8, apply 5 to both sides of F a+k, b = S a F a+k 1, b, and use to eliminate F a 1, b on the right-hand side. a+k+k, b+l+l +m+m. Formula 9 is obtained after several appliations of 5 to F The first intermediate step is: a + k + k, b + l + l F + m + m = Sa k Sl b Sm Pk, l, m a+k, b+l F + m + Sa k SbS l m Qk, l, m a+k+1, b+l F + m. 9

10 a, b Now apply 5 to get terms with F and F a+1, b only, ompare the terms to F a+1, b, and eventually use 8 one to replae Pk, l, m. For a proof of the last two formulas, let us introdue: Pk, l, m Pk + p, l + q, m + r W p,q,r k, l, m := Qk, l, m Qk + p, l + q, m + r. 9 We assume that p, q, r just as k, l, m are integers. Lemma 3.1 The following properties of the W-symbol hold. i W 0,0,0 k, l, m = 0. ii W p,q,r 0, 0, 0 = Qp, q, r. iii For fixed k, l, m the determinants W p,q,r satisfy the ontiguous relations of F 1 funtions. More preisely, if 6 is a ontiguous relation, then 3 n=1 S k+p a S l+q b S m+r f n W p+αn,q+β n,r+γ n k, l, m = 0. iv W 1,0,0 satisfies the first order reurrene relations W 1,0,0 k + 1, l, m = Sa k+1 S m a W 1,0,0 k, l, m, a 1 z W 1,0,0 k, l + 1, m = Sb l S m b + 1 W 1,0,0 k, l, m, b 1 z W 1,0,0 k, l, m + 1 = Sa k Sb l S m z 1 W 1,0,0 k, l, m. a bz Proof. The first two properties are straightforward. If 6 is a ontiguous relation, then Pk, l, m Qk, l, m 3 n=1 S k+p a S l+q b S m+r f n W p+αn,q+β n,r+γ n k, l, m = 3 n=1 S a k+p S l+q b S m+r f n Pk+p+α n, l+q+β n, m+r+γ n 3 n=1 S a k+p S l+q b S m+r f n Qk+p+α n, l+q+β n, m+r+γ n Pk, l, m 0 = Qk, l, m 0 = 0. The reurrene relation for W 1,0,0 with respet to k follows from ontiguous relation : Pk + 1, l, m ÃPk, l, m + W 1,0,0 k+1, l, m = B Pk+1, l, m Qk + 1, l, m ÃQk, l, m + B Qk+1, l, m = ÃW 1,0,0k, l, m, with à = Sk+1 a S m 10 a a 1 z.

11 To prove the reurrene relation for W 1,0,0 with respet to l we use ontiguous relations 4 and [Erd53,.8.36] to derive Pk, l + 1, m = Pk + 1, l + 1, m = APk, l, m + B Pk + 1, l, m, C Pk, l, m + D Pk + 1, l, m, with A = Sa k Sb l b a b, B = Sk a Sb l a b, C = Sk a Sb lsm a 1 b 1 z, D = Sk a Sb lsm There are similar expressions for Qk, l + 1, m and Qk + 1, l + 1, m. Then W 1,0,0 k, l + 1, m = APk, l, m + B Pk+1, l, m C Pk, l, m + DPk+1, l, m AQk, l, m + B Qk+1, l, m C Qk, l, m + DQk+1, l, m = AD B C W 1,0,0 k, l, m, a+b +1 b 1 z. whih is the required relation. The reurrene with respet to m an be proved similarly, using formulas 34 and 38 in [Erd53,.8]. This ompletes the proof of the lemma. QED. The reurrene relations for W 1,0,0 k, l, m imply that W 1,0,0 k, l, m = m m z m z 1 m k l a+1 k b l a m k b m l. 30 The reurrene relations for W p,q,r and an initial ondition i imply that W p,q,r k, l, m = W 1,0,0 k, l, m S k a S l b S m Qp, q, r. 31 Formula 10 follows from these two equations. To prove 11 note that the determinant in 10 an be expressed in two ways as the symbol W. This gives: W k k,l l,m m k, l, m = W k k,l l,m mk, l, m. 3 By setting k = l = m = 0 and hanging the sign of k, l and m one obtains: W k, l, m 0, 0, 0 = W k,l,m k, l, m. 33 Statement ii of lemma 3.1, formulas 31 and 30, and transformation of Pohhammer symbols give the last formula of Theorem 1.1. QED. Referenes [AAR99] G.E. Andrews, R. Askey, and R. Roy. Speial Funtions. Cambridge Univ. Press, Cambridge,

12 [Erd53] A. Erdélyi, editor. Higher Transendental Funtions, volume I. MGraw-Hill Book Company, New-York, [Knu71] D.E. Knuth. The art of omputer programming, volume II, Seminumerial Algorithms. Addison-Wesley, Reading, Mass., [Pau01] P. Paule. Contiguous relations and reative telesoping. Tehnial report, RISC, Austria, 001. [PWZ96] M. Petkovšek, H.S. Wilf, and D. Zeilberger. A=B. A. K. Peters, Wellesley, Massahusetts, [Rai60] [Tak89] [Vid03] [VK00] E.D. Rainville. Speial Funtions. The MaMillan Company, New York, N. Takayama. Gröbner basis and the problem of ontiguous relations. Japan Journal of Applied Mathematis, 6: , R. Vidūnas. A generalization of Kummer s identity. Roky Mount. J., 3, 003. also available at R. Vidūnas and T.H. Koornwinder. Web-page of the NWO projet Algorithmi methods for speial funtions by omputer algebra. thk/spefun/ompalg.html,