Função para Cálculo pelo Método Avalia Hornet para Gretl

function [y,z]=avalia_horner2(a,x)

    // [y,z]=avalia_horner2(a,x)

    // Avalia o polinomio cujos coeficientes encontram-se no vetor "a" em "x",

    // bem como a sua derivada, usando o método de multiplicações aninhadas de Horner.

    // O valor do polinômio é retornado em "y" e o da derivada em "z".

    // Os coeficientes devem ser armazenados em "a" na ordem a_0, a_1, ...

    // Se "x" é um número complexo a+bi, escreva-o na forma a+b*%i

    // (C) Rudnei Dias da Cunha 2009

    n = max(size(a))-1

    y = a(n+1)

    z = 0.

    for i=n:-1:1

        z = y+x*z

        y = a(i)+x*y

    end

endfunction