double
ran1()
{
  /*
    return ( genunf (0, 1)); 
    */
  	double drand48();
	return( drand48() );
}