model;
{
   for( i in 1 : n ) {
      x[i] ~ dpois(rate[i])
   }
   for( i in 1 : n ) {
      lam[i] ~ dexp(gamma)
   }
   for( i in 1 : n ) {
      rate[i] <- lam[i] * t[i]
   }
   gamma ~ dgamma(0.001,0.001)
}
list(x=c(5, 1, 5, 14, 3 ,19, 1, 1, 4, 22),t=c(94.5, 15.7, 62.9, 126, 5.24, 31.4, 1.05, 1.05, 2.1, 10.5),n=10)
list(gamma=1,lam=c(1,1,1,1,1,1,1,1,1,1))