%x = [x(1) ; x(2) ; x(3)]; %J = zeros(3,3); %J(:,1) = (f(x+[h ; 0 ; 0]) - f(x-[h ; 0 ; 0]))./(2*h); %J(:,2) = (f(x+[0 ; h ; 0]) - f(x-[0 ; h ;...

%x = [x(1) ; x(2) ; x(3)];

%J = zeros(3,3);

%J(:,1) = (f(x+[h ; 0 ; 0]) - f(x-[h ; 0 ; 0]))./(2*h);

%J(:,2) = (f(x+[0 ; h ; 0]) - f(x-[0 ; h ; 0]))./(2*h);

%J(:,3) = (f(x+[0 ; 0 ; h]) - f(x-[0 ; 0 ; h]))./(2*h);

matlab, i just need to change this into a loop