--- ray/src/common/hilbert.c	2011/02/18 18:41:04	3.2
+++ ray/src/common/hilbert.c	2011/04/08 23:23:28	3.3
@@ -1,5 +1,5 @@
 #ifndef lint
-static const char RCSid[] = "$Id: hilbert.c,v 3.2 2011/02/18 18:41:04 greg Exp $";
+static const char RCSid[] = "$Id: hilbert.c,v 3.3 2011/04/08 23:23:28 greg Exp $";
 #endif
 /* See LICENSE below for information on rights to use, modify and distribute
    this code. */
@@ -25,7 +25,7 @@ static const char RCSid[] = "$Id: hilbert.c,v 3.2 2011
  * from his, the first two interfaces and the style of some comments are very
  * much derived from his work. */
 
-
+#include <stdlib.h>
 #include "hilbert.h"
 
 /* implementation of the hilbert functions */
@@ -1256,9 +1256,10 @@ int main()
   printf("How many points? ");
   scanf("%d", &n);
   a = (double*) malloc(2*n*sizeof(double));
-  for (i = 0; i < n; ++i)
-    a[2*i] = drand48()-0.5, a[2*i+1] = drand48()-0.5;
-
+  for (i = 0; i < n; ++i) {
+    a[2*i] = rand()*(1./RAND_MAX) - 0.5;
+    a[2*i+1] = rand()*(1./RAND_MAX) - 0.5;
+  }
   qsort(a, n, 2*sizeof(double), cmp);
 
   for (i = 0; i < n; ++i)