function Nevill(x,y,n,xx)
% Nevill插值算法地實現
% x=[3.60 3.65 3.70 3.75];
% y=[36.598 38.475 40.447 42.521];
% n=4;
% xx=3.62;
A=zeros(n,n);h=0.05;
for j=1:n
    A(j,1)=y(j);
end
for j=1:n
    for i=1:j
        A(j,i)=((xx-x(i-1))*A(j-1,i-1)-(xx-x(i))*A(j,i-1))/(x(i)-x(i-1));
    end
    A(j,i)
end

%A
%f=A(1,1);
%for j=2:n
%    tmp=1;
%    for i=1:j-1
%        tmp=tmp*(xx-x(i))/(h*i);%jiecheng xx-x(i)
%    end
%    f
%    f=f+A(j,j)*tmp;
%end
fprintf('x         ');
for j=0:n-1
    fprintf('   %d階Nevill',j);
end
fprintf('\n');
for j=1:n
    fprintf('%-12f',x(j));
    for i=1:j
        fprintf('%-12f',A(j,i));
    end
    fprintf('\n');
end
%fprintf('N(%d)(%.2f)=%f',n-1,xx,f);