nt1.m

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

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