数值计算方法牛顿插值matlab实现

xn=[0,1,2,4]';
yn=[3,6,11,51]';
x=0.5;
n=length(xn);
A(:,1)=yn;
for j=2:n
    for i=j:n
        A(i,j)=(A(i,j-1)-A(i-1,j-1))/(xn(i)-xn(i-j+1));
    end
end
y=0;
for i=1:n
    t=1;
    for j=1:i-1
        t=t*(x-xn(j));
    end
    y=y+A(i,i)*t;
end
y