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

function [matriz,termo_indep]=horner_quadratico(a,x,y)

  // [matriz,termo_indep]=horner_quadratico(a,x,y)

  // Aplica o m�todo quadr�tico de Horner sobre o polinomio cujos coeficientes

  // encontram-se no vetor "a", usando o m�todo de multiplica��es aninhadas de Horner.

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

  // (C) Rudnei Dias da Cunha 2009

  n = max(size(a))-1

  b(n+1) = a(n+1)

  b(n) = a(n)+x*b(n+1)

  for i=n-2:-1:0

      b(i+1) = a(i+1)+x*b(i+2)+y*b(i+3)

  end

  c(n+1) = b(n+1)

  c(n) = b(n)+x*c(n+1)

  for i=n-2:-1:1

      c(i+1) = b(i+1)+x*c(i+2)+y*c(i+3)

  end

  termo_indep = -[ b(1) ; b(2) ]

  matriz = [ c(2) c(3); c(3) c(4) ]

endfunction