%
% Computes forward model derivatives for the EM-38 problem.
%
function J=jac(sigma)
%
J=zeros(9,M);
I=eye(M);
h=1.0e-8;
for i=1:M
  J(:,i)=(fun(sigma+h*I(:,i))-fun(sigma))/h;
end
