Re: How to use GNU MP in octave?

[Doug Stewart]

Javier Arantegui wrote:

> Hello,
>
> El Martes, 24 de Mayo de 2005 05:28, Keith Goodman escribió:
>  
>
> > Here's a quote from google
> >    
> [...]
>  
>
> > Let x = lg(8600!). Since lg(xy)=lg(x)+lg(y) you have
> >
> > x = sum(lg(k),k=1..86000)
> >
> > Define floor(x) as the integer part of x. Then
> >
> > 86000! = 10^(x-floor(x))*10^floor(x).
> >
> > You can evaluate x to the degree of accuracy you want, for example,
> > using Maple I got:
> >
> > 86000! = 7.91222558*10^372239
> >    
>
> Unfortunately, Octave doesn't buy the trick :-(:-(
>
> octave:4> k=1:1:86000;
> octave:5> x=sum(log10(k))
> x =  3.8702e+05
> octave:6> 10^(x-floor(x))*10^floor(x)
> ans = Inf
>
> Javier
> P.S. The truth is that I send this message to thank all the people who have be 
> so kind to answer my question about ^ and .^. Today I feel a little less 
> ignorant.
>
>
>  


I think Maple has a problem

I get
86000!= 2.55683586378090 * 10^387020

using Maxima.

Doug Stewart





-------------------------------------------------------------
Octave is freely available under the terms of the GNU GPL.

Octave's home on the web:  http://www.octave.org
How to fund new projects:  http://www.octave.org/funding.html
Subscription information:  http://www.octave.org/archive.html
-------------------------------------------------------------