cs.m

来自「Lagrange插值、Newton插值函数等插值函数」· M 代码 · 共 8 行

function  d=cs(y,h)

for i=2:size(y,2)-1
    a=(y(i+1)-y(i))/h(i);
    b=(y(i)-y(i-1))/h(i);
    d(i-1)=6*(a-b)/(h(i)+h(i-1));
    
end