Metodo Interpolador de Lagrange

function [y0,P,L] = lagrange_interp(xi, fi, x0)
% xi es el vector que contiene los puntos.
% fi es el vector que contiene los valores de los puntos evaluados en f(x)
% x0 es el punto para el cual se quiere interpolar
% y0 es la solucion por medio del Polinomio de Lagrange
    y0 = 0;
    n = length(xi);
    L=zeros(n,n);
    for j = 1 : n
            t = 1; %Este es para hallar el valor para el x0
            V=1; % Este es para hallar la matriz de Lagrange
            for i = 1 : n
                    if i~=j
                            t = t * (x0-xi(i))/(xi(j)-xi(i));
                            V=conv(V,poly(xi(i)))/(xi(j)-xi(i));
                    endif
            endfor
            y0 = y0 + t*fi(j);
            L(j,:)=V;
    endfor
    P=fi*L;
endfunction