function A = GetJacobianInv2(x,p); k1 = size(x, 1); Id = eye(k1); [r] = Fcn2(x,p); f1 = r.out; x0 = x; for j = 1:k1-1; x = x0 + 1e-8*Id(:, j); [r] = Fcn2(x,p); J(:, j) = (r.out - f1)/1e-8; end A = inv(J);