expand_spline function

Expand feature matrix using basis splines