arxiv:math/ v1 [math.ca] 9 Oct 1995

Transcription

1 arxiv:math/9503v [math.ca] 9 Oct 995 UPWARD EXTENSION OF THE JACOBI MATRIX FOR ORTHOGONAL POLYNOMIALS André Ronveaux and Walter Van Assche Facultés Universitaires Notre Dame de la Paix, Namur Katholieke Universiteit Leuven March 3, 995 Abstract. Orthogonal polynomials on the real line always satisfy a three-term recurrence relation. The recurrence coefficients determine a tridiagonal semi-infinite matrix Jacobi matrix) which uniquely characterizes the orthogonal polynomials. We investigate new orthogonal polynomials by adding to the Jacobi matrix r new rows and columns, so that the original Jacobi matrix is shifted downward. The r new rows and columns contain r new parameters and the newly obtained orthogonal polynomials thus correspond to an upward extension of the Jacobi matrix. We give an explicit expression of the new orthogonal polynomials in terms of the original orthogonal polynomials, their associated polynomials and the r new parameters, and we give a fourth order differential equation for these new polynomials when the original orthogonal polynomials are classical. Furthermore we show how the orthogonalizing measure for these new orthogonal polynomials can be obtained and work out the details for a one-parameter family of Jacobi polynomials for which the associated polynomials are again Jacobi polynomials.. Introduction The construction of families of orthogonal polynomials on the real line from a given system of orthogonal polynomials or from a given orthogonalizing weight µ) has been the subject of various investigations [6], [4], [4], [], [3], [5], [7], [9], [6], [5]. Let P n n = 0,,,...) be a sequence of monic orthogonal polynomials on the real line, with orthogonality measure µ, then these polynomials satisfy a three-term recurrence relation.) P n+ x) = x b n )P n x) a np n x), n 0, with b n R and a n > 0 and initial conditions P 0 =, P = 0. The corresponding orthonormal polynomials are p n x) = a a...a n P n x), 99 Mathematics Subject Classification. 4C05. Key words and phrases. Orthogonal polynomials, Jacobi matrices, perturbations. The second author is a Senior Research Associate of the Belgian National Fund for Scientific Research

2 ANDRÉ RONVEAUX AND WALTER VAN ASSCHE and they satisfy the three-term recurrence relation.) xp n x) = a n+ p n+ x) + b n p n x) + a n p n x), n 0. Putting the recurrence coefficients in an infinite tridiagonal matrix gives the Jacobi matrix b 0 a b J = a b, or the symmetric Jacobi matrix b 0 a J s = a b a a b a Interesting new Jacobi matrices can be obtained by deleting the first r rows and columns, and the corresponding orthogonal polynomials are then the associated polynomials of order r, denoted by p r) n. Instead of deleting rows and columns we will add r new rows and columns at the beginning of the Jacobi matrix, thus introducing r new parameters. We call the corresponding new orthogonal polynomials anti-associated of order r and denote them by p r) n. In Section 4 we will give explicit formulas for these anti-associated polynomials in terms of the original orthogonal polynomials and the associated polynomials. We also give a fourth order differential equation in case the original orthogonal polynomials are classical. The construction of the orthogonality measure for the anti-associated polynomials is done in Section 5. The analysis requires the knowledge of the orthogonality measure of the original system and of the associated polynomials, and also the asymptotic behavior of both P n and P n ) as n. These things are known for a particular one-parameter family of Jacobi polynomials which we will call Grosjean polynomials after Grosjean who gave a number of interesting properties of this family in [9]) and which are considered in Section 3. The construction of new orthogonal polynomials by changing and shifting recurrence coefficients has been under investigation by others. If the original system of orthogonal polynomials has constant recurrence coefficients, then one essentially deals with Chebyshev polynomials of the second kind. Changing a finite number of the recurrence coefficients leads to Bernstein-Szegő weights [,.6], i.e., weights of the form x /ρx), where ρ is a positive polynomial on the orthogonality interval [, ]. Such finite perturbations have been considered among others by Geronimus [6], Grosjean [7], Sansigre and Valent [], and Zabelin [5]. The limiting case, when the number of changed recurrence coefficients tends to infinity with changes becoming smaller) is treated by Geronimo and Case [4], Dombrowski and Nevai [3], and in [4]. Instead of starting with constant recurrence coefficients one can also start with periodic recurrence coefficients. Finite perturbations of periodic recurrence coefficients are considered by Geronimus [6], Grosjean [8] and the limiting case by Geronimo and Van Assche [5]. Changing a finite number of recurrence

3 r) r) r) UPWARD EXTENSION OF A JACOBI MATRIX 3 Marcellán et al. [] and the limiting case by Nevai and Van Assche [5]. Associated polynomials correspond to deleting rows and columns of the Jacobi matrix, or equivalent to a positive shift in the recurrence coefficients. This situation is rather well known since it corresponds to numerator polynomials in Padé approximation. A general study of such a shift in the recurrence coefficients can be found in Belmehdi [] or Van Assche [3], and for Jacobi polynomials we refer to Grosjean [9] and Lewanowicz [0]. The case where both a shift in the recurrence coefficients is made together with changing a finite number of coefficients is treated by Nevai [4] and Peherstorfer [6]. In particular Nevai shows that when the original system is orthogonal on an interval with weight function w, then the orthogonal polynomials with weight function wx)/ f + x) on, where fz) = a + SBw, z) and f + x) = lim ǫ 0+ fx + iǫ), with SBw, z) the Stieltjes transform of the function Bw and B a polynomial of degree l, have recurrence coefficients which can be obtained by shifting the original recurrence coefficients and changing a finite number of the initial recurrence coefficients. Note that this generalizes the Bernstein-Szegő polynomials.. Background and notation Let us consider the family of monic orthogonal polynomials P n, n = 0,,,..., where P n has degree n, defined by the three-term recurrence relation.) P n+ x) = x b n )P n x) a n P n x), n 0, a n 0, with P = 0 and P 0 =. The sequences b n n 0) and a n n ) also generate the associated monic polynoimials of order r r a positive integer) P n r), n = 0,,,..., by the shifted recurrence P r) n+ x) = x b n+r)p n r) x) a n+r P r) x), n 0, with P r) r) = 0 and P 0 =. The recurrence relation.) can also be written in operator form using the non-symmetric Jacobi matrix J given by J = b 0 a b a b a r b r n a r+ b r , as x P = J P, where P = P 0, P, P, ) t. In the same way the recurrence relation defining the associated polynomials of order r is represented by

4 α α, α) 4 ANDRÉ RONVEAUX AND WALTER VAN ASSCHE where P r) = P r) 0, P r), P r), )t and J r) is the Jacobi matrix obtained by deleting the first r rows and columns of J. Observe that J r) = J for every integer r if and only if b n = b for n 0 and a n = a for n > 0. The corresponding family of orthogonal polynomials are, up to standardization of the orthogonality interval to [, ], the monic Chebyshev polynomials of the second kind u n x) = n U n x), for which b = 0 and a = /4. This fundamental family is a particular case α = β = /) of the Jacobi family P α,β) n x), which are orthogonal polynomials with respect to the weight function x) α + x) β on [, ]. Closely related to the family U n x) are the Chebyshev polynomials of the first kind T n x) α = β = /), of the third kind V n x) α = /, β = /) and of the fourth kind W n x) α = /, β = /). These monic Chebyshev polynomials satisfy the same recurrence relation.) P n+ x) = xp n x) 4 P n x), n, but with different initial conditions for P and P. These initial conditions are hidden in the Jacobi matrices J T, J V, J W for respectively the polynomials T n, V n and W n, which are explicitly 0 J T = J U ), J V = ) ), J 4 J W =, U 4 J U where J U is the Jacobi matrix for the Chebyshev polynomials of the second kind U n. These four families can be embedded in the more general two-parameter family represented by the Jacobi matrix ) J ) b =, c and by definition the corresponding monic polynomials P n satisfy the recurrence relation.) with initial conditions P x) = x b, P x) = xp x) cp 0 x). All the associated polynomials for r are the same: J U P r) n x) = u nx), n 0, r. 3. Grosjean polynomials The families {t n, n 0} and {u n, n 0} of Chebyshev polynomials have an interesting extension given by Grosjean. Let us consider first the monic Jacobi family for which the parameters satisfy α + β =, and define G α 0 =, G α nx) = c n P n α, α) x), < α < 0, where c n is a constant making this a monic polynomial, i.e., c n = n / ) n n. These are monic orthogonal polynomials with respect to the weight function 3.) w G x) = sin πα) ) α x, < x <. π + x + x A second family of monic Jacobi polynomials with α + β = is denoted by

5 UPWARD EXTENSION OF A JACOBI MATRIX 5 with d n = n ) n+, n which are monic orthogonal polynomials with weight function ) α sin πα x 3.) w g x) = + x), < x <. α α)π + x Obviously we have the special cases Gn / x) = t n x) and gn / x) = u n x). We will refer to the polynomials G α n and gn α as the Grosjean polynomials of the first and second kind respectively. This terminology is justified since Grosjean [9] showed, using direct verification, that 3.3) G α n) ) = g α n, and Ronveaux [9] gave the generalization of the differential link t n x) = nu n x) between the Chebyshev polynomials 3.4) G α n) x) = ) n ngn α x) = ng+α n. Recall that for α = / we have u n x) = ) n u n x). Grosjean actually shows [9, p. 75] that the only Jacobi polynomials for which the associated polynomials are again Jacobi polynomials are the Chebyshev polynomials of the first, second, third and fourth kind, and the Grosjean polynomials. Property 3.3) can easily be checked from the recurrence coefficients, which are 3.5a) b n = α + 4n, n + α)n α) a n =, for G α n ) n, and 3.5b) b n = α + 4n + ), a n = n + α)n + α) n + ), for g α n, see, e.g., Chihara [, p. 0]) by changing n to n + and α to α in 3.5a), which α,β) gives 3.5b). Property 3.4) follows by the differential property ˆP n ) x) = α+,β+) α,β) n ˆP n x) for the monic Jacobi polynomials ˆP n and the symmetry property p α,β) n x) = ) n p β,α) n x) see, e.g., Szegő [, p. 63]). It is already interesting to note that for the Grosjean polynomials the a n are rational function of n consisting of the ratio of two quadratic polynomials in n as in the Legendre case), whereas in general for the Gegenbauer family one has a ratio of two cubic polynomials in n and for the Jacobi polynomials one deals with quartic polynomials in the degree n. The pair of Grosjean polynomials also give an answer to the following question. Let D denote differentiation D = d dx and let L = σx)d + τx)d + λ n be the second order hypergeometric or degenerate hypergeometric) differential operator for the classical orthogonal polynomials Jacobi, Laguerre, Hermite, Bessel), and {P n, n 0} be the corresponding family of orthogonal polynomials, where P n corresponds with the eigenvalue λ n = n[n )σ + τ ]/, so that L P n ) = 0. If L is the formal adjoint of L,

6 6 ANDRÉ RONVEAUX AND WALTER VAN ASSCHE then for which orthogonal family of polynomials {Pn, n 0} does one have L P n ) = 0? Of course, the Legendre polynomials σ = τ) solve this problem since then L = L. For Grosjean polynomials the operator L is and L G,α,n = x )D + α x)d + n, L G,α,n G α n) = 0, L g,α,n = x )D + α 3x)D + nn + ), L g,α,n g α n ) = 0, and we have L G,α,n = L g, α,n, Therefore, if P n = G α n, then P n = g α n. From 3.4) we also have D L G,α,n = L g,+α,n D. The relative position of the zeros x α,n < x α,n < < x α n,n of G α n compared to the zeros x j,n = cos j )π of T n n is controlled by a classical comparison theorem due to Markov [, Theorem 6..]. The ratio between the Chebyshev weight w T x) =, < x <, π x and the weight w G for the Grosjean polynomials G α n given by 3.) is ) α+/ w T x) + x w G x) = const. x Now + x)/ x) is an increasing function of x on the interval [, ], hence for / < α < 0 the ratio w T /w G is an increasing function on [, ], and consequently we have the following inequalities x α j )π j,n < cos, j =,,..., n, / < α < 0. n For < α < / the ratio w T /w G is decreasing, and thus the inequalities for the zeros are reversed x α j )π j,n > cos, j =,,..., n, < α < /. n Similar conclusions can be made for the zeros yj,n α j =,,..., n) of Grosjean polynomials gn α jπ of the second kind as compared to the zeros cos n+ j =,,..., n) of the Chebyshev polynomials of the second kind. The ratio of the two weights is ) α / w U x) + x w g x) = const., x and thus for / < α < this ratio is an increasing function so that yj,n α jπ < cos, j =,,..., n, / < α <. n + For < α < / the inequalities are reversed yj,n α > cos jπ, j =,,..., n, < α < /. n + Finally we note that when we are dealing with Grosjean polynomials of the first kind, the product of the weight function of the orthogonal polynomials and the weight function of the associated polynomials is constant, as in the case of Chebyshev polynomials of the first kind. This is not true when dealing with the

7 r) r) r) UPWARD EXTENSION OF A JACOBI MATRIX 7 4. Anti-associated orthogonal polynomials The situation described in suggests to construct new families of orthogonal polynomials, which we will denote by P r) n+r, obtained by pushing down a given Jacobi matrix and by introducing in the empty upper left corner new coefficients b i i = r, r,..., ) on the diagonal and new coefficients a i 0 i = r, r,..., 0) on the lower subdiagonal. The new Jacobi matrix is then of the form b r a r+ b r+ a r+ b r+ J r) = a b 0 a b We will call the orthogonal polynomials P n r) for this Jacobi matrix J r) antiassociated polynomials for the family P n. They contain r new parameters and satisfy [P r) n+r x)]k) = P k r) n+r x). For r = we have and for r = we have J ) b = J a 0 ), b 0 J ) = a b. 0 a 0 J This new family of anti-associated polynomials P n r) can easily be represented as a combination of the original family P n and the associated polynomials P ) n. First, denote by Q n the orthogonal polynomials for the finite Jacobi matrix b r a r+ b r+ a r+ b r ,. a b so that they satisfy Q 0 =, Q = 0 and Then, clearly Q n+ x) = x b r+n )Q n x) a r+n Q n x), n r. P r) n x) = Q n x), 0 n r. For n > r the anti-associated polynomials satisfy the three-term recurrence relation

8 r) 8 ANDRÉ RONVEAUX AND WALTER VAN ASSCHE so that P r) n+r x) is a solution of the three-term recurrence relation of the original family P n x). The initial conditions however are P r r) x) = Q r x) and P r) r x) = Q r x), and since every solution of the three-term recurrence relation.) is a linear combination of P n x) and P ) n x), we have P r) n+r x) = AP nx) + BP ) n x), n 0. Using the initial conditions for n = 0 and n = gives A = Q r x) and B = a 0 Q r x), and thus we have 4.) P r) n+r x) = Q rx)p n x) a 0 Q r x)p ) n x), n 0. From this representation it is easy to construct a differential equation satisfied by the family P r) n+r if the original family P n is itself solution of a linear differential equation of second order, for instance when P n are the classical polynomials Jacobi, Laguerre, Hermite, Bessel), then they are a solution of the hypergeometric differential equation ) n 4.) L y σx)y + τx)y + n )σ τ n y = 0. The techniques used in [7], [8], [9], [0] and the fact that 4.3) L P ) n = σ τ )P n easily give a fourth order differential equation satisfied by the anti-associated polynomials P r) n+r x). The general technique is the following. Let Bx) = a 0 Q r x) and put Jx) = Bx)P ) n x). Then we can transform the equation 4.3) to L and introducing the differential operator J B = σ τ )P n, R = σb D + [σ τ)b σbb ]D + σb ) σbb BB σ τ) [ σ n n ) + τ + n)]b this is equivalent to 4.4) R J = σ τ )B 3 P n. From 4.) we then have R P r) n+r = R [Q r P n ] σ τ )a 6 0Q 3 r P n, and eliminating the second derivative using 4.) then leads to

9 UPWARD EXTENSION OF A JACOBI MATRIX 9 where M 0 and N 0 are polynomials. Taking the derivative in 4.5) and using 4.) to eliminate P n also gives σ[r P r) n+r ] = M P n + N P n, where M and N are polynomials, and repeating this also gives This shows that 4.6) det σ[σ[r P r) n+r ] ] = M P n + N P n. R P r) n+r M 0 N 0 σ[r P r) n+r ] M N σ[σ[r P r) n+r ] ] M N = 0, which is the desired fourth order differential equation. When P n x) = G α n x), then σ τ = 0 so that 4.4) simplifies and becomes homogeneous. The differential equation 4.6) however remains one of the fourth order, except when α = /, because then P r) n+r x) = A 0x)T n x) + B 0 x)t n x), where A 0 and B 0 are polynomials. Similar as in the above reasoning we then get a second order differential equation det P r) n+r A 0 B 0 σ[p r) n+r ] A B σ[σ[p r) n+r ] ] A B which is equivalent to the equation given in []. = 0, 5. Construction of the orthogonality measure If we use probability measures throughout the analysis, and if we use lower case p and q for the orthonormal polynomial, then we have p n = γ n P n, where Similarly and p ) n = γ) n P ) n polynomials are q n x) = where γ) n γ n = a a a n ). Q n x) a r+ a r+ a r+n, n r, = a γ n. Thus the orthonormal anti-associated p r) n+r x) = a r+ a 0 ) γ n P r) n+r x), and using this in 4.) gives 5.) p r) n+r x) = q rx)p n x) a 0 a q r x)p ) n x). The orthonormal polynomials are useful in obtaining the weight function for the anti-associated polynomials. Indeed, we can compute the weight function using Christoffel functions n λ n x) = p jx) j=0

10 0 ANDRÉ RONVEAUX AND WALTER VAN ASSCHE Theorem MNT. Suppose p n x) are orthonormal polynomials for a measure µ on [, ] and let µ belong to Szegő s class, i.e., Then log µ x) x dx >. 5.) lim n nλ nx) = πµ x) x holds almost everywhere on [, ]. The Szegő condition can be relaxed and in fact it suffices to assume that µ is a regular measure on [, ] i.e., suppµ) = [, ] and lim n γn /n = ) and b a log µ x) dx > in order that 5.) holds almost everywhere on [a, b], ) [, Thm. 8]. We will assume these conditions for the measure µ and moreover we allow the addition of a finite number of mass points to µ. Then 5.) will still hold almost everywhere on [a, b]. Indeed, if we add a mass point c to µ then by the extremum property λ n x; µ) = min qnt) dµt), q n x)= where the minimum is taken over all polynomials q n of degree at most n which take the value at the point x, we see that for the measure µ c = µ + ǫδ c i.e., the measure µ to which we add a mass point at c with mass ǫ) ) λ n x; µ c ) = min qnt) dµt) + ǫqnc) q n x)= min q n x)= = λ n x; µ), qn t) dµt) so that lim inf n nλ n x; µ) πµ x) x, almost everywhere on [a, b]. On the other hand, consider the polynomial q n t) = t c)r n t)/x c), where r n is the minimizing polynomial for the measure dνt) = t c) dµt), then for x c λ n x; µ c ) x c) = x c) t c) rn t) dµt) rn t) dνt) = x c) λ n x; ν) so that lim sup n nλ n x; µ c ) lim n x c) nλ n x; ν) = πµ x) x almost everywhere on [a, b], since ν is regular and satisfies the Szegő condition on [a, b]. Combined with the previous inequality this gives

11 UPWARD EXTENSION OF A JACOBI MATRIX almost everywhere on [a, b], independent of the mass point c. This procedure can be repeated, so that adding a finite number of mass points does not change the behavior in 5.). In the previous discussion we wanted to allow the measure µ to have a larger support than [, ] by allowing mass points outside [, ]) but in such a way that 5.) still holds. Alternative we can use a weaker result by Nevai [3, Thm. 54, p. 04] in which the support of µ is allowed to be [, ] E, where E contains at most a denumerable number of points which can only accumulate at ±. Theorem N. Suppose µ M0, ), i.e., the recurrence coefficients a n and b n have asymptotic behavior given by Then lim a n = /, n holds for almost every x suppµ). lim b n = 0. n lim supnλ n x; µ) = πµ x) x n With this result, if we can compute the limit of nλ n x; µ) for x, ), then this limit is almost everywhere equal to πµ x) x, which allows us to compute the weight µ on, ) without having to worry about the fact that suppµ) may be larger than [, ]. The existence of such a limit is however not guaranteed by Nevai s theorem, contrary to the theorem of Máté, Nevai and Totik where the existence of the limit is in the conclusion of the theorem but where suppµ) = [, ] is required. If we assume that the original family p n belongs to a measure µ in the class M0, ), then the anti-associated polynomials p r) n will also have a measure µ r) which belongs to the class M0, ), and thus we know that µ r) has support [, ] E, where E is at most denumerable with the only accumulation points at ± [4, Thm., p. 437]. In fact, there can be at most r mass points, since the associated polynomials of order r of the anti-associated polynomials p r) n are again p n and these have all their zeros inside [, ], and the interlacing property of orthogonal polynomials and associated orthogonal polynomials shows that adding one row and one column in the Jacobi matrix to form the anti-associated polynomials p n ) can at most add a mass point to the left of and to the right of. Adding r rows and columns thus can add at most r mass points to the left of and r mass points to the right of. If the original measure µ belongs to the Szegő class which is a subclass of M0, )), then the measure µ r) restricted to [, ] also belongs to the Szegő class. This is so because if we denote by µ r) [,] the restriction of µ r) to the interval [, ], then lim n γ n µ r) )/γ n µ r) [,]) exists and is strictly positive [3, Thm. 5 on p. 36] and since γ n µ r) ) = γ n µ)a r+ a 0 ) and µ belongs to the Szegő class, which is equivalent with the statement that lim n γ n / n exists and is strictly positive, it follows that lim n γ n µ r) [,] )/n exists and is strictly positive, so that µ r) [,] belongs to the Szegő class. In order to use the result in 5.) we observe that for n r n [p r) j x)] = r q j x) + n r [p r) j+r x)],

12 ANDRÉ RONVEAUX AND WALTER VAN ASSCHE and using 5.) this gives n r [p r) j=0 n r j+r x)] = qrx) k=0 p kx) + a0 a k= ) n r qr x) [p ) k x)] k=0 a n r 0 q r x)q r x) p k x)p ) k a x). In order to be able to compute these sums, we will need to be able to compute Christoffel functions of associated polynomials and sums of mixed form containing the product p n x)p ) n x). This is in general not so easy. However, when we use Grosjean polynomials of the first kind, then the associated polynomials are Grosjean polynomials of the second kind, so that we are always dealing with Jacobi polynomials, for which we can compute these sums, at least as n becomes large. For the anti-associated polynomials corresponding to Grosjean polynomials we thus can prove the following result. Theorem. If the original system of orthogonal polynomials consists of Grosjean polynomials of the first kind, then the anti-associated polynomials p r) n are orthogonal with respect to a measure µ r) which is absolutely continuous on [, ] with density w r x) = sin πα) x) α π + x) α+ q rx) a 0 e iαπ x)α q r x) + x) α+. In addition, there may be at most r mass point outside [, ] r to the left of and r to the right of ). These mass points are the roots of the equation q r x) a 0 q r x) sign x x α = 0. x + α+ Proof. If the original family consists of Grosjean polynomials, then by 4.3) and 3.) we have for < x < n r 5.3) lim p n n kx) = k=0 x + x) α+ sin πα) x) α, and this even holds uniformly on closed subintervals of, ), and similarly by using 3.) with α replaced by α) n r 5.4) lim n n k=0 [p ) k x)] = α + α) x sin πα) x) α + x) α+. For the remaining sum, which consists of a mixture of the original polynomials and the associated polynomials, we can use Darboux s extension of the Laplace-Heine formula for Legendre polynomials [, Thm on p. 96] 5.5) nπp α,β) n cos θ) = cos ) α / cos θ ) β / [n + α + β + )/]θ α + / ) π sin θ + O/n),

13 UPWARD EXTENSION OF A JACOBI MATRIX 3 which holds uniformly for x = cos θ on closed intervals of, ). For Grosjean polynomials P n x) = G α n x) we have nπ p n x) = [ + O/n)] sin πα) P n α, α) x), p ) n x) = [ + O/n)] and hence for x = cos θ with ǫ θ π ǫ p n x) = sin πα) and p ) n x) = Using this gives nπα + α) sin πα) sin θ ) α / cos θ α+/ cos ) P α,+α) n x), α + α) sin θ ) α / cos θ α 3/ cos sin πα) ) n r lim p k x)p ) α + α) n k n x) = sin πα) sinθ/ cosθ/ k= n r lim cos kθ α + / ) π cos n n k= Using cosacos b = cosa + b) + cosa b) gives n r n k= cos kθ α + / π ) cos kθ k= ) α + / π nθ α + / π ) + O/n) nθ + α / π kθ ) + O/n). ) α + / π. = n r [cos kθ π/) + cos απ] n and thus uniformly for x on closed subsets of, ) we have n r 5.6) lim p k x)p ) α + α) cos πα n k n x) =. sin θ sin πα) k= Alternatively, 5.5) can also be obtained by using a Christoffel-Darboux type formula [, corollary.] which for orthonormal polynomials is n [ ] pk x)p ) k x) = a n+ p n+ x)p) n x) p n x)p) n x).

14 4 ANDRÉ RONVEAUX AND WALTER VAN ASSCHE When the original system consists of Grosjean polynomials of the first kind G α n one knows by 3.3) that the associated polynomials are Grosjean polynomials of the second kind gn α and by 3.4) we know that the derivative of G α n is ng+α n. For the orthonormal polynomials this gives nπ p n x) = [ + O/n)] sin πα) P n α, α) x), p ) n x) = [ + O/n)] and hence using 5.5) gives lim n n n k= nπα + α) sin πα) P n α,+α) x), p n x) = [ + O/n)] πn n sin πα) P +α, α) n x), p k x)p ) k x) α + α) = sin πα) sin θ lim cos n [ cos nθ α + 3/ π n + )θ α + 3/ ) cos n + )θ ) π cos nθ )] α + / π ) α + 3/ π Simple trigonometry then gives 5.6). Combining the limiting relations 5.3), 5.4), and 5.6) gives uniformly for x on closed intervals of, ) lim n n n j=0 [p r) j x)] = r + x)α+ [q x) x sin πα) x) α Observe now that a = α + α), thus we find lim n n n j=0 [p r) j x)] a 0 q r x)q r x) cosαπ α + α) a ) ] a0 + qr x)α x) α) + α) + x) α+ a + x) α+ = sin πα) x x) α q rx) a 0 e iαπ x)α q r x) + x) α+ and this holds uniformly on closed intervals [a, b], ). From 5.) we can then conclude that the orthogonality measure µ r) for the anti-associated polynomials has an absolutely continuous part on, ) with weight function 5.7) w r x) = sin πα) x)α qr x) a 0 e iαπ q r x) x)α,,

15 UPWARD EXTENSION OF A JACOBI MATRIX 5 which is the weight function of the original system p n x) Grosjean polynomials of the first kind) divided by a positive factor containing the new parameters b r,..., b and a r+,..., a, a 0. This factor cannot vanish on, ), because then both the real part and the imaginary part of q r x) a 0 e iαπ x)α q r x) + x) α+ need to vanish. The imaginary part can only vanish when q r x) = 0, and assuming this, the real part can then only vanish when also q r x) = 0. This is impossible since two consecutive orthogonal polynomials cannot have a common zero. Since the original orthogonal polynomials are Jacobi polynomials, they will belong to the class M0, ). Moreover, the original system satisfies 4a k+ + b k ) <, k=0 and hence also the new system of anti-associated polynomials satisfies this trace class condition. But then we know [4, Theorem 6] that the orthogonality measure µ r) is absolutely continuous on, ) and we have obtained the weight function in 5.7). The mass points outside, ) are those points x R \, ) for which k=0 [p r) k ] <. This means that at a mass point x we have p n x) 0, and from 5.) this implies that lim p nx) n [ n x) q r x) a 0 q r x) p) a p n x) ] = 0. For x / [, ] we know that p n x) increases exponentially fast, hence at a mass point we always have q r x) a 0 p ) n q r x) lim x) a n p n x) = 0. The limit of the ratio p ) n x)/p nx) can be found by using Markov s theorem, which states that this limit is the Stieltjes transform of the measure µ, or it can be obtained from Darboux s generalization of the Laplace-Heine formula for Legendre polynomials for x outside [, ] [, Theorem 8..7]. Both methods give p ) n lim x) n a p n x) = x )α x + ) α+, where the right hand side is to be taken positive if x > and negative if x <. This means that a mass point satisfies x )α q r x) a 0 q r x) = 0. x + ) α+

16 6 ANDRÉ RONVEAUX AND WALTER VAN ASSCHE Remark: In order to determine the measure µ r) for the anti-associated polynomials corresponding to Grosjean polynomials, one can also use the results from Theorem 3.9 in Peherstorfer [6] using the Cauchy principal value and the Stieltjes transform. Indeed, the Grosjean polynomials are one of the seldom cases where a nice explicit expression for the Cauchy principal value and the Stieltjes transform exists and for this reason the results of [6] can be applied without problems. Observe that this technique is basically also the one used by Grosjean in [9]. Our approach using Christoffel functions has the advantage that it advoids taking boundary values of a Stieltjes transform or evaluating a Cauchy principal values and uses only information of the orthogonal polynomials on the interval [, ]. In particular our method also works when appropriate asymptotic information of the orthogonal polynomials and the associated orthogonal polynomials on the real line is available. 6. Examples The simplest examples occur when we take α = /, in which case the original system consists of Chebyshev polynomials of the first kind. The weight function for the anti-associated polynomials then becomes w r x) = π x q r x) ia 0 q r / x = x π x )qr x) + a 0 q r x). Hence this weight function is the weight function of Chebyshev polynomials of the second kind, divided by a polynomial which is positive on [, ]. Such orthogonal polynomials are known as Bernstein-Szegő polynomials [,.6]. In case the original system consists of Chebyshev polynomials of the second kind, for which all the recurrence coefficients are constant, we have p n = U n and p ) n = U n. Relation 5.) then becomes p r) n+r x) = q rx)u n x) a 0 q r x)u n x). Using U n x) = sinn + )θ/ sinθ x = cos θ), one easily shows for < x < and and hence n n j=0 [p r) j x)] = lim n n lim n n n j=0 U j x) = x ) n U j x)u j x) = j= and hence the weight function becomes x x ), [ q x ) r x) 4a 0 xq r x)q r x) + 4a 0qr x) ], wx) = x,

17 UPWARD EXTENSION OF A JACOBI MATRIX 7 i.e., this is again a Bernstein-Szegő weight. Observe that it can be written as wx) = x π q r x) a 0 e iθ q r x). For the mass points we see that they can occur only for x / [, ] when U n x) q r x) a 0 q r x) lim n U n x) = 0, and since U n x)/u n x) /x + x ) for x / [, ], where the limit is to be taken positive for x > and negative for x <, it follows that x is a mass point only when 6.) q r x) a 0 q r x) x + x = 0. If we multiply the left hand side by q r x) a 0 q r x)/x x ), then this implies q rx) 4a 0 xq r x)q r x) + 4a 0q r x) = 0. So the mass points are zeros of the polynomial in the denominator of the weight function, but only those zeros for which 6.) holds. For r = we can consider the following cases: ) a 0 = / and b = 0. In this case the weight function of the antiassociated polynomials is wx) = /π x ) and thus we have the Chebyshev polynomials of the first kind. ) a 0 = / and b = /. Then the weight function is the one for Jacobi polynomials with α = / and β = / and thus we have Chebyshev polynomials of the third kind V n x). 3) Finally when a 0 = / and b = / the weight function becomes the one for Jacobi polynomials with α = / and β = /, so that we have Chebyshev polynomials of the fourth kind W n x). References. S. Belmehdi, On the associated orthogonal polynomials, J. Comput. Appl. Math ), T. S. Chihara, An Introduction to Orthogonal Polynomials, Gordon and Breach, New York, J. Dombrowski and P. Nevai, Orthogonal polynomials, measures and recurrence relations, SIAM J. Math. Anal ), J. S. Geronimo and K. M. Case, Scattering theory and polynomials orthogonal on the real line, Trans. Amer. Math. Soc ), J. S. Geronimo and W. Van Assche, Orthogonal polynomials with asymptotically periodic recurrence coefficients, J. Approx. Theory ), Ya. L. Geronimus, On some finite difference equations and corresponding systems of orthogonal polynomials, Zap. Mat. Otd. Fiz.-Mat. Fak. i Kharkov Mat. Obsc. 5 4), in Russian) 7. C. C. Grosjean, The measure induced by orthogonal polynomials satisfying a recursion formula with either constant or periodic coefficients, part I: constant coefficients, Acad. Analecta,

18 8 ANDRÉ RONVEAUX AND WALTER VAN ASSCHE 8. C. C. Grosjean, The measure induced by orthogonal polynomials satisfying a recursion formula with either constant or periodic coefficients, part II: pure or mixed periodic coefficients, Acad. Analecta, Koninkl. Acad. Wetensch. Lett. Sch. Kunsten België 48 5) 986), C. C. Grosjean, The weight functions, generating functions and miscellaneous properties of the sequences of orthogonal polynomials of the second kind associated with the Jacobi and Gegenbauer polynomials,, J. Comput. Appl. Math. 6, S. Lewanowicz, Properties of the polynomials associated with the Jacobi polynomials, Math. Comp ), F. Marcellán, J. S. Dehesa and A. Ronveaux, On orthogonal polynomials with perturbed recurrence relations, J. Comput Appl. Math ), 03.. A. Máté, P. Nevai and V. Totik, Szegő s extremum problem on the unit circle, Ann. of Math ), P. G. Nevai, Orthogonal Polynomials, Memoirs Amer. Math. Soc. 3, Providence, RI, P. Nevai, A new class of orthogonal polynomials, Proc. Amer. Math. Soc ), P. Nevai and W. Van Assche, Compact perturbations of orthogonal polynomials, Pacific J. Math ), F. Peherstorfer, Finite perturbations of orthogonal polynomials, J. Comput. Appl. Math ), A. Ronveaux, Fourth order differential equations for numerator polynomials, J. Phys. A 988), A. Ronveaux, 4th order differential equations and orthogonal polynomials of the Laguerre- Hahn class, in Orthogonal Polynomials and their Applications C. Brezinski et al.), eds.), IMACS Annals of Computing and Applied Mathematics, vol. 9, J. C. Baltzer AG, Basel, 99, pp A. Ronveaux, On a first integral relating classical orthogonal polynomials and their numerator polynomials, Simon Stevin 66 99), A. Ronveaux, A. Zarzo and E. Godoy, Fourth order differential equations satisfied by the generalized co-recursive of all classical orthogonal polynomials. A study of their distribution of zeros, J. Comput. Appl. Math. to appear).. G. Sansigre and G. Valent, A large family of semi-classical polynomials: the perturbed Tchebichev, J. Comput. Appl. Math. to appear).. G. Szegő, Orthogonal Polynomials, Amer. Math. Soc. Colloq. Publ., vol. 3, Providence, RI, 975 4th edition). 3. W. Van Assche, Orthogonal polynomials, associated polynomials and functions of the second kind, J. Comput. Appl. Math ), W. Van Assche, Asymptotics for orthogonal polynomials and three-term recurrence, in Orthogonal Polynomials: Theory and Practice P. Nevai, ed.), NATO ASI Series C, vol. 94, Kluwer, Dordrecht, 990, pp V. V. Zabelin, Polynomial properties, Vestnik Moskov. Univ. Ser., Mat. Meh. nr. 989), in Russian); Mosc. Univ. Math. Bull ), Département de Physique, Facultés Universitaires N.D.P., Rue de Bruxelles 6, B-5000 Namur BELGIUM) address: aronveaux@cc.fundp.ac.be Department of Mathematics, Katholieke Universiteit Leuven, Celestijnenlaan 00 B, B-300 Heverlee BELGIUM) address: Walter.VanAssche@wis.kuleuven.ac.be