> > > > I tried > > def f(x) : > > return 1 - x**2 > > > This works for me.. fixed_point(f, 0.5, xtol=1e-12) 0.6180339887498949 >>> f(_) 0.61803398874989479 It looks like it found the place where f(x)=x pretty well. -Travis