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.
poly_orth(x, degree = length(unique(x)) - 1, norm = TRUE)
A list of orthogonal polynomials as a polylist object
A numeric vector
The desired maximum degree
Logical: should polynomials be normalised to length one?
x <- c(0:3, 5) P <- poly_orth(x) plot(P, lty = "solid") Pf <- as.function(P) zap(crossprod(Pf(x)))
Run the code above in your browser using DataLab