create_bivariate_design: Creates a design matrix from a bivariate smoothing algorithm
Description
create_bivariate_design accepts two numeric vectors of equal length as inputs. From these inputs, a bivariate smoothing design matrix is produced using thin plate splines.