[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: Problems with gamma, lgamma, and alternatives

From: Mike Miller
Subject: Re: Problems with gamma, lgamma, and alternatives
Date: Tue, 25 Mar 2003 16:36:36 -0600 (CST)

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.



Michael B. Miller, Ph.D.
Assistant Professor
Division of Epidemiology
and Institute of Human Genetics
University of Minnesota

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

Octave's home on the web:
How to fund new projects:
Subscription information:

reply via email to

[Prev in Thread] Current Thread [Next in Thread]