model;
{
   g ~ dgamma(0.01,0.01)I(1.0E-31,1.0E+31)
   x <- round(g*(pow(10,-round(log(g)/log(10)-0.5)))-0.5)
   u ~ dunif(0,1)
   y0 <- x*u
   y <- round(y0*(pow(10,-round(log(y0)/log(10)-0.5)))-0.5)
}

data;
list()

#  Program runs with above (empty) data set.  Data in book below.

data;
list(a=c(28778, 453, 83849, 30513, 110912, 9251, 127869, 43069, 1399, 337032,
    547026, 108178, 248577, 6, 131944, 93030, 103000, 70283, 301225, 157, 2586,
    316, 1, 40844, 324219, 312677, 92082, 237500, 61, 504782, 449964, 41293,
    23623, 130439, 20768, 78772, 14121, 5571000, 0.44, 255804),
f=c(10, 7, 6, 6, 3, 2, 2, 1, 3))