arxiv:math/ v1 [math.ca] 7 May 2002

Transcription

1 Large Parameter Cases of the Gauss Hypergeometric Function Nico M. Temme arxiv:math/565v [math.ca] 7 May ABSTRACT CWI P.O. Box GB Amsterdam The Netherlands nicot@cwi.nl We consider the asymptotic behaviour of the Gauss hypergeometric function when several of the parameters a b c are large. We indicate which cases are of interest for orthogonal polynomials Jacobi but also Krawtchouk Meixner etc. which results are already available and which cases need more attention. We also consider a few examples of 3 F functions of unit argument to explain which difficulties arise in these cases when standard integrals or differential equations are not available. Mathematics Subject Classification: 33C5 33C45 4A6 3C5 4A. Keywords & Phrases: Gauss hypergeometric function asymptotic expansion. Note:. Introduction Work carried out under project MAS. Analysis Asymptotics and Computing. This report has been accepted for publication in the proceedings of The Sixth International Symposium on Orthogonal Polynomials Special Functions and their Applications 8- June Rome Italy. The Gauss hypergeometric function c ; z = + ab c z + aa+bb+ z +... = cc+! where a n = Γa+n Γa is defined for z < and c.... A simplification occurs when b = c: b ; z = n n! n= a n b n c n n! z n. a. n = z a..3 When a or b are non-positive integers the series in. will terminate and F reduces to a polynomial. We have nb n n m b m n n ; z = z m = m bm z m..4 c c m m! m c m k= m=

2 The value at z = is defined when Rc a b > and is given by c ; = ΓcΓc a b Γc aγc b..5 For the polynomial case we have nb ; = c b n n = c c n Generalizing let pq =... with p q +. Then a...a p pf q ; z = b...b q This series converges for all z if p < q + and for z < if p = q +. n= a n a p n z n b n b q n n!..7 Sometimes we know the value of a terminating function at z = Saalschütz s theorem: where n =.... n a b 3F a +a+b c n ; = c a nc b n c n c a b n.8 The behaviour of the Gauss hypergeometric function F;c;z for large z follows from the transformation formula c ; z = ΓcΓb a aa c+ ΓbΓc a z a a b+ ; z + ΓcΓa b ΓaΓc b z b bb c+ b a+ ; z.9 where ph z < π and from the expansion c ; z = + ab c z + aa+bb+ z +... cc+! a n b n = z n c n n! n=. with a n = Γa+n Γa = n n! and where z <. We always assume that c.... a n =.... n The asymptotic behaviour of F;c;z for the case that one or more of the parameters a b or c are large is more complicated except when only c is large. Several contributions in the literature are available for the asymptotic expansions of functions of the type Fa+e λb+e λ;c+e 3 λ;z e j = ± λ..

3 In this paper we give an overview of these results and we show how this set of 6 cases can be reduced by using several types of transformation formulas. In particular we indicate which cases are of interest for orthogonal polynomials that can be expressed in terms of the Gauss hypergeometric function. The Jacobi polynomial P n αβ x is an important example and we mention several cases of this polynomial in which n is a large integer and α and/or β are large. When α or β are negative the zeros of P n αβ x may be outside the interval [ ] and for several cases we give the distribution of the zeros in the complex plane. We mention a few recent papers where certain uniform expansions of the Gauss function are given in terms of Bessel functions and Airy functions from which expansions the distribution of the zeros can be obtained. Many cases need further investigations. Of a completely different nature is the asymptotics of generalised hypergeometric functions. We consider a few cases of 3 F terminating functions of argument and and show for these cases some ad hoc methods. In general no standard methods based on integrals or differential equations exist for these quantities.. Asymptotics: a first example Consider the asymptotics of aβ +λ γ +λ ; z λ.. We use β +λ γ +λ and we try using.3 aβ +λ aβ +λ γ +λ ; z F β +λ ; z = z a?. Observe first that if Rγ a β > aβ +λ γ +λ ; Γγ +λγγ a β = Γγ +λ aγγ β..3 We see that. cannot hold for z if Rγ a β >. However we can use c ; z = z a ac b ; c 3 z..4 z This gives aβ +λ γ +λ ; z [ = z a + = z a aγ β γ +λ ; z z aγ β γ +λ [ = z a aγ β + γ +λ z aa+γ βγ β + + z γ +λγ +λ+! ] z z +O λ z...] z.5

4 4 as λ with z fixed. This is the beginning of a complete asymptotic expansion which converges if z/z <. It is an asymptotic expansion for λ large and all fixed zz. We can also use other transformation formulae: c ; z = z b c z ; = z c a b c ac b c ; z.ens6 z c These give with.4 the three relations aβ +λ γ +λ ; z = z a aγ β γ +λ ; z = z b γ +λ aβ +λ ; γ +λ z z z = z γ a β γ +λ aγ β ; z γ +λ.6 We see in.7 that the large parameter λ can be distributed over other parameter places. In the present case only the second form is suitable for giving an asymptotic expansion when using the Gauss series. The third form gives a useless Gauss series for large λ. In the final form the Gauss series converges at z = if Ra+β γ > but the series does not have an asymptotic property for large values of λ. The transformation formulae in.7 are an important tool for obtaining asymptotic expansions. We will use these and other formulae for investigating all cases of.. So far we mention the Gauss series for obtaining an asymptotic expansion. However to obtain an asymptotic expansion of. that holds uniformly if z is close to we need a different approach as will be explained in later sections. 3. Some history and recent activities [] Watson 98 studied the cases α+λβ +λ ; γ +λ z α+λβ λ ; z c γ +λ ; z. All results are of Poincaré type i.e. negative powers of the large parameter and hold for large domains of complex parameters and argument. Watson began with contour integrals and evaluated them by the method of steepest descent. [] Luke 969 Vol. I p. 35 summarizes these results and gives many other results also for higher p F q functions. In particular Luke investigates extended Jacobi polynomials which are of the type nn+λa...a p p+f q ; z. b...b q For p = q = and integer n this is the Jacobi polynomial. [3] Jones considers a uniform expansion of α+λβ λ ; c z

5 and gives a complete asymptotic expansion in term of I Bessel functions with error bounds. Applications are discussed for Legendre functions. [4] Olde Daalhuis and give new uniform expansions of aβ λ γ +λ ; z in terms of parabolic cylinder functions and of α+λβ +λ ; z c in terms of Airy functions. [5] Other recent contributions on uniform expansions for Gauss functions follow from Ursell 984 Legendre functions Boyd & Dunster 986 Legendre functions Temme 986 see later section Frenzen 99 Legendre functions Dunster 99 conical functions Dunster 999 Jacobi and Gegenbauer polynomials Wong and co-workers 99 Jacobi polynomials 996 Pollaczek 997 Jacobi function 998 and 999 Meixner Krawtchouk Meixner- Pollaczek. [6] On Bessel polynomials which are of F type see Wong & Zhang 997 and Dunster. On Charlier polynomials F type see Bo Rui & Wong 994. On Laguerre polynomials F type see Frenzen & Wong 988 and Temme e e e 3 e e e 3 e e e Table. All 7 cases of F a+e λb+e λ;c+e 3 λ;z e j = ±. In the table we write e j = ±. 4. From 7 cases to only a few Skipping the dummy using the symmetry between the a and b parameters: ba c ; z = F c ; z 4.

6 6 and by using the transformation formulae of.6 we reduce the 7 cases to the 8 cases: e e e This set is obtained by using 4. and.6 that is only 4 of Kummer s 4 solutions of the hypergeometric differential equation. A further reduction of the set of 8 cases can be obtained by other solutions of the differential equation. Using the connection formula c ; z = Γc Γa c+γb c+ z c z c a b a b ΓaΓbΓ c c + Γb c+γa c+ Γa+b c+γ c F a+b c+ ; z we see that the second case in fact reduces to +. Similarly we have ; z 4. c ; z = e aπi ΓcΓb c+ aa c+ Γa+b c+γc a z a + ΓcΓb c+ ΓaΓb a+ za c z c a b a+b c+ ; z ac a b a+ ; z 4.3 which reduces the third case + again to +. Finally we use c ; z = e cπiγc Γb c+γ a ΓbΓc aγ c +e c aπiγ aγb c+ Γ cγb a+ za c z c a b z c z c a b a b c ac a b a+ ; z which reduces the forth case + to the fifth case +. ; z 4.4

7 7 Consequently we have the five remaining cases: e e e These five cannot be reduced further by using relations between the functions a+e λb+e λ ; z e c+e 3 λ j = ±. However Olde Daalhuis using c ; z = e aπiγcγb a ΓbΓc a z a +e bπiγcγa b ΓaΓc b z b aa c+ a b+ ; z bb c+ b a+ ; z 4.5 and c ; z = e c bπiγcγb a ΓbΓc a zb c z c a b bc b a b+ ; z +e c aπiγcγa b ΓaΓc b za c z c a b ac a b a+ ; z 4.6 we see that the cases + + and + both can be handled by the case + +. Hence the smallest set is the following set of four: e e e 3 A + B + C + D + + if we accept the special case D: a+λb+λ ; z. c

8 8 4.. Quadratic transformations For special combinations of the parameters c other types of relations between F functions exist. For example we have the quadratic transformation b ; z For large a the left-hand side is of type = z a/ F a b+ ; = A whereas the right-hand side is of type C = Which special functions are involved? z. 4z 4 The next step is to indicate which special functions are associated with the above four cases. 5.. Legendre functions We have the Legendre function Case A : + P ν µ z = z + µ z µ νν + Γ µ µ ; z. 5. We see that case A occurs when µ. Several relations between the Legendre function and the Gauss function give also the case µ +. Case B : + Because = B see 3.5 we see that case B occurs when in P µ ν z both parameters ν and µ tend to +. Case C : + + This also occurs for P ν µ z with ν large and µ fixed. In the results for Legendre functions in the literature rather flexible conditions on the parameters µ and ν are allowed. For example in Dunster 99 it is assumed that the ratio ν/µ is bounded. Any hypergeometric function for which a quadratic transformation exists can be expressed in terms of Legendre functions. So many special cases and mixed versions of the cases ABCD are possible for Legendre functions. 5.. Hypergeometric polynomials Several orthogonal polynomials of hypergeometric type see the Askey scheme in Koekoek & Swartouw 998 have representations with Gauss hypergeometric functions. Next to the special cases of the Jacobi polynomials Gegenbauer Legendre we have the following cases.

9 9 Pollaczek P n cosθ; = e nθi n F +iφθ ; e iθ φθ = acosθ +b. sinθ For uniform expansions in terms of Airy functions see Bo Rui & Wong 996. The asymptotics is for n θ = t/ n with t bounded and bounded away from zero. Meixner-Pollaczek P λ n x;φ = λ n n! e nφi nλ+ix ; e iφ. λ For uniform expansions in terms of parabolic cylinder functions see Li & Wong. The asymptotics is for n x = αn with α in a compact interval containing the point α =. Meixner n x M n x;βc = F ;. β c In Jin & Wong 998 two uniform expansions in terms of parabolic cylinder functions are given. The asymptotics is for n x = αn with α > in compact intervals containing the point α =. Krawtchouk n x K n x;pn = F N ; n =...N. p A uniform expansion in terms of parabolic cylinder functions is given in Li & Wong. The asymptotics is for n x = λn with λ and ν = n/n in compact intervals of Jacobi polynomials The Jacobi polynomial has the representation: P αβ n x = n+α n nα+β +n+ ; α+ x 5. and we see that case C : + + applies if n. Uniform expansions use Bessel functions for describing the asymptotic behaviour at the points x = ± Hilb-type formula see Szegö 975 or for complete asymptotic expansions Wong & Zhang 996. A well-known limit is x lim α α n/ P n a+αb+α α = H nx n n!

10 where H n x is the Hermite polynomial a special case of the parabolic cylinder function. Approximations of H n x are available in terms of Airy functions. Several other limits are known see Koekoek & Swartouw998. In Chen & Ismail 99 and Dunster 999 asymptotic expansions are given for large positive parameters α and/or β. Of particular interest are asymptotic expansions for large n and large negative values of the parameters α β. See Bosbach & Gawronski 999 and Martínez-Finkelshtein et al.. With α = β = n we have in 5. the case + = A. These are non-classical values for the parameters. We have taking α = β = n P n αβ z = n! n/ nc n n x = n m= n!n! n m x n m. m!m!n m! This representation shows that all zeros are on the imaginary axis. Other representations can be used: β n+α +x P n αβ x = n α n+β x = n n α+n+ β n ; α+ x β +n+ α n ; β + +x 5.3 and we see that case D : + + applies if n and in the first case α is constant and β equals b 3n and similarly in the second case. See Figure for the zero distribution of the Jacobi polynomial P n αβ z n = 3 α = β = 3n+. 5 Figure. Zeros of P n αβ z n = 3α = β = 3n+.

11 The case B : + + applies in the representations P αβ n x = α+β +n x n α+β +n x+ = n n n α n n n α β ; x n β n n α β ; +x 5.4 when α = a nβ = b n first case or β = b nα = a n second case with again non-classical parameter values. Olde Daalhuis gives an expansion of the Gauss function for this case in terms of parabolic cylinder functions. See Figure for the zero distribution of the Jacobi polynomial P n αβ z n = 5 α = n+ β = n+. As n the zeros approach the curve defined by 3 z +z = Figure. Zeros of P n αβ z n = 5α = n+ β = n+. Uniform expansions for Jacobi polynomials are available in terms of J Bessel functions. Details on uniform expansions of Jacobi polynomials in terms of Bessel functions can be found in Wong & Zhang 996 Dunster 999 and in cited references in these papers Gegenbauer polynomials For Gegenbauer ultraspherical polynomials we have several representations. For example C λ n x = λ n n! and for λ = n n we have nn+λ λ+ ; x + = A. We also have C λ n x = nλ n n! nn+λ ; x

12 and for λ = n n we have + = C. Also in this case a sort of quadratic transformation is used Asymptotics of the A case We give a few steps in deriving the expansion for the case A = +. As always we can use the differential equation or an integral representation. We take c ; z Γc t b t c b = ΓbΓc b zt a dt 5.5 where Rc > Rb > ph z < π. A simple transformation gives the Laplace transform representation c+λ ; z Γc+λ = t b e λt ftdt ΓbΓc+λ b with λ and where e t b ft = e ct z +ze t a. t This standard form in asymptotics can be expanded by using Watson s lemma see Olver 974 or Wong 989. For fixed c and z: c+λ ; z Γc+λ Γc+λ b where f z = and other coefficients follow from ft = s= s= f s z b s λ b+s f s zt s z t < minπ t t = ln z. When z is large the singularity t is close to the origin and the expansion becomes useless. We can write c+λ ; z = Γc+λ t b e λt gt ΓbΓc+λ b t+ζ a dt where ζ = t = ln z gt = t+ζ a ft z with g regular at t = t. Expanding gt = g s zt s gives c+λ ; z Γc+λ Γc+λ b g s zb s ζ b a+s Ub+sb a++sζλ s=

13 in which U is the confluent hypergeometric function. This expansion holds if λ uniformly with respect to small values of ζ that is large values of z. This method can be used to obtain an expansion of the Legendre function P ν ±m z for large values of m; see Gil Segura & Temme where we used this type of expansion for computing this Legendre function for large m and z. Other sources with related expansions for this type of Legendre functions are Boyd & Dunster 986 Dunstger 99 Frenzen 99 Olver 974 and Ursell Another A case The case is very important for all kinds of orthogonal polynomials and can be reduced to the previous case A = +. It is of interest to give a direct approach of 3 nb ; z n 5.6 c whether or not n is an integer. First we recall the limit ac lim b F c ; z/b = F ; z which may be used as a definition for the F function the Kummer or confluent hypergeometric function. We are interested in the asymptotics behind this limit and we expect a role of the Kummer function when in 5.6 n becomes large and z is small. We again can use 5.5 and consider I n = t b t c b tz n dt where Rc > Rb > ph z < π. We transform Then we have tz = z u = e uln z. [ ] c ln z I n = z c b J n z J n = fuu b u c b e ωu du [ ] z u b [ ] z u c b fu = uln z uln z ω = n+ln z. The function fu is holomorphic in a neighborhood of [ ]. Singularities occur at u k = kπi ln z v m = + mπi km ZZ\{}. ln z So when z ranges through compacta of C\{} the singularities of f are bounded away from []. If Rω < the dominant point in the integral is u = ; if Rω > then u = is the

14 4 dominant point. To obtain an asymptotic expansion for large n that combines both cases and which will give a uniform expansion for all z z > δ > f should be expanded at both end-points and. More details can be found in Temme 986 where the results have been applied to a class of polynomials biorthogonal on the unit circle Asymptotics of the B case We consider a simple case: n n+ ; z with z near the point. We write and we have = n+ [ t+zt] n dt 5.7 [ t+zt] n = e nφt φt = ln[ t+zt] φ t = tz +z t+zt. We see that the integrand has a peak value at t = z z. So if z = the peak is at t = if z > then t and if z < then t <. See Figure 3. z = /4 z = z = 4 Figure 3. t+zt for a few values of z. The same situation occurs for the integral e n w αw dw 5.8 which has a peak value at w = α. This integral is an error function a special case of the parabolic cylinder functions. In a uniform expansion of 5.7 that holds for large n and z in a neighborhood of we need an error function. This explains why in Case B in a more general form than 5.7 parabolic cylinder functions may occur for special values of z.

15 5 We transform the integral in 5.7 by writing t+zt = e w +αw with the conditions t = w = t = w = t = t w = α. The quantity α follows from satisfying the matching of t with α: t +zt = e w +αw giving where We obtain α = ln [ n n+ ; z = n+ ] z signα = signz. z + fw = dt dw = w α t+zt. t t z e n w αw fwdw 5.9 giving A first approximation follows by replacing fw by n n+ ; z nfα fα = +z z e n w αw dw = πn +z 4z e nα erfc α n/. as n uniformly with respect to z in a neighborhood of z =. This expansion is in agreement with the general case F λ;c+λ; z considered in Olde Daalhuis. The relation with the Jacobi polynomials is: n n+ ; z See Figure 4 for the distribution of the zeros of = n+! n 3 +z n P n n+ n n P n αβ z n = 3 α = n+ β = n. z. +z

16 6.5.5 Figure 4. Zeros of P n αβ z n = 3α = β = n+. As n the zeros approach the curve z =. 6. Asymptotics of some 3 F polynomials In the previous examples we used integral representations of the Gauss function. The differential equation can also be used for obtaining asymptotic expansions. As a rule methods based on differential equations provide more information on the coefficients and remainders in the expansions than methods based on integral representations. For the generalized hypergeometric functions asymptotic expansion are known when the argument z is large and the other parameters are fixed; see Braaksma 963 and Luke 969. For the p F q functions a differential equation is available of order maxpq+. However for higher order equations no methods are available for deriving uniform expansions for large values of parameters. The p F q functions can be written as a Mellin-Barnes contour integral but again no methods are available for deriving uniform expansions for large values of parameters from these integrals. The terminating p F q functions one of the parameters a k is equal to a non-positive integer and of unit argument z = are of great interest in special function theory and in applications. In many cases recursion relations are available. For the examples to be considered here we use ad hoc methods. It is of interest to investigate recursion methods for the same and other problems. 6.. A first example As a first example Larcombe & Koornwinder consider the asymptotics of n fn = 3 F n; n

17 It was conjectured by Larcombe that lim n fn = and Koornwinder gave a proof. By using k n k = k n! πγn+ k =...n we can write fn as an integral fn = n! π n t k t n k dt t t n n t n; dt t and the Gauss function is a Legendre polynomial. We obtain fn = n n!n! π π n sin n θ P n dθ. 6. n sinθ There is no straightforward way to obtain asymptotics out of this. A few manipulations with 6. give the result n!n!n! fn = n n c n n where c n is the coefficient in the Bessel function expansion [ ] e w/ I w/ = and now the asymptotics easily follows. We have c n = [ e w/ I w/ ] dw = πi C w n+ πi C n= c n w n e w [ ] w fwdw fw = e w/ I n+ w/ and the contour C is a circle around the origin. The main contribution comes from the saddle e point of w w that is from w = w n+ = n/. Because fw we have πn c n n πn n! n. This gives finally fn. A complete asymptotic expansion follows by using more terms in the expansion of fw at w = w. 6.. Generalizing Kummer s identity Kummer s identity reads: where +a b.... +a b ; = Γ+a bγ + a Γ+aΓ + a b. 7

18 8 A generalization Vidunas is written in the form a+nb +a b ; = Pn Γa bγ a+ ΓaΓ a+ +QnΓa bγ a b ΓaΓ a b. Vidunas showed that for n =...a n and a b...: Pn = n+3 F n n a b a ; Qn = n+ n+ 3 F n+ n a+ b 3 a+ ; = 3 F n n b n a ; = 3 F n+ nb n a + ;. 6. We are interested in the asymptotic behaviour of Pn and Qn for large values of n. Using if Rd > Rc > c k d k = Γd Γd cγc t d c t c+k dt and a few manipulations of the Gauss functions we obtain the integral representations for Ra > Rb > Pn = n Γ a ΓbΓ a b Qn = n Γ a+ ΓbΓ a+ b sinh a b t cosh a n t coshn+tdt sinh a b t cosh a n t sinhn+tdt 6.3 We can use standard methods for obtaining the asymptotic behaviour of these integrals. First we write the hyperbolic functions coshn + t and sinhn + t as exponential functions. We obtain for n Pn b Γ a Γb+ ΓbΓ a b c k k n b+ k k= Qn b Γ a+ Γb+ ΓbΓ a+ b c k k n b+ k k= where c = c = c = 5b 4a+3. These expansions are not valid for b = a b = a + respectively. In these cases we have Pn = n+ n+ Qn = n+ as follows from 6.. Similar for the cases b = a+m b = a+ +m m =...

19 9 For negative values of n different integral representations should be used. We have P n = n a n 3F n n+ a b b n a n ; Q n = n n a n 3F n+ n+ a+ b 3 b n a+ n ;. with the integral representations P n = n+ Γ b π Γ a bγ a n+ Γ b Q n = Γ a+ bγ a cos n+ a t sin c t cosntdt π cos n+ a t sin c t sinntdt. By integrating on contours i π π π+i similar integrals arise as in 6.3 and again standard methods can be used for obtaining the asymptotic behaviour. 7. Concluding remarks [.] The asymptotic analysis of all 6 cases in can be reduced to the 4 cases a+e λb+e λ ; z e c+e 3 λ j = ± e e e 3 A + B + C + D + + [.] These cases are of interest for orthogonal polynomials and special functions. [3.] New recent results on uniform asymptotic expansions have been published or announced for special cases. [4.] The distribution of the zeros of Jacobi polynomials for non-classical values of the parameters α and β shows interesting features. New research is needed for describing the asymptotics of these distributions. [5.] The asymptotics of 3 F terminating functions with large parameters is quite difficult. Standard methods based on integrals and differential equations are not available. Recursion relations may be explored further.

20 8. Bibliography [ ] M. Abramowitz and I.A. Stegun 964 Handbook of mathematical functions with formulas graphs and mathematical tables Nat. Bur. Standards Appl. Series 55 U.S. Government Printing Office Washington D.C. paperback edition published by Dover New York. [ ] G.E. Andrews R. Askey R. Roy Ranjan 999. Special functions Cambridge University Press Cambridge. [ 3] Christof Bosbach and Wolfgang Gawronski 999. Strong asymptotics for Jacobi polynomials with varying weights. Methods Appl. Anal [ 4] W.G.C. Boyd and T.M. Dunster 986. Uniform asymptotic solutions of a class of second-order linear differential equations having a turning point and a regular singularity with an application to Legendre functions SIAM J. Math. Anal [ 5] B.L. Braaksma 963. Asymptotic expansions and analytic continuation for a class of Barnes-integrals Noordhoff Groningen. [ 6] L.-C. Chen and M.E.H. Ismail 99. On asymptotics of Jacobi polynomials SIAM J. Math. Anal [ 7] T.M. Dunster 99. Conical functions with one or both parameters large Proc. Royal Soc. Edinburgh 9A [ 8] T.M. Dunster 999. Asymptotic approximations for the Jacobi and ultraspherical polynomials and related functions Methods and Applications of Analysis [ 9] T.M. Dunster. Uniform asymptotic expansions for the reverse generalized Bessel polynomials and related functions SIAM J. Math. Anal [ ] C.L. Frenzen and R. Wong 988. Uniform asymptotic expansions of Laguerre polynomials SIAM J. Math. Anal [ ] C.L. Frenzen 99. Error bounds for a uniform asymptotic expansion of the Legendre function Q m n cosh z SIAM J. Math. Anal [ ] Amparo Gil Javier Segura and Nico M. Temme. Computing toroidal functions for wide ranges of the parameters J. Comput. Phys [ 3] X.-S. Jin and R. Wong 998. Uniform asymptotic expansions for Meixner polynomials Constructive Approximation [ 4] X.-S. Jin and R. Wong. Asymptotic formulas for the zeros of the Meixner-Pollaczek polynomials and their zeros Constr. Approx [ 5] D.S. Jones. Asymptotics of the hypergeometric function Math. Methods Appl. Sci [ 6] R. Koekoek & R.F. Swartouw 998. The Askey-scheme of hypergeometric orthogonal polynomials and its q analogue. Technical University Delft. Report [ 7] X.-C. Li and R. Wong. A uniform asymptotic expansion for Krawtchouk polynomials J. Approximation Theory [ 8] X. Li and R. Wong. On the asymptotics of the Meixner polynomials J. Approximation Theory [ 9] Y.L. Luke 969. The special functions and their approximations Vol. I II Academic Press New York. [ ] A. Martínez-Finkelshtein P. Martínez-González and R. Orive. Zeros of Jacobi polynomials with varying non-classical parameters In Special functions Hong Kong 999 Charles Dunkl Mourad Ismail and Roderick Wong eds. World Sci. Publishing River Edge NJ [ ] A. Martínez-Finkelshtein and N.M. Temme 999. Private communications.

21 [ 3] A. Olde Daalhuis. Private communications. [ 4] A. Olde Daalhuis. Uniform asymptotic expansions for hypergeometric functions with large parameters I II. Submitted. [ 5] F.W.J. Olver 974 & 997. Asymptotics and Special Functions. Academic Press New York. Reprinted in 997 by A.K. Peters [ 6] Bo Rui and R. Wong 994. Uniform asymptotic expansion of Charlier polynomials Methods and Applications of Analysis [ 7] Bo Rui and R. Wong 996. Asymptotic behavior of the Pollaczek polynomials and their zeros Stud. Appl. Math [ 8] G. Szegö 975. Orthogonal polynomials 4th edition Amer. Math. Soc. Colloq. Publ. 3 Providence R.I. [ 9] N.M. Temme 986. Uniform asymptotic expansion for a class of polynomials biorthogonal on the unit circle Constr. Approx [ 3] N.M. Temme 99. Asymptotic estimates for Laguerre polynomials ZAMP [ 3] N.M. Temme 996. Special functions: An introduction to the classical functions of mathematical physics Wiley New York. [ 3] F. Ursell 984. Integrals with a large parameter: Legendre functions of large degree and fixed order Math. Proc. Camb. Philos. Soc [ 33] R. Vidunas. A generalization of Kummer s identity. Accepted for publication in Rocky Mountain J. Math.. [ 34] G.N. Watson 98. Asymptotic expansions of hypergeometric functions Trans. Cambridge Philos. Soc [ 35] R. Wong 989. Asymptotic approximations of integrals Academic Press New York. [ 36] R. Wong and Q.-Q. Wang 99. On the asymptotics of the Jacobi function and its zeros SIAM J. Math. Anal [ 37] R. Wong and J.-M. Zhang 996. A uniform asymptotic expansion for the Jacobi polynomials with explicit remainder Appl. Anal [ 38] R. Wong and J.-M. Zhang 997. The asymptotics of a second solution to the Jacobi differential equation Integral Transforms Spec. Funct [ 39] R. Wong and J.-M. Zhang 997. Asymptotic expansions of the generalized Bessel polynomials J. Comput. Appl. Math