pspline function

Smoothing splines using a pspline basis