msci Code

#include <stdio.h>
#include <sys/timeb.h>

#include <gsl/gsl_rng.h>

int main(int argc, char **argv) {
  long n;
  double num;
  gsl_rng *rann;
  long seed;
  timeb now;

  float min, max;

  if (argc != 2) {
    printf("Usage: random n\n");
    return 1;
  }

  sscanf(argv[1], "%d", &n);

  rann=gsl_rng_alloc(gsl_rng_mt19937);
  min=gsl_rng_min(rann);
  max=gsl_rng_max(rann);

  //printf("%g %g\n", min, max);

  ftime(&now);
  seed=(long) now.time + ((long) now.millitm) << 7;

  //printf("%d\n", seed);

  gsl_rng_set(rann, seed);

  for (long i=0; i<n; i++) {
    //num=(gsl_rng_get(rann)-min)/(max-min);
    num=gsl_rng_uniform(rann);
    printf("%lg\n", num);
  }

  gsl_rng_free(rann);

  return 0;

}