Metode numerice - Lab09

function rez = Neville(x, f, y)
    n = length(x);
    m = length(f);
    temp = f;
    for j=2:n
        for i = n:-1:j
            temp(i) = ((y - x(i-j+1))*temp(i) - (y-x(i))*temp(i-1))/ ...
                ( x(i) - x(i-j+1);
        endfor
    endfor
    rez = temp(n);
endfunction