MathGroup Archive 2011

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: FixedPoint[Cos, 1.0]

  • To: mathgroup at smc.vnet.net
  • Subject: [mg119397] Re: FixedPoint[Cos, 1.0]
  • From: Bob Hanlon <hanlonr at cox.net>
  • Date: Wed, 1 Jun 2011 07:34:15 -0400 (EDT)
  • Reply-to: hanlonr at cox.net

What you meant was

Length@Union@Take[FixedPointList[Cos, 1.0, 1000], -10]

2

i.e., oscillating between two values.

Take[FixedPointList[Cos, 1.0, 1000], -3] // InputForm

{0.7390851332151608, 0.7390851332151605, 0.7390851332151608}

$Version

"8.0 for Mac OS X x86 (64-bit) (February 23, 2011)"


Bob Hanlon

---- "Szabolcs Horv=C3=A1t" <szhorvat at gmail.com> wrote:

=============
On 2011.05.31. 13:49, J Siehler wrote:
> My numerical analysis students this past term pointed out to me that
> the command
>
> FixedPoint[Cos, 1.0]
>
> which is the first example in the Information for FixedPoint doesn't
> ever finish running (or, not in any reasonable amount of time).  We
> tried this on serveral versions of Mathematica (6.x, 7.x, and 8.x),
> all on OS X, and got the same nonresult.
>
> It can be forced by specifying something like FixedPoint[Cos, 1.0,
> SameTest ->  (Abs[#1 - #2]<  10^-10&)] or FixedPoint[Cos, 1.0, 40],
> but it seems like that should be unnecessary, and the form given in
> the documentation should work.  Can anybody else confirm this
> behavior?  Am I missing something totally obvious here?
>

On Windows / Pentium M it finishes if I pass a machine number (1.0) to
it.  It does not finish if I pass it a number with fixed precision (e.g.
1.0`10) because the Precision of the result keeps increasing with
successive applications of Cos.

I'm curious, on your machine will the function keep oscillating between
two values?  What (and how many) values do you get if you do e.g.

Union@Take[Length@FixedPointList[Cos, 1.0, 1000], -10]

?

  -- Szabolcs



  • Prev by Date: Re: Mapping to Create Nested Loops
  • Next by Date: Re: Dt@x@1
  • Previous by thread: Re: FixedPoint[Cos, 1.0]
  • Next by thread: Re: FixedPoint[Cos, 1.0]