Re^3: Fastest way to calculate hypergeometric distribution probabilities (i.e. BIG factorials)?

in reply to Re^2: Fastest way to calculate hypergeometric distribution probabilities (i.e. BIG factorials)?
in thread Fastest way to calculate hypergeometric distribution probabilities (i.e. BIG factorials)?

Good idea, but no cigar.

$ perl -MMath::Gsl::Sf=:Gamma -e'print gamma(300)'
gsl: gamma.c:1111: ERROR: overflow
Default GSL error handler invoked.
Aborted (core dumped)
$
[download]

The GSL doesn't handle large numbers.

There's also a math error to correct: gamma($n) == factorial($n - 1)

After Compline,
Zaxo

In Section Seekers of Perl Wisdom