nt2.m

牛顿插值法

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