[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