| #include <stdio.h> |
| #include <gsl/gsl_rng.h> |
| #include <gsl/gsl_randist.h> |
| |
| 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); |
| printf (" %u", k); |
| } |
| |
| printf ("\n"); |
| gsl_rng_free (r); |
| return 0; |
| } |