[SciPy-Dev] Comments on optimize.newton function

[Gökhan Sever]

On Sun, May 22, 2011 at 11:12 AM,  <josef.pktd at gmail.com> wrote:
> I proposed taking logs of left hand side and right hand side
> separately, since you are just looking for a zero,
> with myfunc.log(),  it is not simplified

OK, I have gotten this right this time:

myfunc =  rh - exp(kelvin/x) * (x**3 - rd**3) / (x**3 - rd**3 * (1.0 - kappa))
myfunc2 = log(myfunc) = log(rh) - (kelvin/x) + log(x**3 - rd**3) -
log(x**3 - rd**3*(1.0 - kappa))

myfunc2_prime = -3*x**2/(rd**3 - x**3) - 3*x**2/((kappa - 1.0)*rd**3 +
x**3) + kelvin/x**2

How can I proceed this point onwards?


>
> (I don't have a quick way to do the symbolic derivative, but there
> shouldn'd be any exp left in the expression)