#include <gsl/gsl_rng.h>

#include <gsl/gsl_sort_double.h>

int

main (void)

{

  const gsl_rng_type * T;

  gsl_rng * r;

  size_t i, k = 5, N = 100000;

  double * x = malloc (N * sizeof(double));

  double * small = malloc (k * sizeof(double));

  gsl_rng_env_setup();

  T = gsl_rng_default;

  r = gsl_rng_alloc (T);

  for (i = 0; i < N; i++)

    {

      x[i] = gsl_rng_uniform(r);

    }

  gsl_sort_smallest (small, k, x, 1, N);

  printf ("%zu smallest values from %zu\n", k, N);

  for (i = 0; i < k; i++)

    {

      printf ("%zu: %.18f\n", i, small[i]);

    }

  free (x);

  free (small);

  gsl_rng_free (r);

  return 0;

}