splint function

Cubic spline interpolation