Função para Avaliação Polinomial para Gretl

function y=num_avalia_polinomial(x,a)

    // y=num_avalia_polinomial(x,a)

    // Avalia o ajuste polinomial y(x)=a(0)+a(1)*x+a(2)*x**2+...+a(p)*x**p

    // (C) Rudnei Dias da Cunha 2014

    p = max(size(a))

    y = a(p)

    for i=p-1:-1:1 // Utiliza multiplicações aninhadas

        y = a(i)+x.*y

    end

endfunction