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

function J=jacobiano(Func,x)

    // J=jacobiano(F,x)

    // Calcula uma aproximação numerica para o jacobiano de F, em x=[x_1;...;x_n].

    // F deve ser transmitida como um vetor de "strings", cada uma contendo a função escrita de

    // acordo com a sintaxe do Scilab. P.ex., para o sistema de equações não-lineares

    // {x^2-y^2 = 1,y = e^x}, escreve-se: ["x(1)^2-x(2)^2-1";"exp(x(1))-x(2)"].

    // C.T.Kelley, "Iterative Methods for Linear and Nonlinear Equations", Def. 5.4.1, pp.80

    // (C) Rudnei Dias da Cunha 2008

    n = max(size(x));

    J = zeros(n,n);

    for k=1:n

        J(:,k)=dernum(Func,x,k);

    end

endfunction