/*
returns random double in range (0.0,1.0)
written by R. Geist
*/
#include <stdio.h>
#include <math.h>

double genrand();
long random();

double 
genrand()
{
static flag;
if (flag==0){
	srandom(123456789);
	srand(1);
	flag=1;
	}
return(((double)(random()+1))/2147483649.);
}