arxiv: v3 [math.ca] 26 Jun 2013

Transcription

1 arxiv:0.0039v3 [math.ca] 6 Jun 03 Generalizations and specializations of generating functions for Jacobi, Gegenbauer, Chebyshev and Legendre polynomials with definite integrals Howard S. Cohl and Connor MacKenzie Applied and ComputationalMathematics Division, NationalInstitute ofstandards and Technology, Gaithersburg, Maryland, U.S.A. Department of Mathematics, Westminster College, 39 South Market Street, New Wilmington, Pennsylvania, U.S.A. Abstract In this paper we generalize and specialize generating functions for classical orthogonal polynomials, namely Jacobi, Gegenbauer, Chebyshev and Legendre polynomials. We derive a generalization of the generating function for Gegenbauer polynomials through extension a two element sequence of generating functions for Jacobi polynomials. Specializations of generating functions are accomplished through the re-expression of Gauss hypergeometric functions in terms of less general functions. Definite integrals which correspond to the presented orthogonal polynomial series expansions are also given.

2 Introduction This paper concerns itself with analysis of generating functions for Jacobi, Gegenbauer, Chebyshev and Legendre polynomials. The analysis involves generalization and specialization by re-expression of Gauss hypergeometric generating functions for these orthogonal polynomials. The generalizations that we present here are for two of the most important generating functions for Jacobi polynomials, namely [4, 4.3.-]. In fact, these are the first two generating functions which appear in Section 4.3 of [4]. As we will show, these two generating functions, traditionally expressed in terms of Gauss hypergeometric functions, can be re-expressed in terms of associated Legendre functions and also in terms of Ferrers functions, associated Legendre functions on the real segment,. Our Jacobi polynomial generating function generalizations, Theorem, Corollary and Corollary, generalize the generating function for Gegenbauer polynomials. The presented proofs of these generalizations rely upon the series re-arrangment technique. The motivation for the proofs of our generalizations was purely intuitive. Examination of the two Jacobi polynomial generating functions which we generalize, indicate that these generating functions represent two elements of an infinite sequence of eigenfunction expansions. The resulting formal proofs of our generalizations are then simple consequences. Our generalized expansions and hypergeometric orthogonal polynomial generating functions are given in terms of Gauss hypergeometric functions. The Gauss hypergeometric Jacobi polynomial generating functions which we generalize, as well as their eigenfunction expansion generalizations, are all re-expressible in terms of associated Legendre functions. Associated Legendre functions [9, Chapter 4] are given in terms of Gauss hypergeometric functions which satisfy a quadratic transformation of variable. These have an abundance of applications in Physics, Engineering and Applied Mathematics for solving partial differential equations in a variety of contexts. Recently, efficient numerical evaluation of these functions has been investigated in []. Associated Legendre functions are more elementary than Gauss hypergeometric functions because Gauss hypergeometric functions have three free parameters, whereas associated Legendre functions have only two. One can make this argument of elementarity with all of the functions which can be expressed in terms of Gauss hypergeometric functions, namely inverse trigonometric, inverse hyperbolic, exponential, logarithmic, Jacobi, Gegenbauer, and Chebyshev polynomials, and complete elliptic integrals. We summarize how associated Legendre functions, Gegenbauer, Chebyshev and Legendre polynomials and complete elliptic integrals of the first kind are interrelated. In a one-step process, we obtain definite integrals from our orthogonal polynomial expansions and generating functions. To the best of our knowledge our generalizations, re-expressions of Gauss hypergeometric generating functions for orthogonal polynomials and definite The interesting question of which orthogonal polynomial generating functions are important is addressed in [4, Section 4.3] see also R. A. Askey s MathSciNet review of Srivastava & Manocha 984 []. As well as their companion identities, see Section.

3 integrals are new and have not previously appeared in the literature. Furthermore, the generating functions presented in this paper are some of the most important generating functions for these hypergeometric orthogonal polynomials and any specializations and generalizations will be similarly important. This paper is organized as follows. In Sections, 3, 4, 5, 6, we present generalized and simplified expansions for Jacobi, Gegenbauer, Chebyshev of the second kind, Legendre, and Chebyshev of the first kind polynomials respectively. In A we present definite integrals which correspond to the derived hypergeometric orthgonal polynomial expansions. Unless stated otherwise the domains of convergence given in this paper are those of the original generating function and/or its corresponding definite integral. Throughout this paper we rely on the following definitions. Let a,a,a 3,... C. If i,j Z and j < i, then j n=i a n = 0 and j n=i a n =. The set of natural numbers is given by N := {,,3,...}, the set N 0 := {0,,,...} = N {0}, and Z := {0,±,±,...}. Let D := {z C : z < } be the open unit disk. Expansions over Jacobi polynomials The Jacobi polynomials P n α,β : C C can be defined in terms of the terminating Gauss hypergeometric series as follows [9, 8.5.7] P n α,β z := α+ n n,n+α+β + F n! α+ ; z for n N 0, and α,β > such that if α,β,0 then α+β + 0. The Gauss hypergeometric function F : C C\N 0 D C see Chapter 5 in [9] is defined as F a,b c ;z := a n b n c n z n n!, where the Pochhammer symbol rising factorial n : C C [9, 5..4] is defined by n z n := z +i, i= wheren N 0. Note that the Gausshypergeometricfunction can be analytically continued through for instance the Euler s integral representation for z C \, see for instance [, Theorem..]. Consider the generating function for Gegenbauer polynomials see 3 for their definition given by [9, 8..4], namely +ρ ρx ν = ρ n Cn ν x., 3

4 We attempt to generalize this expansion using the representation of Gegenbauer polynomials in terms of Jacobi polynomials given by [9, 8.7.], namely Cnx ν = ν n ν + n P ν /,ν / n x. By making the replacement ν / to α and β in using, we see that there are two possibilities for generalizing the generating function for Gegenbauer polynomials to a generating function for Jacobi polynomials. These two possibilities are given below, namely 3, 7. The first possibility is given for ρ D\,0] by [9, 8..3] +ρ α+β+ F = ρ+x α+β+, α+β+ β + ; ρ+x +ρ β/ Γβ + R α+ P β α ζ + = where we have used the definitions R = Rρ,x := +ρ ρx, ζ ± = ζ ± ρ,x := α+β + n β + n ρ n P α,β n x, 3 ±ρ +ρ ρx. Note that the restriction given by ρ D\,0] is so that the values of ρ are ensured to remain in the domain of P ν µ, but may otherwise be relaxed to D by analytic continuation if one uses the Gauss hypergeometric representation. The Ferrers function of the first kind representation given below provides the analytic continuation to the segment,0]. Here P ν µ : C\,] C is the associated Legendre function of the first kind see Chapter 4 in [9], which can be defined in terms of the Gauss hypergeometric function as follows [9, 4.3.6, 5.., 4.i] P µ ν z := Γ µ µ/ z + ν,ν + F z µ ; z. 4 The associated Legendre function of the first kind can also be expressed in terms of the Gauss hypergeometric function as see [9, 4.3.8, 4.iii], namely P ν µ µ z ν+µ ν µ z = Γ µz µ/, ν µ+ F ; µ z, 5 where argz < π. We have used 5 to re-express the generating function 3. We will refer to a companion identity as one which is produced by applying the map x x to an expansion over Jacobi, Gegenbauer, Chebyshev, or Legendre polynomials with argument x, in conjunction with the parity relations for those orthogonal polynomials. 4

5 In our first possibility for generalizing the generating function for Gegenbauer polynomials to a generating function for Jacobi polynomials, namely 3, we use the parity relation for Jacobi polynomials see for instance [9, Table 8.6.] P n α,β x = n P n β,α x, 6 and the replacement α,β β,α. This produces a companion identity which is the second possibility for generalizing the generating function for Gegenbauer polynomials to a generating function for Jacobi polynomials for ρ 0, by ρ α+β+ F = ρ x α+β+, α+β+ α+ ; ρ x ρ α/ Γα+ R β+ P α β ζ = α+β + n ρ n P n α,β x. 7 α+ n Note that the restriction given by ρ 0, is so that the values of ρ are ensured to remain in the domain of P µ ν, but may otherwise be relaxed to D by analytic continuation if one uses the Gauss hypergeometric representation. Here P µ ν :, C is the Ferrers function of the first kind associated Legendre function of the first kind on the cut through [9, 4.3.], defined as µ/ P µ +x ν,ν + νx := F ; x. 8 Γ µ x µ The Ferrersfunction ofthe first kind can alsobe expressedin terms ofthe Gauss hypergeometric function as see [7, p. 67], namely P µ ν x = µ x ν+µ ν µ Γ µ x µ/, ν µ+ F ; µ x, 9 for x 0,. We have used 9 to express 7 in terms of the Ferrers function of the first kind. One can easily see that 3 and 7 are generalizations of the generating function for Gegenbauer polynomials by taking α = β = ν /. The right-hand sides easily follow using the identification and the left-hand sides follow using [, ]. There exist natural extensions of 3, 7 in the literature see [4, 4.3.]. The extension corresponding to 7 is given for ρ 0, by α+β ++ρ α+β+, α+β+3 ρ α+β+ F ; ρ x α+ ρ α/ α+β ++ργα+ = ρ x R β+ P α β+ ζ 5

6 = n+α+β + α+β + n ρ n P n α,β x, 0 α+ n and its companion identity corresponding to 3, for ρ D\,0] is α+β + ρ α+β+, α+β+3 +ρ α+β+ F ; ρ+x β + +ρ β/ α+β + ργβ + = ρ+x R α+ P β α+ ζ + = n+α+β + α+β + n ρ n P n α,β x. β + n We have used 5, 9 to re-express these Gauss hypergeometric function generating functions as associated Legendre functions. We have not seen the companion identity in the literature, but it is an obvious consequence of [4, 4.3.] using parity. On the other hand, we have not seen the associated Legendre function representations of 0, in the literature. Upon examination of these two sets of generating functions, we suspected that these were just two examples of an infinite sequence of such expansions. This led us to the proof of the following theorem, which is a Jacobi polynomial expansion which generalizes the generating function for Gegenbauer polynomials. According to Ismail 005 [4, 4.3.], the generating functions 0,, their generalizations Theorem, Corollary and their corresponding definite integrals 45, 45, are closely related to the Poisson kernel for Jacobi polynomials, so our new generalizations will have corresponding applications. Theorem Let m N 0, α,β > such that if α,β,0 then α+β+ 0, x [,], ρ D\,0]. Then +x β/ R α+m+ P β α+mζ + = where a α,β n,m : D\,0] C is defined by a α,β n,m ρp α,β n x, a α,β n,m ρ := n+α+β +Γα+β +n+α+β +m+ n β/ Γβ +n+ +ρ ρ α+/ ρ mp α β n m. ρ Proof. Let ρ 0,ǫ with ǫ sufficiently small. Then using the definition of the following Gauss hypergeometric function α+β+m+, α+β+m+ F ; ρ+x β + +ρ 6

7 = α+β+m+ n α+β+m+ ρ n +x n n n!α+ n +ρ n, 3 the expansion of +x n in terms of Jacobi polynomials is given by +x n = n β + n n k=0 k n k α+β +k +α+β + k α+β + n+k+ β + k P α,β k x, 4 whose coefficients can be determined using orthogonality of Jacobi polynomials see Appendix combined with the Mellin transform given in [9, ]. By inserting 4 in the right-hand side of 3, we obtain an expansion of the Gauss hypergeometric function on the left-hand side of 3 in terms of Jacobi polynomials. By interchanging the two sums with justification by absolute convergence, shifting the n-index by k, and taking advantage of standard properties such as n k k = k n+k!, n! a n a n+k = a k a+k n, a+ = n a n, n n,k N 0, a C, produces a Gauss hypergeometric function as the coefficient of the Jacobi polynomial expansion. The resulting expansion is α+β+m+, α+β+m+ F ; ρ+x β + +ρ = f n,m α,β ρp n α,β x, 5 where f α,β n,m : 0,ǫ R is defined by f α,β n,m ρ := n+α+β +α+β + nα+β +m+ n ρ n β + n α+β + k+ +ρ n α+β+m+n+, α+β+m+n+ F α+β +n+ ; 4ρ +ρ. TheaboveexpansionisactuallyanalyticonD. However,ifweexpressitinterms of associated Legendre functions, then we must necessarily subdivide it into two regions. The Gauss hypergeometric function coefficient of this expansion, as well as the Gauss hypergeometric function on the left-hand of 5 are realized to be associated Legendre functions of the first kind through 5. Both sides of the resulting Jacobi polynomial expansion are analytic functions on ρ D\, 0]. Since we know that is valid for ρ 0,ǫ, then by the identity theorem for analytic functions, the equation holds on this domain. This completes the proof. 7

8 Note that the left-hand side of Theorem can be rewritten as α+β+m+, α+β+m+ ρ/ β/ Γβ ++ρ α+β+m+ F β + ; ρ+x +ρ. We have also derived the companion identity to, which we give in the following corollary. Corollary Let m N 0, α,β > such that if α,β,0 then α+β+ 0, x [,], ρ 0,. Then x α/ R β+m+ P α β+m ζ = where b α,β n,m : 0, R is defined by b α,β n,m ρpα,β n x, 6 b α,β n,m ρ := n+α+β +Γα+β +n+α+β +m+ n α/ Γα+n+ ρ ρ β+/ +ρ mp α β n m. +ρ Proof. We start with 5 and apply the parity relation for Jacobi polynomials 6. Let ρ 0,. The Gauss hypergeometric function coefficient of the Jacobi expansion is seen to be a Ferrers function of the first kind 8. After the application of the parity relation, the left-hand side also reduces to a Ferrers function of the first kind through 9. This completes the proof. Theorem generalizes 3,, while Corollary generalizes 7, 0. Both Theorem and Corollary generalize the generating function for Gegenbauer polynomials, which is its own companion identity.. Expansions and definite integrals from the Szegő transformation If one applies on the complex plane, the Szegő transformation conformal map z = +ρ ρ, 7 which maps a circle with radius less than unity to an ellipse with the foci at ± to the expansion in Theorem, then one obtains a new expansion. By [3, Theorem.7.3], this new Jacobi polynomial expansion is convergent for all x C within the interior of this ellipse. Applying 7 to yields the following corollary. 8

9 Corollary Let m N 0, α,β > such that if α,β,0 then α+β+ 0, x,z C, with z C \,] on any ellipse with the foci at ± and x in the interior of that ellipse. Then +x β/ z x α+m+/p β α+m +z z z = z z x where c α,β n,m : C\,] C is defined by c α,β n,m zpα,β n c α,β n,m z := n+α+β +Γα+β +n+α+β +m+ n β α m / Γβ +n+ z z m/ z + z mp α β n m x, 8 z +. z We would just like to briefly note that one may use the Szegő transformation 7 to obtain new expansion formulae and corresponding definite integrals from all the Jacobi, Gegenbauer, Legendre and Chebyshev polynomial expansions used in this paper. For the sake of brevity, we leave this to the reader. 3 Expansions over Gegenbauer polynomials The Gegenbauer polynomials C n µ : C C can be defined in terms of the terminating Gauss hypergeometric series as follows [9, 8.5.9] C n µ z := µ n n,n+µ F n! µ+ ; z, 9 for n N 0 and µ /, \ {0}. The Gegenbauer polynomials 9 are defined for µ /, \ {0}. However many of the formulae listed below actually make sense in the limit as µ 0. In this case, one should take the limit of the expression as µ 0 with the interpretation of obtaining Chebyshev polynomials of the first kind see 6 for the details of this limiting procedure. Corollary 3 Let m N 0, µ /, \{0}, x [,]. If ρ D\,0], then R µ+mcµ mζ + = Γµ+m +ρ m!ρ µ ρ m n+µµ+m n P µ n m C ρ nx, µ 0 and if ρ 0, then R µ+mcµ m ζ = Γµ+m m!ρ µ +ρ m n+µµ+m n P µ n m ρ C n µ +ρ x. 9

10 Proof. Using, substitute α = β = µ / along with and [9, 4.3.], namely P / µ n+µ / z = µ / Γµn! z µ/ /4 C µ πγµ+n n z. Through 4, we see that the Gauss hypergeometric function in the definition of the associated Legendre function of the first kind on the right-hand side is terminating and therefore defines an analytic function for ρ D. The analytic continuation to the segment ρ 0,] is provided by replacing the associated Legendre function of the first kind with the Ferrers function of the first kind with argument ρ/+ρ. As an example for re-expression using the elementarity of associated Legendre functions which was mentioned in the introduction, we now apply to two generating function results of Koekoek et al. 00 [6, 9.8.3] and Rainville 960 [0, 44.8]. Theorem Let λ C, µ /, \{0}, ρ 0,, x [,]. Then x /4 µ/ P / µ µ λ / R+ρ P/ µ µ λ / R ρ = / µ λ n µ λ n Γµ+ µ n Γµ+ C µ +nρµ /+n n x. Proof. [6, 9.8.3] give a generating function for Gegenbauer polynomials, namely λ,µ λ F ; R ρ λ,µ λ F µ+ ; R+ρ µ+ = λ n µ λ n µ n µ+ ρ n C nx. µ n Using [9, 5.8.7] to do a quadratic transformation on the Gauss hypergeometric functions and then using 5, with degree and order given by µ λ /, / µ, respectively and either z = R+ρ or z = R ρ. Simplification completes the proof. Theorem 3 Let α C, µ /, \{0}, ρ 0,, x [,]. Then x /4 µ/ R /+α µ P / µ µ α / ρx = ρ/µ / R Γµ+ α n µ n ρ n C µ nx. 3 0

11 Proof. On p. 79 of [0, 44.8] there is a generating function for Gegenbauer polynomials, namely α ρx α F, α+ ; ρ x α n µ+ ρx = ρ n C n µ µ x. n Using 5 to rewrite the Gauss hypergeometric function on the left-hand side of the above equation completes the proof. 4 Expansions over Chebyshev polynomials of the second kind The Chebyshev polynomials of the second kind can be obtained from the Gegenbauer polynomials using [9, 8.7.4], namely U n z = C nz, 4 for n N 0. Hence and through 9, the Chebyshev polynomials of the second kind U n : C C can be defined in terms of the terminating Gauss hypergeometric series as follows U n z := n+ F n,n+ 3 Corollary 4 Let m N 0, x [,]. If ρ D\,0] then R m+u mζ + = m+ ρ ρ m and if ρ 0, then R m+u mζ = m+ ρ+ρ m n+m+ n P n ; z. 5 m n+m+ n P n m +ρ ρ U n x, 6 ρ U n x. 7 +ρ Proof. Using 0, with 4 and π m+ U m z = z /4P / m+/ z, 8 which follows from 4, 5, and [9, 5.8.]. Analytically continuing to the segment ρ 0, completes the proof. Note that using [, 8.6.9], namely z Pν / /4 [ z = z + ν+/ z π ν + z + ] ν / z, and 8 one can derive the elementary function representation for the Chebyshev polynomials of the second kind [8,.5].

12 Corollary 5 Let ρ D, x [,]. Then x /4 P / / λ R+ρP / / λ R ρ = 5/ ρ π λ n λ n n ρ n U n x. 9 n+! Proof. Substituting µ = into, and using 4 with simplification, completes the proof. Corollary 6 Let α C, ρ 0,, x [,]. Then R / α x /4P / / α ρx = R ρ π α n n+! ρn U n x. 30 Proof. Using 4, and substituting µ = into 3 with simplification, produces this generating function for Chebyshev polynomials of the second kind. 5 Expansions over Legendre polynomials Legendre polynomials can be obtained from the Gegenbauer polynomials using [9, 8.7.9], namely P n z = Cn / z, 3 for n N 0. Hence and through 9, the Legendre polynomials of the second kind P n : C C can be defined in terms of the terminating Gauss hypergeometric series as follows P n z := F n,n+ ; z. 3 Using 3 we can write the previous expansions over Gegenbauer polynomials in terms of expansions over Legendre polynomials. Corollary 7 Let m N 0, x [,]. If ρ D\,0] then R m+p mζ + = ρ m ρ and if ρ 0, then R m+p mζ = +ρ m ρ n+m+ n P n m n+m+ n P n m +ρ P n x,33 ρ ρ P n x.34 +ρ

13 Proof. Using 0, substitute µ = / with 3 and P m z = Pm 0 z, which follows from 4, 3. Analytic continuation to ρ 0, completes the proof. Corollary 8 Let λ C, ρ {x C : z < }, x [,]. Then P λ R+ρP λ R ρ = λ n λ n n! ρ n P n x. 35 Proof. Substituting µ = / into and using 3 with simplification completes the proof. Note that Corollary 8 is just a restatement of [4, Theorem A], and therefore Theorem is a generalization of Brafman s theorem. Corollary 9 Let α C, ρ D, x [,]. Then ρx R α P α = R α n ρ n P n x. 36 n! Proof. Substituting µ = / in 3 with simplification completes the proof. As a further example of the elementarity of associated Legendre functions mentioned in the introduction, we apply to the recent generating function results of Wan & Zudelin 0 [4]. Theorem 4 Let x,y be in a neighborhood of. Then π n n! P n = π +xy xy +xy x +xy y n x+y xy x y x y+xy +xy x y K K x y x K K y x y K x K y +xyk x K y if x = y =, if x,y, if x and y, if x and y, if x,y. 3

14 Proof. If we start with 0 from [4], namely n n! P n n x+y xy x y x y+xy +xy = +xy F, ; x F, ; y, and use [9, 5.9.] we can express the Gauss hypergeometric functions as Legendre functions. For instance F, ; x = P / x, with x C\,0]. This domain is because the Legendre function of the first kind P ν and the Ferrers function of the first kind P ν, both with order µ = 0, are given by the same Gauss hypergeometric function and are continuous across argument unity cf. 4, 8. So there is no distinction between these two functions, except that the Ferrers function has argument on the real line with modulus less than unity and the Legendre function is defined on C \, both being well defined with argument unity. Hence there really is no need to use two different symbols to denote this function. The proof is completed by noting the two formulae P / z = π z + K z, z + P / x = x π K [, 8.3., 8.3.8], where K : [0, [π/, is the complete elliptic integral of the first kind defined by [9, 9..8] Kk := π F, ;k. Theorem 5 Let x,y be in a neighborhood of. Then = 3 3 n 3 n! np 3n 3n x+y x y x y x y. +4xyx+y +4xyx+y 4

15 = +4xyx+y if x = y = P /3 x 3 P /3 y 3 if x,y > P /3 x 3 P /3 y 3 if x < and y > P /3 x 3 P /3 y 3 if x > and y < P /3 x 3 P /3 y 3 if x,y < Proof. This follows by [4, ] and [9, 5.9.]. 6 Expansions over Chebyshev polynomials of the first kind The Chebyshev polynomials of the first kind T n : C C can be defined in terms of the terminating Gauss hypergeometric series as follows [7, p. 57] n,n T n z := F ; z, 37 for n N 0. The Chebyshev polynomials of the first kind can be obtained from the Gegenbauer polynomials using [, 6.4.3], namely T n z = ǫ n lim µ 0 n+µ µ Cµ nz, 38 where the Neumann factor ǫ n {,}, commonly seen in Fourier cosine series, is defined as ǫ n := δ n,0. Corollary 0 Let m N 0, x [,]. If ρ D\,0] then R mt mζ + = and if ρ 0, then ρ m R mt mζ = +ρ m Proof. Using 0, 38, and T m z = ǫ n m n P n m ǫ n m n P n m +ρ T n x, 39 ρ ρ T n x. 40 +ρ π z /4 P / m / z, 4 which follows from 4, 37, [9, 5.8.]. Analytic continuation to ρ 0, completes the proof. 5

16 Note that using [, 8.6.8], namely Pν / z = z /4[ z + ν+/ z + z + ] ν / z, π and 4 one can derive the elementary function representation for the Chebyshev polynomials of the first kind [5, p. 77]. A Definite integrals As a consequence of the series expansions given above, one may generate corresponding definite integrals in a one-step procedure as an application of the orthogonality relation for these hypergeometric orthogonal polynomials. We now describe this correspondence. Given an expansion over a set of orthogonal polynomials p n such that fx = a n p n x, 4 and the orthogonality relation where w :, [0,, then using 43 one has fxp n xwxdx = p n xp m xwxdx = c n δ n,m, 43 a m p n xp m xwxdx = a n c n, m=0 and therefore a n = fxp n xwxdx. 44 c n The definite integral expression 44 for the coefficient a n is of equal importance to the expansion 4, since one may use it to derive the other. Integrals of such sort are always of interest since they are very likely to find applications in applied mathematics and theoretical physics and could be included in tables of integrals such as [3]. We now give the orthogonality relations for the orthogonal polynomials used in the main text. For Jacobi, Gegenbauer, Chebyshev of the second kind, Legendre, and Chebyshev of the first kind polynomials, the orthogonality relations can be found in [9, 8.., 8..5, Table 8.3.]. Using the above procedure, we obtain the following definite integrals for products of special functions with Jacobi, Gegenbauer, Chebyshev and Legendre polynomials. Let m,n N 0, α,β > such that if α,β,0 then α + β + 0, ρ {z C : 0 < z < }\,0]. Then x α +x β/ R α+m+ P β α+m ζ +P n α,β xdx 6

17 Let ρ 0,. Then = α+β/+ Γα+n+α+β +m+ n P α β n n!ρ α+/ ρ m m x α/ +x β R β+m+ P α β+m ζ P n α,β xdx = α/+β+ Γβ +n+α+β +m+ n P α β n n!ρ β+/ +ρ m m Let µ /, \{0}, ρ {z C : 0 < z < }\,0]. Then x µ / C m µ ζ +C n µ xdx R µ+m = µ πγµ+nγµ+n+m m!n!γ µρ µ ρ m P n µ m +ρ. ρ ρ. +ρ +ρ. ρ A similar integral on ρ 0, can be obtained using. It should be noted that by using4, 3, 38, the previous definite integral over Gegenbauer polynomials can be written as an integral over Chebyshev polynomials of the first and second kind, and well as Legendre polynomials. Let λ C, ρ D. Then x µ/ /4 P / µ µ / λ R+ρP/ µ µ / λ R ρcµ nxdx = λ nµ λ n ρ n+µ / µ /. n+µγµµ+/ n n! Let ρ 0,. Then x /4 µ/ P / µ ρx R / µ+λ µ λ / C µ R nxdx = λ n πρ n+µ / / µ. n + µγµn! The previous two definite integrals over Gegenbauer polynomials can also be written as integrals over Chebyshev polynomials of the second kind and Legendre polynomials using 4, 3. Acknowledgements. The authors thank Hans Volkmer for valuable discussions. C. MacKenzie would like to thank the Summer Undergraduate Research Fellowship program at the National Institute of Standards and Technology for financial support while this research was carried out. We also would like to acknowledge two anonymous referees whose comments greatly improved the paper. This work was conducted while H. S. Cohl was a National Research Council Research Postdoctoral Associate in the Applied and Computational Mathematics Division at the National Institute of Standards and Technology, Gaithersburg, Maryland, U.S.A. 7

18 References [] M. Abramowitz and I. A. Stegun. Handbook of mathematical functions with formulas, graphs, and mathematical tables, volume 55 of National Bureau of Standards Applied Mathematics Series. U.S. Government Printing Office, Washington, D.C., 97. [] G. E. Andrews, R. Askey, and R. Roy. Special functions, volume 7 of Encyclopedia of Mathematics and its Applications. Cambridge University Press, Cambridge, 999. [3] I. S. Gradshteyn and I. M. Ryzhik. Table of integrals, series, and products. Elsevier/Academic Press, Amsterdam, seventh edition, 007. [4] M. E. H. Ismail. Classical and quantum orthogonal polynomials in one variable, volume 98 of Encyclopedia of Mathematics and its Applications. Cambridge University Press, Cambridge, 005. With two chapters by Walter Van Assche, With a foreword by Richard A. Askey. [5] A. Jakimovski, A. Sharma, and J. Szabados. Walsh equiconvergence of complex interpolating polynomials. Springer Monographs in Mathematics. Springer, Dordrecht, 006. [6] R. Koekoek, P. A. Lesky, and R. F. Swarttouw. Hypergeometric orthogonal polynomials and their q-analogues. Springer Monographs in Mathematics. Springer-Verlag, Berlin, 00. With a foreword by Tom H. Koornwinder. [7] W. Magnus, F. Oberhettinger, and R. P. Soni. Formulas and theorems for the special functions of mathematical physics. Third enlarged edition. Die Grundlehren der mathematischen Wissenschaften, Band 5. Springer- Verlag New York, Inc., New York, 966. [8] J. C. Mason and D. C. Handscomb. Chebyshev polynomials. Chapman & Hall/CRC, Boca Raton, FL, 003. [9] F. W. J. Olver, D. W. Lozier, R. F. Boisvert, and C. W. Clark, editors. NIST handbook of mathematical functions. Cambridge University Press, Cambridge, 00. [0] E. D. Rainville. Special functions. The Macmillan Co., New York, 960. [] B. I. Schneider, J. Segura, A. Gil, X. Guan, and K. Bartschat. A new Fortran 90 program to compute regular and irregular associated Legendre functions. Computer Physics Communications. An International Journal and Program Library for Computational Physics and Physical Chemistry, 8:09 097, 00. [] H. M. Srivastava and H. L. Manocha. A treatise on generating functions. Ellis Horwood Series: Mathematics and its Applications. Ellis Horwood Ltd., Chichester,

19 [3] G. Szegő. Orthogonal polynomials. American Mathematical Society Colloquium Publications, Vol. 3. Revised ed. American Mathematical Society, Providence, R.I., 959. [4] J. Wan and W. Zudilin. Generating functions of Legendre polynomials: a tribute to Fred Brafman. Journal of Approximation Theory, 644: , 0. 9