Hi there, It seems like when calling leastsq(), if the Jacobian matrix is organized as rows of partial derivatives (w.r.t. each variable), then no optimization is done at all -- the return value is the same as the initial guess. It only works when the matrix is columns of derivatives and set the parameter col_deriv to 1. --lin