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 mailing list [email protected] http://lists.gnu.org/mailman/listinfo/bug-gsl
