%--------------------------------------------------------------------------
% Data una lista de datos correspondientes a la variable independiente
% xdata y a la variable dependiente fdata, esta función calcula el valor del
% polinomio de Lagrange en xdata, pasando por todos los puntos (xdata,fdata)
%--------------------------------------------------------------------------

function finterpolation  =  LagrangianInterpolationOneD(xinterpolation,xdata,fdata)
  
n             =  length(xdata);
L             =  ones(n,1);
for i=1:n
  for j=1:n
      if abs(i-j)>0
        L(i)  =  L(i)*(xinterpolation-xdata(j))/(xdata(i)-xdata(j));
      end
  end              
end  
 
finterpolation    =  0; 
for i=1:n
  finterpolation  =  finterpolation + fdata(i)*L(i);
end
 
end
