Re: How to use GNU MP in octave?

[Javier Arantegui]

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.


-- 
Javier Arántegui
Dept. Tecnologia de Alimentos / Dept. of Food Technology
Universitat de Lleida / University of Lleida (Spain)
         
Tel. +34 973702595
Fax  +34 973702596
IM: Jabber - javier.arantegui (AT) jabberes.org
http://www.tecal.udl.es



-------------------------------------------------------------
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
-------------------------------------------------------------