MathGroup Archive 2006

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

Search the Archive

Re: Bug in calculating -3^49999.0


x=(-3)^49999.;

x/.z_Complex:>Re[z]*Chop[1+I*Im[z]/Re[z]]

-3.8513654349963041533518106689182`15.954589770191*^23855

x/. z_Complex:>Abs[z]*Chop[x/Abs[z]]

-3.85136543499630415335181066893086`15.47746851547134*^23855


Bob Hanlon
hanlonr at cox.net



On Apr 26, 2006, at 4:37 AM, kowald at molgen.mpg.de wrote:

> Hello everybody,
>
> I'm using Mathematica 5.1 under Win2k and I found a bug when I try  
> to compute
> -3^49999. I get:
>
> -3^49999 //N  =>  -3.85 * 10^23855          okay
> -3^49999.      =>  -3.85 * 10^23855          okay
> (-3)^49999.    =>  -3.85 * 10^23855 + 4.56*10^23844 i       wrong
> (-3)^49999     =>  -3.85 * 10^23855          okay
>
> Is this a known problem?
> Is there a work around ?
> Obviously this is part of a more complicated calculation and so I
> cannot simply leave out the brackets.
>
>  Many thanks,
>
>                 Axel
>


  • Prev by Date: Re: error function
  • Next by Date: Re: Setting up equations
  • Previous by thread: Bug in calculating -3^49999.0
  • Next by thread: Re: Bug in calculating -3^49999.0