create_bivariate_design function

Creates a design matrix from a bivariate smoothing algorithm