powered by
bernstein(v,x,n)
K <- 3 index.b <- matrix(0:2^K) ## Bernstein polynomial B <- apply(index.b,1,bernstein,x,n=2^K)