From: nico on
Is it possible to obtain the Jacobian with the function fminsearchbnd (like fminsearch) ?