[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Bug-gsl] a bug in gsl - gamma distribution function
From: |
Gad Kimmel |
Subject: |
[Bug-gsl] a bug in gsl - gamma distribution function |
Date: |
Sun, 8 May 2005 11:13:01 +0300 (IDT) |
Hello,
I think you have a bug in gamma distribution function.
I am using Pentium 4, Debian Linux. I downloaded GSL from Debian
(testing).
gcc version: (GCC) 3.3.5 (Debian 1:3.3.5-12)
gsl version: 1.6-2
Here is the program:
#include <gsl/gsl_cdf.h>
#include <iostream>
using namespace std;
int main ()
{
cout << gsl_cdf_gamma_P (3,1,2) << endl;
return 0;
}
Now for the results:
$ g++ -o a a.cc -lgsl -lgslcblas
$ a
0.77687
While the true value of gsl_cdf_gamma_P (3,1,2) is 0.99752 and not
0.77687.
Please inform me if I'm correct.
best,
Gad.
- [Bug-gsl] a bug in gsl - gamma distribution function,
Gad Kimmel <=