#include <stdlib.h>
#include <stdio.h>
#include <gsl/gsl_rng.h>
#include <gsl/gsl_cdf.h>
#include <gsl/gsl_randist.h>

#define MEAN 3.0

int main()
{
double mean = MEAN, sigma = 1.0, maxval = 5.0 ;
double v1, v2, v3 ;
gsl_rng *ran ;

ran = gsl_rng_alloc(gsl_rng_cmrg) ;
if(!ran)
    {
    printf("gsl_rng_alloc() failed\n") ;
    exit(0) ;
    }

v1 = gsl_cdf_gaussian_P(-mean,sigma) ;
v2 = gsl_cdf_gaussian_Q(maxval - mean,sigma) ;
v3 = gsl_cdf_poisson_P(5,mean) ;
/*
v3 = gsl_cdf_poisson_Q(5,mean) ;
*/

printf("v1=%f  v2=%f  v3=%f\n",v1,v2,v3) ;

exit(0) ;
}