function [x,y,C]=newdunchashang(n)
% 计算插商数据点 x y
y=zeros(n+1,1);
x=zeros(n+1,1);
for i=0:1:n
x(i+1,1)=-1+(2*i)/n;
y(i+1,1)=1/(1+25*x(i+1,1).^2);
end
plot(x,y,'b');
grid on
% 计算牛顿插值多项式系数
n=length(x);
C=zeros(n,1);
for j=2:n %别写成j=1:n,这样i-1=0,变量指标不能为0
for i=n:(-1):j
y(i,1)=(y(i,1)-y(i-1,1))./(x(i,1)-x(i-j+1,1)); %计算牛顿插商表
end
end
C=y(i,1);
%