Bug in fixed_point?!

C Barrington-Leigh cpblpublic at gmail.com
Tue Dec 21 23:33:01 EST 2010


On Dec 21, 9:36 am, Robert Kern <robert.k... at gmail.com> wrote:
> When you do, please provide the information that Terry Reedy asked for.
>

Sorry; quite right. For completeness I'll post here as well as over on
scipy.

Here's the actual code:
---------------------
from scipy import optimize
from math import exp
xxroot= optimize.fixed_point(lambda xx: exp(-2.0*xx)/2.0, 1.0,
args=(), xtol=1e-12, maxiter=500)
print ' %f solves fixed point, ie f(%f)=%f ?'%
(xxroot,xxroot,exp(-2.0*xxroot)/2.0)


Here is the output
------------------
Python 2.6.6 (r266:84292, Sep 15 2010, 16:22:56)
In [1]: run tmp.py
 0.332058 solves fixed point, ie f(0.332058)=0.257364 ?





More information about the Python-list mailing list