nt2.m

来自「牛顿插值法」· M 代码 · 共 9 行

function r=nt2(x,y)
    m=length(x);
    h=0.1;
for i=2:m
    t=x(i)-h;
        s=spline(x,y,t);
        r(i)=(y(i)-s)/h; 
    end