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

function y=num_avalia_potencial(x,a_0,a_1)

    // y=num_avalia_potencial(x,a_0,a_1)

    // Avalia o ajuste potencial y(x)=a_0*x**a_1

    // (C) Rudnei Dias da Cunha 2014

    if (a_0>0.0) & (x>0.0) then

        y = a_0*(exp(a_1*log(x))) // y=a_0*x**a_1

    end

endfunction