#include <iostream>
#include <gsl/gsl_rng.h>
#include <gsl/gsl_randist.h>

using namespace std;

int main (void)
{
    
    const gsl_rng_type * T;
    gsl_rng * r;
    int i, n = 10;
    double mu = 3.0;
    /* create a generator chosen by the environment variable GSL_RNG_TYPE */
    gsl_rng_env_setup();
    T = gsl_rng_default;
    r = gsl_rng_alloc (T);
    
    /* print n random variates chosen from the poisson distribution with mean
    parameter mu */
    
    for (i = 0; i < n; i++){
        unsigned int k = gsl_ran_poisson (r, mu);
        cout << k << " ";
    }
    
    cout << endl;    

    system("PAUSE");

    return 0;
}