diff -r 756c7a550542 scripts/optimization/fsolve.m --- a/scripts/optimization/fsolve.m Tue Feb 28 16:54:08 2017 -0800 +++ b/scripts/optimization/fsolve.m Sat Mar 04 11:26:59 2017 +0100 @@ -235,6 +235,9 @@ ## Calculate function value and Jacobian (possibly via FD). if (has_jac) [fvec, fjac] = fcn (reshape (x, xsiz)); + if all(size(fjac) ~= [m, n]) + error ("fsolve: jacobian has wrong size"); + endif ## If the Jacobian is sparse, disable Broyden updating. if (issparse (fjac)) updating = false;