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

function [x,k,resultado]=newton_raphson_num(x,eps,del,kmax)

    // [x,k,resultado]=newton_raphson_num(x,eps,del,kmax)

    eps_m = number_properties('eps')

    for k=1:kmax

        f_x = f(x)

        df_x = df_num(x)

        printf('k=%g x=%g f(x)=%g df(x)=%g\n',k,x,f_x,df_x);

        if abs(df_x)<sqrt(eps_m) then

            resultado = -1

            break       

        end

        if abs(f_x)<eps then

            resultado = 0

            break

        end

        x_novo = x-f_x/df_x       

        if abs(x_novo-x)<del then

            resultado = -2

            break

        end

        x = x_novo

    end

    if k>kmax then

        resultado = -3

    end

endfunction