[Bug-gsl] a bug in gsl - gamma distribution function

[Gad Kimmel]

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.