The Stenger conjectures and the A-stability of collocation Runge-Kutta methods

Stenger conjectures are claims about the location of the eigenvalues of matrices whose elements are certain integrals involving basic Lagrange interpolating polynomials supported on the zeros of orthogonal polynomials. In this paper, we show the validity of the extended Stenger conjecture for families of classical orthogonal polynomials. We also show the validity of the restricted Strenger conjecture for a family of Jacobi and generalized Laguerre orthogonal polynomials. A connection with the A-stability of the collocation Runge-Kutta methods is investigated.


Introduction
Given a non-negative weight function ω over an interval [a, b], -∞ ≤ a < b ≤ +∞, denote by π k , k ≥ 0, the orthogonal polynomials with respect to the scalar product f , g ω = b a f (x)g(x)ω(x) dx.

Extended Stenger Conjecture
Given an almost everywhere positive weight function ω on [a, b] such that -∞ ≤ a < b ≤ +∞, or equivalently, the associated orthogonal polynomials π n , the eigenvalues of each of the matrices U n and V n defined in (1.1) lie in the open right half of the complex plane.
The validity of the conjecture is established for Sinc interpolation in [3].The conjecture is shown to hold for Legendre polynomials in [4] and its equivalence to the A-stability of the n-stage Gauss-Runge-Kutta method is proved.Moreover, in [4] the conjecture is shown to hold for U n in the case of Jacobi orthogonal polynomials with parameters (1,0).In [4], the authors provided a counter-example to the conjecture for some piecewise constant weight functions.In this paper, we show the validity of the conjecture for a specific family of weight functions, which enables us in particular to prove the conjecture for a family of Jacobi, generalized Laguerre orthogonal polynomials and Jacobi-Koonwinder orthogonal polynomials.
In the restricted Stenger conjecture, the matrices U n (with a finite) and V n (with b finite) are defined by It has been shown in [4] that the restricted conjecture is wrong in general.For instance, they showed, at least numerically, that the real part of some eigenvalues of the matrix U n with n = 5 and x i , i = 1, 2, . . ., 5, the zeros of Gegenbauer polynomial C (10)  5 , are negative.However, several conjectures on the range of validity of the restricted Stenger conjecture for Jacobi and generalized Laguerre orthogonal polynomials were advanced in [4].In this paper, we show the validity of the restricted Stenger conjecture for a family of classical orthogonal polynomials.

Extended Stenger conjecture for a family of weight functions
In the following, we give a general result on the validity of the extended Stenger conjecture for a specific family of weight functions.The theorem and its proof were inspired by the methodology introduced in [4].
supposed to be positive on (a, b).Assume that, for any complex polynomial P of degree at most (n -1), there exist a complex polynomial Q of degree at most n and a real differentiable strictly increasing positive function on (a, b] meeting the following three requirements: (2.1) Then the extended Stenger conjecture associated with ω is satisfied for the matrix U n defined in (1.1).
Proof Let λ be an eigenvalue of the matrix U n .Thus, there exists a complex vector (y 1 , y 2 , . . ., y n ) = (0, 0, . . ., 0) (2.3) such that Denote by y the unique complex polynomial of degree at most (n -1) such that y(x i ) = y i .
We have According to condition (2.1), there exists a complex polynomial Q of degree at most n such that where is a real differentiable strictly increasing positive function on (a, b] and thus is non-negative on [a, b].Let ω i , i = 1, 2, . . ., n, be the weights of the Gaussian quadrature with respect to the weight function ω.Denote by ȳ the conjugate of y.From (2.4), we have Thus, according to (2.5), we have Applying Gaussian quadrature to the polynomial ȳQ of degree at most (2n -1), we obtain (2.6) The real part of the left-hand side of (2.6) is given by Moreover, using the product derivative formula, we note that Thus, the real part of the left-hand side of (2.6) is given by while the real part of the right-hand side of (2.6) is equal to Thus, from (2.6), to prove that the real part of λ is non-negative, it is sufficient to prove the non-negativity of (2.7).Since the function ), integration by parts shows that (2.7) is equal to Applying similar arguments, the matrices V n leads to the following theorem.

Theorem 2.2
Let ω be weight function on [a, b], -∞ < a < b ≤ +∞, supposed to be positive on (a, b).Assume that, for any complex polynomial P of degree at most (n -1), there exist a complex polynomial Q of degree at most n and a real differentiable strictly decreasing positive function on (a, b] meeting the following three requirements: Then the extended Stenger conjecture associated with ω is satisfied for the matrix V n defined in (1.1).

Application to Jacobi polynomials
In the following, we shall use Theorem 2.1 to prove the validity of the extended Stenger conjecture for a family of Jacobi orthogonal polynomials.Jacobi polynomials P (α,β) n are orthogonal polynomials on [-1, 1] with respect to the weight function They admit the following explicit expressions Lemma 3.1 Let α > -1 be a real number and P a polynomial of degree m.Then, x -1 where R is the polynomial .
Proof We proceed by induction of the degree of the polynomial P. For constant polynomials, the claim is trivial.Let P be a polynomial of degree n.Integration by parts leads to x -1 We conclude the proof using the induction hypothesis.

Corollary 3.1 The extended Stenger conjecture holds true in each of the following case:
-(a) for U n in the case of Jacobi polynomials Proof (a) We shall consider the cases m = 0, 1 separately.
• Case m = 1: According to Lemma 3.1, for any polynomial P of degree at most n -1, we have where Q is the polynomial of degree at most n given by Thus, the weight function ω (1,α) . The proof of (b) is a direct consequence of the fact that the eigenvalues of V n for P (α,β) n coincide with the eigenvalues of U n for P (β,α) n (see [4]).

Application to generalized Laguerre polynomials
In the following, we shall use Theorem 2.2 to prove the validity of the extended Stenger conjecture for a family of generalized Laguerre orthogonal polynomials.Generalized Laguerre polynomials L (α) n are orthogonal polynomials on [0, +∞) with respect to the weight function They admit explicit expressions as Lemma 4.1 Let P be a polynomial of degree m.Then, where R is the polynomial Proof We proceed by induction of the degree of the polynomial P. For constant polynomials, the claim is trivial.Let P be a polynomial of degree m.Integration by parts leads to +∞ x P(t)e -t dt = P(x)e -x + +∞ x P (t)e -t dt.
We conclude the proof using the induction hypothesis.

Corollary 4.1 The extended Stenger conjecture holds true for V n in the case of generalized
Laguerre polynomials L (m) n for all n ≥ 1 and m = 0, 1.
Proof We shall consider the cases m = 0 and 1 separately.
• Case m = 1: According to Lemma 3.1, for any polynomial of degree n -1, we have where Q is the polynomial of degree at most n Thus, the weight function ω (1) (t) = te -t satisfies the conditions of Theorem 2.2 with (x) = e -x .
Remark 4.1 Similar arguments as in the proof of Theorem 2.1 can be applied to show that the extended Stenger conjecture is valid for U n for the Legendre weight ω ≡ 1 and for V n for the Laguerre weight ω = e -x if, in the definition of the matrices U n and V n , we take the real numbers x 1 , x 2 , . . ., x n as the nodes of a (2n, n) positive quadrature.Similarly, the extended Stenger conjecture can be shown to hold when we take x 1 , x 2 , . . ., x n as the zeros of the othogonal polynomials with respect to the Sobolev inner product for U n for Jacobi weights ω (m,β) with m = 0 or 1 and β > -1 (Jacobi-Koornwinder polynomials [5]) and for V n for the Laguerre weights x m e -x with m = 0 or 1 (Laguerre-Koornwinder polynomials [6]) when we replace r 0 f (a)g(a) by r 0 f (b)g(b) in the Sobolev inner product.

Restricted Stenger conjecture
For the rest of the paper, we shall study the restricted Stenger conjecture.It has been shown in [4] that the conjecture is wrong in general.For instance, they showed, at least numerically, that the real part of some eigenvalues of the matrix U n with n = 5 and x i , i = 1, 2, . . ., 5, the zeros of Gegenbauer polynomial C (10)  5 , are negative.The validity of the restricted Stenger conjecture for Legendre polynomial is shown in [4] and its equivalence with the A-stability of the n-stage Gauss-Runge-Kutta method is established.Several conjectures on the range of validity of the restricted Stenger conjecture for Jacobi and generalized Laguerre orthogonal polynomials were advanced in [4].
In the following, we use the techniques introduced in [4] to prove the validity of the restricted Stenger conjecture for monotonous differentiable weight functions.Proof Let λ be an eigenvalue of the matrix U n .Thus there exists a non-zero complex vector y = (y 1 , y 2 , . . ., y n ) such that Denote by y the unique complex polynomial of degree at most (n -1) such that y(x i ) = y i .Thus,

The polynomial
is of degree at most n such that Q(x i ) = 0 for i = 1, 2, . . ., n.Therefore, there exists a complex constant K such that The real part of the left side of (5.2) is given by while the real part of right side of (5.2) is given by Thus, to show the non-negativity of Re(λ), we should prove the non-negativity of (5.In this case, one can use the methodology in [4] (page 5) to conclude that y is orthogonal to all polynomials of degree at most (n -1) and in particular b a y 2 (t) dt = 0 and thus y ≡ 0. This contradiction shows that Re(λ) > 0. Similar arguments can be used to prove the theorem for the matrices V n .
As an application of Theorem 5.1, we give instances for which the restricted Stenger conjecture holds for Jacobi and generalized Laguerre orthogonal polynomials.
The following is a direct consequence of Theorem 5. for U n in the case of generalized Laguerre polynomials L (α)  n for all n ≥ 1 and -1 < α ≤ 0.
Proof The statement (a) is a direct consequence of the fact that for α ≥ 0 and -1 < β ≤ 0 the Jacobi weight function (3.1) is a non-increasing functionon [-1, 1] since Statement (b) is a consequence of the fact proved in [4] that , i.e.; the matrice U n for the Jacobi polynomials with parameters (α, β) coincide with the matrix V n for Jacobi polynomials with parameters (β, α).To prove (c), we simply remark that, for -1 < α ≤ 0, the generalized Laguerre weight function is non-increasing on [0, +∞[, i.e., This concludes the proof of the corollary.

Restricted Stenger conjecture for Jacobi orthogonal polynomials
In this section, we show that the eigenvalues of the restricted matrices U n and V n defined in (1.2) coincide with the zeros of certain polynomials with coefficients expressed in terms of the value of the successive derivatives of the orthogonal polynomials at the endpoints of the interval.This allows us to restate the restricted Stenger conjecture for Jacobi orthogonal polynomials as a result already proved in [7] and permits us to improve the results of Corollary 5.1.To state the results in their full generality, we shall assume that in the definition of the restricted matrices U n and V n , the real numbers x 1 , x 2 , . . ., x n are arbitrary distinct real numbers in the interval [a, b].

Proposition 6.1 The eigenvalues of the restricted matrices U n and V n coincide with the zeros of the polynomials
respectively, where π n (x) = (xx 1 )(xx 2 ) . . .(xx n ) and π (k) n refer to the kth derivative of the polynomial π n .
Proof Let λ be an eigenvalue of the matrix U n .Thus, from (5.1), there exists a non-zero complex vector y = (y 1 , y 2 , . . ., y n ) and a complex number C such that where y the unique complex polynomial of degree n -1 such that y(x i ) = y i .Writing the polynomial y and the polynomial π n as and solving (6.2), we obtain the system Moreover, using the fact that b 0 = λb 1 + Cα 1 , we obtain Inserting this equation into the relation -λb 0 = Cα 0 , we obtain This proves that λ is a zero of the polynomial n .We shall now show that the eigenvalues of U n are simple, so that the spectrum of U n coincides with the zeros of n .Let y = (y 1 , y 2 , . . ., y n ) and z = (z 1 , z 2 , . . ., z n ) be two eigenvectors associated with the eigenvalue λ.Then, there exist two non-zero constants where y (resp.z) is the unique complex polynomial of degree (n -1) such that y(x i ) = y i (resp.z(x i ) = z i ).Thus, denoting F(x) = (y(x)/C 1 ) -(z(x)/C 2 ), we have Since F is a polynomial, this implies that F ≡ 0 and thus the vectors y and z are linearly dependent.The proof for V n follows similar arguments.Proposition 6.2An eigenvector (y 1 , y 2 , . . ., y n ) of the matrix U n (resp.V n ) associated with the eigenvalue λ is given by y i = y(x i ), where where λ is a zero of n (resp.n ).
Proof We should simply show that there is a constant C such that (6.2) holds.We have Integrating, we obtain This concludes the proof.
In connection with the tau approximation for an eigenvalue problem, Csordas et al. [7] proved the Hurwitz stability of the polynomials More precisely, they showed that if -1 < α ≤ 1 and β > -1, then the zeros of the polynomial (α,β) n lie in the half-plane Re(z) < 0; and thus by virtue of Proposition 6.1, the eigenvalues of V n for Jacobi polynomials lie in the half-plane Re(z) > 0.Moreover, since for Jacobi polynomials, we have U , we conclude Theorem 6.1 The restricted Stenger conjecture holds true in each of the following cases: • for U n in the case of Jacobi polynomials In particular, the previous theorem shows that restricted Stenger conjecture holds for Chebyshev polynomials of the first and second kind T n and S n with Also, Theorem 6.1 shows that the restricted Stenger conjecture holds for Legendre polynomials P (0,0) n and for Gegenbauer (or ultra-spherical) polynomials C (μ)  n when -1 2 < μ ≤ 3 2 , where

Restricted Stenger conjecture for generalized Laguerre polynomials
In the following, we show that the restricted Stenger conjecture for U n holds for generalized Laguerre polynomials L (α)  n when the parameter α is restricted to lie in the interval (-1, 1].An explicit expression of generalized Laguerre polynomials L (α)  n is given by Thus, the associated polynomials n according to Proposition 6.1 are given by To study the restricted Stenger conjecture for Laguerre polynomials, we shall study instead the zeros of the degree n truncation of the binomial function (1 + z) n+α , i.e., We shall need the following Eneström-Kakeya theorem [8].
Theorem 7.1 Let P n (z) = n k=0 a k z k be any polynomial with a k > 0 for all 0 ≤ k ≤ n.Setting then all the zeros of P n are contained in the annulus We introduce the polynomials 2) The connection between the polynomials B (α) n and H (α) n is given by the following.

Lemma 7.1
For any z = -1, we have

.3)
Proof We proceed by induction on n ≥ 1.The claim is trivial for n = 1.Let us assume that (7.3) holds for n.We have This concludes the proof.
Proposition 7.1 For 0 < α ≤ 1, all the zeros of the polynomial B (α) n lie in the half-plane Re(z) < 0.
Proof The condition 0 < α ≤ 1 guarantees that the coefficients of the polynomial H (α) n defined in (7.2) are positive.The successive ratios of the coefficients of H (α)  n are which, for 0 < α ≤ 1, are non-increasing as k runs from 0 to n -1.Thus by the Eneström-Kakeya Theorem 7.1, the zeros of H (α) n are contained in the annulus Since -1 is not a zero of B (α) n and in view of Lemma 7.1, the zeros of B (α) n satisfy Thus, we have and for each zero z of B α n , 1/z lies in a ring with center (-1, 0) and radii smaller than 1.Thus, each zero z of B α n satisfies Re(z) < 0. Remark 7.1 One can improve Corollary 7.1 by using the results in [6], where it is shown that for n ≥ 2 and -1 < α ≤ 1, the polynomial is Hurwitz stable, where L (α,M) n is the Laguerre-Koornwinder polynomials defined by Thus, in virtue of Proposition 6.1, the restricted Stenger conjecture for U n holds for Laguerre-Koornwinder polynomials when -1 < α ≤ 1.The advantage of the proof of Corollary 7.1 we have provided is that it can be applied to other orthogonal polynomials using Eneström-Kakeya theorem or the many existing variants of it.

Restricted Stenger conjecture and blossoming
The notion of polar form (or blossom) (Ramshaw, 1989) associated with a polynomial P of degree n is defined as follows: Definition 8.1 Let P be a complex polynomial of degree less or equal to n.There exists a unique multi-affine, symmetric function in n variables p: C n − → C such that for each z in C we have p(z, z, . . ., z) = P(z).The function p is called the polar form or the blossom of the polynomial P.
Consider the following bilinear de Boor-Fix operator on the space of polynomials of degree less than or equal to n defined as follows: for two given polynomials P, Q of degree less than or equal to n, we define The expression [P, Q] n is independent of τ and for any complex numbers u 1 , u 2 , . . ., u n and any degree n polynomial P, we have where p is the polar form of the polynomial P.
In the following, for generality, we assume that the x 1 , x 2 , . . ., x n in the definition of the resticted matrices U n and V n are arbitrary distinct real numbers in the interval [a, b].
Theorem 8. 1 The eigenvalues of the restricted matrices U n and V n coincide with the zeros of the polynomials and respectively, where f n is the polar form of F n ; the degree n truncation of the exponential function, i.e., Proof Denote by R n the polynomial We have R (k) n (a) = (-1) k z n-k .Thus, the polynomial n defined in (6.1) can be written as Therefore, from (8.2), we have n (z) = n!r n (x 1 , x 2 , . . ., x n ), where r n is the blossom of the polynomial R n .Moreover, from the definition of blossom, it is clear that This concludes the proof for U n .Similar arguments provide a proof for V n .
We need the well-known coincidence theorem of Walsh and thus the need to define circular regions in the complex planes.
A circular region of the complex plane is defined as the image of either the closed or the open unit disc under a non-singular Moebius maps γ (z) of the form Theorem 8.2 (Walsh coincidence theorem [9]).Let f be a symmetric multi-affine function of n complex variables and total degree n.Let u 1 , u 2 , . . ., u n be n complex numbers which lie in a circular region C.Then, there exists a ζ in C such that From the previous theorem, we deduce the following.Proof One can verify that the truncation polynomial of the exponential; F n ; for 1 ≤ n ≤ 4 has all its zeros in the circular region H -= {z ∈ C| Re(z) < 0}.Let us assume that the polynomial n has a zero η with Re(η) ≤ 0.Then, we have Since the complex numbers (ax i )/η, i = 1, 2, . . ., n, belong to the circular region H + = {z ∈ C| Re(z) ≥ 0}, by Walsh coincidence theorem, there exists ξ ∈ H + such that f n (ξ , ξ , . . ., ξ ) = F n (ξ ) = 0.This leads to a contradiction when 1 ≤ n ≤ 4. Similar arguments can be applied to prove the theorem for V n with 1 ≤ n ≤ 4.

Restricted Stenger conjecture and collocation Runge-Kutta methods
Let c 1 , c 2 , . . ., c n , be a set of n distinct collocation points in the interval [0, 1].The collocation method of degree n based on the points c i gives as a solution of the differential equation in each integration interval [t 0 , t 0 + h], as y 1 = P(t 0 + h), where the collocation polynomial P is defined by Given b = (b i ), A = (a i,j ) (i, j = 1, . . ., s) with b i and a i,j real numbers and let c i = s j=1 a i,j .An s-stage Runge-Kutta method for solving the initial-value problem is given by In [10], it is shown that a collocation method based on c 1 , c 2 , . . ., c n is equivalent to the Runge-Kutta method with matrix A and coefficient b given by where n i (τ ) are the elementary Lagrange interpolation polynomials based on c i .Moreover, the order of the collocation method coincides with the order of the underlying quadrature formula.The stability function of the Runge-Kutta method is given by and the stability region of the method is defined by A Runge-Kutta is said to be A-stable if its stability region satisfies H -⊂ S.
Theorem 9.1 The stability function of the collocation Runge-Kutta method with data (9.1) is given by where f n is the blossom of the degree n truncation of the exponential function F n .
Proof One should note that the matrix A is nothing but the matrix U n with x i = c i , i = 1, 2, . . .This concludes the proof.and R has no poles in the complex left half-plane, we conclude that H -⊂ S.This concludes the proof.
From the previous corollary and Theorem 6.1, we conclude the following.

Conclusion
In this paper, we have shown the validity of the extended and restricted Stenger conjectures for a family of weight functions with specific properties.As applications, we proved the validity of the conjecture for a family of Jacobi and generalized Laguerre orthogonal polynomials.We related the restricted Stenger conjecture to the A-stability of collocation Runge-Kutta method when the collocation points are symmetric with respect to the interval [0, 1].This enabled us to prove the A-stability of a large family of collocation methods.
dx, i, j = 1, 2, . . ., n, (1.2) and the conjecture states Restricted Stenger Conjecture Given an almost everywhere positive weight function ω on [a, b] or equivalently the associated orthogonal polynomials π n , the eigenvalues of each of the matrices U n and V n defined in (1.2) lie on the open right half of the complex plane.

Theorem 5 . 1
For differentiable non-increasing (resp.non-decreasing) weight functions ω over [a, b], all the eigenvalues of the matrices U n (resp.V n ) defined in (1.2) have positive real part.

4 )
Since the non-negative differentiable weight function ω is a non-increasing function on [a, b], the expression (5.4) is non-negative and therefore the real part of λ is non-negative.To prove that Re(λ) > 0, we should prove that (5.4) is strictly positive.We first note that (5.4) is equal to zero only when the weight function is a constant function and that b a y(t) dt = 0.(5.5)
where a, b, c, and d are complex numbers such that adbc = 0. Moebius maps are 1 -1 mapping of the extended plane into itself with the property of mapping every circle to either a circle or a line, and every line to either a circle or a line.A circular region is one of the following: an open disk, a closed disk, an open half-plane, a closed half-plane including ∞, the open exterior of a circle including ∞, or a closed exterior of a circle including ∞.

Corollary 8 . 1
The restricted Stenger conjecture holds for U n and V n , with 1 ≤ n ≤ 4, for any weight function ω.

Remark 8 . 1
In the proof of Corollary 8.1, the only used property of the zeros of the orthogonal polynomials is that they are real distinct and lie in the interval [a, b].Thus, the conclusion of Corollary 8.1 remains true if we define the restricted matrices U n and V n with arbitrary distinct real numbers x i , i = 1, 2, . . ., n, in the interval ]a, b[ with 1 ≤ n ≤ 4.

Figure 1 Corollary 9 . 1
Figure 1 Order-star configuration of Legendre (top left), Chebyshev first kind (top right), Chebyshev second kind (bottom left), and Gegenbauer with μ = 1.4 (bottom right) collocation Runge-Kutta methods for n = 5.The configurations show the A-stability of the methods