spl1D function

Fit P-splines