Natural cubic spline interpolation