Re: RealDigits bug?
- To: mathgroup at smc.vnet.net
- Subject: [mg45871] Re: RealDigits bug?
- From: adam.smith at hillsdale.edu (Adam Smith)
- Date: Wed, 28 Jan 2004 05:19:09 -0500 (EST)
- References: <bv00g1$je7$1@smc.vnet.net> <bv2g1t$68$1@smc.vnet.net>
- Sender: owner-wri-mathgroup at wolfram.com
Bobby,
What version of Mathematica are you using? I am running Version 4.1
and it does not include 10 in the Frequencies list
In[3]:=
Frequencies[First[RealDigits[ pisot]]]
Out[3]=
{{2089,0},{422,1},{427,2},{434,3},{467,4},{444,5},{452,6},{447,7},{437,
8},{381,9}}
In response to Enrique. I suspect that the failure is related to a
misunderstanding of the N[] function. In particular look in the Help
for N[] and study the $MaxMachineNumber and $MaxPrecision. Since your
pisot is on the order of 10^3340, you may be running into truncation
errors with N[]. On my machine $MaxMachineNumber is around 10^308.
Since pisot is clearly much greater than this, I would expect
truncation errors.
Adam Smith
drbob at bigfoot.com (Bobby R. Treat) wrote in message news:<bv2g1t$68$1 at smc.vnet.net>...
> I don't know whether the equality should be true, but it's very clear
> that 10 isn't a base-10 digit!
>
> Frequencies[First[RealDigits[
> pisot]]]
> {{2086, 0}, {422, 1},
> {427, 2}, {434, 3},
> {467, 4}, {444, 5},
> {452, 6}, {447, 7},
> {437, 8}, {381, 9}, {3, 10}}
>
> Bobby
>
> Enrique Zeleny <ezeleny at fismat1.fcfm.buap.mx> wrote in message news:<bv00g1$je7$1 at smc.vnet.net>...
> > This is wrong:
> >
> > In[1]:= pisot = N[((1/3)*((1/2)*(27 + 3*Sqrt[69]))^(1/3) + (2/(27 +
> > 3*Sqrt[69]))^(1/3))^27351, 6000];
> >
> > In[2]:= pisot == FromDigits[RealDigits[pisot]]
> >
> > Out[2]:= False
> >
> > If you look at the value returned by RealDigits[pisot] contains three
> > 10's!!!
> >
> >
> > Enrique Zeleny
> > Universidad Autonoma de Puebla
> > Mexico