tp function

Generate basis functions for penalized spline smoothing.