function retval = lsqjac(x)
% usage:  y = lsqjac(x)
% description: returns least squares
% vector function J(x), jacobian of lsqfn

% no bulletproofing
retval = [x(2)^3-7, 3*(x(1)+3)*x(2)^2; ...
          cos(x(2)*exp(x(1))-1)*x(2)*exp(x(1)), cos(x(2)*exp(x(1))-1)*exp(x(1))];