poly_orth {PolynomF}

Simpl orthogonal polynomials

Description

Generate a list of polynomials up to a specified degree, orthogonal with respect to the natural inner product on a discrete, finite set of x-values with equal weights.

Usage

poly_orth(x, degree = length(unique(x)) - 1, norm = TRUE)

Arguments

x

A numeric vector

degree

The desired maximum degree

norm

Logical: should polynomials be normalised to length one?

Value

A list of orthogonal polynomials as a polylist object

Examples

x <- c(0:3, 5)
P <- poly_orth(x)
plot(P, lty = "solid")
Pf <- as.function(P)
zap(crossprod(Pf(x)))