powered by
Construct the J basis functions for the splines evaluated on a grid.
B_basis(xgrid, knots)
n x (K+1) - matrix representing the value of each basis function evaluated on xgrid
grid where we want to evaluate the spline functions (vector of length n)
vector of knots for the splines (vector of length K)