#include <stdio.h>
#include <math.h>

#include "random.h"

/* ****************************************************************
*  ************************************************************* */

double random_number()
{/* this gives a REAL number in the range of [0,1) */


  if (dseed == 0.0) printf("dseed = %f \n",dseed) ; 
  dseed = fmod(16807*(dseed),2147483647.0);
  return(dseed / 2147483655.0);

/*


 m = 2147483647.0 ; 
 a = 16807.0  ; 
 b = 127773.0 ; 
 c = 2836.0 ; 

 kl = dseed / b ; 
 dseed = a * fmod(dseed,b) - kl * c ; 
 if (dseed <= 0.0) dseed = dseed + m ; 
 return(dseed/2147483655.0) ; 
*/

 
}