Função para Cálculo pelo Coeficiente de Ajuste Linear para Gretl

function [a_0,a_1,rmq]=num_ajuste_linear(x,y)

    // [a_0,a_1,rmq]=num_ajuste_linear(x,y)

    // Calcula os coeficientes de ajuste linear g(x)=a_0+a_1*x

    // 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))

    H = [ m sum(x) ; sum(x) sum(x.^2) ]

    b = [ sum(y) ; sum(x.*y) ]

    a = H\b

    a_0 = a(1)

    a_1 = a(2)

    rmq = sqrt(sum((num_avalia_linear(x,a_0,a_1)-y).^2))

endfunction