#The model:
model{
   for(i in 1:n){
 y[i]~dnorm(mu[i],taue)
      mu[i] <- alpha + beta*x[i]

   }
   
taue~dgamma(0.01,0.1)
   alpha~dnorm(0,0.01)
   beta~dnorm(0,0.01)
}