我的一个ＭＡＴＬＡＢ牛顿插值程序 我还很菜　请多多指教　批准加入　谢谢！

				%Newton插值法
				function y=Newton(x1,y1,x);
				m=length(x);
				n=length(x1);
				for i=2:n
				    for j=n:-1:i
				        y1(j)=(y1(j)-y1(j-1))/(x1(j)-x1(j-i+1));
				    end
				end
				for j=1:m
				    p(j)=y1(1);
				    for k=2:n
				        l=1;   
				        for i=1:k-1
				            l=l*(x(j)-x1(i));
				       end
				    p(j)=p(j)+y1(k)*l;
				   end
				end
				re=[x'  p']