/* random numbers using gsl routine
   gsl_rng_uniform
   and seeding it so you get different
   random numbers _every time_ you run.
   omitting the function call to seed 
   gsl_rng_uniform results in the same
   sequence of random numbers generated.
*/


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


using namespace::std;

int main()
{
  int N = 10;  // number of random numbers to generate
  
  // random nmbr declarations

  const gsl_rng_type * T;
  gsl_rng * r;

  T = gsl_rng_default;
  r = gsl_rng_alloc (T);


  // prepare for seeding of gsl random nmbr generator:

  srand((time(0)));                // srand & time are built-in
  unsigned long int s = random();  // gsl_rng_uniform will eventually
                                   // want a non-negative "long" integer
 

  gsl_rng_env_setup();
  gsl_rng_set(r,s); // seed the random number generator;

  // generate N random numbers:

  for(int i = 0; i < N; ++i) {
    cout << " random number = " << gsl_rng_uniform (r) << endl;;
    }

  gsl_rng_free (r);

  return 0;
}