myinterp1

function yy= myinterp1(x,y,xx)
%本次实验我写的是牛顿法差值多项式 
% x 是节点的横坐标向量， y 是纵坐标向量 得到在x=xx处的，牛顿差值多项的yy
n=length(x);  
%定义商差表
A=zeros(n);  
%差商表第一列为y   
for i=1:n  
    A(i,1)=y(i);
end
 %j为列标
for j=2:n  
%i为行标                                            
    for i=j:n
 %计算差商表
         A(i,j)=(A(i,j-1)-A(i-1,j-1))/(x(i)-x(i-j+1));  
    end
end
c=diag(A); 
%根据差商表,求对应的牛顿插值多项式在x=xx处的值yy
%将x=xx带入牛顿插值多项式，得到的yy的值
yy=nest(c(end:-1:1),xx,x(end:-1:1));   
end