Modelo de Derivada da Função para Cálculo no Gretl

function dy=df_num(x)

    h = sqrt(number_properties('eps'))

    dy = (f(x+h)-f(x))/h

endfunction