polynomial.roots {orthopolynom}R Documentation

Create a list of polynomial roots

Description

This function returns a list with $n$ elements containing the roots of the order $k$ monic orthogonal polynomials for orders $k$ = 1, ... , $n$ using a data frame with the monic polynomial recurrence parameter vectors bf{a} and bf{b}

Usage

polynomial.roots(m.r)

Arguments

m.r monic recurrence data frame with parameters a and b

Details

The parameter m.r is a data frame with $n$+1 rows and two names columns. The columns which are names a and b correspond to the above referenced vectors. Function jacobi.matrices is used to create a list of symmetric, tridiagonal Jacobi matrices from these named columns. The eigenvalues of the k times k Jacobi matrix are the roots or zeros of the order $k$ monic orthogonal polynomial.

Value

A list with n elements each of which is a vector of polynomial roots

1 roots of the order 1 monic polynomial
2 roots of the order 2 monic polynomial
n roots of the order $n$ monic polynomial

Author(s)

Frederick Novomestky fnovomes@poly.edu

References

Press, et. al. (1992)

See Also

monic.polynomial.recurrences jacobi.matrices, eigen

Examples

r <- chebyshev.t.recurrences( 10, normalized=TRUE )
m.r <- monic.polynomial.recurrences( r )
p.roots <- polynomial.roots( m.r )

[Package orthopolynom version 1.0-1 Index]