三次样条函数及MATLAB

function [s0,s1,s2,s3]=cubic_spline(x,y) 
 
if any(size(x) ~= size(y)) || size(x,2) ~= 1
   error('inputs x and y must be column vectors of equal length');
end
 
n = length(x);
 
h = x(2:n) - x(1:n-1);
d = (y<