Blame doc/examples/randwalk.c

Packit 67cb25
#include <stdio.h>
Packit 67cb25
#include <gsl/gsl_rng.h>
Packit 67cb25
#include <gsl/gsl_randist.h>
Packit 67cb25
Packit 67cb25
int
Packit 67cb25
main (void)
Packit 67cb25
{
Packit 67cb25
  int i;
Packit 67cb25
  double x = 0, y = 0, dx, dy;
Packit 67cb25
Packit 67cb25
  const gsl_rng_type * T;
Packit 67cb25
  gsl_rng * r;
Packit 67cb25
Packit 67cb25
  gsl_rng_env_setup();
Packit 67cb25
  T = gsl_rng_default;
Packit 67cb25
  r = gsl_rng_alloc (T);
Packit 67cb25
Packit 67cb25
  printf ("%g %g\n", x, y);
Packit 67cb25
Packit 67cb25
  for (i = 0; i < 10; i++)
Packit 67cb25
    {
Packit 67cb25
      gsl_ran_dir_2d (r, &dx, &dy;;
Packit 67cb25
      x += dx; y += dy; 
Packit 67cb25
      printf ("%g %g\n", x, y);
Packit 67cb25
    }
Packit 67cb25
Packit 67cb25
  gsl_rng_free (r);
Packit 67cb25
  return 0;
Packit 67cb25
}