Mathematica 9 is now available
Services & Resources / Wolfram Forums / MathGroup Archive
-----

MathGroup Archive 2011

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

Search the Archive

Re: FixedPoint[Cos, 1.0]

  • To: mathgroup at smc.vnet.net
  • Subject: [mg119450] Re: FixedPoint[Cos, 1.0]
  • From: Dana DeLouis <dana.del at gmail.com>
  • Date: Sat, 4 Jun 2011 06:21:09 -0400 (EDT)

Hi.  Just to mention...
On your two tests that you used:
FixedPoint[Cos, 1.0, SameTest -> (Abs[#1 - #2] < 10^-10 &)]
FixedPoint[Cos, 1.0, 40],

One of many alternatives to get a little better precision might be:

FixedPoint[Cos,1.0,SameTest->(Equal[##]&)]
0.7390851332151572`

Not great, but it=92s a little better.
We know from other posts that the value cycles between:
0.7390851332151605`
0.7390851332151608`

Another SameTest example similar to yours:

FixedPoint[Cos,1.0,SameTest->(Abs[Subtract[##]]< Power[10.5,-15]&)]
0.7390851332151605`

The two numbers above end in 605,and 608.
Interesting is that the numbers ending in 606 and 607 stop at 607.

n=0.7390851332151606`  ;
FixedPointList[Cos,n,10] 
{0.7390851332151606`,0.7390851332151607`}

n=0.7390851332151607`  ;
FixedPointList[Cos,n,10] 
{0.7390851332151607`,0.7390851332151607`}

The number ending in 607 seems to be the answer:
Cos[0.7390851332151607`]
    0.7390851332151607`

Higher precision seems to indicate 607 also.

FixedPoint[Cos,1`30,100]//N
0.7390851332151607`

= = = = = = = = = =
HTH  : >)
Dana DeLouis



On May 31, 7:49 am, J Siehler <jsieh... at gmail.com> 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?




  • Prev by Date: Running P&L
  • Next by Date: Re: ContourPlot ColorFunction Question
  • Previous by thread: Re: FixedPoint[Cos, 1.0]
  • Next by thread: Re: FixedPoint[Cos, 1.0]