TAYLOR SERIES FOR THE ASKEY-WILSON OPERATOR AND CLASSICAL SUMMATION FORMULAS arxiv:math/ v2 [math.ca] 26 Aug 2005

Transcription

1 TAYLOR SERIES FOR THE ASKEY-WILSON OPERATOR AND CLASSICAL SUMMATION FORMULAS arxiv:math/0348v [math.ca] 6 Aug 005 BERNARDO LÓPEZ, JOSÉ MANUEL MARCO AND JAVIER PARCET Abstract. An analogue of Taylor s formula, which arises by substituting the classical derivative by a divided difference operator of Askey-Wilson type, is developed here. We study the convergence of the associated Taylor series. Our results complement a recent work by Ismail and Stanton. Quite surprisingly, in some cases the Taylor polynomials converge to a function which differs from the original one. We provide explicit expressions for the integral remainder. As application, we obtain some summation formulas for basic hypergeometric series. As far as we know, one of them is new. We conclude by studying the different forms of the binomial theorem in this context.. Introduction and definitions The problem of expanding a function with respect to a given polynomial basis has many implications in analysis. The simplest example of this kind is the Taylor s expansion theorem. In this paper, we replace the classical derivative by a difference operator of Askey-Wilson type. Our results complement the paper [3] of Ismail and Stanton and are a natural continuation of the point of view presented in [5], where a new approach to the theory of classical hypergeometric polynomials is given. In contrast with [3], our aim is to find sufficient conditions for the Taylor series to converge, but not necessarily to the original function. In this more general setting, we may consider non-necessarily entire functions and we give an explicit expression for the limit of the remainders in terms of a contour integral. Using this and a new estimate for the q-shifted factorials, which might be of independent interest, we obtain a summation formula which is new as far as we know. As we explain below, it can be regarded as a non-symmetrized version of the non-terminating q-saalschütz sum. As applications, we also provide a new proof of the q-gauss summation formula and a list of binomial type summation formulas in the same line than Ismail s paper []. Now we give some definitions which will be used in what follows. The notions we are presenting were already introduced in [5] with the aim of studying some aspects of the theory of hypergeometric polynomials. The relevance of this approach is justified in [5], where a more detailed exposition is given. Definition.. Let P be the set of quadratic symmetric polynomials P(x, y), normalized so that P(x, y) x y axy b(x y) c with a, b, c complex numbers. A sequence of complex numbers (x t ), with the index t running over the set Z {k/ : k Z}, will be called a P-sequence if P(x t, y) ( t )( t ) for all t Z. 000 Mathematics Subject Classification. Primary: 33D5. Key words and phrases: q-taylor series, Askey-Wilson operator, Basic hypergeometric function.

2 B. LÓPEZ, J.M. MARCO AND J. PARCET Recalling the natural action of the affine group of the complex plane Aff(C) on the set P, we can consider the orbits of this action. As it becomes clear in [5], these orbits provide a very natural classification in the theory. In particular, each symmetric polynomial in P can be rewritten, under affine transformations, in a canonical form. In the following table, we summarize the canonical forms and the corresponding P-sequences which we shall use in our study. We shall also use the parameter λ 0 defined by the relation a (λ λ ). Here we recall that the basis q will be given by λ. Canonical form P-sequence T x y axy a (λt u λ t u ) G x y axy λ t u Q (x y) (x y) 6 t tu u A (x y) 4 t u C (x y) u Table I. Canonical P-sequences. In the first two rows we assume a ±. The capital letters in the left column are acronyms of the names we adopted in [5] for the canonical forms: Trigonometric, Geometric, Quadratic, Arithmetic and Continuous. We shall say that x 0 is the base point of the P-sequence (x t ). Obviously, P(x, y) y A(x)y B(x), where A(x) ax b and B(x) x bx c. This allows us to consider the discriminant of P as a function of the variable x which, up to a constant factor, is given by δ(x) (a )x b(a )x b c. Remark.. P-sequences arise from the recurrence x t± A(x t) ± δ(x t ). In particular, given a complex number ξ, there are at most two P-sequences with base point ξ, one for each choice of the sign for the square root δ(ξ). Definition.3. For any complex number x, we consider a P-sequence (x t ) with base point x. Then we define Φ 0 (x, y) and the polynomials Φ k (x, y), for any positive integer k, as follows k Φ k (x, y) ( j k ). Remark.4. Recall that Φ k (x, ) does not depend on the chosen P-sequence (x t ).. Taylor series Given an open subset Ω of the complex plane, we denote the space of analytic functions in Ω by H(Ω). Also, by γ 0 (mod Ω) we mean that γ is a cycle in Ω homologous to zero with respect to Ω. Finally, given z Ω, Ind(γ, z) denotes the index of z with respect to γ. Lemma.. Given an open subset Ω of the complex plane, let us denote by Ω the diagonal of Ω Ω. Then, for any f H(Ω), the function f d (u, v) f(u) f(v) u v for (u, v) (Ω Ω) \ Ω

3 TAYLOR SERIES FOR THE ASKEY-WILSON OPERATOR 3 can be continuously extended to an analytic function f d : Ω Ω C. Moreover, if γ 0 (mod Ω) and Ind(γ, u) Ind(γ, v), then we have f d (u, v) πi (y u)(y v) dy. Proof. It is a simple consequence of Cauchy s integral formula. γ Let us consider the set Ω {x C : A(x)± δ(x) Ω}. By Lemma., given a quadratic symmetric polynomial P P, we define the Askey-Wilson operator D : H(Ω) H(Ω ) as follows () Df(x) f d ( A(x) δ(x), A(x) δ(x) ) πi γ P(x, y) dy, with γ 0 (mod Ω) and Ind(γ, A(x) ± δ(x)). On the other hand, if we define recursively the sets Ω k by Ω k {x C : A(x) ± δ(x) Ω k } with Ω 0 Ω, we may consider the iterated operators D k : H(Ω) H(Ω k ). Lemma.. Let Ω be an open set of the complex plane and γ 0 (mod Ω). Then, given f H(Ω) and x Ω k, we have the following identity (with the obvious limits for λ ±) D k f(x) ( k λ k j λ j k ) λ λ πi γ dy for k 0 Φ k (x, y) if x is the base point of a P-sequence (x t ) with Ind(γ, x j k ) for j 0,,..., k. In particular, we shall consider the operator k : H(Ω) H(Ω k ) defined by k f(x) πi γ Φ k (x, y) dy. Proof. The cases k 0 and k follow from Cauchy s integral formula and () respectively. Therefore, the general case follows by induction from the following relation [ D ](x) Φ k(x /, y) Φ k (x /, y) Φ k (, y) x / x / which is not difficult to check with the aid of the identity x k/ x k/ x / x / λk λ k λ λ. λk λ k λ λ Φ k(x, y), This identity was proved in [5] for any P-sequence. This completes the proof. The following result provides the sequence of Taylor polynomials associated to a given function f H(Ω) and the corresponding remainder term, with respect to the Askey-Wilson operator. Theorem.3. Let Ω be an open subset of the complex plane and γ 0 (mod Ω). Then, given an analytic function f H(Ω) and x Ω with Ind(γ, x), we can recover f(x) as f(x) n k k f(z k/ ) (x z j ) R n f(x)

4 4 B. LÓPEZ, J.M. MARCO AND J. PARCET where R n f(x) is given by R n f(x) πi γ n x z j y z j dy, with (z t ) a P-sequence such that z j Ω and Ind(γ, z j ) for j 0,,...,n. Proof. The relation below follows by induction on n Φ n (z n/, x) n ()Φ n (z n/, y) Φ k (z (k )/, x) Φ k (z k/, y). Multiplying by (πi), we are done by Lemma. and Cauchy s formula. Remark.4. When dealing with the continuous form C, it turns out that the Askey-Wilson operator can be regarded as the classical derivative, see [5] for the details. In particular, Theorem.3 reduces to the classical Taylor series where f(x) n f (k) (z) (x z) k f n (x)(x z) n, k! f n (x) πi γ (y z) n () dy. Remark.5. If the points z j are pairwise distinct for j 0,,...,n, we can use Lemma. to express k f(z k/ ) as a sum of residues and Theorem.3 gives f(x) R n f(x) n ( k f(z j ) i j (z j z i ) ) k (x z j ), which is Newton s divided difference formula for the interpolation polynomial. In particular, Theorem.3 holds for any collection z 0, z,..., z n of pairwise distinct points. However, as we shall see below, the relevance of Theorem.3 lies in the established connection with the Askey-Wilson operator. Now we study the convergence of the Taylor series for P-sequences (z t ) with z 0, z, z,... bounded. This covers the geometric canonical form for q < and the continuous form, whose associated P-sequences are constant. Given r > 0, we shall denote by D r {z C : z < r} the open disk of radius r. As it was pointed out by the referee, the result below is closely related to Wallisser s paper [6]. Theorem.6. Let (z t ) be a P-sequence satisfying that the subsequence z 0, z, z,.. is bounded. Then, given z limsup k 0 z k, δ > 0 and f H(D rδ ) with r > z, the Taylor polynomials S n f(x) n k k f(z k/ ) (x z j ) converge uniformly to f(x) as n on the disk D x for any x < r z. Proof. By Theorem.3, it suffices to check R n f 0 uniformly on the disk D x. Let us take γ to be a circumference centered at 0 with radius r > x z. It turns out that, if M r (f) denotes the supremum of when y runs over γ, we have R n f(x) r M r(f) r x ( x z r z ) n 0 as n.

5 TAYLOR SERIES FOR THE ASKEY-WILSON OPERATOR 5 Since the given bound holds for every x D x, the proof is complete. The following is a convergence theorem with integral remainder. This case is applicable to the trigonometric and quadratic canonical forms as well as to the geometric canonical form with q >. Let (z t ) be a P-sequence satisfying k 0 / z k <. Then, the product ( ) H(x) x/zj is an entire function of x with zeros at x z j for j 0,,,... Theorem.7. Let (z t ) be a P-sequence satisfying that the subsequence z 0, z, z,... lies in R {x R : x < 0} and k 0 / z k <. Assume that f is analytic in an open neighborhood of the left half-plane Ω { z C : Re z 0 }, and is controlled by () f(z) C ( z ) M for all z Ω and some absolute constants C and M. In that case the pointwise limit exists when Re x < 0 and we have with R f(x) i πi i S f(x) lim n S nf(x) f(x) S f(x) R f(x) H(x) H(y) dy i πi i Proof. According to Theorem.3 we have R n f(x) ( ir n x z j dy πi ir y z j γ r n x z j y z j dy. x z ) j dy, y z j whenever Re x < 0 and where r and γ r : [0, π] C are determined by { } max x, z 0, z,..., z n < r and γ r (t) ire it (0 t π). On the other hand, it easily follows from () that for any y γ r [0, π] we have n x z j C ( r ) M n, y z j for some constant C depending on x, z 0, z,...,z n but not on y. Therefore, since we may take the parameter r arbitrary large, we deduce that for any integer n M we have lim r γ r n x z j y z j dy 0.

6 6 B. LÓPEZ, J.M. MARCO AND J. PARCET In other words, for n M we have R n f(x) n ( ) i x/zj πi i Then, estimating as above we find for n M n ( ) y/zj C ( y ) M n n ( ) y/zj dy. ( y z j ) / C ( y ). Therefore, the assertion follows from the dominated convergence theorem. Remark.8. Similar arguments show that the convergence in Theorem.7 is uniform over the compact sets of {x C : Re x < 0}. On the other hand, it is clear that Theorem.7 can be restated by taking z 0, z, z,... in R, Ω { } z C : Rez 0 and R f(x) πi i i x z j y z j dy. This shows that the quadratic P-sequences are contemplated by Theorem.7. We have emphasized this point following a suggestion from the referee. However, when dealing with the trigonometric form, we might consider a growth restriction weaker than (). We omit the proof for lack of space and refer the reader to [3]. Remark.9. The main result in Ismail/Stanton s paper [3, Theorem 3.] considers Taylor series of entire functions with respect to the trigonometric canonical form. Our main contribution in Theorem.7 is that we do not require f to be entire. Note that the regularity and growth restrictions in [3] allowed the authors to obtain the Taylor series of products of q-shifted factorials. Theorem.7 and Lemma 3. below allow us to consider in the next section quotients of q-shifted factorials. 3. Summation formulas Now we compute the Taylor coefficients of a couple of functions made up of q-shifted factorials. Then we analyze the convergence of the corresponding Taylor series. We obtain a new proof of the q-gauss summation formula and, as far as we know, a new summation formula. 3.. Geometric case. We begin by applying Theorem.6 to a particular case. This will provide a new proof of the q-gauss summation formula. We shall work with the Askey-Wilson operator which arises from the geometric canonical form Df(z) f(q/ z) f(q / z). q / z q / z Given α, β C and q <, we consider the function f(x) (αx; q) q j αx (βx; q) q j βx. By induction, it is not difficult to check that k f(z) ( k β q j α ) (αzq k/ ; q) q j. (βzq k/ ; q)

7 TAYLOR SERIES FOR THE ASKEY-WILSON OPERATOR 7 Therefore, taking ξ C and z t q t ξ in Theorem.6, we obtain (αx; q) (βξ; q) (ξ/x; q) ( k k ) (β q j α) x k. (βx; q) (αξ; q) (q; q) k (αξ; q) k In particular, with the usual notation for basic hypergeometric series [ ] ξ/x, α/β φ ; q, βx (αx, βξ; q). αξ (βx, αξ; q) In other words, we have obtained the q-gauss classical summation formula. 3.. An estimate for q-shifted factorials. In this paragraph we provide an estimate for q-shifted factorials that will be needed below. Although we just need a weaker estimate, we have included the sharpest result we know since it might be of independent interest. Lemma 3.. The following assertions hold for 0 < q < : (a) There exists some C q > 0 such that every x C satisfies ( q j x ) Cq x log x log q. (b) Let A be any closed set in C with qa A and q j / A for j 0,,,... Then, given any δ > 0 there exists a constant C q (A, δ) > 0 such that for all x A \ D δ, we have Proof. Let us write h(x) ( q j x ) Cq (A, δ) x log x log q. ( q j x ) and h n (x) n ( q j x ). We begin by proving the upper estimate. Since the function x log x log q explotes at x 0, we may assume that x >. Then we fix the only integer m satisfying the estimates q m x < q m x and observe that (3) h(x) h(q m x) m ( q j x ) x m q (m ) h(a)hm (b), where a q m x and b q m x. Our choice of m implies m Moreover, a simple computation gives log x ρ for some 0 ρ <. log q (4) x m q ) (m ρ(ρ ) q x log x log q. Therefore, since a and q b <, we have h(a)h n (b) h( ) and we deduce the assertion with constant C q q /8 h( ). Now we prove the lower estimate. Let us consider the set { } B z C : z A, q z.

8 8 B. LÓPEZ, J.M. MARCO AND J. PARCET The function h does not vanish on A B. Therefore, since h n h uniformly over compact sets, we deduce that min { h n (z), h(z) } C q (A) for some C q (A), for all n and all z A B with z. Thus, if x > we deduce from (3) and (4) (recall that a A and b B) that ( q j x ) Cq (A) x log x log q for all x A \ D. The assertion of x A with δ x follows by replacing C q (A) by C q (A, δ). Remark 3.. Let us observe how to construct sets of type A. Since q j / A for j 0,,,... and A is closed, there must exists an open neighborhood U of such that A U. Now, recalling that qa A if and only if q A c A c, we must have that the neighborhood of q j defined by q j U is contained in A c. Thus, the biggest A s are sets generated by a neighborhood U of in the following way ( A q j U ) c Trigonometric case. In this paragraph we apply Theorem.7 to ( u u ) f αβ (x) f αβ (αu; q) (α/u; q) with x (βu; q) (β/u; q) (u u ) and (α, β) R R. As a function of x, f α,β is well-defined since (γu, γ/u; q) is symmetric under the mapping u u. Moreover, it is analytic in a neighborhood of the left half-plane (as a function of x) since (γu, γ/u; q) is analytic in C \ {0} for any γ C and has zeros in R when γ R. We shall use the Askey-Wilson operator which arises from the trigonometric canonical form ( u u ) Df (λ λ )(u u ) [ f ( λu λ u ) f ( λ u λu ) ] with q λ and 0 < λ <. Grouping the common factors, it can be checked that ( u u ) Df αβ α β ( u u ) q f λα,λ β. Moreover, by induction on k we obtain ( u u ) k f αβ k (q; q) k k ( u u (β q j ) α) f λk α,λ k β. In particular, taking z t ( q t ξ q t ξ ) with ξ < 0, we can write k f αβ (z k/ ) (β)k (α/β; q) ( k ξ ξ ) (q, αξ, β/q k f αβ. ξ; q) k According to Theorem.3 and ζ ζ (5) η ( η (ζη) / (ζη) /)( (ζ/η) / (ζ/η) /), the Taylor polynomials of f α,β are ( u u ) ( ξ ξ ) n (α/β, ξu, ξ/u; q) k S n f αβ f αβ q k. (q, αξ, qξ/β; q) k

9 TAYLOR SERIES FOR THE ASKEY-WILSON OPERATOR 9 In other words, in terms of Theorem.7 we have ( u u ) (6) S f αβ (αξ, α/ξ; q) [ α/β, ξu, ξ/u 3φ (βξ, β/ξ; q) αξ, qξ/β ] ; q, q. On the other hand, the subsequence z 0, z, z,... lies in R with k 0 / z k <. Thus, in order to apply Theorem.7, it suffices to check whether the bound () holds for some M > 0. To that aim we shall assume that u and { x max α, α, β, β }. This implies 4 u x u so that α/u, β/ u / and the functions (α/u; q), (β/u; q) are bounded from above and below. Moreover, according to Lemma 3., we obtain the following estimate (αu; q) C q (α) u log u log u log α log q α log q u log q C q (α) u log u log α log q log q, with C q (α) C q α log α log q. Similarly, when Re u 0 (equivalently Re x 0) (βu; q) Cq (β) u log u In summary, we have seen that log u log β log q β log q u fα,β (x) Cq (α, β) u M 4C q (α, β) ( x ) M log q C q (β) u log u log β log q log q. with M log α/β log q, whenever Re x 0 and x max( α, α, β, β ). However, the last restriction on x can be dropped by continuity and we are in the hypotheses of Theorem.7. This gives the identity f α,β (x) S f α,β (x) R f α,β (x) for Re x < 0. Letting f γ (x) (γu, γ/u; q) (so that f α,β f α /f β ), we have H(x) ( u ) u q j ξ q j ξ Therefore, we have (7) R f α,β (x) πi i i ( q j ξu)( q j ξ/u) q j ξ ( ξ ; q ) f ξ(x). f α,β (y) H(x) H(y) dy i f α (y)f ξ (x) πi i f β (y)f ξ (y) Taking y ( v v ), it follows from (6) and (7) that (8) (αu, α/u; q) (αξ, α/ξ; q) [ α/β, ξu, ξ/u 3φ (βu, β/u; q) (βξ, β/ξ; q) αξ, qξ/β πi i i ] ; q, q dy. (αv, α/v, ξu, ξ/u; q) ( v )u (βv, β/v, ξv, ξ/v; q) (u v)(uv )v dv, for (α, β) R R, ξ < 0 and Re u < 0. This summation formula is new as far as we know. It can be regarded as a non-symmetrized version of the non-terminating q-saalschütz sum. To see this, we observe that an obvious reformulation of Theorem.7 gives (8) when Re x 0 with R f α,β (x) replaced by R (x; α, β, ξ) ( x δi πi i γ x i xδi ) fα (y)f ξ (x) f β (y)f ξ (y) dy,

10 0 B. LÓPEZ, J.M. MARCO AND J. PARCET with γ x (t) x iδeit for 0 t π. Transposing (β, ξ) and changing signs ( x δi i ) fα ( y)f β (x) R ( x; α, ξ, β) πi f ξ ( y)f β ( y) πi i ( xδi i γx γ x xδi i x δi ) fα (y)f β (x) f β (y)f ξ (y) dy, where γ x (t) x iδe it for 0 t π. Calculating a residue we conclude (9) f α (x) f β (x)f γ (x) R (x; α, β, ξ) f ξ (x) R ( x; α, ξ, β). f β (x) dy y x The non-terminating q-saalschütz sum follows easily from (8) and (9) in the form (αu, α/u; q) (αξ, α/ξ; q) [ ] α/β, ξu, ξ/u 3φ ; q, q (βu, β/u, ξu, ξ/u; q) (βξ, β/ξ, ξu, ξ/u; q) αξ, qξ/β (αβ, α/β; q) [ α/ξ, βu, β/u ] 3φ ; q, q. (βu, β/u, ξβ, ξ/β; q) αβ, qβ/ξ Note that our restrictions on α, β and ξ disappear after analytic extension. Remark 3.3. Going back to the example considered in Paragraph 3., we may consider the P-sequence z t q t ξ and apply Theorem.7 instead of Theorem.6. In that case, similar arguments give rise to a non-symmetrized version of the non-terminating q-vandermonde sum (αx; q) (βx; q) (αξ; q) (βξ; q) φ [ α/β, x/ξ q/βξ, ] ; q, q i (αy, x/ξ; q) πi i (βy, y/ξ; q) 4. Binomial theorem dy. As announced in the Introduction, we finish with a list of binomial type formulas adapted to our classification into canonical forms. We must recall that the same techniques we are using here were already employed by Ismail in [] to obtain the binomial type sum for the trigonometric canonical form, see below. Let f P n [x] and (z t ) a P-sequence with z 0, z,...z n pairwise distinct. Then, it follows from Remark.5 that R m f 0 for any m n. Moreover, by a simple continuity argument, we may drop the assumption that z 0, z,... z n are pairwise distinct. We need to use the q-binomial coefficients [ r k ] q n (q; q) r where (x; q) n ( q k x). (q; q) k (q; q) r k Theorem 4.. Given two P-sequences (y t ) and (z t ), we have n (x z k ) n [ n ] k q q k(n k)/ k (x y j ) Proof. We claim that the following identity holds k Φ n (x, ) ( k n k (y k/ z j k λ n j λ j n λ k j λ j k ) Φ n k (x, ), 0 k n. ).

11 TAYLOR SERIES FOR THE ASKEY-WILSON OPERATOR It can be checked by induction on k and applying the last identity in the proof of Lemma.. In particular, since Φ n (z (n )/, ) P n [x], Theorem.3 gives Φ n (z (n )/, x) n ( k λ n j λ j n ) k λ k j λ j k Φ n k (z (n )/, y k/ ) By the definition of Φ k and q-binomial coefficient, the proof is completed. (x y j ). Theorem 4. is symmetric under the transformation q q. Moreover, when q we must take the obvious limits. Clearly, for the continuous canonical form, Theorem 4. is nothing but Newton s binomial theorem. Other explicit formulas of binomial type arise from the remaining canonical forms. Moreover, them can be rewritten in terms of the basic hypergeometric function. In the following table we summarize the binomial type expressions and we mention the corresponding hypergeometric formula. Binomial form Sum C n ( n (x z) n (x y) k) (y z) n k Newton A Q G T ( x z ) n ( x y )( y z ) Gauss n k n k ( u w ) n ( u v ) ( v w ) (u w) n (u v) k (v w k) n k Pfaff n k n k n n [ n ] k n k (x q k z) (x q j y) (y q j z) q-gauss k q (wu, w/u; q) n n (vu, v/u; q) k (wvq k, w/v; q) ( n k w ) k q-saalschütz (q; q) n (q; q) k (q; q) n k v Table II. Binomial type summation formulas. The reader is referred to Gasper and Rahman s book [] to see these summation formulas. The symbol means that, if we rewrite the binomial type formula in terms of hypergeometric summation formulas, what we obtain is a terminating form of the corresponding non-terminating summation formula. Although we leave the details for the interested reader, we should point out the changes of variables employed in the quadratic and trigonometric canonical forms. (Q) x u y k (v k) z k (w k) (T) x (u u ) y k (qk v q k v ) z k (qk w q k w ). In the trigonometric case, we also need to use the factorization formula (5). Acknowledgments The authors wish to thank the referee for several comments that gave rise to a significant improvement of this paper. This research was supported in part by the Project MTM , Spain. References. G. Gasper and M. Rahman, Basic hypergeometric series. With a foreword by Richard Askey., Encyclopedia of Mathematics and its Applications 35, Cambridge Univ. Press, M.E.H. Ismail, The Askey-Wilson operator and summation theorems, M. Ismail, M.Z. Nashed, A. Zayed, A. Ghaleb (Eds.), Mathematical Analysis, Wavelets and Signal Processing, Contemp. Math. 90 (995), 7-78.

12 B. LÓPEZ, J.M. MARCO AND J. PARCET 3. M.E.H. Ismail and D. Stanton, q-taylor theorems, polynomial expansions, and interpolation of entire functions, J. Approx. Th. 3 (003), M.E.H. Ismail and D. Stanton, Applications of q-taylor theorems, J. Comp. Appl. Math. 53 (003), J.M. Marco and J. Parcet, A new approach to the theory of classical hypergeometric polynomials. To appear in Trans. Amer. Math. Soc. 6. R.V. Wallisser, On entire functions assuming integer values in a geometric sequence, Théorie des nombres (Quebec, PQ, 987), , de Gruyter, Berlin, 989. Department of Mathematics, Universidad Autónoma de Madrid. address: bernardo.lopez@uam.es address: javier.parcet@uam.es