model{
  theta ~ dgamma(alpha,beta)
  for(i in 1:n){
    x[i] ~ dinterval(y[i],L)
    y[i] ~ dexp(theta)
  }  
}