nt1.m

牛顿插值法

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