htswanalysis/MACS/lib/gsl/gsl-1.11/doc/examples/randwalk.c

   1 #include <stdio.h>
   2 #include <gsl/gsl_rng.h>
   3 #include <gsl/gsl_randist.h>
   4
   5 int
   6 main (void)
   7 {
   8   int i;
   9   double x = 0, y = 0, dx, dy;
  10
  11   const gsl_rng_type * T;
  12   gsl_rng * r;
  13
  14   gsl_rng_env_setup();
  15   T = gsl_rng_default;
  16   r = gsl_rng_alloc (T);
  17
  18   printf ("%g %g\n", x, y);
  19
  20   for (i = 0; i < 10; i++)
  21     {
  22       gsl_ran_dir_2d (r, &dx, &dy);
  23       x += dx; y += dy;
  24       printf ("%g %g\n", x, y);
  25     }
  26
  27   gsl_rng_free (r);
  28   return 0;
  29 }