Re: Problems with gamma, lgamma, and alternatives

[Mike Miller]

On Tue, 25 Mar 2003, Mike Miller wrote:

> To use llgamma efficiently with giant numbers, you'd have to first
> transform to log(n):
>
> x = log(n)
> f(x) = x + log(x + 1)

ERROR!  That should have been...

f(x) = x + log(x - 1)

...and this...

> octave:1> x=log(6.384)+4356*log(10);
> octave:2> x+log(x+1)
> ans =  1.0041e+04

...should have had x+log(x-1), but the answer was the same.

Best,

Mike

-- 
Michael B. Miller, Ph.D.
Assistant Professor
Division of Epidemiology
and Institute of Human Genetics
University of Minnesota
http://taxa.epi.umn.edu/~mbmiller/



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