function[c,l]=lagran(x,y)
%x為n個節點的橫坐標所組成的向量,y為縱坐標組成的向量
%c為所得插值函數的系數組成的向量
%conv為多項式乘多項式函數
w=length(x);
n=w-1;
l=zeros(w,w);
for k=1:n+1
    v=1;
    for j=1:n+1
        if k~=j
            v=conv(v,poly(x(j)))/(x(k)-x(j));
        end
    end
    l(k,:)=v;
end
c=y*l;