function y=LagrangePolynomial(x,xp,yp)

nn=1:length(xp);
y=zeros(1,length(x));
for i=nn
    nni=find(nn~=i);
    for j=1:length(x)
        y(j)=y(j)+yp(i)*prod((x(j)-xp(nni))./(xp(i)-xp(nni)));
    end
end