#include <stdio.h>
#include <stdlib.h>

#include <gsl_rng.h>
#include <gsl_randist.h>
#include <gsl_cdf.h>

int main( int avgc, char **argv ) {

     double nu = 1.5, P, v;

     P = 1.93238e-01;
     v = gsl_cdf_chisq_Pinv(  P, nu );
     printf( "P(%.17e) = %.17e\n", P, v );

     P = 1.932383e-01;
     v = gsl_cdf_chisq_Pinv( P, nu );
     printf( "P(%.17e) = %.17e\n", P, v );

     P = 1.93238145206123590e-01;
     v = gsl_cdf_chisq_Pinv( P, nu);
     printf( "P(%.17e) = %.17e\n", P, v );

     exit( EXIT_SUCCESS );
}
