double random_double_ab ( double a, double b )
{
  double r;
  double value;
  
  r = ( double ) rand ( ) / ( double ) RAND_MAX;

  value = a + ( b - a ) * r;

  return value;
}