#include <math.h>
#include "gsl/gsl_rng.h"

int main()
{
  // return code
  int retval = 1;

  // create a generator
  gsl_rng *generator;
  generator = gsl_rng_alloc(gsl_rng_mt19937);

  // Read a value.
  double const Result = gsl_rng_uniform(generator);

  // Check value
  double const expectedResult( 0.999741748906672 );
  if( fabs( expectedResult - Result ) < 1.0e-6 )
    retval = 0;

  // free allocated memory
  gsl_rng_free(generator);
  return retval;
}
