smooth_spline function

Spline interpolation