Função para Ajuste Polinomial para Gretl

function [a,rmq]=num_ajuste_polinomial(p,x,y)

    // [a,rmq]=num_ajuste_polinomial(p,x,y)

    // Calcula os coeficientes de ajuste polinomial

    // g(x)=a_0+a_1*w*x+...+a_p*w*x**p de grau "p" para os pontos

    // (x_i,y_i) e o resíduo do ajuste de mínimos quadrados

    // (C) Rudnei Dias da Cunha 2014

    m = max(size(x))

    for j=0:p

        for i=0:p

            H(i+1,j+1) = sum(x.^(i+j))

        end

        b(j+1) = sum((x.^j).*y)

    end

    a = H\b

    rmq = sqrt(sum((num_avalia_polinomial(x,a)-y).^2))

endfunction