function [yp] = Lagrange(x,y,xp)
% Lagrange interpolation
n=length(x)-1; % Βαθμός πολυωνύμου
% Εύρεση της τιμής των Li(x)στο xp
for i=1:n+1 % n+1 σημεία παρεμβολής
 L(i)=1; 
 for j=1:n+1 
 if i~=j 
 L(i)=L(i)*(xp-x(j))/(x(i)-x(j)); 
 end
 end
end
% Εύρεση της τιμής του πολυωνύμου για x=xp
yp=0;
for i=1:n+1
 yp=yp+L(i)*y(i);
end