[Bug-gsl] Incomplete Gamma Functions flipped?

[Vladimir Sakharuk]

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.