algoritmo y codigo en matlab para lagrange

% Algoritmo para calcular la interpolacion de lagrange
% n: numero de datos
% f: arreglo. Contiene las imagenes de de f(xi)
% x: arreglo. Comtiene los xi
% xa: punto(s) a evaluar en el polinomio de interpolacion
% i, j: Contadores o variables de control
% yres: reultado(s) de interpolacion
% w: variable para un calculo parcial
%
%
% 1.  INICIO
% 2.  DEFINIR n, y, x, i, j, xi, yi, w
% 3.  LEER n
% 4.  HACER DESDE i=1 HASTA i<=n
% 4.1.     HACER z= 1
% 4.2.     HACER DESDE j=1 HASTA j<=n
% 4.2.1.      SI i <> j ENTONCES
% 4.2.1.1.       HACER w = w * (xa - x(j)) / (x(i) - x(j))
%          FIN DE SI
% 4.2.2.      j = j + 1
%
%       FIN DE HACER DESDE
% 4.3.    yi = yi + w * f(i)
% 4.4.    i = i + 1
%     FIN DE REPETIR DESDE
% 5.  FIN
%
%
%   Funcion en MATLAB

function yi = lagrange_(x, y, xi)

yi = zeros(size(xi));
n= length(y);

for i=1:n
  w=ones(size(xi));
  for j=1:n
    if i ~= j,
       w = w .* (xi - x(j)) / (x(i) - x(j));
    end
  end
  yi = yi + w * y(i);
end