matlab拉格朗日插值多项式代码

个人用，大家也可以直接复制，直接调用即可

function yy = lag(x,y,xx)
    m=length(x);
    n=length(y);
 if m~=n,
     error('')
 end
    s=0;
 for i=1:n
    t=ones(1,length(xx))  ;
        for j=1:n
           if j~=i,
               t=t.*(xx-x(j))/(x(i)-x(j));
           end
        end
    s=s+t*y(i);  
 end
 yy=s;