clear; 
N=input; 
x=-1:2/N:1; 
y=1/(1+25*x^2); 
Ln=0
for i=0:n 
for j=0:n
if j~=i lx=1,else lx=1;   
lx=lx.*(x-x(j))/(x(i)-x(j)); 
end 
Ln=y(i).*lx+Ln; 
end 
fplot('1/(25+x^2)',[-1,1]) 
hold on 
plot(x,Ln) 
xlabel('x') 
ylabel('f(x)') 
title('拉格朗日插值') 
这是关于多项式插值振荡的程序，为什么不行，我的程序是什么问题？非常感谢