mpoly (version 1.1.1)

bernstein: Bernstein polynomials

Description

Bernstein polynomials

Usage

bernstein(k, n, indeterminate = "x")

Arguments

k

Bernstein polynomial k

n

Bernstein polynomial degree

indeterminate

indeterminate

Value

a mpoly object

Examples

# NOT RUN {
bernstein(0, 0)

bernstein(0, 1)
bernstein(1, 1)

bernstein(0, 1, "t")

bernstein(0:2, 2)
bernstein(0:3, 3)
bernstein(0:3, 3, "t")


bernstein(0:4, 4)
bernstein(0:10, 10)
bernstein(0:10, 10, "t")
bernstein(0:20, 20, "t")

# }
# NOT RUN {
  # visualize the bernstein polynomials

library(ggplot2); theme_set(theme_classic())
library(tidyr)

s <- seq(0, 1, length.out = 101)
N <- 10 # number of bernstein polynomials to plot
(bernPolys <- bernstein(0:N, N))

df <- data.frame(s, as.function(bernPolys)(s))
names(df) <- c("x", paste0("B_", 0:N))
head(df)

mdf <- gather(df, degree, value, -x)
head(mdf)

qplot(x, value, data = mdf, geom = "line", color = degree)

# }
# NOT RUN {

# }