Hello all,
Running 32 bit gsl on windows.
Got flipped results for incomplete gamma Functions.
#include <iomanip>
#include <iostream>
#include <gsl\gsl_specfunc.h>
using namespace std;
int main( int argc, char* argv[])
{
cout<<"GSL incomplete gamma=" << gsl_sf_gamma_inc(1.000000, 3.000000)<<
endl;
cout<<"GSL incomplete gamma P=" << gsl_sf_gamma_inc_P(1.000000, 3.000000)<<
endl;
return 0;
}
Output:
GSL incomplete gamma=0.0497871
GSL incomplete gamma P=0.950213
IMHO, GSL incomplete gamma should be equal to 0.950213
This is output of incomplete gamma P
Regards V.
