natural_cspline function

Natural cubic spline