[SciPy-User] Simple root solving issues
Gökhan Sever
gokhansever at gmail.com
Sun Oct 24 17:38:15 EDT 2010
On Sun, Oct 24, 2010 at 4:29 PM, Warren Weckesser
<warren.weckesser at enthought.com> wrote:
>
> Your function is basically c1 - c2/log(gsd), so yeah, it will have a problem
> if the initial guess is 1.0. It works fine if the initial guess is, say,
> 1.01.
Yes, catch this better now. Interestingly it suffers when x0 < 1.00
too. Probably log(gsd<1) explodes the function by making the right
term positive and restricting the function to find a zero.
--
Gökhan
More information about the SciPy-User
mailing list