[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