Função para Cálculo pelo Método Avalia Polinômio para Gretl

function y=avalia_polinomio(coef,x,a)

    // y=avalia_polinomio(coef,x,a)

    // Avalia o polinomio interpolador de Newton, com

    // coeficientes "coef" e nós "x", no ponto "a",

    // usando multiplicações aninhadas (método de Horner)

    k = max(size(coef))-1;

    y = coef(k+1);

    for i=k-1:-1:0

        y = (a-x(i+1))*y+coef(i+1);

    end

endfunction