4 #include <gsl/gsl_sf_bessel.h>
5 #include <gsl/gsl_rng.h>
6 #include <gsl/gsl_randist.h>
8 gsl_rng *rng; /* global generator */
11 const gsl_rng_type * T;
15 double y = gsl_sf_bessel_J0 (x);
16 printf ("J0(%g) = %.18e\n", x, y);
20 rng = gsl_rng_alloc (T);
21 seed = (unsigned long)getpid();
22 gsl_rng_set(rng, seed);
23 printf("generator type: %s\n", gsl_rng_name(rng));
24 printf("seed = %lu\n", gsl_rng_default_seed);
25 printf("first value = %lu\n", gsl_rng_get(rng));
27 /* number of first element */
28 r = gsl_ran_hypergeometric(rng, 5, 5, 3);
29 printf("rhg = %u\n", r);