在MATLAB编程实现拉格朗日插值多项式

function Ys=langr(X,Y,Xs)
%拉格朗日插值
n=length(X);
nl=length(Xs);

for i=1:nl
x=Xs(i);
s=0;
for j=1:n
p=1;
for k=1:n
if k~=j
p=p*(x-X(k))/(X(j)-X(k));
end
end
s=s+p*Y(j);
end
Ys(i)=s;
end

这个应该行吧,对给定的一组数据,求插值点的值