poly_orth: Simpl orthogonal polynomials

View source: R/Orthogonal_polynomials.R

poly_orthR Documentation

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)))

PolynomF documentation built on May 29, 2024, 12:06 p.m.