function [yy,c]=lagrange(x,y,z)
nx=length(x);ny=length(y);
n=min(nx,ny);
yy=0;
for i=1:n
    p=1.0;
    for j=1:n
        if i～=j
            if abs(x(i)-x(j))<eps
                error('error');
                return;
            end
         ll=[1.0-x(j)]/(x(i)-x(j));
         p=conv(p,ll);
     end
 end
 yy=pplus(yy,p*y(i));
end
if nargout=2
    c=polyval(yy,2);
end
x=[1 2 3 4 5]
y=1.30 1.80 2.20 2.9 3.5]
p=lagrange(x,y) 




运行不了。据说是少个主程序去运行它。。谁帮我看看，怎么整阿。。