arxiv:math/ v1 [math.ca] 9 Oct 1995
|
|
- Alisha Hodge
- 5 years ago
- Views:
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
ON LINEAR COMBINATIONS OF
Física Teórica, Julio Abad, 1 7 (2008) ON LINEAR COMBINATIONS OF ORTHOGONAL POLYNOMIALS Manuel Alfaro, Ana Peña and María Luisa Rezola Departamento de Matemáticas and IUMA, Facultad de Ciencias, Universidad
More informationLAMÉ DIFFERENTIAL EQUATIONS AND ELECTROSTATICS
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 00, Number 0, Xxxx XXXX, Pages 000 000 S 0002-9939(XX)0000-0 LAMÉ DIFFERENTIAL EQUATIONS AND ELECTROSTATICS DIMITAR K. DIMITROV AND WALTER VAN ASSCHE
More informationPositivity of Turán determinants for orthogonal polynomials
Positivity of Turán determinants for orthogonal polynomials Ryszard Szwarc Abstract The orthogonal polynomials p n satisfy Turán s inequality if p 2 n (x) p n 1 (x)p n+1 (x) 0 for n 1 and for all x in
More informationEntropy of Hermite polynomials with application to the harmonic oscillator
Entropy of Hermite polynomials with application to the harmonic oscillator Walter Van Assche DedicatedtoJeanMeinguet Abstract We analyse the entropy of Hermite polynomials and orthogonal polynomials for
More informationORTHOGONAL POLYNOMIALS WITH EXPONENTIALLY DECAYING RECURSION COEFFICIENTS
ORTHOGONAL POLYNOMIALS WITH EXPONENTIALLY DECAYING RECURSION COEFFICIENTS BARRY SIMON* Dedicated to S. Molchanov on his 65th birthday Abstract. We review recent results on necessary and sufficient conditions
More informationOrthogonal polynomials
Orthogonal polynomials Gérard MEURANT October, 2008 1 Definition 2 Moments 3 Existence 4 Three-term recurrences 5 Jacobi matrices 6 Christoffel-Darboux relation 7 Examples of orthogonal polynomials 8 Variable-signed
More informationBounds on Turán determinants
Bounds on Turán determinants Christian Berg Ryszard Szwarc August 6, 008 Abstract Let µ denote a symmetric probability measure on [ 1, 1] and let (p n ) be the corresponding orthogonal polynomials normalized
More informationInterpolation and Cubature at Geronimus Nodes Generated by Different Geronimus Polynomials
Interpolation and Cubature at Geronimus Nodes Generated by Different Geronimus Polynomials Lawrence A. Harris Abstract. We extend the definition of Geronimus nodes to include pairs of real numbers where
More informationSimultaneous Gaussian quadrature for Angelesco systems
for Angelesco systems 1 KU Leuven, Belgium SANUM March 22, 2016 1 Joint work with Doron Lubinsky Introduced by C.F. Borges in 1994 Introduced by C.F. Borges in 1994 (goes back to Angelesco 1918). Introduced
More informationSZEGÖ ASYMPTOTICS OF EXTREMAL POLYNOMIALS ON THE SEGMENT [ 1, +1]: THE CASE OF A MEASURE WITH FINITE DISCRETE PART
Georgian Mathematical Journal Volume 4 (27), Number 4, 673 68 SZEGÖ ASYMPOICS OF EXREMAL POLYNOMIALS ON HE SEGMEN [, +]: HE CASE OF A MEASURE WIH FINIE DISCREE PAR RABAH KHALDI Abstract. he strong asymptotics
More informationHermite Interpolation and Sobolev Orthogonality
Acta Applicandae Mathematicae 61: 87 99, 2000 2000 Kluwer Academic Publishers Printed in the Netherlands 87 Hermite Interpolation and Sobolev Orthogonality ESTHER M GARCÍA-CABALLERO 1,, TERESA E PÉREZ
More informationORTHOGONAL POLYNOMIALS FOR THE OSCILLATORY-GEGENBAUER WEIGHT. Gradimir V. Milovanović, Aleksandar S. Cvetković, and Zvezdan M.
PUBLICATIONS DE L INSTITUT MATHÉMATIQUE Nouvelle série, tome 8498 2008, 49 60 DOI: 02298/PIM0898049M ORTHOGONAL POLYNOMIALS FOR THE OSCILLATORY-GEGENBAUER WEIGHT Gradimir V Milovanović, Aleksandar S Cvetković,
More informationA Birth and Death Process Related to the Rogers-Ramanujan Continued Fraction
A Birth and Death Process Related to the Rogers-Ramanujan Continued Fraction P. R. Parthasarathy, R. B. Lenin Department of Mathematics Indian Institute of Technology, Madras Chennai - 600 036, INDIA W.
More informationOrthogonal polynomials with respect to generalized Jacobi measures. Tivadar Danka
Orthogonal polynomials with respect to generalized Jacobi measures Tivadar Danka A thesis submitted for the degree of Doctor of Philosophy Supervisor: Vilmos Totik Doctoral School in Mathematics and Computer
More informationUNIFORM BOUNDS FOR BESSEL FUNCTIONS
Journal of Applied Analysis Vol. 1, No. 1 (006), pp. 83 91 UNIFORM BOUNDS FOR BESSEL FUNCTIONS I. KRASIKOV Received October 8, 001 and, in revised form, July 6, 004 Abstract. For ν > 1/ and x real we shall
More informationDifference Equations for Multiple Charlier and Meixner Polynomials 1
Difference Equations for Multiple Charlier and Meixner Polynomials 1 WALTER VAN ASSCHE Department of Mathematics Katholieke Universiteit Leuven B-3001 Leuven, Belgium E-mail: walter@wis.kuleuven.ac.be
More informationJacobi-Angelesco multiple orthogonal polynomials on an r-star
M. Leurs Jacobi-Angelesco m.o.p. 1/19 Jacobi-Angelesco multiple orthogonal polynomials on an r-star Marjolein Leurs, (joint work with Walter Van Assche) Conference on Orthogonal Polynomials and Holomorphic
More informationWeighted Sums of Orthogonal Polynomials Related to Birth-Death Processes with Killing
Advances in Dynamical Systems and Applications ISSN 0973-5321, Volume 8, Number 2, pp. 401 412 (2013) http://campus.mst.edu/adsa Weighted Sums of Orthogonal Polynomials Related to Birth-Death Processes
More informationQuadratures and integral transforms arising from generating functions
Quadratures and integral transforms arising from generating functions Rafael G. Campos Facultad de Ciencias Físico-Matemáticas, Universidad Michoacana, 58060, Morelia, México. rcampos@umich.mx Francisco
More informationCOMPACT PERTURBATIONS OF ORTHOGONAL POLYNOMIALS
PACIFIC JOURNAL OF MATHEMATICS Vol. 153, No. 1, 1992 COMPACT PERTURBATIONS OF ORTHOGONAL POLYNOMIALS PAUL NEVAI AND WALTER VAN ASSCHE We investigate orthogonal polynomials on the real line defined by a
More informationOn some tridiagonal k Toeplitz matrices: algebraic and analytical aspects. Applications
On some tridiagonal k Toeplitz matrices: algebraic and analytical aspects Applications R Álvarez-Nodarse1,, J Petronilho 3, and NR Quintero 4, 1 Dpto Análisis Matemático, Universidad de Sevilla, Apdo 1160,
More informationMultiple Orthogonal Polynomials
Summer school on OPSF, University of Kent 26 30 June, 2017 Plan of the course Plan of the course lecture 1: Definitions + basic properties Plan of the course lecture 1: Definitions + basic properties lecture
More informationMultiple Orthogonal Polynomials
Summer school on OPSF, University of Kent 26 30 June, 2017 Introduction For this course I assume everybody is familiar with the basic theory of orthogonal polynomials: Introduction For this course I assume
More informationMultiple orthogonal polynomials. Bessel weights
for modified Bessel weights KU Leuven, Belgium Madison WI, December 7, 2013 Classical orthogonal polynomials The (very) classical orthogonal polynomials are those of Jacobi, Laguerre and Hermite. Classical
More informationIntroduction to orthogonal polynomials. Michael Anshelevich
Introduction to orthogonal polynomials Michael Anshelevich November 6, 2003 µ = probability measure on R with finite moments m n (µ) = R xn dµ(x)
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis Volume 14, pp 127-141, 22 Copyright 22, ISSN 168-9613 ETNA etna@mcskentedu RECENT TRENDS ON ANALYTIC PROPERTIES OF MATRIX ORTHONORMAL POLYNOMIALS F MARCELLÁN
More informationReferences 167. dx n x2 =2
References 1. G. Akemann, J. Baik, P. Di Francesco (editors), The Oxford Handbook of Random Matrix Theory, Oxford University Press, Oxford, 2011. 2. G. Anderson, A. Guionnet, O. Zeitouni, An Introduction
More informationORTHOGONAL POLYNOMIALS OF DISCRETE VARIABLE AND BOUNDEDNESS OF DIRICHLET KERNEL
ORTHOGONAL POLYNOMIALS OF DISCRETE VARIABLE AND BOUNDEDNESS OF DIRICHLET KERNEL JOSEF OBERMAIER AND RYSZARD SZWARC Abstract. For orthogonal polynomials defined by compact Jacobi matrix with exponential
More informationKatholieke Universiteit Leuven Department of Computer Science
Separation of zeros of para-orthogonal rational functions A. Bultheel, P. González-Vera, E. Hendriksen, O. Njåstad Report TW 402, September 2004 Katholieke Universiteit Leuven Department of Computer Science
More informationOPSF, Random Matrices and Riemann-Hilbert problems
OPSF, Random Matrices and Riemann-Hilbert problems School on Orthogonal Polynomials in Approximation Theory and Mathematical Physics, ICMAT 23 27 October, 2017 Plan of the course lecture 1: Orthogonal
More informationTwo special equations: Bessel s and Legendre s equations. p Fourier-Bessel and Fourier-Legendre series. p
LECTURE 1 Table of Contents Two special equations: Bessel s and Legendre s equations. p. 259-268. Fourier-Bessel and Fourier-Legendre series. p. 453-460. Boundary value problems in other coordinate system.
More informationOberwolfach Preprints
Oberwolfach Preprints OWP 2013-23 CLEONICE F. BRACCIALI AND JUAN JOSÉ MORENO- BALCÁZAR Mehler-Heine Asymptotics of a Class of Generalized Hypergeometric Polynomials Mathematisches Forschungsinstitut Oberwolfach
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 9, 1999, pp. 39-52. Copyright 1999,. ISSN 1068-9613. ETNA QUADRATURE FORMULAS FOR RATIONAL FUNCTIONS F. CALA RODRIGUEZ, P. GONZALEZ VERA, AND M. JIMENEZ
More informationPainlevé VI and Hankel determinants for the generalized Jacobi weight
arxiv:98.558v2 [math.ca] 3 Nov 29 Painlevé VI and Hankel determinants for the generalized Jacobi weight D. Dai Department of Mathematics, City University of Hong Kong, Tat Chee Avenue, Kowloon, Hong Kong
More informationDepartment of Applied Mathematics Faculty of EEMCS. University of Twente. Memorandum No Birth-death processes with killing
Department of Applied Mathematics Faculty of EEMCS t University of Twente The Netherlands P.O. Box 27 75 AE Enschede The Netherlands Phone: +3-53-48934 Fax: +3-53-48934 Email: memo@math.utwente.nl www.math.utwente.nl/publications
More informationReconstruction of sparse Legendre and Gegenbauer expansions
Reconstruction of sparse Legendre and Gegenbauer expansions Daniel Potts Manfred Tasche We present a new deterministic algorithm for the reconstruction of sparse Legendre expansions from a small number
More informationGenerating Function: Multiple Orthogonal Polynomials
Applied Mathematical Sciences, Vol. 0, 06, no. 6, 76-77 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/0.988/ams.06.50669 Generating Function: Multiple Orthogonal Polynomials W. Carballosa Universidad
More informationULTRASPHERICAL TYPE GENERATING FUNCTIONS FOR ORTHOGONAL POLYNOMIALS
ULTRASPHERICAL TYPE GENERATING FUNCTIONS FOR ORTHOGONAL POLYNOMIALS arxiv:083666v [mathpr] 8 Jan 009 Abstract We characterize, up to a conjecture, probability distributions of finite all order moments
More informationRatio Asymptotics for General Orthogonal Polynomials
Ratio Asymptotics for General Orthogonal Polynomials Brian Simanek 1 (Caltech, USA) Arizona Spring School of Analysis and Mathematical Physics Tucson, AZ March 13, 2012 1 This material is based upon work
More informationGAUSS-KRONROD QUADRATURE FORMULAE A SURVEY OF FIFTY YEARS OF RESEARCH
Electronic Transactions on Numerical Analysis. Volume 45, pp. 371 404, 2016. Copyright c 2016,. ISSN 1068 9613. ETNA GAUSS-KRONROD QUADRATURE FORMULAE A SURVEY OF FIFTY YEARS OF RESEARCH SOTIRIOS E. NOTARIS
More informationINTEGER POWERS OF ANTI-BIDIAGONAL HANKEL MATRICES
Indian J Pure Appl Math, 49: 87-98, March 08 c Indian National Science Academy DOI: 0007/s36-08-056-9 INTEGER POWERS OF ANTI-BIDIAGONAL HANKEL MATRICES João Lita da Silva Department of Mathematics and
More informationIntroduction to Orthogonal Polynomials: Definition and basic properties
Introduction to Orthogonal Polynomials: Definition and basic properties Prof. Dr. Mama Foupouagnigni African Institute for Mathematical Sciences, Limbe, Cameroon and Department of Mathematics, Higher Teachers
More informationSpectral Theory of Orthogonal Polynomials
Spectral Theory of Orthogonal Polynomials Barry Simon IBM Professor of Mathematics and Theoretical Physics California Institute of Technology Pasadena, CA, U.S.A. Lectures 11 & 12: Selected Additional,
More informationA Right Inverse Of The Askey-Wilson Operator arxiv:math/ v1 [math.ca] 5 Oct 1993 B. Malcolm Brown and Mourad E. H. Ismail
A Right Inverse Of The Askey-Wilson Operator arxiv:math/9310219v1 [math.ca] 5 Oct 1993 B. Malcolm Brown and Mourad E. H. Ismail Abstract We establish an integral representation of a right inverse of the
More informationRational approximations to π and some other numbers
ACTA ARITHMETICA LXIII.4 (993) Rational approximations to π and some other numbers by Masayoshi Hata (Kyoto). Introduction. In 953 K. Mahler [2] gave a lower bound for rational approximations to π by showing
More informationSpecial Functions of Mathematical Physics
Arnold F. Nikiforov Vasilii B. Uvarov Special Functions of Mathematical Physics A Unified Introduction with Applications Translated from the Russian by Ralph P. Boas 1988 Birkhäuser Basel Boston Table
More informationNonlinear Integral Equation Formulation of Orthogonal Polynomials
Nonlinear Integral Equation Formulation of Orthogonal Polynomials Eli Ben-Naim Theory Division, Los Alamos National Laboratory with: Carl Bender (Washington University, St. Louis) C.M. Bender and E. Ben-Naim,
More informationFinite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product
Chapter 4 Hilbert Spaces 4.1 Inner Product Spaces Inner Product Space. A complex vector space E is called an inner product space (or a pre-hilbert space, or a unitary space) if there is a mapping (, )
More informationMultiple Meixner polynomials and non-hermitian oscillator Hamiltonians arxiv: v2 [math.ca] 8 Nov 2013
Multiple Meixner polynomials and non-hermitian oscillator Hamiltonians arxiv:1310.0982v2 [math.ca] 8 Nov 2013 F Ndayiragie 1 and W Van Assche 2 1 Département de Mathématiques, Université du Burundi, Campus
More informationINFINITE-DIMENSIONAL JACOBI MATRICES ASSOCIATED WITH JULIA SETS
proceedings of the american mathematical society Volume 88. Number 4. August 1983 INFINITE-DIMENSIONAL JACOBI MATRICES ASSOCIATED WITH JULIA SETS M. F. BARNSLEY1. J. S. GERÓNIMO2 AND A. N. HARRINGTON Abstract.
More informationHOMEWORK PROBLEMS FROM STRANG S LINEAR ALGEBRA AND ITS APPLICATIONS (4TH EDITION)
HOMEWORK PROBLEMS FROM STRANG S LINEAR ALGEBRA AND ITS APPLICATIONS (4TH EDITION) PROFESSOR STEVEN MILLER: BROWN UNIVERSITY: SPRING 2007 1. CHAPTER 1: MATRICES AND GAUSSIAN ELIMINATION Page 9, # 3: Describe
More informationInverse Polynomial Images which Consists of Two Jordan Arcs An Algebraic Solution
Inverse Polynomial Images which Consists of Two Jordan Arcs An Algebraic Solution Klaus Schiefermayr Abstract Inverse polynomial images of [ 1, 1], which consists of two Jordan arcs, are characterised
More informationOn the remainder term of Gauss Radau quadratures for analytic functions
Journal of Computational and Applied Mathematics 218 2008) 281 289 www.elsevier.com/locate/cam On the remainder term of Gauss Radau quadratures for analytic functions Gradimir V. Milovanović a,1, Miodrag
More informationCONVERGENCE OF APPROXIMATING FIXED POINTS FOR MULTIVALUED NONSELF-MAPPINGS IN BANACH SPACES. Jong Soo Jung. 1. Introduction
Korean J. Math. 16 (2008), No. 2, pp. 215 231 CONVERGENCE OF APPROXIMATING FIXED POINTS FOR MULTIVALUED NONSELF-MAPPINGS IN BANACH SPACES Jong Soo Jung Abstract. Let E be a uniformly convex Banach space
More informationml then either k + m < n or k + n < m and one of the above scalar
SIAM J. MATH. ANAL. Vol. 23, No. 4, pp. 959-964, July 1992 1992 Society for Industrial and Applied Mathematics 010 ORTHOGONAL POLYNOMIALS AND A DISCRETE BOUNDARY VALUE PROBLEM I* RYSZARD SZWARC" Abstract.
More informationORTHOGONAL POLYNOMIALS
ORTHOGONAL POLYNOMIALS 1. PRELUDE: THE VAN DER MONDE DETERMINANT The link between random matrix theory and the classical theory of orthogonal polynomials is van der Monde s determinant: 1 1 1 (1) n :=
More informationRecurrence Coef f icients of a New Generalization of the Meixner Polynomials
Symmetry, Integrability and Geometry: Methods and Applications SIGMA 7 (011), 068, 11 pages Recurrence Coef f icients of a New Generalization of the Meixner Polynomials Galina FILIPUK and Walter VAN ASSCHE
More informationJACOBI TYPE AND GEGENBAUER TYPE GENERALIZATION OF CERTAIN POLYNOMIALS. Mumtaz Ahmad Khan and Mohammad Asif. 1. Introduction
MATEMATIQKI VESNIK 64 (0) 47 58 June 0 originalni nauqni rad research paper JACOBI TYPE AND GEGENBAUER TYPE GENERALIZATION OF CERTAIN POLYNOMIALS Mumtaz Ahmad Khan and Mohammad Asif Abstract. This paper
More informationAN EXAMPLE OF SPECTRAL PHASE TRANSITION PHENOMENON IN A CLASS OF JACOBI MATRICES WITH PERIODICALLY MODULATED WEIGHTS
AN EXAMPLE OF SPECTRAL PHASE TRANSITION PHENOMENON IN A CLASS OF JACOBI MATRICES WITH PERIODICALLY MODULATED WEIGHTS SERGEY SIMONOV Abstract We consider self-adjoint unbounded Jacobi matrices with diagonal
More informationPUBLICATIONS DE L'INSTITUT MATHÉMATIQUE Nouvelle série, tome 58 (72), 1995, Slobodan Aljan»cić memorial volume AN ESTIMATE FOR COEFFICIENTS OF
PUBLICATIONS DE L'INSTITUT MATHÉMATIQUE Nouvelle série, tome 58 (72), 1995, 137 142 Slobodan Aljan»cić memorial volume AN ESTIMATE FOR COEFFICIENTS OF POLYNOMIALS IN L 2 NORM. II G. V. Milovanović and
More informationAdvanced Computational Fluid Dynamics AA215A Lecture 2 Approximation Theory. Antony Jameson
Advanced Computational Fluid Dynamics AA5A Lecture Approximation Theory Antony Jameson Winter Quarter, 6, Stanford, CA Last revised on January 7, 6 Contents Approximation Theory. Least Squares Approximation
More informationLaura Chihara* and Dennis Stanton**
ZEROS OF GENERALIZED KRAWTCHOUK POLYNOMIALS Laura Chihara* and Dennis Stanton** Abstract. The zeros of generalized Krawtchouk polynomials are studied. Some interlacing theorems for the zeros are given.
More informationChapter 1. Preliminaries. The purpose of this chapter is to provide some basic background information. Linear Space. Hilbert Space.
Chapter 1 Preliminaries The purpose of this chapter is to provide some basic background information. Linear Space Hilbert Space Basic Principles 1 2 Preliminaries Linear Space The notion of linear space
More informationMarkov s Inequality for Polynomials on Normed Linear Spaces Lawrence A. Harris
New Series Vol. 16, 2002, Fasc. 1-4 Markov s Inequality for Polynomials on Normed Linear Spaces Lawrence A. Harris This article is dedicated to the 70th anniversary of Acad. Bl. Sendov It is a longstanding
More informationINTRODUCTION TO PADÉ APPROXIMANTS. E. B. Saff Center for Constructive Approximation
INTRODUCTION TO PADÉ APPROXIMANTS E. B. Saff Center for Constructive Approximation H. Padé (1863-1953) Student of Hermite Histhesiswon French Academy of Sciences Prize C. Hermite (1822-1901) Used Padé
More informationINTERLACING PROPERTIES FOR HERMITIAN MATRICES WHOSE GRAPH IS A GIVEN TREE
INTERLACING PROPERTIES FOR HERMITIAN MATRICES WHOSE GRAPH IS A GIVEN TREE C M DA FONSECA Abstract We extend some interlacing properties of the eigenvalues of tridiagonal matrices to Hermitian matrices
More informationOn Gauss-type quadrature formulas with prescribed nodes anywhere on the real line
On Gauss-type quadrature formulas with prescribed nodes anywhere on the real line Adhemar Bultheel, Ruymán Cruz-Barroso,, Marc Van Barel Department of Computer Science, K.U.Leuven, Celestijnenlaan 2 A,
More informationSpectral analysis of two doubly infinite Jacobi operators
Spectral analysis of two doubly infinite Jacobi operators František Štampach jointly with Mourad E. H. Ismail Stockholm University Spectral Theory and Applications conference in memory of Boris Pavlov
More informationMeasures, orthogonal polynomials, and continued fractions. Michael Anshelevich
Measures, orthogonal polynomials, and continued fractions Michael Anshelevich November 7, 2008 MEASURES AND ORTHOGONAL POLYNOMIALS. MEASURES AND ORTHOGONAL POLYNOMIALS. µ a positive measure on R. 2 MEASURES
More informationAN INEQUALITY OF TURAN TYPE FOR JACOBI POLYNOMIALS
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 32, Number 2, April 1972 AN INEQUALITY OF TURAN TYPE FOR JACOBI POLYNOMIALS GEORGE GASPER Abstract. For Jacobi polynomials P^^Hx), a, ß> 1, let RAX)
More informationA trigonometric orthogonality with respect to a nonnegative Borel measure
Filomat 6:4 01), 689 696 DOI 10.98/FIL104689M Published by Faculty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat A trigonometric orthogonality with
More informationP AC COMMUTATORS AND THE R TRANSFORM
Communications on Stochastic Analysis Vol. 3, No. 1 (2009) 15-31 Serials Publications www.serialspublications.com P AC COMMUTATORS AND THE R TRANSFORM AUREL I. STAN Abstract. We develop an algorithmic
More informationMeasures, orthogonal polynomials, and continued fractions. Michael Anshelevich
Measures, orthogonal polynomials, and continued fractions Michael Anshelevich November 7, 2008 MEASURES AND ORTHOGONAL POLYNOMIALS. µ a positive measure on R. A linear functional µ[p ] = P (x) dµ(x), µ
More informationZeros and ratio asymptotics for matrix orthogonal polynomials
Zeros and ratio asymptotics for matrix orthogonal polynomials Steven Delvaux, Holger Dette August 25, 2011 Abstract Ratio asymptotics for matrix orthogonal polynomials with recurrence coefficients A n
More informationPainlevé equations and orthogonal polynomials
KU Leuven, Belgium Kapaev workshop, Ann Arbor MI, 28 August 2017 Contents Painlevé equations (discrete and continuous) appear at various places in the theory of orthogonal polynomials: Discrete Painlevé
More informationAsymptotics of Integrals of. Hermite Polynomials
Applied Mathematical Sciences, Vol. 4, 010, no. 61, 04-056 Asymptotics of Integrals of Hermite Polynomials R. B. Paris Division of Complex Systems University of Abertay Dundee Dundee DD1 1HG, UK R.Paris@abertay.ac.uk
More informationHankel Operators plus Orthogonal Polynomials. Yield Combinatorial Identities
Hanel Operators plus Orthogonal Polynomials Yield Combinatorial Identities E. A. Herman, Grinnell College Abstract: A Hanel operator H [h i+j ] can be factored as H MM, where M maps a space of L functions
More informationPower series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) a n z n. n=0
Lecture 22 Power series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) Recall a few facts about power series: a n z n This series in z is centered at z 0. Here z can
More informationOn Tricomi and Hermite-Tricomi Matrix Functions of Complex Variable
Communications in Mathematics and Applications Volume (0), Numbers -3, pp. 97 09 RGN Publications http://www.rgnpublications.com On Tricomi and Hermite-Tricomi Matrix Functions of Complex Variable A. Shehata
More informationOn some tridiagonal k-toeplitz matrices: Algebraic and analytical aspects. Applications
Journal of Computational and Applied Mathematics 184 (005) 518 537 wwwelseviercom/locate/cam On some tridiagonal k-toeplitz matrices: Algebraic and analytical aspects Applications R Álvarez-Nodarse a,b,,
More informationRecurrence Relations and Fast Algorithms
Recurrence Relations and Fast Algorithms Mark Tygert Research Report YALEU/DCS/RR-343 December 29, 2005 Abstract We construct fast algorithms for decomposing into and reconstructing from linear combinations
More informationHydrogen atom and orthogonal polynomials
AZA Preprint: 1990/7C. Hydrogen atom and orthogonal polynomials J.S. Dehesa, F. Domínguez Adame, E.R. Arriola, A. Zarzo Updated preprint version: October 1992 of Ref.: C. Brezinski, L. Gori and A. Ronveaux
More informationDeterminantal point processes and random matrix theory in a nutshell
Determinantal point processes and random matrix theory in a nutshell part I Manuela Girotti based on M. Girotti s PhD thesis and A. Kuijlaars notes from Les Houches Winter School 202 Contents Point Processes
More informationMarkov operators, classical orthogonal polynomial ensembles, and random matrices
Markov operators, classical orthogonal polynomial ensembles, and random matrices M. Ledoux, Institut de Mathématiques de Toulouse, France 5ecm Amsterdam, July 2008 recent study of random matrix and random
More informationSome Open Problems Concerning Orthogonal Polynomials on Fractals and Related Questions
Special issue dedicated to Annie Cuyt on the occasion of her 60th birthday, Volume 10 2017 Pages 13 17 Some Open Problems Concerning Orthogonal Polynomials on Fractals and Related Questions Gökalp Alpan
More informationJournal of Mathematical Analysis and Applications. Matrix differential equations and scalar polynomials satisfying higher order recursions
J Math Anal Appl 354 009 1 11 Contents lists available at ScienceDirect Journal of Mathematical Analysis and Applications wwwelseviercom/locate/jmaa Matrix differential equations and scalar polynomials
More informationDiscriminants of Polynomials Related to Chebyshev Polynomials: The Mutt and Jeff Syndrome
Discriminants of Polynomials Related to Chebyshev Polynomials: The Mutt and Jeff Syndrome Khang Tran University of Illinois at Urbana-Champaign Abstract The discriminants of certain polynomials related
More informationON SUM OF SQUARES DECOMPOSITION FOR A BIQUADRATIC MATRIX FUNCTION
Annales Univ. Sci. Budapest., Sect. Comp. 33 (2010) 273-284 ON SUM OF SQUARES DECOMPOSITION FOR A BIQUADRATIC MATRIX FUNCTION L. László (Budapest, Hungary) Dedicated to Professor Ferenc Schipp on his 70th
More informationClass notes: Approximation
Class notes: Approximation Introduction Vector spaces, linear independence, subspace The goal of Numerical Analysis is to compute approximations We want to approximate eg numbers in R or C vectors in R
More informationScattered Data Interpolation with Polynomial Precision and Conditionally Positive Definite Functions
Chapter 3 Scattered Data Interpolation with Polynomial Precision and Conditionally Positive Definite Functions 3.1 Scattered Data Interpolation with Polynomial Precision Sometimes the assumption on the
More informationTHE ZERO-DISTRIBUTION AND THE ASYMPTOTIC BEHAVIOR OF A FOURIER INTEGRAL. Haseo Ki and Young One Kim
THE ZERO-DISTRIBUTION AND THE ASYMPTOTIC BEHAVIOR OF A FOURIER INTEGRAL Haseo Ki Young One Kim Abstract. The zero-distribution of the Fourier integral Q(u)eP (u)+izu du, where P is a polynomial with leading
More informationSpectral Theory of Orthogonal Polynomials
Spectral Theory of Orthogonal Polynomials Barry Simon IBM Professor of Mathematics and Theoretical Physics California Institute of Technology Pasadena, CA, U.S.A. Lecture 1: Introduction and Overview Spectral
More informationTHE CIRCLE THEOREM AND RELATED THEOREMS FOR GAUSS-TYPE QUADRATURE RULES. Dedicated to Ed Saff on the occasion of his 60th birthday
THE CIRCLE THEOREM AND RELATED THEOREMS FOR GAUSS-TYPE QUADRATURE RULES WALTER GAUTSCHI Dedicated to Ed Saff on the occasion of his 6th birthday Abstract. In 96, P.J. Davis and P. Rabinowitz established
More informationON A WEIGHTED INTERPOLATION OF FUNCTIONS WITH CIRCULAR MAJORANT
ON A WEIGHTED INTERPOLATION OF FUNCTIONS WITH CIRCULAR MAJORANT Received: 31 July, 2008 Accepted: 06 February, 2009 Communicated by: SIMON J SMITH Department of Mathematics and Statistics La Trobe University,
More informationRELATION BETWEEN SMALL FUNCTIONS WITH DIFFERENTIAL POLYNOMIALS GENERATED BY MEROMORPHIC SOLUTIONS OF HIGHER ORDER LINEAR DIFFERENTIAL EQUATIONS
Kragujevac Journal of Mathematics Volume 38(1) (2014), Pages 147 161. RELATION BETWEEN SMALL FUNCTIONS WITH DIFFERENTIAL POLYNOMIALS GENERATED BY MEROMORPHIC SOLUTIONS OF HIGHER ORDER LINEAR DIFFERENTIAL
More informationPHYS 404 Lecture 1: Legendre Functions
PHYS 404 Lecture 1: Legendre Functions Dr. Vasileios Lempesis PHYS 404 - LECTURE 1 DR. V. LEMPESIS 1 Legendre Functions physical justification Legendre functions or Legendre polynomials are the solutions
More informationMultiple orthogonal polynomials associated with an exponential cubic weight
Multiple orthogonal polynomials associated with an exponential cubic weight Walter Van Assche Galina Filipu Lun Zhang Abstract We consider multiple orthogonal polynomials associated with the exponential
More informationGegenbauer Matrix Polynomials and Second Order Matrix. differential equations.
Gegenbauer Matrix Polynomials and Second Order Matrix Differential Equations Polinomios Matriciales de Gegenbauer y Ecuaciones Diferenciales Matriciales de Segundo Orden K. A. M. Sayyed, M. S. Metwally,
More informationarxiv:math/ v1 [math.ca] 6 Sep 1994
NUMERICAL COMPUTATION OF REAL OR COMPLEX arxiv:math/909227v1 [math.ca] 6 Sep 199 ELLIPTIC INTEGRALS B. C. Carlson Ames Laboratory and Department of Mathematics, Iowa State University, Ames, Iowa 50011-020,
More informationBernstein-Szegö Inequalities in Reproducing Kernel Hilbert Spaces ABSTRACT 1. INTRODUCTION
Malaysian Journal of Mathematical Sciences 6(2): 25-36 (202) Bernstein-Szegö Inequalities in Reproducing Kernel Hilbert Spaces Noli N. Reyes and Rosalio G. Artes Institute of Mathematics, University of
More information