/* @(#)Copyright (c), 1987, 1992 StatSci, Inc.  All rights reserved. */
static char whatssi[] = "@(#)randu.c version 3.10 created 10/31/92 ";

randu(seed,n,output)
long *seed, *n; float *output;
{
	long i = *n, j = *seed;

	while(i--){
		j = j*65539 & 0x7fffffff;
		*output++ = j/2147483648.;
	}
}

/*
randu <- function(n, seed=12345)
	.C("randu", as.integer(seed), as.integer(n), result=real(n))$result
*/