[seqfan] Re: Extension of A097486, Calculating Pi in the Mandelbrot Set

Hans Havermann pxp at rogers.com
Sat Jan 9 03:38:58 CET 2010


Robert Munafo:

> I discovered that the last term in the current entry, A(7)=31415928,  
> is off
> by one...

> Can anyone with Mathematica or Maple or Matlab verify?

> First 11 terms: 3, 33, 315, 3143, 31417, 314160, 3141593, 31415927,
> 314159266, 3141592655, 31415926537

It took me a while to figure out how to squeeze the extra precision  
out of the coding in Mathematica but I think/hope that the following,  
which utilizes 128-digit accuracy, does the job:

$MinPrecision=128; Do[c=SetPrecision[.1^n*I-.75,128]; z=c; a=0;  
While[Abs[z]<4, z=z^2+c; a++]; Print[a], {n,0,8}]

3
33
315
3143
31417
314160
3141593
31415927
314159267

So, yes and no.




More information about the SeqFan mailing list