[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