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

Sturla Molden sturla at molden.no
Thu Oct 1 09:01:08 EDT 2009


Pauli Virtanen skrev:
> 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.
>
>   

Thanks Pauli.

I differentiated the Michaelis-Menten function instead of the residuals. 
That goes for "incredibly stupid" then.

Sturla
 



More information about the SciPy-User mailing list