[SciPy-User] optimize.leastsq does not converge with full Jacobian?

Pauli Virtanen pav+sp@iki...
Thu Oct 1 02:41:15 CDT 2009


Thu, 01 Oct 2009 09:12:19 +0200, Sturla Molden wrote:
[clip]
> def residuals(t, x, y):
>     return y - michaelis_menten(t, x)
[clip]
> def jacobian(t, x, y):
>     j = np.zeros((2,x.shape[0]))
>     Vmax, Km = t
>     j[0,:] = x/(Km + x)
>     j[1,:] = -Vmax*x/((Km + x)**2)
>     return j

Sign error.

-- 
Pauli Virtanen



More information about the SciPy-User mailing list