arXiv:math/9411226v1 [math.CA] 21 Nov 1994 M ujc classification. Subject AMS assoc functions, weight functions, generating fractions, ewrsadphrases. and Keywords * eerhprilyspotdb S rn M 235 n NS and 9203659 DMS grant NSF by supported partially Research soitdcontinuous Al-Salam- associated associated Al-Salam-Chihara, associated Wall, fteseta esr sotie.Eee ii ae r d are big cases associated limit include Eleven These c obtained. continuous detail. is some absolutely measure in The spectral the fo derived. expression of is explicit an fraction and continued identified is solution minimal yascae otnosdual continuous associated by Abstract. I.Ascae otnosDual Continuous Associated III. otgosRltos ai yegoercFunctions, Hypergeometric Basic Relations, Contiguous xlctsltosfrtetretr eurnesatisfied recurrence three-term the for solutions Explicit n rhgnlPolynomials. Orthogonal and ai yegoercsre,cniuu eain,conti relations, contiguous series, hypergeometric Basic eateto Mathematics of Department Mathematics of Department eateto Mathematics of Department ap,Foia33620-5700 Florida Tampa, oot,MS11Canada 1A1 M5S Toronto, oot,MS11Canada 1A1 M5S Toronto, q Hriepolynomials. -Hermite nvriyo Toronto of University nvriyo Toronto of University nvriyo Florida of University 31,3D5 07,39A10. 30B70, 33D45, 33D15, ...Ismail M.E.H. ..Masson D.R. ..Gupta D.P. q U.S.A. Hh oyoil r band A obtained. are polynomials -Hahn q Hh oyoil * Polynomials -Hahn ae rhgnlpolynomials. orthogonal iated q Lger associated , -Laguerre alt ,and I, Carlitz h related the r omponent R rn A5384. grant ERC iscussed nued 1. Introduction. The discrete q -Hahn and the discrete dual q - were introduced by Hahn in [10] and [11] respectively. These are particular cases of polynomials introduced by Askey and Wilson and called q -Racah polynomials [3]. Askey and Wilson [3] have also considered continuous dual q -Hahn polynomials as a particular case of the Askey- Wilson polynomials. The objective of the present study is to generalize the continuous dual q -Hahn polynomials to the associated continuous dual q -Hahn polynomials. The q 1 limit gives the case of associated continuous Hahn polynomials which have have been→ studied by Ismail et al [13]. It may be mentioned that in two earlier papers [7], [8], we have discussed the associated continuous Hahn (for continuous Hahn polynomials see Askey [1]) and the associated big q -Jacobi polynomials respectively. In both of these associated cases we made extensive use of contiguous relations for hypergeometric and q -hypergeometric functions. Although the use of contiguous relations in connection with continued fractions goes back to Gauss [23], the importance of contiguous relations in relation to the theory of was first stressed by Wilson [25]. In Section 2, we obtain six solutions to the three-term recurrence relation satisfied by associated continuous dual q -Hahn polynomials. This is done with the help of three- term contiguous relations satisfied by balanced 3φ2 ’s. It is also demonstrated how an existing three-term transformation formula for balanced 3φ2 ’s connects any three of these solutions. By examining the large n asymptotics of the solutions and the associated second order difference equation we show that one of the solutions is a minimal [6],[16] solution. When one of the four parameters is equal to ‘ q ’, another of our solutions reduces to the continuous dual q -Hahn polynomial solution [3]. In Section 3, the related infinite continued fraction is obtained. Following the proce- dure employed in several other cases (see [21], [22], [23]) we then derive the explicit weight function for the absolutely continuous component of the spectrum. Section 4 is devoted to obtaining a generating function and hence an explicit expression for the associated continuous dual q -Hahn polynomials. The method employed is the same as by Ismail and Libis [14] for big q -Laguerre polynomials. In Section 5, we examine four limiting cases of the original recurrence relation together with their solutions, related continued fractions and explicit polynomials. The first two limits are associated big q -Laguerre and associated Wall polynomials. These are on the 2φ1 and 1φ1 levels respectively. Two further limits are found at the 0φ1 level. In Section 6, we consider seven additional limiting cases. These include the associated cases for Al-Salam-Chihara, Al-Salam-Carlitz I, and continuous q -Hermite polynomials. In Section 7, we give the connection between solutions to the associated Askey-Wilson [15],[9] and the associated continuous dual q -Hahn polynomial recurrence relations.

2 2. Three-term contiguous relations and solutions. The recurrence relation satisfied by associated continuous dual q -Hahn polynomials can be expressed as

2 (2.1) X (z a )X + b X − = 0 n+1 − − n n n n 1 an := an(z; A,B,C,D) 1 1 1 1 = + + + qn (1 + q)q2n−1 A B C D −   2 2 bn := bn(z; A,B,C,D) q = (1 Aqn−1)(1 Bqn−1)(1 Cqn−1)(1 Dqn−1). ABCD − − − − The symmetry with respect to the parameters A,B,C,D is obvious. With this form of the recurrence, it is easy to take successive limits A,B,C,D . The limit D gives associated big q -Laguerre polynomials (see Ismail and Libis→ ∞ [14]) and a subsequent→ ∞ limit C gives associated Wall polynomials (see Chihara [4, p. 198]). Finally B and then→∞A give additional cases. →∞ →∞ Note that in (2.1) the an can also be expressed as 1 q a = (λ + µ ) + + n − n n AB CD λ = (1 Aqn)(1 Bqn)/AB n − − µ = q(1 Cqn−1)(1 Dqn−1)/CD. n − − This means that with a renormalization and a translation of the coordinate z we may re-express (2.1) and the aforementioned limits as birth and death processes with birth and death rates λn and µn respectively [12]. A second family with µ0 := 0 should also be investigated. This has already been done in the more general case of associated Askey-Wilson polynomials [15]. The solutions to (2.1) and its limit cases will be expressed in terms of the basic hypergeometric functions ∞ a1,a2, ,ar (a1,a2, ,ar)k k k(k−1)/2 1+s−r k rφs ··· ; z = ··· [( 1) q ] z , z < 1, b1,b2, ,bs (b1,b2, ,bs, q)k − | |  ···  kX=0 ··· where ∞ j−1 n (a)∞ = (1 aq ), (a) = (a)∞/(aq )∞,n integer, − n j=1 Y m (a ,a ,...,a ) = (a ) , n integer or . 1 2 m n k n ∞ kY=1 We will use the notation a,b,c de de (2.2) φ := φ(a,b,c,d,e) = φ ; , < 1, 3 2 d, e abc abc  

3 de and its analytic continuation for abc 1 given by the following transformations [5, (III.9), (III.10), p. 241] ≥

a,b,c de (e/a, de/bc)∞ a,d/b,d/c e (2.2a) 3φ2 ; = 3φ2 ; d, e abc (e, de/abc)∞ d, de/bc a     a,b,c de (b, de/ab, de/bc)∞ d/b, e/b, de/abc (2.2b) 3φ2 ; = 3φ2 ; b . d, e abc (d, e, de/abc)∞ de/ab, de/bc    

To obtain solutions to (2.1) we use 3φ2 contiguous relations and the usual notation

φ(a ) = φ(aq±1,b,c,d,e), ± ±1 ±1 ±1 ±1 ±1 φ± = φ(aq , bq , cq , dq , eq ).

Two such contiguous relations are [18]:

(1 b)(1 c) de (2.3) φ φ(a+) + − − φ = 0, − (1 d)(1 e) abcq + − −

e de (2.4) (1 d)(1 e)φ +(d a) 1 φ (a ) (1 a) 1 φ = 0. − − − − a + − − − − abcq +   Changing (a,b,c,d,e) (a,b/q,c/q,d/q,e/q) in (2.4) and then eliminating φ(a+) and −→ φ+ from the resulting equation together with (2.3) and (2.4), we obtain

(1 b)(1 c)(1 d )(1 e ) de (2.5) − − − a − a φ (a ) (1 d)(1 e) bcq + − − − de d e de (1 a) 1 + a 1 1 + (1 b)(1 c) φ − − − abcq − aq − aq abcq − −   d e    + 1 1 φ−(a+) = 0. − q − q   With the replacements b = Bqn, c = Cqn, d/a = Dqn, e/a = Aqn , this becomes

(1 Bqn)(1 Cqn)(1 Dqn)(1 Aqn) DAa2 (2.6) − − − − φ (1 aDqn)(1 aAqn) BCq n+1 − − DAa DAa (1 a)(1 ) + a(1 Dqn−1)(1 Aqn−1) + (1 Bqn)(1 Cqn) φ − − − BCq − − BCq − − n   n−1 n−1 + (1 Daq )(1 Aaq )φ − = 0, − − n 1 where a,Bqn, Cqn DAa φ = φ ; . n 3 2 aDqn,aAqn BC   4 We write z = q/aDA + a/BC so that

BC 4q (2.7) a = z z2 2 ± r − ABCD ! = BCλ±, say.

After renormalization, (2.6) becomes (2.1) with a solution

(1),± (1),± (2.8) Xn = Xn (z; A,B,C,D) n n (A,B,C,D)n n BCλ±, Bq , Cq = (λ±) 3φ2 n n ; ADλ± . (BCDλ±,ABCλ±) BCDλ±q ,ABCλ±q n   (1),− We shall show later that, with a suitable choice of square root branch, Xn is a minimal solution of (2.1). Because of the symmetry in (2.1) and the fact that the minimal solution is unique up to a constant multiple, the parameter interchanges A C or B D in the (1)↔,− ↔ above solution must yield only an n independent multiple of Xn . Let us verify this. (1),± If we make the interchange A C,Xn changes to ↔ n n (1),± (A,B,C,D)n n BAλ±, Bq , Aq (2.9) ξn = (λ±) 3φ2 n n ; CDλ± . (BADλ±,ABCλ±) BADλ±q ,ABCλ±q n   If we apply the transformation (2.2a) to (2.9) we obtain

n (1),± (A,B,C,D)n n (ADλ±,BCDλ±q )∞ ± ξn = (λ ) n (BADλ±,ABCλ±)n (BADλ±q ,CDλ±)∞

n n Bq , Cq ,BCλ± 3φ2 n n ; ADλ± × ABCλ±q ,BCDλ±q   (1),± which is clearly just a constant multiple of Xn . In order to obtain a second solution of (2.1) we start from a different three-term contiguous relation [18];

(2.10) [de(a b c) + abc(d + e + q a aq)]φ − − − − + (1 a)(de abcq)φ(a+) + bc(d a)(e a)φ(a ) = 0. − − − − − Changing a aqn , we can write the above as → n n n n+1 [de(aq b c) + abcq (d + e + q aq aq )]Yn+ − − −n − n+1 n n (1 aq )(de abcq )Y + bc(d aq )(e aq )Y − =0, − − n+1 − − n 1 where n aq ,b,c de − Y = φ ; q n . n 3 2 d, e abc   5 Writing Aq Aq q 1 1 a = A, b = ABλ , c = ABλ−, d = , e = and ( + ) = z, + C D CD b c and renormalizing we again arrive at equation (2.1) with a new solution

n (A,B)n Aq ,ABλ ,ABλ− 1 − (2.11) X(2) = φ + ; q n+1 . n (AB)n 3 2 Aq/C, Aq/D B   From (2.11) we obtain additional solutions by parameter interchanges due to the symmetry of (2.1). However, we find that this is not true for the interchanges B C or B D . (2) ↔ ↔ This can be seen by applying transformation (2.2a) to Xn . We get

−n+1 n −n+1 (2) (A,B)n (q /D, Aq/B)∞ Aq ,ADλ−,ADλ+ q Xn = − 3φ2 ; (AB)n (Aq/D,q n+1/B)∞ Aq/C, Aq/B D   n −n+1 (A, D)n (Aq/B, q/D)∞ Aq ,ADλ−,ADλ+ q (2.12) = 3φ2 ; (AD)n (Aq/D, q/B)∞ Aq/C, Aq/B D   (2) where the right side is a constant multiple of the B D interchange applied to Xn . It is similarly seen that B C does not yield a new↔ solution. However the interchanges A B, A C and A ↔D do yield the following new solutions: ↔ ↔ ↔ n −n+1 (A,B)n Bq ,ABλ ,ABλ− q (2.13) X(3) = φ + ; , n (AB)n 3 2 Bq/C,Bq/D A   n −n+1 (B, C)n Cq ,BCλ ,BCλ− q (2.14) X(4) = φ + ; , n (BC)n 3 2 Cq/A, Cq/D B   n −n+1 (B, D)n Dq ,BDλ ,BDλ− q (2.15) X(5) = φ + ; . n (BD)n 3 2 Dq/C, Dq/A B   It can be shown that the three-term transformation formula [5, (III.33), p. 245] connects (1) (2) (3) (4) (5) Xn with any two of the solutions Xn ,Xn ,Xn and Xn . One such relation works out to be

q A (1),+ Cq (4) (2.16) (ABCλ ,ACλ−, , )∞X (A,Aλ−,ABλ , )∞X + D C n − + D n A Aq (C,Cλ−, C , D ,BCλ+,CDλ+,ABDλ+)∞ (2) = C Xn . ( A ,ADλ+,BCDλ+)∞

Another three-term contiguous relation satisfied by balanced 3φ2 ’s also yields solu- tions to (2.1). The required contiguous relation, which can be deduced from (2.3), (2.4) and (2.10), is (1 a)(1 b)(1 c) de (2.17) − − − (de abcq)φ (1 d)(1 e) abcq − + − − d e +[abc(d + e q) + de(1 + q a b c)]φ + abcq(1 )(1 )φ− = 0. − − − − − q − q

6 Replacing (a,b,c,d,e) by (aq−n, bq−n, cq−n, dq−n, eq−n) and writing

−n −n −n aq , bq , cq de n Z = φ − − ; q , n 3 2 dq n, eq n abc   we have from (2.17) qn+1 qn+1 1 1 qn+1 qn q−n−1(1 )(1 )Z + + + (qn + qn+1 a b c) Z − d − e n+1 d e − de abc − − − n   n 1 n 1 n 1 n de n q (1 a q )(1 b q )(1 c q )(1 abcq q ) (2.18) + − − − − Zn−1 = 0. de (1 1 qn)(1 1 qn) − d − e Choosing the parameters a = q/B, b = q/A, c = q/D, d = Cqλ+ and e = Cqλ− and renormalizing we again obtain the equation (2.1) with a solution

(6) n q n −n(n−1)/2 ABD ABD (2.19) X = ( 1) q λ , λ− n − ABD q + q n −n+1 −n+1 −n+1 q  /B,q /A, q /D n  3φ2 −n+1 −n+1 ; Cq × Cλ+q , Cλ−q   and three similar solutions obtained by parameter interchanges C A or C B or (6) ↔ ↔ C D . However, the solution Xn and its C A, C B or C D interchanges ↔ ↔ ↔ (2) (3)↔ (4) (5) do not give new solutions. They are related to the solutions Xn ,Xn ,Xn and Xn by the transformation formula (2.2b). For example, with the help of this formula we find that (6) (Cq/A,Cq/D,Cq/B)∞ (4) Xn = Xn (C,Cλ+q,Cλ−q)∞ (6) (4) and thus Xn is the same solution as Xn except for a constant factor. We next show that continuous dual q -Hahn polynomials ([3, pp. 3, 28]) are obtained (2) as a particular case of the above solutions. This is true for the solution Xn if C = q or (6) Xn if A,B or D = q . That is why our general case represents associated continuous (2) dual q -Hahn polynomials. In order to show this for Xn we first apply the transformation formula [5, (III.34), p. 245] connecting 3φ2 ’s of type I and II

a,b,c de (e/b, e/c)∞ d/a,b,c 3φ2 ; = 3φ2 ; q d, e abc (e, e/bc)∞ d, bcq/e     (d/a,b,c,de/bc)∞ e/b, e/c, de/abc + 3φ2 ; q . (d, e, bc/e, de/abc)∞ de/bc, eq/bc   We have from (2.11)

−n+1 (2) (A,B)n (ACλ−,ACλ+)∞ q /C,ABλ+,ABλ− (2.20) Xn = 3φ2 ; q (AB)n (Aq/D, C/B)∞ Aq/C,Bq/C    −n+1 −n+1 (q /C,ABλ+,ABλ−,Aq/B)∞ q /B,ACλ−,ACλ+ + − 3φ2 ; q . (Aq/C,Aq/D,B/C,q n+1/B)∞ Aq/B,Cq/B   7 When we write C = q in (2.20), the right side becomes −n (A,B)n (Aqλ−,Aqλ+)∞ q ,ABλ+,ABλ− (2.21) n Aq q 3φ2 ; q . (AB) ( , )∞ A,B D B   In order to compare this with continuous dual q -Hahn polynomials we apply a transfor- mation which will ultimately change the interval of orthogonality from [ 1/α, 1/α] when α = 1 ABCD/q is real, to [ 1, 1] . Take x = αz = cos θ, u = e−iθ , which means 2 − u =2αλ =1/2αλ− . Thus, omitting constant factors we can write (2.21) as p+ −n (A,B)n q ,ABu/2α,AB/2αu (2.22) φ ; q (AB)n 3 2 A,B   which except for a normalization factor of (2α)−n , is the same as the continuous dual q -Hahn polynomials of Askey and Wilson [3]. (1),− We now proceed to show that Xn is a minimal solution of (2.1) for a particular branch in the complex plane. This is done by evaluating the large- n asymptotics of (1),± (1),± Xn and the equation (2.1). Applying the transformation (2.2b) to (2.8) i.e. to Xn we obtain for n →∞ n n 2 (1),± (A,B,C,D)n n (Bq ,ACDλ±q ,ABCDλ±)∞ ± Xn = (λ ) n n (BCDλ±,ABCλ±)n (BCDλ±q ,ABCλ±q ,ADλ±)∞

CDλ±,ACλ±,ADλ± n (2.23) 3φ2 n 2 ; Bq × ACDλ±q ,ABCDλ±  n  const (λ±) . ≈ On the other hand, asymptotics of the second order difference equation (2.1) is given by q (2.24) X zX + X − = 0, n+1 − n ABCD n 1 from which we have for large n n (2.25) X const (λ±) . n ≈ Choosing the square root branch for which λ−/λ+ < 1, the minimal (or the sub- dominant) solution in terms of the variable x =| αz is| given by (min) (1),− (2.26) Xn (x) = Xn n n n (A,B,C,D)n Bq , Cq ,BCλ− = (λ−) 3φ2 n n ; ADλ− (BCDλ−,ABCλ−) BCDλ−q ,ABCλ−q n   1 ABCD 1 valid for zα = x C [ 1, 1], α = , λ± = (x √x2 1) . Summarizing ∈ \ − 2 q 2α ± − what we have done so far, we have q (1),± (k) Theorem 1. The functions Xn , Xn , k =2, 3, 4, 5 of (2.8), (2.11), (2.13), (2.14) and (2.15) respectively are solutions to the recurrence relation (2.1) for associated continuous dual q -Hahn polynomials. These solutions are pairwise linearly independent. The minimal solution of (2.1) is given, up to a multiplicative factor, by (2.26) with the square root branch n 1 chosen so that (λ−/λ ) 0 as n with zα = x C [ 1, 1], α = ABCD/q . + → →∞ ∈ \ − 2 p 8 3. The continued fraction and measure.

The infinite continued fraction associated with (2.1) is

b2 b2 (3.1) CF (z) = z a 1 2 , b2 =0, n > 0. − 0 − z a z a n 6 − 1 − − 2 −··· Pincherle’s theorem [16] connects the minimal solution of (2.1) with the continued fraction (3.1) by the formula

1 X(min)(x) 1 (3.2) = 0 , zα = x, α = ABCD/q. CF (z) 2 (min) 2 b0X−1 (x) p Therefore from Theorem 1, we obtain the continued fraction representation

BCλ−,B,C 3φ2 ; ADλ− 1 ABCDλ− BCDλ−,ABCλ− (3.3) =   . BCDλ− ABCλ− CF (z) q(1 )(1 ) BCλ−,B/q,C/q q q 3φ2 ; ADλ− − − BCDλ−/q,ABCλ−/q   We can also write (3.3) as

BCλ−,B,C 3φ2 ; ADλ− 1 1 q/Aλ+,q/Dλ+ (3.3a) = 1 1   . CF (z) λ+(1 )(1 ) BCλ−,B/q,C/q Aλ+ Dλ+ φ ; ADλ− − − 3 2 1/Aλ , 1/Dλ  + +  (3.3) or the alternative form (3.3a) are valid for

1 zα = x C [ 1, 1], λ−/λ < 1, α = ABCD/q. ∈ \ − | +| 2 p In the particular case C = q (the case of continuous dual q -Hahn polynomials), (3.3) reduces to

1 ABDλ− Bqλ−,B,q (3.4) = 3φ2 ; ADλ− CF (z) (1 BDλ−)(1 ABλ−) DBqλ−,ABqλ− − −   which can also be written with the help of (2.2b) in the form

2 1 ABDλ−(q,ABDλ−,ABDqλ−)∞ BDλ−,ABλ−,ADλ− (3.5) = 3φ2 2 ; q . CF (z) (BDλ−,ABλ−,ADλ−)∞ ABDλ−,ABDqλ−  

with explicit pole terms given by the zeros of the denominator (BDλ−,ABλ−,ADλ−)∞ . These pole singularities and their residues determine the discrete component of the spectral measure of orthogonality for continuous dual q -Hahn polynomials.

9 We now determine the absolutely continuous part of the spectrum for the general associated case. If x = zα ( 1, 1) , then (2.1) has linearly independent solutions given by the boundary values∈ of (2.26)− as x approaches ( 1, 1) from above and below. With now 1 − λ± = [x i√1 x2], x = zα , we have the large n asymptotics 2α ± − 2 (min) n (A,B,C,D)n(ABCDλ−)∞ (3.6) Xn (x + i0) (λ−) , ≈ (BCDλ−,ABCλ−,ADλ−)∞ and 2 (min) n (A,B,C,D)n(ABCDλ+)∞ (3.7) Xn (x i0) (λ+) . − ≈ (BCDλ+,ABCλ+,ADλ+)∞ Since the minimal solution changes as we cross the line segment zα = x ( 1, 1) , we have the representation ∈ − 1 ω(t)dt (3.8) = + possible pole terms, CF (z) R z t/α Z − Also ω(x), x ( 1, 1) can be obtained by using the formula [21] ∈ − (min) (min) 1 W (X− (x + i0),X− (x i0)) (3.9) ω(x) = 1 1 − 2 (min) (min) 2πiα b X− (x + i0)X− (x i0) 0 1 1 − (min) (min) 1 W (Xn (x + i0),Xn (x i0)) = lim − , n→∞ 2 2 2 4 (min) (min) 2πiα b b ...b b X− (x + i0)X− (x i0) 1 2 n 0 1 1 − where W (X ,Y ) = X Y X Y . n n n n+1 − n+1 n Using (2.1), (2.26), (3.6) and (3.7) and simplifying we have

1 2 1 (A,B,C,D)∞( 2 ,u )∞ (3.10) ω(x) = u √ 2 2α 2α 2α 2α 2αq 2αq 2π 1 x ( Au , A u, Du , D u, BCu , BC u)∞ − − B C B C 1 − BCλ , q , q BCλ+, q , q − 3φ2 BCDλ− ABCλ− ; ADλ 3φ2 BCDλ+ ABCλ+ ; ADλ+ . × " q , q ! q , q !# In the particular case C = q this reduces to

1 2 1 (A,B,q,D)∞( 2 ,u )∞ (3.11) ω(x) = u . 2π√1 x2 BD BD AB AB AD AD − ( A /u, A u, D /u, D u, B /u, B u)∞ q q q q q q Taking the appropriate value of the parameters from (2.22) we find that this weight function is the same as the one obtained by Askey and Wilson ([3, p. 11, Theorem 2.2 in the special case d = 0 ]) for continuous dual q -Hahn polynomials. Summarizing the above we have

10 Theorem 2. The associated continuous dual q -Hahn polynomials Pn(x/α) given by (4.12) of the next section are orthogonal with respect to a measure with absolutely con- tinuous component given by the weight function (3.10) on ( 1, 1) . In the particular case C = q this absolutely continuous component reduces to (3.11)−and the discrete spectrum is given by the zeros of (BDλ−,ABλ−,ADλ−)∞ where λ− =(x √x2 1)/2α for x > 1 and x < 1 respectively. ∓ − − The q 1 limit of continuous dual q -Hahn polynomials yields the case of continuous dual Hahn polynomials.→ The corresponding results of this Section for associated continuous dual Hahn polynomials with q = 1 are given by Ismail, Letessier and Valent [13].

11 4. Generating function.

The associated continuous dual q -Hahn polynomials Pn(z; A,B,C,D) satisfy the second order difference equation (2.1) i.e. the equation

1 1 1 1 (4.1) X z + + + qn +(1+ q)q2n−1 X n+1 − − A B C D n     1 n−1 n−1 n−1 n−1 + (1 Aq )(1 Bq )(1 Cq )(1 Dq )X − = 0. 4α2 − − − − n 1

1 A renormalized form of (4.1) is, with x = αz, α = 2 ABCD/q , p 1 1 1 1 (4.2) (1 Aqn)(1 Dqn)ζ 2x 2α + + + qn +2α(1 + q)q2n−1 ζ − − n+1 − − A B C D n     n−1 n−1 + (1 Bq )(1 Cq )ζ − = 0. − − n 1 This is satisfied by the polynomials

n (2α) Pn(z; A,B,C,D) ζn(x; A,B,C,D) = . (A)n(D)n

Let the generating function of the polynomials ζn be

∞ (2α)nP (z; A,B,C,D) (4.3) G(x,t) = n tn (A) (D) n=0 n n X∞ n = ζn(x; A,B,C,D)t . n=0 X An explicit form for the generating function may be obtained by employing the procedure given in [14]. First multiply (4.2) by tn and sum up the resulting equations from n =0 through n = . Using the initial conditions for the polynomial solutions of the first kind, we obtain ∞

t 2αtq A D 2αt 2αt (4.4) (1 ut)(1 )G(x,t) (1 )( + )G(x,tq) − − u − − AD q q − B − C 2αtq 2αtq2 AD A D + (1 )(1 ) G(x,tq2) = (1 )(1 ). − AD − AD q2 − q − q

In (4.4) we put

2αtq (4.5) G(x,t) = AD ∞ f(t) (tu)∞ 12 to obtain

t A D 2αt 2αt (4.6) (1 )f(t) ( + )f(tq) − u − q q − B − C AD A D (tuq)∞ + (1 tuq)f(tq2) = (1 )(1 ) . q2 q q 2αtq − − − AD ∞

∞ n  In the left side of (4.6) we write f(t) = n=0 fnt and on the right side we use q -binomial 2αtq theorem to replace (tuq)∞/ by AD ∞ P ∞ ADu 2αtq n 2α n . (q) AD n=0 n X   If we now equate coefficients of tn on both sides, we obtain the first order difference equation

n−1 n−1 1 2αu n−1 2αu n−1 (4.7) (1 Aq )(1 Dq )f = (1 q )(1 q )f − − − n u − B − C n 1 2αq n A D ADu + (1 )(1 ) 2α n . AD − q − q (q)   n Rewriting (4.7) as

n (A)n(D)n n−1 (A)n−1(D)n−1 − (4.8) u 2αu 2αu fn = u 2αu 2αu fn 1 B n C n B n−1 C n−1 n A D ADu   2αq  ( )n( )n( )n + un q q 2α , AD ( 2αu ) ( 2αu ) (q)   B n C n n the general solution of (4.7) is

2αu 2αu n ( A ) ( D ) ( ADu ) j ( B )n( C )n q j q j 2α j 2αq j (4.9) fn = E + u , un(A) (D)  ( 2αu ) ( 2αu ) (q) AD  n n j=0 B j C j j X     where E is a constant which by the boundary conditions may be taken as 0. Consequently we have the generating function

∞ (2α)nP (z; A,B,C,D) (4.10) n tn (A) (D) n=0 n n X 2αtq ∞ n 2αu 2αu A D ADu j ( )∞ ( )n( )n ( )j( )j( )j 2αq = AD B C q q 2α u−n+jtn, (tu)∞ (A) (D) ( 2αu ) ( 2αu ) (q) AD n=0 j=0 n n B j C j j X X   13 where z = x/α and x = cos θ, u = eiθ . If we interchange C D in (4.10) and write C = q we obtain, with a =2α/B,b =2α/D,c =2α/A , ↔ ∞ (2α)nP (z; A,B,C) n tn (A) (q) n=0 n n X ∞ n (ct)∞ (au) (bu) t (4.11) = n n (tu)∞ (q) (ab) u n=0 n n X   (ct)∞ au, bu t = 2φ1 ; . (tu)∞ ab u   This gives the generating function result for the continuous dual q -Hahn polynomials (see [13, (3.3.7), p. 55]). By comparing coefficients of tn on the left and right sides of (4.10) we obtain an explicit expression for our monic associated continuous dual q -Hahn polynomials.

Pn(z) = Pn(z; A,B,C,D) 2αq n −n 2αu 2αu ℓ u n (A,D, )n (q , , )ℓ AD (4.12) = ADu B C ADu −n 2α (q)n (( 2α q ,A,D)ℓ 2αu   Xℓ=0   ℓ A D ADu j ( , , )j 2αuq q q 2α , × (q, 2αu , 2αu ) AD  j=0 B C j X    1 1 α = ABCD/q, x = (u +1/u), z = x/α . 2 2 p Note that Pn(z) is of course symmetric in the parameters A,B,C,D . The symmetry under the interchanges A D or B C is obvious from (4.12). However, the symmetry ↔ ↔ under the interchanges A B,A C, D B or D C is hidden. Pn(z) is also symmetric under the interchange↔ u↔ u−1↔. Again this↔ is not apparent from (4.12). Applying any one of these hidden symmetry↔ interchanges to (4.12) gives us a type of transformation formula. A different expression for Pn(z) is obtained in Section 7. It is derived from the associated Askey-Wilson polynomial formula of Ismail and Rahman [15]. In order to contrast (7.3) with (4.12) we repeat it here as

n (q−n, BCq u, BCq 1 ) qk (B, C)n AD AD u k (4.13) Pn(z; A,B,C,D) = n (BC)  q (q,B,Cq)k kX=0 n−k A D k+1 k−n j ( q, q , q , q )j BCq k k −n . × (q,Cq , Bq , q )j AD  Xj=0    This formula also does not reveal the full symmetry with respect to A,B,C,D . However, it does make explicit the u u−1 symmetry. ↔ 14 5. Four limiting cases.

We now take successive limits D , C ,B , A . →∞ →∞ →∞ →∞ 5.1 Associated big q -Laguerre ( D ). →∞ The recurrence relation (2.1) becomes 1 1 1 (5.1) L z + + qn +(1+ q)q2n−1 L n+1− − A B C n   n   q n−1 n−1 n−1 (1 Aq )(1 Bq )(1 Cq )L − = 0. − ABC − − − n 1 If we write A = aq, B = bq, C = abq/s and z = t/abq , and renormalize with a factor (abq)n , (5.1) gives the equation− (4.3) of Ismail and Libis [14] for associated big q -Laguerre polynomials in monic form. The solutions of (5.1) are obtained as D limits of solutions of (2.1). We have → ∞ λ z and λ− q/ABCDz and consequently from (2.8) + ≈ ≈ (1) (1),− (5.2) Ln (z; A,B,C) = lim Xn (z; A,B,C,D) D→∞ n(n+1)/2 n n q (A,B,C)n Bq , Cq q = ( 1)n φ ; , − (ABCz)n ( q ) 2 1 qn+1/Az BCz Az n   and from (2.12)

(2) (q/B)∞ (2) (5.3) Ln (z; A,B,C) = lim Xn (z; A,B,C,D) (Aq/B)∞ D→∞ Aqn, q = ( 1)nqn(n−1)/2(A) A−n φ BCz ; Azq−n+1 . − n 2 2 Aq , Aq  B C  (3) (4) Interchanging A B and A C in (5.3) yields the limits of solutions Xn ,Xn . We have ↔ ↔ Bqn, q (5.4) L(3)(z; A,B,C)= ( 1)nqn(n−1)/2(B) B−n φ ACz ; Bzq−n+1 , n − n 2 2 Bq , Bq  C A  and Cqn, q (5.5) L(4)(z; A,B,C) = ( 1)nqn(n−1)/2(C) C−n φ ABz ; Czq−n+1 . n − n 2 2 Cq , Cq  A B  (6) Limits of Xn given by (2.19) and its parameter interchanges will be in terms of 2φ1 and (2) (3) (4) these three limits will simply be transforms of Ln ,Ln and Ln . We write below one of these limits

(5) (6) Ln (z; A,B,C) = lim Xn (z; A,B,C,D) D→∞ q−n+1/A, q−n+1/B (5.6) = zn(1/Cz) φ ; Cqn . n 2 1 Czq−n+1   15 Next we take the D limit of the explicit expression (4.12) for the polynomial solution to get → ∞

(1) (5.7) Pn (z; A,B,C) = lim Pn(z; A,D,C,B) D→∞ q n −n ABCz ℓ ℓ(ℓ−1)/2 ℓ (A,B, ) (q , )ℓ( 1) q AB n ABz n q − = z −n (q)n ( (ABzq ,A,B)ℓ C Xℓ=0   ℓ A B j ( , ,ABz)j Cq q q ( 1)jq−j(j−1)/2 . × ABCz − AB  j=0 ( q , q)j X    Note that before taking limit D of (4.12) we have made the parameter interchange B D in (4.12). →∞ ↔ (1) The minimal solution of (5.1) is Ln and we therefore have the continued fraction representation

B, C φ ; q 1 b2 b2 2 1 q/Az BCz (5.8) 1 2 =   , z a z a z a B/q,C/q − 0 − − 1 − − 2 −··· z(1 1/Az) φ ; q − 2 1 1/Az BCz   1 1 1 a = ( + + )qn (1 + q)q2n−1 n A B C − qn b2 = (1 Aqn−1)(1 Bqn−1)(1 Cqn−1). n −ABC − − − The associated orthogonality is discrete and only explicit in the case A,B or C = q .

16 5.2. Associated Wall (C, D ) . →∞ With C the recurrence relation (5.1) changes to (5.9) →∞ 2n−1 1 1 n 2n−1 q n−1 n−1 W z + q +(1+ q)q W + (1 Aq )(1 Bq )W − = 0. n+1 − − A B n AB − − n 1     Solutions are

(1) (1) (5.10) Wn (z; A,B) = lim Ln (z; A,B,C) C→∞ n n n+1 q n(n−1) (A,B)n Bq q = q 1φ1 n+1 ; , ABz (q/Az)n q /Az Bz    

(2) (2) (5.11) Wn (z; A,B) = lim Ln (z; A,B,C) C→∞ n − − Aq − = ( 1)nqn(n 1)/2(A) A n φ ; Azq n+1 , − n 1 1 Aq/B   and

(3) (3) (5.12) Wn (z; A,B) = lim Ln (z; A,B,C) C→∞ n − − Bq − = ( 1)nqn(n 1)/2(B) B n φ ; Bzq n+1 , − n 1 1 Bq/A   which is just A B interchange of (5.11). ↔ Also from (5.6)

(4) (5) (5.13) Wn (z; A,B) = lim Ln (z; A,B,C) C→∞ − q−n+1/A, q−n+1/B q2n 1 = zn φ ; . 2 0 z  − 

The series representing the above 2φ0 converges only when it terminates. the relevant terminating cases are when A = q or B = q . A limit C of (5.7) gives explicit expression for associated Wall polynomials →∞ (2) (1) (5.14) Pn (z; A,B) = lim Pn (z; A,B,C) C→∞ n −n ℓ(ℓ−1) 2 2 ℓ n (q/ABz,A,B)n (q )ℓq A B z = z −n (q)n ( (q ABz,A,B)ℓ q Xℓ=0   ℓ (A/q,B/q,ABz) q 2j j z−j q−j(j−1) . × (q) AB  j=0 j X    17  The minimal solution of (5.9) is given by (5.10) and therefore we have the associated continued fraction representation

B φ ; q 1 b2 b2 1 1 q/Az Bz (5.15) 1 2 =   , z a z a z a B/q − 0 − − 1 − − 2 −··· z(1 1/Az) φ ; 1 − 1 1 1/Az Bz   where 1 1 a = ( + )qn (1 + q)q2n−1 n A B − 2n−1 q − − b2 = (1 Aqn 1)(1 Bqn 1). n AB − − (2) When A or B = q , (5.15) has associated with it a discrete orthogonality for Pn (z). We note here that our associated Wall polynomials given by (5.9) reduce to the Wall polynomials when B = q . In fact, if we make the substitutions A = aq, B = q, z = x/aq in (5.9) and renormalize, the equation can be written in the form (see [4])

n n+1 n n+1 n n n n q (1 aq )p (x) q (1 aq ) + aq (1 q ) x p (x) + aq (1 q )p − (x)=0 − n+1 − − − − n − n 1 with   n n −n(n−1)/2 (aq) (2) x pn(x) = pn(x; a)=( 1) q Wn ( ; aq,q) − (aq)n aq (2) where we use solution Wn from (5.11). From (5.28) we obtain the standard expression for Wall polynomials (see [17, (3.20.1), p. 83] and [4, p. 198])

q−n, 0 p (x; a) = φ ; qx . n 2 1 aq  

5.3. Limit Wall (B,C,D ) . →∞ The three-term recurrence now becomes n 3n−2 q 2n−1 q n−1 (5.16) U z +(1+ q)q U (1 Aq )U − = 0. n+1 − − A n − A − n 1   Using (5.10), (5.11) and (5.13), we have the solutions

(1) (1) (5.17) Un (z; A) = lim Wn (z; A,B) B→∞ q (A) q2n+1 n n 3n(n−1)/2 n n = ( 1) ( ) q 0φ1 q −+1 ; , − Az ( q ) z Az n  Az  (2) (2) (5.18) Un (z; A) = lim Wn (z; A,B) B→∞ n − (A) Aq − = ( 1)nqn(n 1)/2 n φ ; Azq n+1 , − An 1 1 0   18 and

−n q +1 q2n−1 (3) (4) n A , 0 (5.19) Un (z; A) = lim Wn (z; A,B) = z 2φ0 ; , B→∞ z  −  (3) where Un converges when it terminates with say A = q . (2) A direct limit of (5.14) i.e. Pn (z; A,B) as B leads to an indeterminate form. However we can obtain the explicit form of the polynomials→∞ by applying the method of Section 4 abinitio to equation (5.16). The result is

n2 n −n (3) q (A)n (q )ℓ ℓ −ℓ(ℓ−1)/2 ℓ (5.20) Pn (z; A) = n ( 1) q (Az) A (q)n (A)ℓ − Xℓ=0  ℓ (A/q) − − j qj(j 1)/2( Az) j . × (q) −  j=0 j X  The minimal solution of (5.16) being (5.17) we have the related continued fraction

q 0φ1 − ; 1 b2 b2 1 q/Az z (5.21) 1 2 =   , z a0 z a1 z a2 z 1 − − − − − −··· (1 1/Az)0φ1 − ; − 1/Az qz   A φ ; q 1 1 1 0 Az =   , z A/q φ ; 1 1 1 0 Az   where qn a = (1 + q)q2n−1 n A − q3n−2 b2 = (1 Aqn−1). n − A − The second expression on the right side of (5.21) comes from the transformation

1 z 0φ1 − ; cz = 1φ1 ; c c (c)∞ 0     which can be derived from (5.30) by letting a = 0. Note that when A = q , the right side of (5.21) is a meromorphic function of z . When A q the singularities6 coalesce at z = 0 to produce an essential singularity. →

5.4 A fourth limit (A,B,C,D ) . →∞ 19 The three-term recurrence is now

2n−1 4n−3 (5.22) V [z +(1+ q)q ]V + q V − = 0. n+1 − n n 1 Using (5.17) and (5.19) we have the solutions

2n+1 (1) (1) 2n(n−1) q n q (5.23) Vn (z) = lim Un (z; A) = q ( ) 0φ1 − ; , A→∞ z 0 z  

2n−1 (2) (3) n 0, 0 q (5.24) Vn (z) = lim Un (z; A) = z 2φ0 ; . A→∞ z  − 

(2) The solution Vn is divergent and is thus only a formal solution. The associated polyno- mials are given by

(4) (3) (5.25) Pn (z) = lim Pn (z; A) A→∞

n n2 n(n−1)/2 n ℓ j(j−1) ( 1) q q − − − q − = − (q n) q ℓ(ℓ 1)zℓ (qz) j . (q)  ℓ (q)  n j=0 j Xℓ=0  X 

The minimal solution of (5.22), being given by (5.23), yields the associated continued fraction

0φ1 − ; q/z 1 b2 b2 1 0 (5.26) 1 2 =   z a z a z a z − 0 − − 1 − − 2 −··· φ ;1/qz 0 1 −0   where

2n−1 an = (1 + q)q 2 −4n−3 bn = q .

We can also write (5.26) more explicitly as

2 ∞ n − 5 q n 1 q q 1 n=0 n z (5.27) = (q) . −1 3 ∞ n2−2n z +(1+ q)q z +(1+ q)q z +(1+ q)q z P q z−n − − −··· n=0 (q)n P Using (5.22), (5.23) and (5.26) we have the following

20 Corollary. If 0 < q < 1 and n is an integer, then

q2n+1 fn(z) = 0φ1 − ; 0 z  

has only real simple negative zeros which interlace those of fn+1(z) . Proof. If 0 < q < 1 then (5.26) is a completely convergent positive definite J -fraction which can be represented as a Stieltjes transform of a unique positive discrete measure [21]. This means that (5.26) can have only simple pole sigularities on the real axis with positive residues. Hence we must have simple real intertwining zeros for f0(z) and f−1(z). From k k− ∞ q ( 1) 2n+1 k − the series representation fn(z) = k=0 (q)k (q /z) we see that the zeros of f 1(z) and f (z) must be negative. This establishes the result for n = 1. The proof for other 0 P values of n is the same if one starts from the continued fraction −

1 b2 1 f (z) n+2 = n+1 . z a z a z f (z) − n+1 − − n+2 −··· n

Note 1. There are similar corollaries associated with the positive definite cases for the continued fractions (5.8), (5.15) and (5.21). These require special parameter conditions. See also the remarks after (6.18), (6.28), (6.33) and (6.63).

Note 2. The identities

a,b (b,az)∞ c/b, z (5.28) 2φ1 ; z = 2φ1 ; b c (c, z)∞ az     a,b (az)∞ a, c/b (5.29) 2φ1 ; z = 2φ2 ; bz c (z)∞ c, az     (z)∞ a, 0 a (5.30) 2φ1 ; z = 1φ2 ; cz (az)∞ c c, az     1 z = 1φ1 ; c (c)∞ az   can be used to relate some of the above solutions. For (5.28) and (5.29) see [5, (III.1), (III.4), p. 241]. (5.30) follows from b 0 in (5.28) and (5.29). →

21 6. Additional Limits. There are other less obvious limiting cases which we can obtain from (2.1) and its solutions. These may also be re-expressed as birth and death processes. In the cases of Sections 6.4 and 6.5 there are process based seperately on the even and odd approximants with z2 replaced by z [7, Section 4]. We begin with 6.1 Associated Al-Salam-Chihara. In (2.1) we put D = δC , multiply by C , replace zC by z and renormalize and let C 0 to get → −1 n q n−1 n−1 (6.1) Q (z (1 + δ )q )Q + (1 Aq )(1 Bq )Q − =0 n+1 − − n ABδ − − n 1 with

n (6.2) Qn = lim C Xn (z/C; A,B,C,Cδ) . C→0 This can be recognized as the recurrence relation for associated Al-Salam-Chihara poly- nomials given in [2, (3.54) with a = (1+ γ2b)/γ ] if we make the replacements

(6.3) (A,B,δ−1, z) (bγ/c,γq,bγ2,x) → and renormalize (see also [20]). If A = q or B = q then (6.1) becomes the recurrence for monic Al-Salam-Chihara polynomials [2]. For other references to this case see [17, p. 63]. We record here the solutions to (6.1) based on (6.2) and the solutions we obtained for Xn in Sections 2 and 4. Using (2.8) we have

(1),± n (1),± z (6.4) Qn (z; A,B,δ) = lim C Xn ; A,B,C,δC C→0 C  n (A,B)n n BΛ±, Bq = (Λ±) 2φ1 n ; AδΛ± , (ABΛ±) ABΛ±q n   1 (6.5) Λ± = (z z2 γ2), 2 ± − 1 γ = 2(q/ABδp) 2 .

Also from (2.14), (2.15) and (2.19) respectively we similarly obtain

− BΛ ,BΛ− − (6.6) Q(2)(z; A,B,δ)= (B) B n φ + ; q n+1/B , n n 2 1 q/δ   − BδΛ , BδΛ− − (6.7) Q(3)(z; A,B,δ)= (B) (δB) n φ + ; q n+1/B , n n 2 1 qδ   (6.8) Q(4)(z; A,B,δ)= ( q/ABδ)nq−n(n−1)/2 n − q−n+1/A, q−n+1/B (ABδΛ+/q,ABδΛ−/q)n 2φ2 −n+1 −n+1 ; q/δ , × Λ+q , Λ−q   22 (4) (2) with Qn proportional to Qn via (5.28) and (5.29). Using (4.12) we obtain the explicit polynomial formula (first make the interchange B D ) ↔ n (6.9) Qn(z; A,B,δ) = lim C Pn(z/C; A,δC,C,B) C→0 (A,B) = (γu/2)n n (q)n n (q−n, 2u/γδ, 2u/γ) ℓ ( 1)ℓu−2ℓqnℓq−ℓ(ℓ−1)/2 × ( (A,B)ℓ − Xℓ=0 ℓ (A/q, B/q) ( 1)ju2jqj(j+1)/2 j − , × (q, 2u/γδ, 2u/γ)  j=0 j X  −1 1 z = γ(u + u )/2, γ = 2(q/ABδ) 2 . 

With the choice of square root branch chosen so that Λ−/Λ+ < 1 for x = z/γ C [ 1, 1], we have the minimal solution to (6.1) given by| (6.5). As| a consequence we may∈ give\ − an explicit expression for the corresponding continued fraction and the absolutely con- tinuous component of the measure which gives its representation as a Stieltjes transform. The calculations proceed as in Section 3 and yield the following. For z/γ = x C [ 1, 1] ∈ \ − and Λ−/Λ < 1, | +|

23 BΛ−,B 2 2 2φ1 ; AδΛ− 1 b b ABδΛ− ABΛ− (6.10) 1 2 =   z a0 z a1 z a2 q(1 ABΛ−/q) BΛ−,B/q − − − − − −··· − 2φ1 ; AδΛ− ABΛ−/q   1 ω(t)dt = + possible pole terms, − z γt Z 1 − with

2 2 1 (A,B, 1/u ,u )∞ (6.11) ω(x) = 2π√1 x2 (Aδγu/2,Aδγ/2u,ABγu/2q,ABγ/2qu)∞ − −1 BΛ−,B/q BΛ+,B/q 2φ1 ; AδΛ− 2φ1 ; AδΛ+ . × ABΛ−/q ABΛ+/q      Note that (6.11) agrees with the weight function derived in [2, (3.64) with a = q(1 + δ−1)/B,b = q2/B2δ, c = q/ABδ, γ = B/q and x2 replaced by x2/4c ].

6.2 Associated Al-Salam-Carlitz I. We take the B limit of (6.1) and its corresponding solutions to obtain the case of associated Al-Salam-Carlitz→∞ I. The recurrence becomes

n −1 n q n−1 (6.12) R [z (1 + δ )q ]R (1 Aq )R − = 0 n+1 − − n − Aδ − n 1 with solutions from (6.5) – (6.8) given by

(1) (1),− (6.13) Rn (z; A, δ) = lim Qn (z; A,B,δ) B→∞

(A)n − q/Azδ = ( q/Aδz)nqn(n 1)/2 φ ; qn+1/z , (q/δz) − 1 1 qn+1/zδ n   (2) (2) (6.14) Rn (z; A, δ) = lim Qn B→∞ − q/Azδ − = ( 1)nqn(n 1)/2 φ ; zq n+1 , − 1 1 q/δ   (3) (3) (6.15) Rn (z; A, δ) = lim Qn B→∞ − − q/Az − = ( δ) nqn(n 1)/2 φ ; δzq n+1 , − 1 1 qδ   (4) (4) (6.16) Rn (z; A, δ) = lim Qn B→∞ q−n+1/A = (1/z) zn φ ; q/δ , n 1 1 zq−n+1   24 (4) (2) with Rn proportional to Rn via (5.29) and (5.30); namely c/b (bz)∞ z (6.17) 1φ1 ; bz = 1φ1 ; c . c (c)∞ bz     The minimal solution to (6.12) is given by (6.13). Using Pincherle’s theorem we obtain the continued fraction representation 1 q(1 A)/Aδ q2(1 Aq)/Aδ (6.18) − − z (1 + δ−1) + z (1 + δ−1)q + z (1 + δ−1)q2 + − − − ··· q/Azδ φ ; q/z 1 1 1 q/zδ =   . z(1 1/δz) q/Azδ − φ ;1/z 1 1 1/zδ   This is a positive definite J -fraction in the case when 0 < q < 1 and A < 1,Aδ < 0. We may then deduce that the zeros of the 1φ1 ’s in the numerator and denominator on the right side of (6.18) are real and simple and interlace. See for example the Corollary in Section 5.4. When A = q , the pole singularities in (6.18) become explicit, since we may then use (6.17) to obtain 1 (1 q)/δ q(1 q2)/δ (6.18′) − − z (1 + δ−1) + z (1 + δ−1)q + z (1 + δ−1)q2 + − − − ··· (q/δz)∞ 1/zδ = 1φ1 ; q/z z(1/z)∞(1/zδ)∞ q/zδ ∞   qn qn = − + . (z qn)(q) (δq) (δ 1)∞ (z qn/δ)(q) (q/δ) (δ)∞ n=0 n n n n X  − −  In the last equality we have assumed that δ = q−m, m an integer. The explicit polynomial solution to (6.12) can be obtained from the 6 B limit of (6.9). We get →∞ (6.19) Rn(z; A, δ)= lim Qn(z; A,B,δ) B→∞ (A) n (q−n, 1/zδ, 1/z) = ( q/Aδz)n n qn(n−1)/2 ℓ − (q)n × ( (A)ℓ Xℓ=0 2 ℓ (A/q) qj (Aδz2)−j q−ℓ(ℓ−1)(Aδz2/q)ℓqnℓ j . × (q, 1/zδ, 1/z)  j=0 j X  (4) When A = q , the expression for Rn(z; q, δ) must be equal to Rn (z; q, δ) given by (6.16), since they are both monic polynomial solutions. Equating these two expressions we obtain −n − − q , 1/zδ, 1/z (6.20) R (z; q, δ) =( δz) nqn(n 1)/2 φ ; δz2qn n − 3 0  −  q−n = zn(1/z) φ ; q/δ . n1 1 zq−n+1   25 The above connection between a terminating 3φ0 and a terminating 1φ1 appears to be new. Note that (6.20) differs also from the standard expression in [17, (3.24.1) with a = δ−1 and x replaced by z ] which gives

−n −1 − − q , z (6.21) R (z; q, δ) =( δ) nqn(n 1)/2 φ ; qzδ n − 2 1 0   −n n q = z (1/zδ) φ − ; qδ . n 1 1 q n+1zδ   For the last equality we have used (5.29) with c =0. 6.3 Limit Al-Salam-Carlitz I. We now take the A limit of (6.12) and its corresponding solutions. The recur- rence becomes →∞

2n−1 −1 n q (6.22) S [z (1 + δ )q ]S + S − = 0 n+1 − − n δ n 1 with solutions

(1) (1) (6.23) Sn (z; δ) = lim Rn (z; A, δ) A→∞ 2 qn 0 = φ ; qn+1/z , (q/δz) (δz)n 1 1 qn+1/zδ n   (2) (2) (6.24) Sn (z; δ) = lim Rn (z; A, δ) A→∞ − 0 − = ( 1)nqn(n 1)/2 φ ; zq n+1 , − 1 1 q/δ   (3) (3) (6.25) Sn (z; δ) = lim Rn (z; A, δ) A→∞ − − 0 − = ( δ) nqn(n 1)/2 φ ; δzq n+1 , − 1 1 qδ   (4) (4) (6.26) Sn (z; δ) = lim Rn (z; A, δ) A→∞

n 0 = (1/z) z φ − ; q/δ . n 1 1 zq n+1   (4) (2) Note that Sn is proportional to Sn via the transformation (6.17) which yields the identity

0 (z)∞ 0 (6.27) 1φ1 ; z = 1φ1 ; c . c (c)∞ z     (1) The minimal solution to (6.22) is Sn (z; δ) . Using Pincherle’s theorem we then obtain the evaluation of the continued fraction associated with (6.22). Namely

26 1 q/δ q3/δ (6.28) z (1 + δ−1) z (1 + δ−1)q z (1 + δ−1)q2 − − − − − −··· 0 φ ; q/z 1 1 1 q/zδ =   . z(1 1/δz) 0 − φ ; q/z 1 1 1/zδ   This is a positive definite J -fraction if 0 < q < 1 and δ > 0 and we may then deduce that the zeros of the 1φ1 ’s on the right of (6.28) are real, simple and interlacing. See the Corollary in Section 5.4. The monic polynomials are given explicitly by

(6.29) Sn(z; δ) = lim Rn(z; A, δ) A→∞ −n n2 n (zδ) q −n ℓ −3ℓ(ℓ−1)/2 = (q , 1/zδ, 1/z)ℓ( δ) q (q)n ( − Xℓ=0 ℓ −2j −j 3j(j−1)/2 − z ( δ) q z2ℓqℓ(n 1) − . × (q, 1/zδ, 1/z)  j=0 j X  A simpler expression is obtained if one applies the generating function method of Section 4 directly to (6.22). This results in

−n n(n+1)/2 n δ q −n ℓ −ℓ(ℓ−1)/2 (6.30) Sn(z; δ) = (q )ℓ(1/z)ℓ( δz) q (q)n ( − Xℓ=0 ℓ ( zδ)−j qj(j−1)/2 − . × (1/z) (q)  j=0 j j X  6.4 Associated continuous q -Hermite.  1 1 If we multiply (6.1) by B 2 , replace zB 2 by z , renormalize and let B 0, we get the recurrence → q n−1 (6.31) H zH + (1 Aq )H − = 0. n+1 − n Aδ − n 1 This will become the continuous q -Hermite case if we set A = q . The solutions to (6.31) are obtained as limits of the solutions to (6.1). They are

± ± 1 (1), n/2 (1), 2 (6.32) Hn (z; A, δ) = lim B Qn (z/B ; A,B,δ) B→0 n n Aq 2 = (A) (µ±) φ ; Aδµ , n 1 1 0 ±   −n+1 (2) n q /A, 0 n 2 (6.33) Hn (z; A, δ) =(µ−) 2φ0 ; q /δµ− ,  −  27 where

′ ′ (6.34) µ± = (z z2 γ 2)/2, γ = 2 q/Aδ. ± − q To derive (6.32) we first made the transformation (5.28) befop re taking the limit and the transformation (5.30) in the form 1 z (6.35) 0φ1 −; cz = 1φ1 ; c c (c)∞ 0     after the limit was taken. (2) (3) To derive (6.33) we transformed either Qn or Qn using an iterate of (5.28), dis- carded factors which were n -independent, multiplied by Bn/2 and then let B 0. Note that (6.33) is only a formal solution unless it terminates by having A = q . → For the general polynomial solution we take the limit of (6.9) to obtain n/2 (6.36) Hn(z; A, δ) = lim B Qn(z; A,B,δ) B→0 n −n ′ (A) (q ) − = (γ u/2)n n ℓ ( 1)ℓu 2ℓqnℓ (q)n ( (A)ℓ − Xℓ=0 ℓ j 2j j(j+1)/2 − − (A/q) ( 1) u q q ℓ(ℓ 1)/2 j − , × (q)  j=0 j X  ′ γ −1 ′ 1 z = (u + u ), γ = 2(q/Aδ) 2 .  2 Note that when A = q , the expressions (6.33) and (6.36) become equal. ′ For z = γ x, x C [ 1, 1] and with the square root branch chosen so that µ−/µ+ < (1),− ∈ \ − | | 1 , we have Hn (z; A, δ) as a minimal solution to (6.31) and hence the continued fraction representation 1 q(1 A)/Aδ q(1 Aq)/Aδ q(1 Aq2)/Aδ (6.37) − − − z z z z − − − −··· A 2 1φ1 ; Aδµ− Aδµ− 0 =   q A/q φ ; Aδµ2 1 1 0 −   1 ω(t)dt = ′ + possible pole terms. − z γ t Z 1 − Repeating the method of Section 3, we find that 2 −2 1 (A,u ,u )∞ (6.38) ω(x) = , 2π√1 x2 A/q A/q φ ; u2/q φ ;1/u2q − 1 1 0 1 1 0     ′ 1 − z/γ = x = (u + u 1). 2 28 When A = q the denominator 1φ1 ’s in (6.37) and (6.38) become equal to 1 and there are no pole terms. This is then the continuous q -Hermite case. See [17, p.88] for a list of references.

6.5 Limit q -Hermite. We take the A limit of the results in the previous Section 6.4. This gives the recurrence relation → ∞ qn (6.39) T zT T − = 0 n+1 − n − δ n 1 with minimal solution

(1) (1),− (6.40) Tn (z; δ) = lim Hn (z; A, δ) A→∞ − = ( 1)nqn(n 1)/2(q/δz)n φ −; qn+2/δz2 − 0 1 0   and the polynomial solution

(6.41) Tn(z; δ) = lim Hn(z; A, δ) A→∞ −n n(n−1)/2 n ( z) q n −n = − (q/δ) (q )ℓ (q)n ( Xℓ=0 2 ℓ z−2j δ−jqj z2ℓqnℓ(δ/q)ℓq−ℓ(ℓ−1) . × (q)  j=0 j X  Using the minimal solution we get the continued fraction result 

2 − 2 2 1 q/δ q /δ 1 0φ1 0 ; q /δz (6.32) = − . z + z + z + z φ ; q/δz2 ··· 0 1 0  This is a positive definite J -fraction if 0 < q < 1 and δ < 0 and one may then has simple real interlacing zeros for the 0φ1 ’s on the right of (6.32) (see the Corollary in Section 5.4).

6.6 Associated continuous big q -Hermite.

In (6.1) let δ =1/aB and then let B 0 to get the recurrence → n aq n−1 (6.43) C (z q )C + (1 Aq )C − = 0. n+1 − − n A − n 1 The solutions to (6.43) may be obtained from the solutions of (6.1) by using

(6.44) Cn(z; A,a) = lim Qn(z; A,B, 1/aB). B→0

29 First making the parameter interchange A B in (6.5), then writing δ = 1/aB and taking limit B 0 we obtain the solution ↔ → ′ n ′ (1),± (A,B)n ′ n AΛ±, Aq Λ± (6.45) Cn (z; A,a) = lim ′ (Λ±) 2φ1 ′ n ; , B→0 (ABΛ ) ABΛ±q a ± n   ′ n ′ ′ n AΛ±, Aq Λ± = (A) (Λ±) φ ; , n 2 1 0 a   ′ 1 2 2 Λ± = (z z γ ) 2 ± − 1 aq q γ = 2 . 1 A r Also from (6.2) we similarly obtain

(2) (2) (6.46) Cn (z; A,a) = lim Qn (z; B,A, 1/aB) B→0 −n+1 ′ ′ q /A, 0 ′ = const (Λ )n(AΛ /aq) φ ; AΛ + − n 2 1 Λ′ q−n+1 −  +  which is a polynomial solution when A = q . In (6.7) we first apply the transformation [5,(III.2), p. 241], write δ =1/aB, B = qm and let m . We obtain the solution →∞ ′ −n+1 ′ 2 2 n−2 ′ q /A, Λ /a A Λ− q (6.47) C(3)(z; A,a) = const (Λ )n φ + ; , n + 2 0 a − ! which is again a polynomial solution for A = q . Using (6.8), transformation (5.30) and taking limits we have the solution

−n+1 (4) ′ n ′ q /A, 0 ′ (6.48) C (z; A,a) = const (Λ ) (AΛ−/aq) φ ′ − ; AΛ− . n + n × 2 1 Λ q n+1  +  Also from (6.9) we obtain the explicit polynomial solution

1 (6.49) Cn(z; A,a) = lim Qn(z; A,B, ) B→0 aB

n n n (A)n (q , 2u/γ1)ℓ ℓ −2ℓ nℓ −ℓ(ℓ−1)/2 = (γ1u/2) ( 1) u q q (q)n ( (A)ℓ − Xℓ=0 ℓ (A/q) ( 1)ju2jqj(j+1)/2 j − , × (q, 2u/γ )  j=0 1 j X  γ1 = 2 aq/A.  p 30 1 If in (6.43) we write z = xγ1 = 2(aq/A) 2 x , and renormalize, we obtain the recur- rence for the associated continuous big q -Hermite polynomials (see [17, (3.18.4)]) viz., the relation

n n−1 (6.50) ξ (2x bq )ξ + (1 Aq )ξ − = 0 n+1 − − n − n 1

1 with b = (A/aq) 2 . Writing A = q in the renormalized solution given by (6.47) yields continuous big q -Hermite polynomials [17, (3.18.1)]. C ′ ′ Choosing, for x = z/γ1 [ 1, 1] , the square root branch for which Λ−/Λ+ < 1, the minimal solution of (6.43)∈ is\ given− by (6.45). Consequently we have| the continued| fraction representation

′ A,AΛ− ′ ′ 2φ1 ; Λ−/a 1 b1 b2 AΛ− 0 (6.51) =  ′  z a0 z a1 z a2 aq A,AΛ− ′ − − − − − −··· φ ; Λ−/a 2 1 0   1 ω(t)dt = + possible pole terms, − z γ t Z 1 − 1 with (6.52) 2 2 ′ ′ ′ ′ −1 1 (A)∞(1/u ,u )∞ A/q,AΛ− Λ− A/q,AΛ+ Λ+ ω(x) = 2φ1 ; 2φ1 ; 2π√1 x2 (γ1u/2a,γ1/2au)∞ 0 a 0 a −      where 1 x = (u + u−1) 2 1 γ1 = 2(aq/A) 2 n an = q qa b2 = (1 Aqn−1). n A − In the particular case A = q ,

2 2 1 (q)∞(1/u ,u )∞ (6.53) ω(x) = . 2π√1 x2 (u/√a, 1/u√a)∞ − This weight function agrees with the form given in [17, (3.18.2] with b in (6.50) replaced by 1/√a and taking into consideration the normalization factor (A/aq)n/2 = a−n/2 = bn used in (6.50) and also that x is replaced by bx .

31 6.7 q -Bessel order.

The limit A of the recurrence (6.43) gives →∞ n n (6.54) B (z q )B aq B − = 0. n+1 − − n − n 1 It is clear that (6.54) will give real orthogonal polynomials only for 0 < q < 1, a < 0 . The solutions of (6.54) can be obtained by taking A limits of (6.45)- (6.48). We have →∞

(1) (1),− n n(n−1)/2 aq/z n+1 (6.55) Bn (z; a) = lim Cn (z; A,a)= ( aq/z) q 1φ1 ; q /z A→∞ − 0   n n(n−1)/2 n+1 n+2 2 (6.56) = ( aq/z) q (q /z)∞ 0φ1 − ; aq /z , − qn+1/z   using the transformation (6.35). Also

(2) (2) n 0, 0 (6.57) Bn (z; a) = lim Cn (z; A,a) = const z (1/z)n2φ1 −n+1 ; aq/z , A→∞ zq  

(3) (3) n 0,z/a n 2 (6.58) Bn (z; a) = lim Cn (z; A,a) = const z 2φ0 ; aq /z , A→∞  − 

(4) (2) (6.59) lim Cn (z; A,a) = const Bn (z; a). A→∞

Using (6.49) we obtain the explicit polynomial

(6.60) Bn(z; a) = lim Cn(z; A,a) A→∞ n n ( a/z) n(n+1)/2 −n = − q (q , 1/z)ℓ (q)n ( Xℓ=0 2 ℓ ℓ j2 2 j − 2 z q (a/z ) q ℓ qnℓ . × a (q) (1/z)  j=0 j j   X  We now demonstrate the relationship of the above solutions with Jackson’s q -analogues of Bessel functions (see [5, Exercise 1.24, p. 25]). Using the notation of [5] and writing

qν = z−1, x2/4 = aq/z, − we find from (6.57) and (6.56) that

2 (1) 1 − (q /z)n 2 n n(n+4)/2 n+ν (2) (6.61) J−ν−n(2i(aq/z) ; q) = const ( 1) q √z( i z/a) Bn (z; a) − − (1/z)n p 32 and

−n(n−1)/2 (2) 1 n q n n+ν (1) (6.62) J (2i(aq/z) 2 ; q) = const ( 1) (z/a√q) (i a/z) B (z; a). ν+n − √z n p This shows that (6.54) is connected with the recurrence for q -Bessel functions with z appearing both in the arguement and the order of the q -Bessel function. Thus we choose to call this the q -Bessel order case. The continued fraction representation obtained with the help of minimal solution of (6.54) given by (6.55) is

2 2 0φ1 − ; aq /z 1 aq aq2 1 q/z (6.63) =  . z 1 + z q + z q2 + (z 1) − − − ··· − φ ; aq/z2 0 1 1−/z  

This is a positive definite J -fraction if 0 < q < 1 and a < 0 and the zeros of the 0φ1 ’s on the right are then real simple and interlacing. A q 1 limiting case is due to Maki [19] (see also [21]). → 7. Limit Askey-Wilson.

In this section we give the connection between solutions to the associated Askey- Wilson and the associated continuous dual q -Hahn polynomial recurrence relations. The associated Askey-Wilson polynomial recurrence relation in monic form is given by

2 (7.1) p (x) (x a )p (x) + b p − (x) = 0, n+1 − − n n n n 1 a 1 a = A B + + , n − n − n 2 2a 2 bn = An−1Bn,

(1 abcdqn+α−1)(1 abqn+α)(1 acqn+α)(1 adqn+α) A = − − − − , n 2a(1 abcdq2n+2α−1)(1 abcdq2n+2α) − − a(1 qn+α)(1 bcqn+α−1)(1 bdqn+α−1)(1 cdqn+α−1) B = − − − − . n 2(1 abcdq2n+2α−2)(1 abcdq2n+2α−1) − − When α = 0 , (7.1) reduces to the non-associated monic Askey-Wilson case [3]. If one further puts d = 0 , one has the monic continuous dual q -Hahn case [3, pp.3, 28]. Our associated dual q -Hahn case may be obtained from (7.1) by multiplying (7.1) by k−1 =2 q/ABCD, replacing (a,b,c,d,qα) by (2kq/AD, 2kq/AC, 2kq/AB, 0,A/q) and renormalizing to monic form. Note that ‘ k ’ is now used to denote the ‘ α ’ of previous p Sections so as to avoid confusion with the ‘ α ’ in (7.1). There are two papers [15], [9], which deal with the associated Askey-Wilson polynomial case. We now indicate how our associated dual q -Hahn recurrence solutions in Section 2 are connected with the solutions in [15] and [9].

33 Using the solutions to (7.1) which Ismail and Rahman have given in [15], we have the following limiting cases.

(B, C)n 2kq 2kq 2kq (7.2) Pn(z; A,B,C,D) = lim pn(x; , , ,d,A/q) d→0 kn AD AC AB where 1 p (x) = p (x; a,b,c,d,qα) = p(α)(x) n n (2a)n n

(α) and pn (x) is given explicitly in [15, (4.15)]. The calculation yields

n −n −1 k (B, C)n (q , BCq/ADu, BCq/ADu )k q (7.3) Pn(z; A,B,C,D) = n (BC) (q,B,C)k (k=0 p p X − n k (A/q,D/q,qk+1, qk−n) BCq j j . × (q,Cqk, Bqk, q−n) AD  j=0 j X   

Ismail and Rahman have also obtained non-polynomial solutions rn+α(u) and sn+α(u) (α) which correspond to pn (x). Using [15, (1.12), (1.13)] we obtain

(1),± ± (B, C)n −α ±1 (7.4) Xn (z; A,B,C,D) = lim C1 (BCλ±) sn+α(u ) d→0 (BC)n  where ± (A,D,BCλ±)∞ C1 = , (BCDλ±,ABCλ±,ADλ±)∞ and

(5) (B, C)n −α (7.5) Xn (z; A,B,C,D) = lim C5 (BCλ−) rn+α(u) d→0 (BC)n  where (C,ADλ+,Dqλ−,ABCλ+/q)∞ C5 = . (B,C,q/B,Dq/C)∞ In [9], the parameters (a,b,c,d,qα) of[15]arereplaced by (α,β,γ,δ,ǫ) . The solutions (r) (r) Xn to (7.1) in [9] will now be denoted by Xn . The δ 0 limits of some of these (1) →(2) (4) (5) (6) solutions are given below in terms of the solutions Xn ,Xn ,Xn ,Xn and Xn of Section 2. We have e

(6) ±1 ± n (1),± (7.6) lim Xn (u ) = D6 k Xn (z; A,B,C,D), δ→0

e ± (BCDλ±,ADλ±,ABCλ±)∞ D6 = ± (A,B,C,D,q/u 2)∞

34 (5) n (5) ′ n (2) (7.7) lim Xn (u) = D5k Xn (z; A,B,C,D) + D5k Xn (z; A,B,C,D), δ→0 A 2 q2 q Dq q e ( q u , Au2 , D ,ACλ−,ABλ−, C , B )∞ D5 = ABC ABD ACD q ( q λ+, q λ+, q λ+,Cqλ−,Bqλ−, u2 )∞ 1 A , × ( D ,Dqλ−,ADλ+)∞ A 2 q2 q Dq ′ ( q u , Au2 , A ,Aqλ−,BDλ−,CDλ−, 2k )∞ D5 = ABC ABD ACD q ( q λ+, q λ+, q λ+,Cqλ−,Bqλ−, u2 )∞ q2 q Aq 2k Du ( , , , , )∞ Du B C D q . D q2 Aq Au 2k × − − ∞ (Dqλ , A , Au , 2k , q , A ,Aqλ ,ADλ+)

(2) n (6) (7.8) lim Xn (u) = D2k Xn (z; C,B,A,D), δ→0 2 2 (A qλ ,Aqλ−)∞ e D = + . 2 Aq Aq Aq BCDλ− u ∞ ( B , C , D , q )

(1) n (5) ′ n (1) (7.9) lim Xn (u) = D1k Xn (z; A,B,C,D) + D1k Xn (z; A,B,C,D), δ→0

(q/B,BCDqλ+/A, BCλ−,Dq/C)∞ e D1 = (D,Cq/A,Bq/A,Dqλ+,ABCλ−/q,ADλ−)∞ q BCq 2k Dq BCDq 2k ( , , , ,BDλ+,CDλ+, λ+,BCDλ−, )∞ D′ = A 2k BC 2k A D . 1 Dq Cq Bq q uq2 A ∞ (B,C,D, A , A , A ,u, u , A , uq , Dλ+)

Parameter interchanges in [9] in some of the cases yield new solutions of (7.1). For ≈ (1) (1) example an interchange of δ β in solution Xn (u) gives a new solution Xn (u) for which we have ↔ e ≈ (1) ′′ n (4) (7.10) lim Xn (u) = D1 k Xn (z; A,B,C,D), δ→0 ′′ (q/B,Cq/D)∞ D1 = . (C,Cqλ+,ACλ−,ABDλ−/q)∞

(1) If we instead make the interchange α δ in solution Xn (u) , then the δ 0 limit n (5) ↔ → would be proportional to k Xn (z; A,B,C,D). e

35 REFERENCES 1. R. Askey, Continuous Hahn polynomials, J. Phys. A. 18 (1985), L1017-L1019. 2. R. Askey and M.E.H. Ismail, Recurrence relations, continued fractions and orthogonal polynomials, Mem. Amer. Math. Soc., 300 (1984). 3. R. Askey and J. Wilson, Some basic hypergeometric orthogonal polynomials, Mem. Amer. Math. Soc., 319 (1985). 4. T.S. Chihara, An Introduction to Orthogonal Polynomials, Gordon and Breach, New York, 1978. 5. G. Gasper and M. Rahman, Basic Hypergeometric Series, Cambridge University Press, Cambridge, 1990. 6. W. Gautschi, Computational aspects of three-term recurrence relations, SIAM Review 9 (1967), 24-82. 7. D.P. Gupta, M.E.H. Ismail and D.R. Masson, Associated continuous Hahn Polynomi- als, Canadian J. of Math. 43 (1991), 1263-1280. 8. D.P. Gupta, M.E.H. Ismail and D.R. Masson, Contiguous relations, basic hypergeo- metric functions and orthogonal polynomials II. Associated big q -Jacobi polynomials, J. Math. Anal. Appl. 171 (1992), 477-497. 9. D.P. Gupta and D.R. Masson, Solutions to the associated q-Askey-Wilson polynomial recurrence relation in: Approximation and Computation 158, Birkhauser, Cambridge, 1994, to appear. 10. W. Hahn, Uber¨ Orthogonalpolynome, die q -Differenzengleichungen gen¨ugen, Math. Nachr. 2 (1949), 4-34. 11. W. Hahn, Uber¨ Polynome, die gleichzeitig zwei verschiedenen Orthogonalsystemen angeh¨oren, Math. Nachr. 2 (1949), 263-278. 12. M.E.H. Ismail, J. Letessier, D.R. Masson and G. Valent, Birth and death processes and orthogonal polynomials in: Orthogonal Polynomials, P. Nevai, ed., Kluwer Academic Publishers, Dordrecht, 1990, 229-255. 13. M.E.H. Ismail, J. Letessier and G. Valent, Quadratic birth and death models and associated continuous dual Hahn polynomials, SIAM J. Math. Anal. 20 (1989), 727- 737. 14. M.E.H. Ismail and C.A. Libis, Contiguous relations, basic hypergeometric functions and orthogonal polynomials I, J. Math. Anal. Appl. 141 (1989), 349-372. 15. M.E.H. Ismail and M. Rahman, The associated Askey-Wilson polynomials, Trans. Amer. Math. Soc. 328 (1991), 201-237. 16. W.B. Jones and W.J. Thron, Continued Fractions: Analytic Theory and Applications, Addison-Wesley, Reading, Mass., 1980. 17. R. Koekoek and R.R. Swarttouw, The Askey-scheme of hypergeometric orthogonal polynomials and its q -analogue, Reports of the Faculty of Technical Mathematics and Informatics no. 94-05, Delft, 1994. 18. C.A. Libis, Private communication. 19. D. Maki, On constructing distribution functions with applications to Lommel polyno- mials and Bessel functions, Trans. Amer. Math. Soc. 130 (1968), 281-297. 20. D.R. Masson, Difference equations revisited in: Canadian Mathematics Society Con- ference Proceedings Vol. 9, J.S. Feldman and L.M. Rosen, eds., AMS, Providence RI,

36 1988, 73-82. 21. , Difference equations, continued fractions, Jacobi matrices and orthogonal polynomials in: Non-linear Numerical Methods and Rational Approximation, 1987, A. Cuyt, ed., D. Reidel, Dordrecht, 1988, 239-257. 22. , Associated Wilson polynomials, Constructive Approximation 7 (1991), 521- 534. 23. , Wilson polynomials and some continued fractions of Ramanujan, Rocky Mountain J. Math. 21 (1991), 489-499. 24. H. S. Wall, Analytic Theory of Continued Fractions, D. Van Nostrad, Princeton, NJ, 1948. 25. J. Wilson, Hypergeometric series, recurrence relations and some orthogonal functions, Doctoral Dissertation, University of Wisconsin, Madison, 1978.

37