title: Can I get Jacobian from fminsearch
body: <p>I wonder if there is an easy way to get Jacobian out from fminsearch in Matlab ? like in 
  <code>
  [OptimizedParamters,residualsNorm,residual,exitflag,output,lambda,jacobian] = 
           lsqnonlin(@function,
               intialparamtersguess,lb,ub,options);
  </code></p>
  
  <p>I've tried 
  <code>
  options = optimset('MaxFunEvals',100,'Jacobian','on');
  [x,fval,exitflag,output] = fminsearch(fun,x0,options)
  </code></p>
  
  <p>but there is no Jacobin in the output</p>
  
  <p>any ideas please</p>