function yint = Interpolacion_Lagrange(xint, x, y) 

suma = 0;
N = length(x);
for n=1:N
    producto = 1;
    for m = 1:N
        if n~=m 
            producto = producto*(xint -x(m))/(x(n) - x(m));
            %disp([ xint, x(m), x(n), x(m), y]);
        end;
    end;
    disp(producto);
    suma = suma + producto*y(n);
end;
yint = suma;
end